Текст подпрограммы и версий ( Фортран ) zp47r.zip zp48r.zip zp49r.zip zp50r.zip zp51r.zip zp52r.zip zp47d.zip zp48d.zip zp49d.zip zp50d.zip zp51d.zip zp52d.zip |
Тексты тестовых примеров ( Фортран ) tzp47r.zip tzp48r.zip tzp49r.zip tzp50r.zip tzp51r.zip tzp52r.zip tzp47d.zip tzp48d.zip tzp49d.zip tzp50d.zip tzp51d.zip tzp52d.zip |
Текст подпрограммы и версий ( Си ) zp47r_c.zip zp48r_c.zip zp49r_c.zip zp50r_c.zip zp51r_c.zip zp52r_c.zip |
Тексты тестовых примеров ( Си ) tzp47r_c.zip tzp48r_c.zip tzp49r_c.zip tzp50r_c.zip tzp51r_c.zip tzp52r_c.zip |
Текст подпрограммы и версий ( Паскаль ) zp47r_p.zip zp48r_p.zip zp49r_p.zip zp50r_p.zip zp51r_p.zip zp52r_p.zip zp47e_p.zip zp48e_p.zip zp49e_p.zip zp50e_p.zip zp51e_p.zip zp52e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp47r_p.zip tzp48r_p.zip tzp49r_p.zip tzp50r_p.zip tzp51r_p.zip tzp52r_p.zip tzp47e_p.zip tzp48e_p.zip tzp49e_p.zip tzp50e_p.zip tzp51e_p.zip tzp52e_p.zip |
Вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева.
Для заданных коэффициентов С (К) вычисляется частичная сумма ряда
N -1 SUM = ∑ C(K+1) * T(K, X) K =0
где Т (К, Х) = cos (К * АRССОS (Х)) - значение ортогонального полинома Чебышева порядка К в точке Х.
М.Аbramowitz, I.А.Stegun, Нandbook of Мathematical Functions, New York, Dover Рublications, Inc., 1965.
SUBROUTINE ZP47R (X, C, N, SUM)
Параметры
X - | заданное значение аргумента (тип: вещественный); |
C - | вещественный вектор длины N, содержащий заданные коэффициенты разложения; при этом полиному порядка К соответствует коэффициент С (К + 1); |
N - | длина вектора С (тип: целый); |
SUM - | вычисленное значение частичной суммы (тип: вещественный). |
Версии
ZP47D - | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева для данных двойной точности. |
ZP48R - ZP48D | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева для данных одинарной и двойной точности. |
ZP49R - ZP49D | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева четного порядка для данных одинарной и двойной точности. |
ZP50R - ZP50D | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева четного порядка для данных одинарной и двойной точности. |
ZP51R - ZP51D | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева нечетного порядка для данных одинарной и двойной точности. |
ZP52R - ZP52D | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева нечетного порядка для данных одинарной и двойной точности. |
Вызываемые подпрограммы: нет
Замечания по использованию
Смещенные ортогональные многочлены Чебышева ТS (К, Х) определяются по формуле ТS (К, Х) = Т (К, 2 * Х - 1), где К - порядок полинома, Х - аргумент и Т (К, Х) - значение обычного полинома Чебышева порядка К в точке Х: Т (К, Х) = cos (К * АRССОS (Х)). |
DIMENSION C(3) DATA C /1., 2., 3./ X = 0.5 N = 3 CALL ZP47R (X, C, N, SUM) Результат: SUМ = 0.5
При тех же самых значениях входных параметров как и в примере выше, при обращении к другим версиям, будут выданы следующие результаты.
ZP48R SUМ = -2.0 ; ZP49R SUМ = -1.5 ; ZP50R SUМ = 2.0 ; ZP51R SUМ = 0.0 ; ZP52R SUМ = 0.0 .