Линейная интерполяция Опции

[volvo]

Что значит "не хочет"?

The Data :
-1.0000000 1.009000000000000E+000
-0.9000000 4.740000000000000E-001
-0.8000000 2.475000000000000E-001
-0.7000000 4.091000000000000E-001
-0.6000000 6.512000000000000E-001
-0.5000000 6.007000000000000E-001
-0.4000000 3.610000000000000E-002
-0.3000000 -7.662000000000000E-001
-0.2000000 -1.381400000000000E+000
-0.1000000 -1.442900000000000E+000
The polynomial :
Poly[ 9]= 7.502480158738636E+003
Poly[ 8]= 3.788740079369168E+004
Poly[ 7]= 8.179012896833815E+004
Poly[ 6]= 9.828562500009664E+004
Poly[ 5]= 7.179005381951237E+004
Poly[ 4]= 3.279548854169682E+004
Poly[ 3]= 9.368500327389358E+003
Poly[ 2]= 1.635103664683947E+003
Poly[ 1]= 1.559646261906029E+002
Poly[ 0]= 4.518900000004579E+000
X Interpolated Y value

Ты же получаешь полином... Что тебе не нравится?

[Jill]

спасибо, volvo, уже разобралась что-то у меня проглючило (...или это меня проглючило? )

СПАСИБО!!!

[Jill]

настал черед коеффициентов полинома Ньютона... вычисляются они по формуле, из которой я никак не могу просчитать дельту...
http://www.eva.ru/pictures/album_photos/79...f?1136569231119 - тут изображение формул
остальное нацарапала вроде:

Program Polyn;
const
TNArraySize = 50;            { Size of vectors }
h = 0.1;
N = 3;
type
  TNvector = array[0..TNArraySize] of Real;
var
  NumPoints : integer;              { Number of data points }
  YData : TNvector;          { Data points (Y) }
  Poly : TNvector;                  { The constructed polynomial }
  k : Integer;

Function Delta(Y:TNvector):Real;     {ВОТ ЗДЕСЬ ПРОБЛЕМА}
var k:integer;
begin
     for k:=0 to N-1 do begin
     Delta:=Y[k+1]-Y[k];
     end;
end;

Function Factorial(k:Integer):Integer;
  Var m,l:Integer;
Begin;
  If k<>0 then Begin
    m:=k;
    For l:=(k-1) downto 2 do m:=m*l;
    Factorial:=m;
  End
  Else Factorial:=1;
End; {Factorial}

Function Step(k:Integer):real;
Begin
     Step:=exp(k*ln(h));
End;

Begin
for k:=0 to n do begin
   write('Y[',k,'] = ');
   readln(YData[k]);
   end;
for k:=0 to n do begin
Poly[k]:=Delta(YData[n-k])/(Factorial(k)*Step(k));  {НУ И ЗДЕСЬ СООТВЕТСТВЕННО...}
end;
end.

как ее, эту дельту, написать?

Сообщение отредактировано: volvo - 5.11.2006 20:20

[volvo]

Вот тут глянь, по-моему, там более доступно, чем в том архиве, на который я раньше давал ссылку...