数值分析课件15样条插值30.pptxVIP

;一、三次样条插值问题的提法;;;;;（2）构造三弯矩方程;（2）构造三弯矩方程;;四、样条插值的收敛性;Chzh1.m;三次样条插值常用命令 (1)y0=interp1(x,y,x0,spline) y0=interp2(x,y,x0,spline) （2）y0=spline(x,y,x0) (一维) （3）pp=spline(x,y);（结构数组） y0=fnval(pp,x0); poly(x0,y0), pp.coefs （系数） drp=fnder(pp,n), drv=fnval(drp,a) 求导数 itp=fnint(pp), itv=fnval(itp,[a,b])*[-1,1] 求积分 fnplt(pp),fnplt(drp,n),fnplt(itp) 绘图 ?（4） pp=csape(x,y,边界条件选项); （5） pp=csapi(x,y); 演示程序 sptest.m ;【例】根据连续时间函数的采样数据，利用spline重构该连续函数，并检查重构误差。 t=-5:0.5:5;w=exp(-abs(t)); N0=length(t);tt=linspace(t(1),t(end),10*N0); ww=spline(t,w,tt); error=max(abs(ww-exp(-abs(tt)))) plot(tt,ww,b);hold on stem(t,w,filled,r);hold off error = 0.0840;【例】对于函数 ，很容易求得 。本例将借此演示样条函数求数值不定积分、导函数的能力。;（2）不定积分样条函数、导数样条函数的使用 % 计算y(x)在区间[1,2]上的定积分 DefiniteIntegral.bySpline=ppval(int_pp,[1,2])*[-1;1]; DefiniteIntegral.byTheory=(1-cos(2))-(1-cos(1)); %计算dy(3)/dx Derivative.bySpline=fnval(der_pp,3); Derivative.byTheory=cos(3); Derivative.byDiference=(sin(3.01)-sin(3))/0.01; DefiniteIntegral,Derivative DefiniteIntegral = bySpline: 0.9563 byTheory: 0.9564 Derivative = bySpline: -0.9895 byTheory: -0.9900 byDiference: -0.9907 ;（3）绘制三个样条函数的图形 fnplt(pp,b-);hold on fnplt(int_pp,m:),fnplt(der_pp,r--);hold off legend(y(x),S(x),dy/dx) ;习题五 P208----- 14 数值实验题五 P210----- 4 ginput(10)