|
Текст подпрограммы и версий 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. )