filter

Измените частотное содержимое timeseries объекты

Описание

tsout = filter(tsin,b,a) применяет рациональную передаточную функцию фильтра b (z−1) / a (z−1) к равномерно разнесенным данным в timeseries tsin объекта. Числитель b и знаменатель a являются векторами, содержащими коэффициенты передаточной функции.

tsout = filter(tsin,b,a,ind) задает индексы столбцов или строк для фильтрации. ind - вектор из целых чисел, представляющих индексы столбцов для ориентированных на столбцы данных (tsin.IsTimeFirst является true) и представление индексов строк для ориентированных на строки данных (tsin.IsTimeFirst является false).

Примеры

свернуть все

Этот пример применяет следующую передаточную функцию к набору данных:

H(z-1)=b(z-1)a(z-1)=2+3z-11+0.2z-1

Создайте timeseries объект из матрицы count в count.dat.

load count.dat
tsin = timeseries(count(:,1),[1:24]);

Введите коэффициенты для знаменателя и числителя передаточной функции. Упорядочьте коэффициенты в возрастающих степенях z-1 представлять 1+0.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')

Figure contains an axes. The axes with title Time Series Plot:unnamed contains 2 objects of type line. These objects represent Original Data, Filtered Data.

Входные параметры

свернуть все

Входные timeseries, заданный как скаляр. tsin должны быть равномерно отобраны.

Типы данных: timeseries

Коэффициенты числителя передаточной функции, заданные в виде скаляра или вектора.

Коэффициенты знаменателя передаточной функции, заданные в виде скаляра или вектора.

Индексы строка или столбец, заданные как положительный целый числовой скаляр или вектор. ind представляет индексы столбцов для ориентированных на столбцы данных (tsin.IsTimeFirst является true) и представляет индексы строк для ориентированных на строки данных (tsin.IsTimeFirst является false).

Типы данных: double | single | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Подробнее о

свернуть все

Рациональная передаточная функция

Область описания операции ввода-вывода filter операция над вектором в области Z-преобразования является рациональной передаточной функцией. Рациональная передаточная функция имеет вид,

Y(z)=b(1)+b(2)z1+...+b(nb+1)znb1+a(2)z1+...+a(na+1)znaX(z),

который обрабатывает как конечные импульсные характеристики, так и БИХ [1]. na - порядок фильтра с обратной связью, и nb - порядок фильтра с feedforward.

Можно также выразить рациональную передаточную функцию как следующее разностное уравнение,

a(1)y(n)=b(1)x(n)+b(2)x(n1)+...+b(nb+1)x(nnb)a(2)y(n1)...a(na+1)y(nna).

Кроме того, вы можете представлять рациональную передаточную функцию с помощью ее прямой реализации в форме II, как на следующей схеме. Из-за нормализации примите a(1) = 1.

Область операции filter на выборке m задается разностными уравнениями во временной области

y(m)=b(1)x(m)+w1(m1)w1(m)=b(2)x(m)+w2(m1)a(2)y(m)       =                 wn2(m)=b(n1)x(m)+wn1(m1)a(n1)y(m)wn1(m)=b(n)x(m)a(n)y(m).

См. также

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте