Лаба 2 Вариант 7

Вариант 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) f0

T '-----------------------------------------------------------------' '-------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*Xt

K=(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); % вектор частот figure

title('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); % вектор частот figure

title('Frequency analysis');

2

Рекурсивный фильтр (IIR).

Нерекурсивный фильтр (FIR)

Клавиша #, частоты 941 и 1477 Герц.

Имитационное моделирование с использованием пакета Simulink.

Схема модели

10

1.4 U

1.2

1

0.8

0.6

0.4

0.2

t

Осциллограмма выходного сигнала

График выходного сигнала

11