魏凤英的五点三次平滑

wmlydia

纠结了好久，魏凤英那本书47页的五点三次平滑，是三次多项式拟合吧？为毛我编出来的结果几乎没有平滑呢？用EXCEL做的三次多项式拟合，却得到了想要的结果，一条光滑的峰值在右侧的类似正弦函数的曲线。。。还请大家来帮忙解释下呗
下面是用的程序代码
implicit none
integer,parameter::nx=164
real  sstapj(nx+7)!!!!另外的7个数据是缺测值!!!!!!!!!!1
real  nh(nx+7)
integer ix,j
!!!!read data!!!!!!!!!!!
open(10,file='I:\ep3.19\hc.grd',form='binary')
read(10) (sstapj(ix),ix=1,nx+7)
close(10)
!!!!!五点三次平滑!!!!!!!!!!!
nh(1)=(69*sstapj(1)+4*sstapj(2)-6*sstapj(3)+4*sstapj(4)-sstapj(5))/70
nh(2)=(2*sstapj(1)+27*sstapj(2)+12*sstapj(3)-8*sstapj(4)-sstapj(5))/35
do ix=3,nx-2
nh(ix)=(-3*sstapj(ix-2)+12*sstapj(ix-1)+17*sstapj(ix)+12*sstapj(ix+1)-3*sstapj(ix+2))/35
enddo
nh(nx-1)=(2*sstapj(nx-4)-8*sstapj(nx-3)+12*sstapj(nx-2)+27*sstapj(nx-1)+2*sstapj(nx))/35
nh(nx)=(-sstapj(nx-4)+4*sstapj(nx-3)-6*sstapj(nx-2)+4*sstapj(nx-1)+69*sstapj(nx))/70
do ix=nx+1,nx+7
nh(ix)=-1e+30
enddo
open(20,file='I:\ep3.19\nhzf.grd',form='binary')
write(20) (nh(ix),ix=1,nx+7)
close(20)
end
错误在哪里呢？{:soso_e149:}

wmlydia

绿色的是平滑后的线，白色的原来的线，两者都差不多重合了。。

lqouc

wmlydia 发表于 2014-4-3 22:24
绿色的是平滑后的线，白色的原来的线，两者都差不多重合了。。

个人感觉程序没什么问题，画的图我觉得也没问题，要是特别光滑了我反而觉得不对。

wmlydia

lqouc 发表于 2014-4-4 08:52
个人感觉程序没什么问题，画的图我觉得也没问题，要是特别光滑了我反而觉得不对。

就是觉得拟合的太好了，想要的平滑效果没得到，还是决定用多项式拟合了，谢谢！

dyhtt29

nickbsb

看看            

短腿的反击

非常有用！