|
Текст подпрограммы и версий 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.
procedure ZP47R(X :Real; var C :Array of Real; N :Integer;
var SUM :Real);
Параметры
| X - | заданное значение аргумента (тип: вещественный); |
| C - | вещественный вектор длины N, содержащий заданные коэффициенты разложения; при этом полиному порядка К соответствует коэффициент С (К + 1); |
| N - | длина вектора С (тип: целый); |
| SUM - | вычисленное значение частичной суммы (тип: вещественный). |
Версии
| ZP47E - | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева для данных расширенной (Extended) точности. |
|
ZP48R - ZP48E | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева для данных одинарной и расширенной (Extended) точности. |
|
ZP49R - ZP49E | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева четного порядка для данных одинарной и расширенной (Extended) точности. |
|
ZP50R - ZP50E | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева четного порядка для данных одинарной и расширенной (Extended) точности. |
|
ZP51R - ZP51E | вычисление частичной суммы ряда Фурье по ортогональным полиномам Чебышева нечетного порядка для данных одинарной и расширенной (Extended) точности. |
|
ZP52R - ZP52E | вычисление частичной суммы ряда Фурье по смещенным ортогональным полиномам Чебышева нечетного порядка для данных одинарной и расширенной (Extended) точности. |
Вызываемые подпрограммы: нет
Замечания по использованию
| Смещенные ортогональные многочлены Чебышева ТS (К, Х) определяются по формуле ТS (К, Х) = Т (К, 2 * Х - 1), где К - порядок полинома, Х - аргумент и Т (К, Х) - значение обычного полинома Чебышева порядка К в точке Х: Т (К, Х) = cos (К * АRССОS (Х)). |
Unit tzp47r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, ZP47R_p;
function tzp47r: String;
implementation
function tzp47r: String;
var
N :Integer;
X,SUM :Real;
const
C :Array [0..2] of Real = ( 1.0,2.0,3.0 );
begin
Result := ''; { результат функции }
X := 0.5;
N := 3;
ZP47R(X,C,N,SUM);
Result := Result + Format(' %3d ',[N]) + #$0D#$0A;
Result := Result + Format(' %20.16f ',[SUM]) + #$0D#$0A;
UtRes('tzp47r',Result); { вывод результатов в файл tzp47r.res }
end;
end.
Результат: SUМ = 0.5
При тех же самых значениях входных параметров как и в примере выше, при обращении к другим версиям, будут выданы следующие результаты.
ZP48R SUМ = -2.0 ;
ZP49R SUМ = -1.5 ;
ZP50R SUМ = 2.0 ;
ZP51R SUМ = 0.0 ;
ZP52R SUМ = 0.0 .