章绍辉版数学建模第五章作业.doc

第五章作业 第一题： 第三种边界条件 第三种边界条件要求给定三次样条s(x)在区间[x0,xn]的左右端点的一阶导数和。我们参考了例5.1.4在210页计算出了左右端点的一阶导数，其中= -3.3667、= 2.3333 MATLAB程序如下： x=[0,1,3,6,8,9];y=[-3.3667,3,1,2,0,2,4,2.33333]; 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,v1,k-.,u,v2,k--,u,v3,k:),hold on end title(第三种边界条件及其一、二、三阶导函数的图像) legend(三次样条（第三种边界条件）,样条的一阶导函数,样条的二阶导函数,样条的三阶导函数) plot([0,1,3,6,8,9],[3,1,2,0,2,4],ko),hold off 运行结果为： pp = form: pp breaks: [0 1 3 6 8 9] coefs: [5x4 double] pieces: 5 order: 4 dim: 1 ans = -0.0389 1.4056 -3.3667 3.0000 -0.3514 1.2890 -0.6722 1.0000 0.1696 -0.8197 0.2664 2.0000 -0.0848 0.7064 -0.0736 0 0.0678 0.1977 1.7345 2.0000 所绘制的图形如下： 结果说明： 计算结果说明该三次样条的分段多项式为： S(x)= 执行以下命令可以验算该三次样条在区间[0，9]的左端点x=0和右端点x=9的一阶导数分别为-3.3667和2.3333： 在MATLAB的command window运行：[1.*pp.coefs(1,3),d1s(9,8,pp.coefs(5,:))] 运行结果为： ans = -3.3667 2.3333 第四种边界条件 第四种边界条件要求给定三次样条s(x)在区间[x0,xn]的左右端点的二阶导数和。我们分别取左右端点的二阶导数为=-2,=3 MATLAB程序如下： x=[0,1,3,6,8,9];y=[-2,3,1,2,0,2,4,3]; pp=csape(x,y,second) 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,v1,k-.,u,v2,k--,u,v3,k:),hold on end title(第四种边界条件及其一、二、三阶导函数的图像) legend(三次样条（第四种边界条件）,样条的一阶导函数,样条的二阶导函数,样条的三阶导函数) plot([0,1,3,6,8,9],[3,1,2,0,2,4],ko),hold off 运行结果为： pp = form: pp breaks: [0 1 3 6 8 9] coefs: [5x4 double] pieces: 5 order: 4 dim: 1