prony

Метод Прони для проектирования фильтров.

Библиотека

EngeeDSP

Синтаксис

Вызов функции

b,a = prony(h,bord,aord) — возвращает коэффициенты числителя и знаменателя для причинной рациональной передаточной функции с импульсной характеристикой h, порядком числителя bord и порядком знаменателя aord.

Аргументы

Входные аргументы

# h — импульсная характеристика
вектор

Details

Импульсная характеристика, заданная как вектор.

Типы данных	`Float32`, `Float64`
Поддержка комплексных чисел	Да

# bord,aord — порядки числителя и знаменателя
положительный целочисленный скаляр

Details

Порядки числителя и знаменателя, заданные как целые положительные скаляры. Если длина h меньше max(bord,aord), то функция заполняет импульсную характеристику нулями.

Если вам нужна передаточная функция только с полюсами, укажите для bord значение 0.
Если вам нужна передаточная функция только с нулями, укажите для aord значение 0.

Типы данных

Float32, Float64

Выходные аргументы

# b,a — коэффициенты передаточной функции
вектор

Details

Коэффициенты передаточной функции, возвращаемые в виде векторов. Вектор b имеет длину bord+1 и a имеет длину aord+1.

Дополнительно

Передаточная функция

Details

Передаточная функция — это Z-преобразование импульсной характеристики :

Рациональная передаточная функция — это отношение многочленов по отрицательным степеням . Это уравнение описывает причинную рациональную передаточную функцию с порядком числителя и порядком знаменателя :

где .

Литература

Parks, Thomas W., and C. Sidney Burrus. Digital Filter Design. New York, NY, USA: Wiley-Interscience, 1987.