МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
Національний аерокосмічний університет ім. М.Є. Жуковського
«Харківський авіаційний інститут»
Кафедра «Проектування радіоелектронних систем літальних апаратів»
ЛАБОРАТОРНА РОБОТА
з дисципліни Цифрова обробка сигналів
на тему: Дослідження цифрових фільтрів з кінцевою імпульсною характеристикою
Харків +2014
Завдання
Розрахувати КИХ-фільтр четвертого порядку методом найменших квадратів.
№ В.Від фільтраw1w2w3w4B1B2B3B4B5D1D2D3D4D51ФНЧ0.20.30.30.310.5000100.7070.7
Таблиця 1
Потрібно розрахувати по заданих параметрах коефіцієнти b0, b1, b2, b3, b4 при N=4, K=2.
Переглянути і замалювати характеристики: АЧХ, ФЧХ, ІПХ, ГВЗ, спектр сигналу, результати фільтрації.
Виконання
. Код програми розрахунку коеффіціентовall; clear all; clc;
w1=0.2; w2=0.3; w3=0.3; w4=0.3; b1=1; b2=0.5; b3=0; b4=0; b5=0; d1=1; d2=0; d3=0.7; d4=0.7; d5=0.7;=[b1 b2 b3 b4 b5]; d=[d1 d2 d3 d4 d5]; w=[w1 w2 w3 w4];
%% t03=b1 * d1 * w1; tqt=0; i=2: 4
tqt=b (i) * d (i) * (w (i) - w (i - 1)) + b5 * d5 * (0.5-w4);
end=t03 + tqt; (t03);
%% t13=1; t13=b1 * d1 * (sin (m * 2 * pi * w1)/(2 * m * pi)); tqt=0; i=2: 4
tqt=b (i) * d (i) * (sin (m * 2 * pi * w (i)) - sin (m * 2 * pi * w (i - 1)))/(2 * m * pi) - b5 * d5 * (sin (m * 2 * pi * w4)/(2 * m * pi));
end=t13 + tqt; display (t13);
%% t23=2; t23=b1 * d1 * (sin (m * 2 * pi * w1)/(2 * m * pi)); tqt=0; i=2: 4
tqt=b (i) * d (i) * (sin (m * 2 * pi * w (i)) - sin (m * 2 * pi * w (i - 1)))/(2 * m * pi) - b5 * d5 * (sin (m * 2 * pi * w4)/(2 * m * pi));
end=t23 + tqt; display (t23);
%% t00=d1 * w1; tqt=0; i=2: 4
tqt=d (i) * (w (i) -w (i - 1)) + d5 * (0.5-w4);
end=t00 + tqt; display (t00);
%% t11=1; t11=d1 * (w1/2 + sin (m * 4 * pi * w1)/(8 * m * pi)); tqt=0; i=2: 4
tqt=d (i) * ((w (i) -w (i - 1))/2 + (sin (m * 4 * pi * w (i)) - sin (m * 4 * pi * w (i - 1)))/(8 * m * pi)) + d5 * ((0.5-w4)/2- (sin (m * 4 * pi * w4))/(8 * m * pi)) ;
end=t11 + tqt; display (t11);
%% t22=2; t22=d1 * (w1/2 + sin (m * 4 * pi * w1)/(8 * m * pi)); tqt=0; i=2: 4
tqt=d (i) * ((w (i) -w (i - 1))/2 + (sin (m * 4 * pi * w (i)) - sin (m * 4 * pi * w (i - 1)))/(8 * m * pi)) + d5 * ((0.5-w4)/2- (sin (m * 4 * pi * w4))/(8 * m * pi)) ;
end=t22 + tqt; display (t22);
%% t01=0; l=1; t01=d1 * ((sin (2 * pi * (ml) * w1)/(4 * (ml) * pi) + sin (2 * pi * (m + l) * w1)/(4 * (m + l) * pi))); tqt=0;
for i=2: 4=d (i) * ((sin (2 * pi * (ml) * w (i)) - sin (2 * pi * (ml) * w (i - 1)) )/(4 * (ml) * pi) + (sin (2 * pi * (m + l) * w (i)) - sin (2 * pi * (m + l) * w (i - 1)) )/(4 * (m + l) * pi)) - d5 * (((sin (2 * pi * (ml) * w4))/(4 * (ml) * pi)) + ((sin (2 * pi * (m + l) * w4))/(4 * (m + l) * pi)));
end=t01 + tqt; display (t01);
%% t02=0; l=2; t02=d1 * ((sin (2 * pi * (ml) * w1)/(4 * (ml) * pi) + sin (2 * pi * (m + l) * w1)/(4 * (m + l) * pi))); tqt=0;
for i=2: 4=d (i) * ((sin (2 * pi * (ml) * w (i)) - sin (2 * pi * (ml) * w (i - 1)) )/(4 * (ml) * pi) + (sin (2 * pi * (m + l) * w (i)) - sin (2 * pi * (m + l) * w (i - 1)) )/(4 * (m + l) * pi)) - d5 * (((sin (2 * pi * (ml) * w4))/(4 * (ml) * pi)) + ((sin (2 * pi * (m + l) * w4))/(4 * (m + l) * pi)));
end=t02 + tqt; display (t02);
%% t10=1; l=0; t10=d1 * ((sin (2 * pi * (ml) * w1)/(4 * (ml) * pi) + sin (2 * pi * (m + l) * w1)/(4 * (m + l) * pi))); tqt=0;
for i=2: 4=d (i) * ((sin (2 * pi * (ml) * w (i)) - sin (2 * pi * (ml) * w (i - 1)) )/(4 * (ml) * pi) + (sin (2 * pi * (m + l) * w (i)) - sin (2 * pi * (m + l) * w (i - 1)) )/(4 * (m + l) * pi)) - d5 * (((sin (2 * pi * (ml) * w4))/(4 * (ml) * pi)) + ((sin (2 * pi * (m + l) * w4))/(4 * (m + l) * pi)));
end=t10 + tqt; display (t10);
%% t20=2; l=0; t20=d1 * ((sin (2 * pi * (ml) * w1)/(4 * (ml) * pi) + sin (2 * pi * (m + l) * w1)/(4 * (m + l) * pi))); tqt=0;
for i=2: 4=d (i) * ((sin (2 * pi * (ml) * w (i)) - sin (2 * pi * (ml) * w (i - 1)) )/(4 * (ml) * pi) + (sin (2 * pi * (m + l) * w (i)) - sin (2 * pi * (m + l) * w (i - 1)) )/(4 * (m + l) * pi)) - d5 * (((sin (2 * pi * (ml) * w4))/(4 * (ml) * pi)) + ((sin (2 * pi * (m + l) * w4))/(4 * (m + l) * pi)));
end=t20 + tqt; display (t20);
%% t12=1; l=2; t12=d1 * ((sin (2 * pi * (ml) * w1)/(4 * (ml) * pi)) + (s...
