Идентифицируйте параметры фильтра дискретного времени из данных о частотной характеристике
По умолчанию, invfreqz
использует ошибочный метод уравнения, чтобы идентифицировать лучшую модель из данных. Это находит b
и a
\in
путем создания системы линейных уравнений и решения их с MATLAB®
\
оператор. Здесь A (ω (k)) и B (ω (k)) является преобразованиями Фурье полиномов a
и b
, соответственно, на частоте ω (k) и n количество точек частоты (длина h
и w
). Этот алгоритм на основе Леви [1].
Начальник (“ошибка на выходе”) алгоритм использует ослабленный метод Ньютона Гаусса в итеративном поиске [2] с выходом первого алгоритма как первоначальная оценка. Это решает прямую задачу минимизации взвешенной суммы квадратичной невязки между фактическим и желаемыми точками частотной характеристики.
[1] Леви, E. C. “Комплексный Curve Fitting”. Транзакции IRE на Автоматическом управлении. Издание AC-4, 1959, стр 37–44.
[2] Деннис, J. E. младший, и Р. Б. Шнабель. Численные методы для оптимизации без ограничений и нелинейных уравнений. Englewood Cliffs, NJ: Prentice Hall, 1983.