Текст подпрограммы и версий ( Фортран ) zp45r.zip zp45d.zip zp46r.zip zp46d.zip |
Тексты тестовых примеров ( Фортран ) tzp45r.zip tzp45d.zip tzp46r.zip tzp46d.zip |
Текст подпрограммы и версий ( Си ) zp45r_c.zip zp46r_c.zip |
Тексты тестовых примеров ( Си ) tzp45r_c.zip tzp46r_c.zip |
Текст подпрограммы и версий ( Паскаль ) zp45r_p.zip zp45e_p.zip zp46r_p.zip zp46e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp45r_p.zip tzp45e_p.zip tzp46r_p.zip tzp46e_p.zip |
Вычисление значений ортогональных полиномов Чебышева.
Подпрограмма вычисляет значения ортогональных полиномов Чебышева для заданного аргумента Х по рекуррентной формуле:
T(K) = 2*X*T(K-1) - T(K-2) , T(1) = 1 , T(2) = X ,
где Т (К) - значение полинома порядка К в точке Х.
SUBROUTINE ZP45R (X, N, T)
Параметры
X - | заданное значение аргумента (тип: вещественный); |
N - | длина вектора Т (тип: целый); |
T - | вещественный вектор длины N, в котором содержатся значения полиномов в порядке возрастания степеней, начиная с полинома степени ноль. |
Версии
ZP45D - | вычисление значений ортогональных полиномов Чебышева для данных двойной точности. |
ZP46R - ZP46D | вычисление значений смещенных ортогональных полиномов Чебышева. для данных одинарной и двойной точности. |
Вызываемые подпрограммы: нет
Замечания по использованию
Смещенные ортогональные многочлены Чебышева ТS (К, Х) определяются по формуле ТS (К, X) = Т (К, 2 * Х - 1), где К - порядок полинома, Х - аргумент и Т (К, Х) - значение обычного полинома Чебышева порядка К в точке Х: Т (К, Х) = cos (К * АRССОS (Х)). |
DIMENSION T(3) X = 0.5 N = 3 CALL ZP45R (X, N, T) Результаты: T = ( 1., 0.5, - 0.5 )