Решетка и лестничная структура решетки реализация фильтра
[f,g] = latcfilt(k,x)
[f,g] = latcfilt(k,v,x)
[f,g] = latcfilt(k,1,x)
[f,g,zf] = latcfilt(...,'ic',zi)
[f,g,zf] = latcfilt(...,dim)
При фильтрации данных для представления могут использоваться решетчатые коэффициенты
Конечная импульсная характеристика
Полнополюсные БИХ
Фильтры БИХ Allpass
Общие БИХ
[f,g] = latcfilt(k,x) фильтры x с конечная импульсная характеристика коэффициентами в векторе k. Результат фильтра прямой решетки f и g - результат обратного фильтра. Если , f соответствует минимально-фазовому выходу, и g соответствует максимальному выходу фазы.
Если k и x являются векторами, результатом является вектор (сигнал). Матричные аргументы разрешены в соответствии со следующими правилами:
Если x является матрицей и k является вектором, каждый столбец x обрабатывается через решетчатый фильтр, заданный как k.
Если x является вектором и k является матрицей, каждый столбец k используется для фильтрации x, и возвращается сигнальная матрица.
Если x и k являются обеими матрицами с одинаковым числом столбцов, затем i-м столбцом k используется для фильтрации i-го столбца x. Возвращается сигнальная матрица.
[f,g] = latcfilt(k,v,x) фильтры x с БИХ коэффициентами k и лестничные коэффициенты v. Оба k и v должны быть векторами, в то время как x может быть сигнальной матрицей.
[f,g] = latcfilt(k,1,x) фильтры x с БИХ решеткой, заданной как k, где k и x могут быть векторами или матрицами. f - результат и g полнополюсного решетчатого фильтра - результат фильтра allpass.
[f,g,zf] = latcfilt(...,'ic',zi) принимает метод length- k векторная zi определение начального условия состояний решетки. Выходные zf представляет собой k вектор, задающий окончательное условие состояний решетки.
[f,g,zf] = latcfilt(...,dim) фильтры x вдоль размерности dim. Как задать dim значение, коэффициенты конечной импульсной характеристики решетки k должен быть вектором, и вы должны задать все предыдущие входные параметры по порядку. Используйте пустой вектор [] для любых параметров, которые вы не хотите задавать. zf возвращает окончательные условия в столбцах, независимо от формы x.
Сгенерируйте сигнал с 512 выборками белого Гауссова шума.
x = randn(512,1);
Фильтрация данных с помощью конечная импульсная характеристика фильтра. Задайте коэффициенты отражения так, чтобы решетчатый фильтр был эквивалентен фильтру скользящего среднего 3-го порядка.
[f,g] = latcfilt([1/2 1],x);
Постройте графики максимального и минимального выходов решетчатого фильтра
subplot(2,1,1) plot(f) title('Maximum-Phase Output') subplot(2,1,2) plot(g) title('Minimum-Phase Output')
