Цифровая фильтрация (стр. 2 из 2)
(1.4)В формуле (1.4) было учтено, что в интервал интегрирования (-T/2,T/2) попадает только одна дельта-функция, соответствующая k=0.
Таким образом, периодическая последовательность дельта-функций может быть представлена в виде комплексного ряда Фурье:
(1.5)где
Умножение сигнала на
соответствует сдвигу спектральной функции на 
поэтому спектр дискретизированного сигнала можно записать следующим образом: 
(1.6)Таким образом спектр дискретизированного сигнала представляет собой бесконечный ряд сдвинутых копий спектра исходного непрерывного сигнала s(t). Расстояние по частоте между соседними копиями спектра равно частоте дискретизации
Следует также отметить, что из-за наличия в формуле (1.6) множителя 1/T спектр дискретизированного сигнала имеет размерность, совпадающую с размерностью сигнала.
Спектр дискретизированного сигнала.
Вывод: Таким образом спектр дискретизированного сигнала представляет собой бесконечный ряд сдвинутых копий спектра исходного непрерывного сигнала s(t).
Дискретное преобразование Фурье.
Пусть последовательность отсчетов {x(k)} является периодической с периодом N:
x(k+N)=x(k) для любого k.
Такая последовательность полностью описывается конечным набором чисел, в качестве которого можно взять произвольный фрагмент длиной N, например {x(k), k=0, 1, …, N-1}. Поставленный в соответствии этой последовательности сигнал из смещенных по времени дельта-функций:
(1.1)Также, разумеется, будет периодическим с минимальным периодом NT.
Так как сигнал (1.1) является дискретным, его спектр должен быть периодическим с периодом
. Так как этот сигнал является также и периодическим, его спектр должен быть дискретным с расстоянием между гармониками, равным 
.Итак, периодический дискретный сигнал имеет периодический дискретный спектр, который также описывается конечным набором из N чисел (один период спектра содержит N гармоник).
Рассмотрим процедуру вычисления спектра периодического дискретного сигнала. Так как сигнал периодический, будем раскладывать его в ряд Фурье. Коэффициенты
этого ряда равны 
(1.2)В выражении (1.2) реальный масштаб времени фигурирует только в множителе 1/T перед оператором суммирования. При рассмотрении дискретных последовательностей обычно оперируют номерами отсчетов и спектральных гармоник без привязки к действительному масштабу времени и частоты. Поэтому множитель 1/T (1.2) удаляют, то есть считают частоту дискретизации равной единице. Удаляют обычно и множитель 1/N. Получившееся выражение называется дискретным преобразованием Фурье.
С помощью программы, написанной на языке MatLab строим график модуля ДПФ.
x1=[10; 10; 10; 10; 10; 10; 10; 10; 10; 10; 10; zeros(8,1)];
y1=fft(x1);
x2=[x1; zeros(18,1)];
y2=fft(x2);
subplot(2, 2, 1),stem(0:18, x1)
xlim([0 37])
subplot(2, 2, 2),stem((0:18)/19,abs(y1))
subplot(2, 2, 3),stem(0:36,x2)
xlim([0 37])
subplot(2, 2, 4)
stem((0:36)/38, abs(y2))
Рис.3 Повышение спектрального разрешения ДПФ при дополнении сигнала нулями: сверху - исходный сигнал и модуль его ДПФ, снизу - сигнал, дополненный 18 нулями, и модуль его ДПФ.
Анализ фильтра.
Графическое и аналитическое представление фильтра.
Расчет и преобразование аналоговых фильтров-прототипов.
Одной из часто возникающих задач является создание фильтров, пропускающих сигналы в определенной полосе частот и задерживающих остальные частоты.
Фильтры высокой частоты (ФВЧ, highpass filter), пропускающие частоты в диапазоне от wн до бесконечности.
Идеальная форма АЧХ не может быть физически реализована. Поэтому в теории аналоговых фильтров разработан ряд методов аппроксимации прямоугольных АЧХ. Расчет аналогового фильтра начинается с расчета так называемого фильтра-прототипа, представляющего ФНЧ.
Через программу SPTool рассчитываем аналоговый фильтр-прототип с помощью фильтра Чебышева первого рода.
Функция передачи фильтра-прототипа Чебышева не имеет нулей, а ее полюсы равномерно расположены на S-плоскости в левой половине эллипса.
Формула АЧХ фильтра Чебышева:
,Где
- частота среза, n-порядок фильтра, 
- полином Чебышева n-порядка, 
- параметр, определяющий величину АЧХ в полосе пропускания.Входной сигнал.
Сигнал = [10(*50раз), 0(*270раз)]
Характеристики фильтра.
Расчет ФВЧ фильтра Чебышева второго рода 2-го порядка с помощью функций MatLab.
[z,p,k]=cheb1ap(5,3);
[b, a]= zp2tf(z,p,k);
w0=2*pi*1e3;
[b,a]=lp2hp(b, a, w0);
f=0:1:20e3;
h=freqs(b,a,2*pi*f);
plot(f/1000,abs(h)),grid
axis tight
figure
plot(f/1000, unwrap(angle(h))), grid
АЧХ фильтра.
Параметры фильтра
Тип: Chebyshev I IIR HighPass(ФВЧ)
Sampling Frequency(Частота Дискрктизации): 1000
Passband:
Fp: 80
Rp: 3
Stopband:
Fs: 50
Rs: 20
ФЧХ фильтра.
Вывод: графики АЧХ и ФЧХ фильтра, построенные через SPTool и MatLab практически сходны.
Фильтрация сигнала.
Для пропускания сигнала через фильтр необходимо выбрать сигнал и фильтр, полученные в предыдущем пункте.
Выходной сигнал.
Входной и выходной сигналы.
Нахождение выходного сигнала с помощью функций MatLab:
s=[10(*50раз), 0(*270раз)]
[b,a]=cheby1(5, 3, 0.75);
s1=filter(b,a,s);figure
plot(s1)
Вывод: графики, построенные с помощью функций MatLab и SPTool практически одинаковы.
