func3p2.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 13 行

%func3p2飞机降落曲线
yy=linspace(0,1000,121)
xx=linspace(0,12000,121)
for i=0:120
  y=1000*sub3p2(i*100)
  xx(i+1)=i*100
  yy(i+1)=y
end
plot(xx,yy,'-')
%figure enclosed
s=dsolve('D4y=0','y(0)=0','y(12000)=1000','Dy(0)=0','Dy(12000)=0')
%解得等价多项式为:-1/864000000*t^3+1/48000*t^2
%图形基本满足实际情况，Hermite插值精度较高