Моделирование работы рекурсивного цифрового фильтра

Федеральное агентство по образованию

Государственное образовательное учреждение

высшего профессионального образования

«Братский государственный университет»

                                                                                                                   Кафедра УТС

КСММ

Лабораторная работа №2

«Моделирование работы рекурсивного цифрового фильтра»

Выполнил:

 Ст-т гр. АТП-03-2                                                                                  Ельников С.С

Проверил:                                                                                               Дьяконица С.А.

Братск 2007

Задание:

1.  Определение структуры спроектированного рекурсивного цифрового фильтра.

2.  Моделирование работы рекурсивного цифрового фильтра

f_n=650

f_з=500

F_max=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))