Текст подпрограммы и версий zp41r_p.zip zp41e_p.zip |
Тексты тестовых примеров tzp41r_p.zip tzp41e_p.zip |
Восстановление коэффициентов многочлена, заданного разложением по многочленам Чебышева.
Пусть известны коэффициенты А = ( a1, a2, ..., aN ) разложения многочлена Р (Х) по многочленам Чебышева { Тk (Х) } :
N P(X) = ∑ ak Tk -1(X) k =1
Подпрограмма вычисляет коэффициенты В = ( b1, b2, ..., bN ) разложения этого многочлена по степеням Х:
N P(X) = ∑ bk X k -1 k =1
procedure ZP41R(var A :Array of Real; var B :Array of Real; N :Integer);
Параметры
A - | вектор заданных коэффициентов разложения многочлена по многочленам Чебышева (тип: вещественный); |
B - | вектор искомых коэффициентов разложения многочлена по степеням аргумента (тип: вещественный); |
N - | размерность векторов А и В (тип: целый). |
Версии
ZP41E - | восстановление коэффициентов многочлена, заданного разложением по многочленам Чебышева, для данных расширенной (Extended) точности. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
Unit TZP41R_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc, UtRes_p, ZP41R_p; function TZP41R: String; implementation function TZP41R: String; var N,_i :Integer; A :Array [0..2] of Real; B :Array [0..2] of Real; begin Result := ''; { результат функции } { прототип оператора DАТА на FORTRANе } A[0] := 1.0; A[1] := 1.0; A[2] := 1.0; N := 3; ZP41R(A,B,N); Result := Result + #$0D#$0A; for _i:=0 to 2 do begin Result := Result + Format('%20.16f ',[B[_i]]); if ( ((_i+1) mod 3)=0 ) then Result := Result + #$0D#$0A; end; Result := Result + #$0D#$0A; UtRes('TZP41R',Result); { вывод результатов в файл TZP41R.res } exit; end; end. Результат: В = ( 0., 1., 2. )