* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
1162
Аналитическое и спектральное моделирование
Построим график АЧХ фильтра:
> p := [seq([j*fs/(T+1),abs(ro[j+1]+io[j+1]*I)],j=0..T/2)]: > plot(p, frequency=0..fs/2, labels=[frequency,gain], title=`A?O oeeuo?a`,color=black);
Он представлен на рис. 15.34. Нетрудно заметить, что и впрямь АЧХ фильтра напоминает АЧХ резонансной цепи – она имеет вид узкого пика. Вы можете легко проверить, что раздвижением частот fl и fh можно получить АЧХ с довольно плоской вершиной и резкими спадами (говорят, что такая характеристика при ближается к прямоугольной).
Рис. 15.34. АЧХ цифрового полосового фильтра
Теперь приступим к тестированию фильтра. Зададим входной сигнал в виде зашумленного меандра с частотой 500 Гц и размахом напряжения 2 В:
> l := round(fs/2/500): > for n from 0 by 2*l to T do > for n2 from 0 to l-1 do > if n+n2 <= T then > x[n+n2] := evalf(-1+rand()/10^12-0.5); > fi; > if n+n2+l <= T then > x[n+n2+l] := evalf(1+rand()/10^12-0.5); > fi; > od: > od:
Временная зависимость синтезированного входного сигнала представлена на рис. 15.35. Вычислим реакцию фильтра на входной сигнал:
> for n from 0 to T do > y[n] := sum(h[k]*x[n-k],k=0..N); > od:
Построим график выходного сигнала:
> p := [seq([j/fs,x[j]],j=0..T)]:q := [seq([j/fs,y[j]],j=0..T)]: > plot(p,time=0..T/fs/4,labels=[time,volts],title=`Aoiaiie neaiae`,color=black);