Фильтрация данных

Фильтрация разностного уравнения

Фильтры являются методами обработки данных, которые могут сгладить высокочастотные колебания данных или удалить периодические тренды определенной частоты из данных. В MATLAB^®, filter функционируйте фильтрует вектор из данных x согласно следующему разностному уравнению, которое описывает фильтр на линии задержки с отводами.

$\begin{array}{l} a (1) y (n) = b (1) x (n) + b (2) x (n - 1) + \dots + b (N_{b}) x (n - N_{b} + 1) \\ - a (2) y (n - 1) - \dots - a (N_{a}) y (n - N_{a} + 1) \end{array}$

В этом уравнении a и b являются векторами из коэффициентов фильтра, _Na является порядком фильтра обратной связи, и _Nb является порядком фильтра прямого распространения. n является индексом текущего элемента x. Выход y (n) является линейной комбинацией текущих и предыдущих элементов x и y.

filter функционируйте использование заданные векторы коэффициентов a и b, чтобы отфильтровать входные данные x. Для получения дополнительной информации о разностных уравнениях, описывающих фильтры, см. [1].

Фильтр скользящего среднего значения данных трафика

Скрипт Open Live Script

filter функция является одним способом реализовать фильтр скользящего среднего значения, который является общими данными, сглаживающими метод.

Следующее разностное уравнение описывает фильтр, который усредняет зависящие от времени данные относительно текущего часа и трех предыдущих часов данных.

$y (n) = \frac{1}{4} x (n) + \frac{1}{4} x (n - 1) + \frac{1}{4} x (n - 2) + \frac{1}{4} x (n - 3)$

Импортируйте данные, которые описывают поток трафика в зависимости от времени, и присвойте первый столбец количеств транспортного средства к векторному x.

load count.dat
x = count(:,1);

Создайте векторы коэффициентов фильтра.

a = 1;
b = [1/4 1/4 1/4 1/4];

Вычислите 4-часовое скользящее среднее значение данных и отобразите на графике и исходные данные и отфильтрованные данные.

y = filter(b,a,x);

t = 1:length(x);
plot(t,x,'--',t,y,'-')
legend('Original Data','Filtered Data')

Изменение амплитуды данных

Скрипт Open Live Script

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

В цифровой обработке сигналов фильтры часто представляются передаточной функцией. Z-преобразование разностного уравнения

$\begin{array}{rcl} a (1) y (n) & = & b (1) x (n) + b (2) x (n - 1) + . . . + b (N_{b}) x (n - N_{b} + 1) \\ - a (2) y (n - 1) - . . . - a (N_{a}) y (n - N_{a} + 1) \end{array}$

следующая передаточная функция.

$Y (z) = H (z^{- 1}) X (z) = \frac{b (1) + b (2) z^{- 1} + . . . + b (N_{b}) z^{- N_{b} + 1}}{a (1) + a (2) z^{- 1} + . . . + a (N_{a}) z^{- N_{a} + 1}} X (z)$

Используйте передаточную функцию

$H (z^{- 1}) = \frac{b (z^{- 1})}{a (z^{- 1})} = \frac{2 + 3 z^{- 1}}{1 + 0.2 z^{- 1}}$

изменить амплитуду данных в count.dat.

Загрузите данные и присвойте первый столбец векторному x.

load count.dat
x = count(:,1);

Создайте векторы коэффициентов фильтра согласно передаточной функции $H (z^{- 1})$ .

a = [1 0.2];
b = [2 3];

Вычислите отфильтрованные данные и отобразите на графике и исходные данные и отфильтрованные данные. Этот фильтр, в основном, изменяет амплитуду исходных данных.

y = filter(b,a,x);

t = 1:length(x);
plot(t,x,'--',t,y,'-')
legend('Original Data','Filtered Data')

Ссылки

[1] Оппенхейм, Алан V, Рональд В. Шафер и Джон Р. Бак. Обработка сигналов дискретного времени. Верхний Сэддл-Ривер, NJ: Prentice Hall, 1999.

Документация