Лаба №2, Вариант 7
Описание файла
PDF-файл из архива "Лаба №2, Вариант 7", который расположен в категории "". Всё это находится в предмете "радиопередающие устройства" из 6 семестр, которые можно найти в файловом архиве РТУ МИРЭА. Не смотря на прямую связь этого архива с РТУ МИРЭА, его также можно найти и в других разделах. Архив можно найти в разделе "лабораторные работы", в предмете "радиопередающие устройства" в общих файлах.
Просмотр PDF-файла онлайн
Текст из PDF
Вариант 7Для решения задач 1.1; 1.2 и 3.2 был создан скрипт (программа) позволяющий решитьзадание согласно варианту с одновременным выводом результатов в «Command Window»Матлаба и построением нужных графиков.Листинг программы:clc %очищаем экран'---1.1 Расчет резонансной частоты колебательного контура----'variant=7 %выбор вариантаL=[1E-4 4E-3 0.5 5E-6 0.02 1E-5 2.5 0.01]; %массив индуктивностейC=[1E-10 3E-9 2E-6 1.5E-11 6.8E-9 3.3E-11 1E-5 2E-9]; %массив емкостейf0=1/(2*pi*sqrt(L(variant)*C(variant)));T=1/f0;'L=', L(variant)'C=', C(variant)f0T'-----------------------------------------------------------------''-------1.2 Построение резонансной кривой----------'U=1; R=10; %входное напряжение и сопротивление потерьstart=(f0/1.1); %задаем начальную точкуstep=(f0/400); % разрешение графика (шаг)finish=(1.1*f0); %задаем конечную точкуf= start:step:finish ; %вектор частот нач.знач.:шаг:конеч.значw=2*pi*f; %вектор круговых частотI =U./sqrt(R^2+(w*L(variant)-1./(w*C(variant))).^2);figure; plot(f,I); grid; xlabel('\omega'); ylabel('\it I');title('Resonance curve');end'---------1.3 Операции с матрицами.--------------''Зададим матрицы случайных чисел'B=rand(3,4) %матрица 3 строки 4 столбцаX=rand(4,3) %матрица 4 строки 3 столбцаA=rand(4,4) %матрица 4*4 (т.к.
в результате операций с В и Х получится матрица аналогичногоразмера)Bt=B'Xt=X'BtXt=Bt*XtK=(A+BtXt)K=inv(K)'------3.2 Спектральный анализ сигналов.------'k=199 %число точекt=0:1/k:1; %вектор времениx=3*sin(2*pi*10*t)+2*sin(2*pi*50*t); %вектор сигналовy=fft(x); %вектор БПФ сигналовm=abs(y); %вектор амплитудf=(0:length(y)-1)*k/length(y); % вектор частотfiguretitle('Frequency analysis');subplot(2,1,1); plot(t,x); grid; xlabel('\it t');ylabel('\it I'); title('Frequency analysis'); %график сигнала (временная характ)subplot(2,1,2); plot(f(1:k/2),m(1:k/2));grid; xlabel('\omega'); ylabel('\it I'); %график сигнала (частотная хар,спектр)'---Дополнительное задание частот 30, 40, 80 Гц----'k=199 %число точекt=0:1/k:1; %вектор времениx=sin(2*pi*30*t)+sin(2*pi*40*t)+sin(2*pi*80*t); %вектор сигналовy=fft(x); %вектор БПФ сигналовm=abs(y); %вектор амплитудf=(0:length(y)-1)*k/length(y); % вектор частотfiguretitle('Frequency analysis');2subplot(2,1,1); plot(t,x); grid; xlabel('\it t');ylabel('\it I'); title('Frequency analysis'); %график сигнала (временная характ)subplot(2,1,2); plot(f(1:k/2),m(1:k/2));grid; xlabel('\omega'); ylabel('\it I'); %график сигнала (частотная хар,спектр)Результаты работы программы в «Command Window»:---1.1 Расчет резонансной частоты колебательного контура---variant =7L=2.5000C=1.0000e-005f0 =31.8310T=0.0314-----------------------------------------------------------------------1.2 Построение резонансной кривой------------------1.3 Операции с матрицами.-------------Зададим матрицы случайных чиселB=0.9501 0.4860 0.4565 0.44470.2311 0.8913 0.0185 0.61540.6068 0.7621 0.8214 0.7919X=0.9218 0.9355 0.05790.7382 0.9169 0.35290.1763 0.4103 0.81320.4057 0.8936 0.0099A=0.1389 0.2722 0.4451 0.84620.2028 0.1988 0.9318 0.52520.1987 0.0153 0.4660 0.20260.6038 0.7468 0.4186 0.6721Bt =0.9501 0.2311 0.60680.4860 0.8913 0.76210.4565 0.0185 0.82140.4447 0.6154 0.7919Xt =0.9218 0.7382 0.1763 0.40570.9355 0.9169 0.4103 0.89360.0579 0.3529 0.8132 0.0099BtXt =1.1272 1.1275 0.7558 0.59801.3259 1.4449 1.0710 1.00120.4856 0.6438 0.7560 0.20981.0315 1.1720 0.9749 0.7382K=1.2661 1.3996 1.2009 1.44421.5287 1.6437 2.0029 1.52630.6844 0.6591 1.2220 0.41251.6353 1.9188 1.3935 1.4103K=20.3706 -22.2631 19.1799 -2.3754-15.1566 14.5336 -13.0369 3.60473-3.3335 3.8873 -2.0591 -0.19120.2950 2.1997 -2.4676 -1.2520------3.2 Спектральный анализ сигналов.-----k=199---Дополнительное задание частот 30, 40, 80 Гц---k=199Графики:Resonance curve0.10.090.080.07I0.060.050.040.030.020.012829303132333435ωРезонансная криваяFrequency analysisI50-500.10.20.30.40.5t0.60.70.80.910102030405060708090100300250I200150100500ωСпектральный анализ (задание из методики)4Frequency analysis32I10-1-2-300.10.20.30.40.5t0.60.70.80.91010203040506070809010010080I6040200ωСпектральный анализ (задание преподавателя)Синтез дискретных фильтров.0 80040040800899,940899,99003730252015дБ10001100120001199,9500037f3551100 9001199,95Рекурсивный фильтр (IIR).АФЧХ фильтра:Magnitude (dB) and Phase Responses800-50480-150160-250-160-350-480Magnitude (dB)Phase (degrees)50-45000.5Numerator:0.709197-9.19819257.867025-233.622836675.980938-1484.1575232555.044543-3515.7683473906.290690-3515.7683472555.044543-1484.157523675.980938-233.62283657.867025-9.1981920.70919711.522.5Frequency (kHz)33.544.5Denominator:1.000000-12.41335974.746440-288.858865800.139718-1682.0134002772.858409-3654.2121113889.093414-3353.3653172335.080620-1299.842526567.433023-187.98498644.639512-6.8032330.5029610.71- 9.2⋅ z-1 + 57.87⋅ z-2 − 233.62⋅ z-3 + ...+ 0.71⋅ z-16H(z) =1-12.41⋅ z-1 + 74.75⋅ z-2 − 288.86⋅ z-3 + ...+ 0.5 ⋅ z-16-800Нерекурсивный фильтр (FIR)АФЧХ фильтра:Magnitude (dB) and Phase Responses0-8-2400-36-4800-64-7200-92-9600Magnitude (dB)Phase (degrees)20-12000.511.522.5Frequency (kHz)33.544.5Порядок фильтра велик, поэтому краткая запись передаточной характеристики:H ( z ) = 0.003 − 0.06 ⋅ z -1 + 0.006 ⋅ z -2 + 0.003 ⋅ z -3 + ...
+ 0.003 ⋅ z -1223.1. Просмотр демонстрационного примера phone.Графики частот для цифр 1, 5, 9, #:Spectrum5Signal Power10010-510-10100.0200.0697770852 94112091336 1477Frequency (Hz)Клавиша 1, частоты 697 и 1209 Герц1800-12000Spectrum510Signal Power010-510-10100.0200.0697770852 94112091336 1477Frequency (Hz)1800Клавиша 5, частоты 770 и 1336 ГерцSpectrum5Signal Power10010-510-10100.0200.0697770852 94112091336 1477Frequency (Hz)1800Клавиша 9, частоты 852 и 1477 ГерцSpectrum5Signal Power10010-510-10100.0200.0697770852 94112091336 1477Frequency (Hz)1800Клавиша #, частоты 941 и 1477 Герц.Имитационное моделирование с использованием пакета Simulink.Схема модели101.4U1.210.80.60.40.2t01234567Осциллограмма выходного сигналаГрафик выходного сигнала118910.