Иллюстрированный самоучитель по Maple 6/7

Проектирование цифрового фильтра

Теперь приступим к тестированию фильтра. Зададим входной сигнал в виде зашумленного меандра с частотой 500 Гц и размахом напряжения 2 В:

> 1: = round(fs/2/500):
> for n from 0 by 2*1 to Т do
> for n2 from 0 to 1-1 do
> if n+n2 <= Т then
> x[n+n2]: = evalf(-l+rand()/10^12-0.5);
> fi:
> if n+n2+1 <= Т then
> x[n+n2-H]:-=evalf(l+ranoX)/10^12-0.5);
> fi;
> od:
> od:

Временная зависимость синтезированного входного сигнала представлена на рис. 17.21.

Иллюстрированный самоучитель по Maple 6/7 › Примеры решения научно-технических задач › Проектирование цифрового фильтра
Рис. 17.21. Синтезированный входной сигнал

Вычислим реакцию фильтра на входной сигнал:

> for n from 0 to T do
> y[n]: = sum(h[k]*x[n-k],k=0..N);
> od:

Построим график выходного сигнала:

> р: = [seq([j/fs, x[j]], j=0..T)]:q: = [seq([j/fs, y[j]], j =0..Т)]:
> plot(p,time=0..T/fs/4.labels=[time,volts],title='Входной сигнал",сolor=black);
> plot(q,tine=0..T/fs/4.labels=[tirae,volts], titlе='Выходной сигнал",color=black);

Временная зависимость выходного сигнала показана на рис. 17.22. Нетрудно заметить, что в конце концов выходной сигнал вырождается в пятую гармонику входного сигнала, но этому предшествует довольно заметный переходной процесс. Он связан с узкополосностью данного фильтра.

Иллюстрированный самоучитель по Maple 6/7 › Примеры решения научно-технических задач › Проектирование цифрового фильтра
Рис. 17.22. Временная зависимость выходного сигнала цифрового фильтра

Если Вы заметили ошибку, выделите, пожалуйста, необходимый текст и нажмите CTRL + Enter, чтобы сообщить об этом редактору.