examp8_13.m

来自「东北大学 薛定谔的高等应用数学问题MATLAB求解的源代码」· M 代码 · 共 8 行

x0=[0,0.4,1 2,pi]; y0=sin(x0); ezplot('sin(t)',[0,pi]); hold on
sp1=csapi(x0,y0); fnplt(sp1,'--');  % 三次分段多项式样条插值
sp2=spapi(5,x0,y0); fnplt(sp2,':')  % 5 次 B 样条插值

x=0:.12:1; y=(x.^2-3*x+5).*exp(-5*x).*sin(x);
ezplot('(x^2-3*x+5)*exp(-5*x)*sin(x)',[0,1]), hold on
sp1=csapi(x,y); fnplt(sp1,'--'); sp2=spapi(5,x,y); fnplt(sp2,':')