Реализация решетчатого и решетчато-лестничного фильтра
[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) принимает длину -k вектор zi задание начального условия состояний решетки. Продукция zf является длинно-k вектор, задающий конечное условие состояний решетки.
[f,g,zf] = latcfilt(...,dim) фильтры x вдоль размера dim. Чтобы указать dim значение, коэффициенты решётки КИХ k должен быть вектором и по порядку указывать все предыдущие входные параметры. Используйте пустой вектор [] для любых параметров, которые не требуется указывать. zf возвращает окончательные условия в столбцах, независимо от формы x.
Создайте сигнал с 512 выборками белого гауссова шума.
x = randn(512,1);
Фильтрация данных с помощью фильтра решетки FIR. Задайте коэффициенты отражения так, чтобы решетчатый фильтр был эквивалентен фильтру скользящего среднего 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')
