Лаба №2, Вариант 7 (1084965)
Текст из файла
Вариант 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.
Характеристики
Тип файла PDF
PDF-формат наиболее широко используется для просмотра любого типа файлов на любом устройстве. В него можно сохранить документ, таблицы, презентацию, текст, чертежи, вычисления, графики и всё остальное, что можно показать на экране любого устройства. Именно его лучше всего использовать для печати.
Например, если Вам нужно распечатать чертёж из автокада, Вы сохраните чертёж на флешку, но будет ли автокад в пункте печати? А если будет, то нужная версия с нужными библиотеками? Именно для этого и нужен формат PDF - в нём точно будет показано верно вне зависимости от того, в какой программе создали PDF-файл и есть ли нужная программа для его просмотра.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
