Лаба 1-8 Лабы и ПЗ [Вариант 1] / 1 / Цель работы

Цель работы:

Условие:

ФНЧ КИХ-фильтр. Частота дискретизации Fs=48000 Гц.

Частота среза Fc 500 Гц.

Вектор воздействия:

0.0031 0.0083 0.016 0.026 0.037 0.05 0.06 0.07 0.08 0.09 0.098 0.09 0.08 0.07

0.06 0.05 0.037 0.026 0.016 0.0083 0.0031

w1T 0.039, w2T 0.131

Ход работы:

1.Определить импульсную характеристику фильтра, подав на вход дельта-функцию.

a = [1];

N = length(b);

delta = [1; zeros(N-1, 1)]; h = filter(b, a, delta); subplot(3,2,1); stem(0:N-1,h, 'fill')

grid on

2. Определить импульсную характеристику фильтра, по отсчетам РУ.

a = delta;

[h, nT] = impz(b, a, N, Fs); subplot(3,2,2)

stem(nT, h, 'fill'); grid on

3.Вычислить реакции ЛДС на входное воздействие на основе разностного уравнения. В качестве входного воздействия использовать сумму двух гармонических колебаний sin(w1*T*n)+sin(w2*T*n).Амплитуда двух колебаний принять равными

1В.

n = 0:20;

x = sin(w1*n) + sin(w2*n); y = filter(b,a,x); subplot(3,2,3) plot(n,x,'-or',n,y,'-ob') grid on

hold on

stem(n,x,'-or', 'fill') stem(n,y, 'fill')

4.Вычислить реакции ЛДС на основе уравнения свертки.В качестве входного воздействия использовать сумму двух гармонических колебаний sin(w1*T*n)+sin(w2*T*n).Амплитуда двух колебаний принять равными 1В.

y = conv(b,x); k = length(y); subplot(3,2,4) hold on plot(n,x,'-or') grid on

stem(n,x,'-or', 'fill') nc = 0:(k-1); plot(nc,y,'-ob') stem(nc,y,'fill')

5.Рассчитать АЧХ и ФЧХ фильтра и нарисовать графики.

N1 = N*10;

h = impz(b,a,N); Ah = fft(h,N);

H = abs(Ah);

ph = angle(Ah); dl = round(N/2); H = H(1:dl);

ph = ph(1:dl); n = 0:dl-1; subplot(3,2,5)

plot(n,H,'-or',n,ph,'-ob') grid on;