|
Текст подпрограммы и версий ( Фортран ) zp44r.zip zp44d.zip |
Тексты тестовых примеров ( Фортран ) tzp44r.zip tzp44d.zip |
|
Текст подпрограммы и версий ( Си ) zp44r_c.zip |
Тексты тестовых примеров ( Си ) tzp44r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp44r_p.zip zp44e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp44r_p.zip tzp44e_p.zip |
Восстановление коэффициентов многочлена, заданного разложением по системе ортогональных многочленов.
Пусть известны коэффициенты C = ( c1, c2, ..., cN ) разложения многочлена Р (Х) по системе { Qk (X) } ортогональных многочленов
N
P(X) = ∑ Ck Qk -1(X) ,
k =1
удовлетворяющих трехчленному рекуррентному соотношению
Q0(X) = 1 , Q1(X) = X - a1 ,
Qk +1(X) = ( X - ak +1 ) Qk(X) - bk Qk -1(X) , k = 1, 2, ..., N-2
с известными коэффициентами А = ( a1, a2, ..., aN - 1 ), B = ( b1, b2, ..., bN - 2 ). Подпрограмма вычисляет коэффициенты D = ( d1, d2, ..., dN ) разложения этого многочлена по степеням Х
N
P(X) = ∑ dk X k -1
k =1
SUBROUTINE ZP44R (A, B, C, D, N)
Параметры
| A - | вещественный вектор длины N - 1, содержащий коэффициенты ak рекуррентного соотношения (тип: вещественный); |
| B - | вещественный вектор длины N - 2, содержащий коэффициенты bk рекуррентного соотношения (тип: вещественный); |
| C - | вектор заданных коэффициентов разложения многочлена по системе ортогональных многочленов (тип: вещественный); |
| D - | вектор искомых коэффициентов разложения многочлена по степеням аргумента (тип: вещественный); |
| N - | размерность векторов С и D (тип: целый). |
Версии
| ZP44D - | вычисление значения многочлена, заданного разложением по ортогональным многочленам, для данных двойной точности. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(2), B(1), C(3), D(3)
DATA A(1), A(2), B(1), C(1), C(2), C(3) /0., 5*1./
N = 3
CALL ZP44R (A, B, C, D, N)
Результат: D = ( 0., 0., 1. )