Интерполяция табличной функции. Вычисление коэф-тов.

[ФАНТОМНЫЙ ЮЗЕР]

Подскажите, пожалуйста, где можно найти алгоритм или прогу на VB6 или FORTRAN интероляции таблично заданной комплексной функции в виде отношения двух полиномов числителя и знаменателя:

Код: Выделить всё

    b0 + b1 * z^-1 + b2 * z^-2 + ... + bm * z^-m
W = --------------------------------------------
    a0 + a1 * z^-1 + a2 * z^-2 + ... + an * z^-n


где m<=n - степень полиномов; а, b - постоянные вещественные или комплексные коэффициенты; z - корни.

Аналог функции на MATLAB - функция invfreqz:
[b, a] = invfreqz(h, w, na, nb, wt, iter, tol, 'trace')
где w - вектор частот; h - вектор значений функции; na, nb - степень интерполирумого полинома знаменателя и числителя.

Очень нужен этот алгоритм!!!
Спасибо.

[ФАНТОМНЫЙ ЮЗЕР]

Мне подсказали посмотреть помощь по MATLAB по данной функции, сказав, что...

"данная функция предназначена для решения задачи идентификации систем, она позволяет определить коэффициенты числителя и знаменателя функции передачи дискретной системы по набору значений этой функции передачи на различных частотах.
Решается оптимизационная задача, которая в одном варианте вызова функции приводит к системе линейных уравнений по алгоритму [1], а в другом решается численно ньютоновским методом [2]. И для того и для другого можно найти готовые библиотеки под разные языки программирования. Остается только написать целевую функцию."

[1] Levi, E.C., "Complex-Curve Fitting," IRE Trans. on Automatic Control, Vol. AC-4 (1959), pp. 37-44.
[2] Dennis, J.E., Jr., and R.B. Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, Prentice-Hall, 1983.

Подскажите, пожалуйста, где можной найти эти алгоритмы

[Faust]

http://www.allmath.ru тебе в помощь!
А ещё посмотри на http://alglib.manual.ru/interpolation/.