Одномерный цифровой фильтр fi объекты
y = filter(b,1,x)
[y,zf] = filter(b,1,x,zi)
y = filter(b,1,x,zi,dim)
y = filter(b,1,x)x используя фильтр, описанный вектором с фиксированной точкой b. Функция возвращает отфильтрованные данные в выходных данных fi объект y. Исходные данные b и x должно быть fi объекты. filter всегда работает по первому несигнлтонному измерению. Таким образом, фильтр работает вдоль первого размера для векторов столбцов и нетривиальных матриц, и вдоль второго размера для векторов строк.
[ предоставляет доступ к начальным и окончательным условиям задержек, y,zf] = filter(b,1,x,zi) zi, и zf. zi - вектор длины length(или массив с первым размером размера b)-1length( и с остальными размерами, соответствующими размерам b)-1x. zi должно быть fi объект с тем же типом данных, что и y и zf. Если не указано значение для zi, по умолчанию используется массив с фиксированной точкой со значением 0 и соответствующих numerictype и размер.
y = filter(b,1,x,zi,dim)[] для входного аргумента zi.
|
Вектор с фиксированной точкой коэффициентов фильтра. |
|
Вектор с фиксированной точкой, содержащий данные для фильтруемой функции. |
|
Вектор с фиксированной точкой, содержащий начальные условия задержек. Если начальные условия задержек равны нулю, можно указать ноль или, если не известен соответствующий размер и Если не указано значение для |
|
Размер, вдоль которого выполняется операция фильтрации. |
|
Выходной вектор, содержащий отфильтрованные данные с фиксированной точкой. |
|
Выходной вектор с фиксированной точкой, содержащий конечные условия задержек. |
Следующий пример фильтрует высокочастотную синусоиду с фиксированной точкой из сигнала, который содержит как низкочастотную, так и высокочастотную синусоиду с фиксированной точкой.
w1 = .1*pi; w2 = .6*pi; n = 0:999; xd = sin(w1*n) + sin(w2*n); x = sfi(xd,12); b = ufi([.1:.1:1,1-.1:-.1:.1]/4,10); gd = (length(b)-1)/2; y = filter(b,1,x); % Plot results, accommodate for group-delay of filter plot(n(1:end-gd),x(1:end-gd)) hold on plot(n(1:end-gd),y(gd+1:end),'r--') axis([0 50 -2 2]) legend('Unfiltered signal','Filtered signal') xlabel('Sample index (n)') ylabel('Signal value')
Полученный график показывает как нефильтрованные, так и отфильтрованные сигналы.
filter функция поддерживает только фильтры FIR. В общем представлении фильтра b/a знаменателем a фильтра FIR является скаляр 1, который является вторым входом этой функции.
numerictype из b может отличаться от numerictype из x.
Если вы хотите указать начальные условия, но не знаете, что numerictype для использования сначала попробуйте фильтровать данные без начальных условий. Это можно сделать, указав [] для ввода zi. После выполнения операции фильтрации у вас есть numerictype из y и zf (по запросу). Потому что numerictype из zi должно соответствовать y и zf, теперь вы знаете numerictype для использования в исходных условиях.
filter функция использует реализацию прямой транспонированной FIR формы следующего уравнения разности:
+... + bL * xn − N
где L - длина фильтра и N - порядок фильтрации.
Следующая диаграмма показывает структуру транспонированного FIR-фильтра прямой формы, используемую filter функция:
