А можно подробнее про интерполяционные формулы, которые Вы используете? Где у них разности интерполируемой функции в знаменателе?
Вот код на фортране:
Код:
SUBROUTINE polint(xa,ya,n,x,y,dy)
INTEGER n,NMAX
REAL dy,x,y,xa(n),ya(n)
PARAMETER (NMAX=10000)
INTEGER i,m,ns
REAL den,dif,dift,ho,hp,w,c(NMAX),d(NMAX)
ns=1
dif=abs(x-xa(1))
do 11 i=1,n
dift=abs(x-xa(i))
if (dift.lt.dif) then
ns=i
dif=dift
endif
c(i)=ya(i)
d(i)=ya(i)
11 continue
y=ya(ns)
ns=ns-1
do 13 m=1,n-1
do 12 i=1,n-m
ho=xa(i)-x
hp=xa(i+m)-x
w=c(i+1)-d(i)
den=ho-hp
if(den.eq.0.)pause 'failure in polint'
den=w/den
d(i)=hp*den
c(i)=ho*den
12 continue
if (2*ns.lt.n-m)then
dy=c(ns+1)
else
dy=d(ns)
ns=ns-1
endif
y=y+dy
13 continue
return
END
xa,ya — входные массивы аргумента и функции, x — интересующее значение аргумента,y — сынтерполированное значение функции.
-- Чт окт 21, 2010 14:32:17 --У Вас по координате Z три идеальных горизонтальных участка (2 ступеньки).
Это лишь начало данных по координате Z.
Вот графики по координатам x, y, z соответственно:
Несмотря на разрывы в начале данных по X и Y, интерполяция происходит хорошо, а вот с медленно изменяющимися внечале данными по Z — проблема.
