Измените частотное содержимое 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).
Этот пример применяет следующую передаточную функцию к набору данных:
Создайте timeseries объект из матрицы count в count.dat.
load count.dat
tsin = timeseries(count(:,1),[1:24]);Введите коэффициенты для знаменателя и числителя передаточной функции. Упорядочьте коэффициенты в возрастающих степенях представлять и .
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]. na - порядок фильтра с обратной связью, и nb - порядок фильтра с feedforward.
Можно также выразить рациональную передаточную функцию как следующее разностное уравнение,
Кроме того, вы можете представлять рациональную передаточную функцию с помощью ее прямой реализации в форме II, как на следующей схеме. Из-за нормализации примите a(1) = 1.
Область операции filter на выборке m задается разностными уравнениями во временной области