第五章数模作业.doc

数模作业 朱镇东 20112201079 1、参照第（2）种边界条件的三次样条结论，不妨给定样条在左右端点的一阶导数，其中s(x)在区间[，]的左端点的一阶导数为-3.37，右端点的一阶导数为2.33。则其三次样条（见下图） 绘制上图的MATLAB脚本如下： x=[0,1,3,6,8,9];y=[-3.37,3,1,2,0,2,4,2.33]; pp=csape(x,y,complete),pp.coefs s=@(t,tj,c)c(1).*(t-tj).^3+c(2).*(t-tj).^2+c(3).*(t-tj)+c(4); d1s=@(t,tj,c)3.*c(1).*(t-tj).^2+2.*c(2).*(t-tj)+c(3); d2s=@(t,tj,c)6.*c(1).*(t-tj)+2.*c(2); d3s=@(t,tj,c)6.*c(1).*ones(size(t)); for k=1:pp.pieces c=pp.coefs(k,:);u=x(k):.01:x(k+1);v=s(u,x(k),c); v1=d1s(u,x(k),c);v2=d2s(u,x(k),c);v3=d3s(u,x(k),c); plot(u,v,k,u,v,k-,