Федеральное агентство по образованию
Государственное образовательное учреждение
высшего профессионального образования
«Братский государственный университет»
Кафедра УТС
КСММ
Лабораторная работа №2
«Моделирование работы рекурсивного цифрового фильтра»
Выполнил:
Ст-т гр. АТП-03-2 Ельников С.С
Проверил: Дьяконица С.А.
Братск 2007
Задание:
1. Определение структуры спроектированного рекурсивного цифрового фильтра.
2. Моделирование работы рекурсивного цифрового фильтра
fn=650
fз=500
Fmax=1000
A=1
N=400
1.Необходимо представить рекурсивный цифровой фильтр в прямой, канонической, каскадной и параллельной форме.
2. Необходимо получить рекурсивный цифровой фильтр на следующие тестовые реакции.
1)
а) прямая форма:
б) каноническая форма:
в) каскадная форма:
[s,q]=tf2sos(a,b)
s = 1.0000 -0.8250 1.0003 1.0000 -0.6438 0.7719
1.0000 -0.8249 0.9997 1.0000 -0.8182 0.7814
1.0000 -0.8246 0.9999 1.0000 -0.5533 0.8970
q =0.7007
в) параллельная форма:
[r,p,k]=residuez(a,b)
r = -0.0499 - 0.0514i
-0.0499 + 0.0514i
0.2545 + 0.0810i
0.2545 - 0.0810i
-0.6607 + 0.4273i
-0.6607 - 0.4273i
p = -0.0671 + 0.9979i
-0.0671 - 0.9979i
-0.1834 + 0.7305i
-0.1834 - 0.7305i
-0.2505 + 0.2674i
-0.2505 - 0.2674i
k = 0.9120
[a1,b1]=residuez(r(1:2),p(1:2),[])
a1 = -0.0997 - 0.0000i 0.0958 - 0.0000i
b1 =1.0000 0.1343 1.0003
[a1,b1]=residuez(r(3:4),p(3:4),[])
a1 =0.5091 - 0.0000i -0.0250 - 0.0000i
b1 =1.0000 0.3668 0.5673
[a1,b1]=residuez(r(5:6),p(5:6),[])
a1 =-1.3214 -0.5596
b1 = 1.0000 0.5011 0.1343
2)
t=0:0.0005:0.16
а) реакция на прямоугольный импульс
s=rectpuls(t-0.05,0.1)
plot(t,s)
Рис.1 прямоугольный импульс
y=filter(az,bz,s)
plot(t,y)
Рис.2 выходной сигнал соответствующий разностному уравнению
h=impz(az,bz,401)
y=conv(s,h)
plot(t,y(1:401))
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.