--实验代数插值实验.docVIP

实验7 代数插值实验（2） 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握spline(样条)插值算法，提高编程能力和解决插值问题的实践技能。 实验内容：编写程序用计算机求解三次压紧样条曲线，经过点（0，0.0），（1，0.5），（2，2.0）和（3，1.5），而且一阶导数边界条件S’(0)=0.2和S’(3)=-1 实验原理 若已知N+1个点的及其一阶导数的边界条件S’(a)=和S’(b)=,则存在唯一的三次样条曲线。 构造并求解下列线性方程组 构造分段函数S(x): ; ; ; 实验说明 需要建立两个数组X，Y分别存放（，…，，（，…，，建立两个变量dxo、 dxn分别存放左右端点的一阶导数 dxo = S(x0) , dxn = S(xn) 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验总结 实验报告: 根据实验情况和结果撰写并递交实验报告。 参考程序（MATLAB程序） 如下：主函数： x=[0 1 2 3] y=[0 0.5 2.0 1.5] dxo=0.2 dxn=-1 s=sanci(x,y,dxo,dxn)（调用） x1=0:0.001:1 y1=polyval(s(1,:),x1-x(1),3); x2=1:0.001:2 y2=polyval(s(2,:),x2-x(2),3); x3=2:0.001:3 y3=polyval(s(3,:),x3-x(3),3); plot(x1,y1,x2,y2,*,x3,y3,-) 调用函数：function s=sanci(x,y,dxo,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)-dxo); 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); end m(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); end m(1)=3*(d(1)-dxo)/h(1)-m(2)/2; m(n+1)=3*(dxn-d(n))/h(n)-m(n)/2; for k=0:n-1 s(k+1,1)=(m(k+2)-m(k+1))/(6*h(k+1)); s(k+1,2)=m(k+1)/2; s(k+1,3)=d(k+1)-h(k+1)*(2*m(k+1)+m(k+2))/6; s(k+1,4)=y(k+1); end 输入后结果如下： x = 0 1 2 3 y = 0 0.5000 2.0000 1.5000 dxo = 0.2000 dxn = -1 s = 0.4800 -0.1800 0.2000 0 -1.0400 1.2600 1.2800 0.5000 0.6800 -1.8600 0.6800 2.0000 然后绘图：