三次样条差值实验报告..docxVIP

数值计算方法实验报告专业：姓名：学号：实验二 三次样条插值1、根据给定的数据点集，求（1）自然三次样条插值，（2）第一类边界三次样条插值，（3）第二类边界（非自然）三次样条插值。分别编写出程序。2、使用上题中的程序，根据点（0，0），（1，4），（2，8），（3，9），（4，9），（5，8 和（ 6 6 ），求 3 种不同的三次样条插值，其中S’(0)=1,S’(6)=2;S’’(0)=1,S’’(6)=1 。在同一坐标系中，画出这 3 个三次样条插值和这些数据点。3、下面的表中给出了在洛杉矶的郊区 12 小时内每个小时的温度（华氏温度）。根据这些数据求自然三次样条插值。在同一个坐标系中，画出自然三次样条插值和这些数据。时间（am）度数时间（am）度数158757258858358960458106455711676571268实验程序：自然三次样条插值代码：function S3=ZiRanSanCi(X,Y)%定义自然三次样条插值函数ZiRanSanCi%X Y为已知数据N=length(X)-1;H=diff(X);D=diff(Y)./H;A=H(2:N-1);B=2*(H(1:N-1)+H(2:N));C=H(2:N);U=6*diff(D);for k=2:N-1 temp=A(k-1)/B(k-1); B(k)=B(k)-temp*C(k-1); U(k)=U(k)-temp*U(k-1);endM(N)=U(N-1)/B(N-1);for k=N-2:-1:1 M(k+1)=(U(k)-C(k)*M(k+2))/B(k);endM(1)=0;M(N+1)=0; %自然三次样条的端点约束起点二阶导数值等于终点二阶导数值等于0for k=0:N-1 S3(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1)); S3(k+1,2)=M(k+1)/2; S3(k+1,3)=D(k+1)-H(k+1)*(2*M(k+1)+M(k+2))/6; S3(k+1,4)=Y(k+1);End第一类边界三次样条插值代码function S1=DiYiSanCi(X,Y,dx0,dxn)N=length(X)-1;H=diff(X);D=diff(Y)./H;A=H(2:N-1);B=2*(H(1:N-1)+H(2:N));C=H(2:N);U=6*diff(D);B(1)=B(1)-H(1)/2;U(1)=U(1)-3*(D(1)-dx0);B(N-1)=B(N-1)-H(N)/2;U(N-1)=U(N-1)-3*(dxn-D(N));for k=2:N-1 temp=A(k-1)/B(k-1); B(k)=B(k)-temp*C(k-1); U(k)=U(k)-temp*U(k-1);endM(N)=U(N-1)/B(N-1);for k=N-2:-1:1 M(k+1)=(U(k)-C(k)*M(k+2))/B(k);endM(1)=3*(D(1)-dx0)/H(1)-M(2)/2;M(N+1)=3*(dxn-D(N))/H(N)-M(N)/2;for k=0:N-1 S1(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1)); S1(k+1,2)=M(k+1)/2; S1(k+1,3)=D(k+1)-H(k+1)*(2*M (k+1)+M(k+2))/6; S1(k+1,4)=Y(k+1);end第二类边界约束（非自然）三次样条代码：function S2=DiErSanCi(X,Y,ddx0,ddxn)%定义第二类边界（非自然）三次样条插值函数%X Y为已知数据，ddx0，ddxn分别为起点终点的二阶导数值N=length(X)-1;H=diff(X);D=diff(Y)./H;A=H(2:N-1);B=2*(H(1:N-1)+H(2:N));C=H(2:N);U=6*diff(D);U(1)=U(1)-H(1)*ddx0;U(N-1)=U(N-1)-H(N)*ddxn;for k=2:N-1 temp=A(k-1)/B(k-1); B(k)=B(k)-temp*C(k-1); U(k)=U(k)-temp*U(k-1);endM(N)=U(N-1)/B(N-1);for k=N-2:-1:1 M(k+1)=(U(k)-C(k)*M(k+2))/B(k);endM(1)=ddx0; %起点值的二阶导数M(N+1)=ddxn; %终点值的二阶导数for k=0:N-1 S2(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1)); S2(k+1