Изменение частотного содержания timeseries объекты
tsout = filter( применяет фильтр рациональной передаточной функции b (z − 1 )/a (z − 1) к равномерно разнесенным данным вtsin,b,a)timeseries объект tsin. Нумератор b и знаменатель a - векторы, содержащие коэффициенты передаточной функции.
tsout = filter( указывает индексы столбцов или строк для фильтрации. tsin,b,a,ind)ind - вектор целых чисел, представляющий индексы столбцов для данных, ориентированных на столбцы (tsin.IsTimeFirst является true) и представляющих индексы строк для данных, ориентированных на строки (tsin.IsTimeFirst является false).
В этом примере к набору данных применяется следующая передаточная функция:
3z-11 + 0 .2z-1
Создать timeseries объект из матрицы count в count.dat.
load count.dat
tsin = timeseries(count(:,1),[1:24]);Введите коэффициенты для знаменателя и числителя передаточной функции. Упорядочить коэффициенты в возрастающих степенях для представления .2x 2-3z-1.
a = [1 0.2]; b = [2 3];
Применение функции переноса с помощью filterи сравнить исходные данные с отфильтрованными данными.
tsout = filter(tsin,b,a); plot(tsin) hold on plot(tsout) legend('Original Data','Filtered Data','Location','NorthWest')
Описание ввода-вывода filter операция над вектором в области Z-преобразования является рациональной передаточной функцией. Рациональная передаточная функция имеет вид,
− 1 +... + a (na + 1) z − naX (z),
который обрабатывает фильтры FIR и IIR [1]. na - порядок фильтра обратной связи, а nb - порядок фильтра прямой связи.
Можно также выразить рациональную передаточную функцию как следующее уравнение разности,
2) y (n − 1) −... − a (na + 1) y (n − na).
Кроме того, можно представить рациональную передаточную функцию, используя ее прямую форму II транспонированной реализации, как на следующей диаграмме. Из-за нормализации, предположим a(1) = 1.
Работа filter при выборке m задается уравнениями разностей во временной области
(m) wn − 1 (m) =