Текст подпрограммы и версий zp42r_p.zip zp42e_p.zip |
Тексты тестовых примеров tzp42r_p.zip tzp42e_p.zip |
Отыскание коэффициентов разложения заданного многочлена по многочленам Чебышева.
Пусть известны коэффициенты А = ( a1, a2, ..., aN ) разложения многочлена Р (Х) по степеням Х
N P(X) = ∑ ak X k -1 k =1
Подпрограмма вычисляет коэффициенты В = ( b1, b2, ..., bN ) разложения этого многочлена по многочленам Чебышева { Тk (Х) } :
N P(X) = ∑ bk Tk -1(X) k =1
procedure ZP42R(var A :Array of Real; var B :Array of Real; var C :Array of Real; N :Integer);
Параметры
A - | вектор заданных коэффициентов разложения многочлена по степеням аргумента (тип: вещественный); |
B - | вектор искомых коэффициентов разложения многочлена по многочленам Чебышева (тип: вещественный); |
C - | рабочий вектор размерности N на (N + 1)/2 (тип: вещественный); |
N - | размерность векторов А и В (тип: целый). |
Версии
ZP42E - | отыскание коэффициентов разложения заданного многочлена по многочленам Чебышева для данных расширенной (Extended) точности. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
Unit TZP42R_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc, UtRes_p, ZP42R_p; function TZP42R: String; implementation function TZP42R: String; var N,_i :Integer; A :Array [0..2] of Real; B :Array [0..2] of Real; C :Array [0..5] of Real; begin Result := ''; { результат функции } { прототип оператора DАТА на FORTRANе } A[0] := 0.0; A[1] := 1.0; A[2] := 2.0; N := 3; ZP42R(A,B,C,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('TZP42R',Result); { вывод результатов в файл TZP42R.res } exit; end; end. Результат: В = ( 1., 1., 1. )