三次样条插值实验报告

合集下载

数值分析实验报告-插值、三次样条

实验报告：牛顿差值多项式&三次样条问题：在区间[-1,1]上分别取n=10、20用两组等距节点对龙格函数21()25f x x作多项式插值及三次样条插值，对每个n 值，分别画出插值函数及()f x 的图形。

实验目的：通过编程实现牛顿插值方法和三次样条方法，加深对多项式插值的理解。

应用所编程序解决实际算例。

实验要求：1．认真分析问题，深刻理解相关理论知识并能熟练应用；2．编写相关程序并进行实验；3．调试程序，得到最终结果；4．分析解释实验结果；5．按照要求完成实验报告。

实验原理：详见《数值分析第5版》第二章相关内容。

实验内容：（1）牛顿插值多项式1.1 当n=10时：在Matlab 下编写代码完成计算和画图。

结果如下：代码：clear allclcx1=-1:0.2:1;y1=1./(1+25.*x1.^2);n=length(x1);f=y1(:);for j=2:nfor i=n:-1:jf(i)=(f(i)-f(i-1))/(x1(i)-x1(i-j+1));endendsyms F x p ;F(1)=1;p(1)=y1(1);for i=2:nF(i)=F(i-1)*(x-x1(i-1));p(i)=f(i)*F(i);endsyms PP=sum(p);P10=vpa(expand(P),5);x0=-1:0.001:1;y0=subs(P,x,x0);y2=subs(1/(1+25*x^2),x,x0);plot(x0,y0,x0,y2)grid onxlabel('x')ylabel('y')P10即我们所求的牛顿插值多项式，其结果为：P10（x）=-220.94*x^10+494.91*x^8-9.5065e-14*x^7-381.43*x^6-8.504e-14*x^5+123.36*x^4+2.0202e-1 4*x^3-16.855*x^2-6.6594e-16*x+1.0并且这里也能得到该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.1）。

关于用追赶法求解系数矩阵的三次样条插值实验报告。

x(n)=y(n); for i=(n-1):-1:1 x(i)=y(i)-u(i)*x(i+1); end 用 matlab 编程如下: function[s,y0]=spline3 (x,y,x0) %x,y 为数表 x0 为插值点 s 表示插值函数 y0 为 x0 对应的插值函数值 syms t n=length(x); %得出 n for i=1:n-1; h(i)=x(i+1)-x(i); end for i=2:n-1; lamda(i)=h(i)/(h(i-1)+h(i)); miu(i)=1-lamda(i); g(i)=3*(lamda(i)*((y(i)-y(i-1))/h(i-1))+miu(i)*((y(i+1)-y( i))/h(i))); end g(1)=3*((y(2)-y(1))/h(1)); g(n)=3*((y(n)-y(n-1))/h(n-1)); %前边求出 lamda miu 和 g 从而可以确定系数矩阵 miu(1)=1; miu(4)=0;
实验 for i=2:(n-1) 内容
L(i)=b(i)-a(i)*u(i-1); y(i)=(d(i)-y(i-1)*a(i))/L(i); u(i)=c(i)/L(i); end L(n)=b(n)-a(n)*u(n-1); y(n)=(d(n)-y(n-1)*a(n))/L(n); %“赶”的过程
数值分析实验报告
实验名称
关于用追赶法求解系数矩阵的三次样条插值实验报告。
实验结合求解线性方程组的解。目的算法如下：
会使用 matlab 语言编程使用追赶法 llowup(a,b,c,d) n=length(d); a(1)=0; %“追”的过程 L(1)=b(1); y(1)=d(1)/L(1); u(1)=c(1)/L(1);

学会运用三次样条函数插值

实验目的：学会运用三次样条函数进行函数插值实验内容：给定插值条件如下：i 0 1 2 3 4 5 6 7X i8.125 8.4 9.0 9.485 9.6 9.959 10.166 10.2 Y i0.0774 0.099 0.280 0.60 0.708 1.200 1.800 2.177 作三次样条函数插值，取第一类边界条件Y0’=0.01087 Y7’=1001.画出插值曲线的图像。

程序如下：%将插值点的坐标值输入matlabx=[8.125 8.4 9.0 9.485 9.6 9.959 10.166 10.2];y=[0.0774 0.099 0.280 0.60 0.708 1.200 1.800 2.177];%调用三次样条函数，并给出第一类边界条件即 Y0’=0.01087 Y7’=100 cs = spline(x,[0.01087 y 100]);%给定x（0）到x(7)之间的作图点数，即1000xx=linspace(8.125,10.2,1000);%绘制图像plot(x,y,'o',xx,ppval(cs,xx),'-')grid on插值曲线：2．逆时针旋转座标轴45o保持（1．）中结点和边界条件的几何关系不变，再次作三次样条函数插值，画出插值曲线的图像。

程序如下：x=[8.125 8.4 9.0 9.485 9.6 9.959 10.166 10.2];y=[0.0774 0.099 0.280 0.60 0.708 1.200 1.800 2.177];%用坐标转换求出坐标轴旋转45度后插值点新的坐标值x1=x*cos(pi/4)+y*sin(pi/4);y1=-x*cos(pi/4)+y*sin(pi/4);%求出坐标轴转换后新的边界值Y0=(0.01087-1)/(0.01087+1);Y7=(100-1)/(100+1);cs = spline(x1,[Y0 y1 Y7]);xx=linspace(x1(1),x1(8),1000);plot(x1,y1,'o',xx,ppval(cs,xx),'-')grid on其中在求解边界条件用转换关系tan(45)θ-，θ为在原图像中边界点的切线与原x坐标轴的夹角差值曲线：。

Matlab实验报告六(三次样条与分段线性插值)范文

1.分析问题
本题是给出粗略等分点让你插入更多点用双线性插值法来作出更清晰的山区地貌图。
2.问题求解
x=0:400:2800;
y=0:400:2400;
z=[1430 1450 1470 1320 1280 1200 1080 940;
1450 1480 1500 1550 1510 1430 1300 1200;
2.分段线性插值与计算量与n无关;n越大，误差越小.
3.三次样条插值比分段线性插值更光滑。
4.‘linear’：分段线性插值；‘spline’：三次样条值。
【实验环境】
MatlabR2010b
二、实验内容
问题1对函数，x[-5,5]，分别用分段线性插值和三次样条插值作插值（其中插值节点不少于20），并分别作出每种插值方法的误差曲线.
本次实验因为是我们课本没有的内容，心理上给了我很大的压力，幸好我们还能根据老师的课件以及例题去掌握这次实验所需要的各种插值法，但结果还好，两道题都做出来了。
plot(x,y,'*',x1,yl,'r',x1,y2,'b')
y0=1./(1+x1.^2);
y3=yl-y0;
y4=ys-y0;
holdon
plot(x1,y3,'y',x1,y4,'g')
3.结果
4误。
问题2山区地貌图在某山区（平面区域（0，2800）（0，2400）内，单位：米）测得一些地点的高程（单位：米）如表1，试作出该山区的地貌图.
1.分析问题
本题先取出少量的插值节点并作出图形，再用分段线性插值法和三次样条插值法做出更精确的图形，最后在作出误差曲线。

样条插值实验报告

四、三次样条插值1. 样条函数插值的原理给定区间［a,b］上划分A:a=x<x<<x<x=b，若分段函数S(x)满足:01n-1n1.S(x)在各个子区间［x,x］,i=0,1,,n-1上均为x的三次多项式；ii+12.S(x)在整个区间［a,b］上有直至二阶的连续导数。

则称S(x)为［a,b］上依次划分的三次样条函数，简称样条函数。

具体地有分段表达式：ax3+bx2+cx+d,x G[x,x]000001ax3+bx2+cx+d,x G[x,x]111112S(x)=\ax3+bx2+cx+d,x G[x,x](1)222223ax3+bx2+cx+d,x G[x,x]、°*n-1n—T•••n-1n-1n-1n共有4n个参数a,b,c,d,i=0,1,,n，它们在内节点处满足iiii'S(x)=S(x),…i-0i+0<S'(x)=S'(x),i=1,2,,n-1.(2)i-0i-0S''(x)=S''(x),Ji-0i+0满足样条函数定义的函数集合称为分划A上的三次样条函数空间，记为S(3,A),可以证明S(3,A)为线性空间。

若S(x)G S(3,A)，且进一步满足插值条件S(x)=y=f(x),i=0,1,,n(3)iii其中y为节点x处的给定函数值(若被插函数了(x)已知;••则用了(x)代替之)，iii则称S(x)为以x,x,,x,x为节点的三次样条函数。

01n-1n其中式(3)插值节点提供了n+1个约束条件；加上式(2)的3n-3个，合起来共有4n-2个；欲求4n个待定参数的唯一解；尚缺两个条件。

这两个条件一般由样条函数的边界条件提供。

常用三类边界条件;他们分别与三次样条函数;构成不同边界条件的样条函数插值问题。

2. 三类样条函数插值问题2.1第二类边界条件给定边界条件两端的一阶导数值：S'(x)=y'=m,S'(x)=y'=m000nnn这相当于样条两短处的方向给定(压铁在两端点的压力方向确定)，对应的插值问题如下：对于分划A:a=x<x<<x<x=b,给定节点对应的函数值01n—1ny,y,y,,y,以及两端点处的一阶导数值y'=m,y'=m，求三次样条函数012n00nnS(x)，使…f S(x)=y,i=0,1,,n2iiI S'(x)=m,S'(x)=mJ00n…n2.2第一类边界条件给定边界两端的二阶导数值：S''(x)=y''=M,S''(x)=y''=M000nnn这相当于在样条两端处外加一个力矩，使梁两端点处有相应的曲率。

三次样条插值方法的应用

CENTRAL SOUTH UNIVERSITY数值分析实验报告三次样条插值方法的应用一、问题背景分段低次插值函数往往具有很好的收敛性，计算过程简单，稳定性好，并且易于在在电子计算机上实现，但其光滑性较差，对于像高速飞机的机翼形线船体放样等型值线往往要求具有二阶光滑度，即有二阶连续导数，早期工程师制图时，把富有弹性的细长木条（即所谓的样条）用压铁固定在样点上，在其他地方让他自由弯曲，然后沿木条画下曲线，称为样条曲线。

样条曲线实际上是由分段三次曲线并接而成，在连接点即样点上要求二阶导数连续，从数学上加以概括就得到数学样条这一概念。

下面我们讨论最常用的三次样条函数及其应用。

二、数学模型样条函数可以给出光滑的插值曲线（面），因此在数值逼近、常微分方程和偏微分方程的数值解及科学和工程的计算中起着重要的作用。

设区间[]b ,a 上给定有关划分b x x n =<<<=Λ10x a ，S 为[]b ,a 上满足下面条件的函数。

● )(b a C S ,2∈；● S 在每个子区间[]1,+i i x x 上是三次多项式。

则称S 为关于划分的三次样条函数。

常用的三次样条函数的边界条件有三种类型：● Ⅰ型 ()()n n n f x S f x S ''0'',==。

● Ⅱ型 ()()n n n f x S f x S ''''0'''',==，其特殊情况为()()0''''==n n x S x S 。

● Ⅲ型 ()()Λ3,2,1,0,0==j x S x S n j j ，此条件称为周期样条函数。

鉴于Ⅱ型三次样条插值函数在实际应用中的重要地位，在此主要对它进行详细介绍。

三、算法及流程按照传统的编程方法，可将公式直接转换为MATLAB可是别的语言即可；另一种是运用矩阵运算，发挥MATLAB在矩阵运算上的优势。

LAB03_三次样条插值实验

Lab03．三次样条插值实验【实验目的和要求】1．使学生深入理解三次样条插值法，深入进行程序设计能力训练；2．对第一与第二种边界条件，按三弯矩法，通过用Matlab 语言设计计算三次样条插值的程序，以提高学生程序设计的能力。

【实验内容】1．根据Matlab 语言特点，描述三次样条插值法。

2．对第一与第二种边界条件，按三弯矩法，用Matlab 语言设计计算三次样条插值的程序。

3对(1) 自然边界条件0)0.1()2.0(=''=''S S ；(2) 第一种边界条件55741.1)0.1( ,20271.0)2.0(='='S S . 输出用追赶法解出的弯矩向量),,(521M M M 和)1.02.0(i S + (i =0，1，…，8）的值，并画出)(x S y =的图形。

4．完成教材P45例8的计算，并将计算结果与Langrage 插值法计算的结果进行比较，由此说明三次样条插值的优越性。

【实验仪器与软件】1．CPU 主频在1GHz 以上，内存在128Mb 以上的PC ；2．Matlab 6.0及以上版本。

实验讲评：实验成绩：评阅教师：2011 年月日Lab03．三次样条插值实验一、算法描述1．根据Matlab语言特点，描述三次样条插值法.答：S(x) 在[x j, x j+1]（j=1,2,⋯,n-1）上是三次多项式，于是S"(x)在[x j, x j+1] 上是一次多项式，如果S"(x) 在[x j,x j+1]（j=1,2,⋯,n-1）两端点上的值已知，设S"(x j)=M j，S"(x j+1)=M j+1，则S"(x) 的表达式为：= ，其中h j ＝x j+1-x j，对S"(x) 进行两次积分，则得到1 个具有2二、程序设计2．对第一与第二种边界条件，按三弯矩法，用Matlab语言设计计算三次样条插值的程序。

Matlab实验报告六(三次样条与分段线性插值)

数学与信息科学系实验报告
实验名称插值与拟合
所属课程数学软件与实验
实验类型综合型实验
专业信息与计算科学
班级
学号
姓名
指导教师
一、实验概述
【实验目的】
学会在matlab环境下使用几种不同的插值法和拟合两种方法构造函数依据已经知道的某些特殊点来推测实际问题中需要知道但又不便于测量出来的量。
【实验原理】
1.z=interp2(x0,y0,z0,x,y,’method’): 要求x0,y0单调；x, y可取为矩阵, 或x取行向量, y取为列向量, x,y的值分别不能超出x0,y0的范围。
2.分段线性插值与计算量与n无关;n越大, 误差越小.
3.三次样条插值比分段线性插值更光滑。
4.‘linear’ : 分段线性插值；‘spline’ : 三次样条
二、实验内容
问题1 对函数, x([-5,5], 分别用分段线性插值和三次样条插值作插值（其中插值节点不少于20）, 并分别作出每种插值方法的误差曲线.
1180 1320 1450 1420 1400 1300 700 900];
mesh(x,y,z)
xi=0:20:2800;
yi=0:20:2400;
zi=interp2(x,y,z,xi',yi,'cubic');
mesh(xi,yi,zi)
3.结果
4.结论及分析
通过实验,结果正确,分析无误。
三、实验小结
1270 1500 1200 1100 1350 1450 1200 1150
1230 1390 1500 1500 1400 900 1100 1060
1180 1320 1450 1420 1400 1300 700 900

三次样条插值

三次样条插值的数值实验姓名：王维滨学号：0842011157 姓名：李佳乐学号：0842011034 姓名：谢朝学号：0842011062 姓名：杨其荣学号：08420110721.实验项目的性质和任务对三次样条插值进一步理解，并编写matlab 程序，实现这些功能2.算法设计和matlab 编程。

总前提：x i 第i 点的横坐标i y 第i 点的纵坐标i M ,,s 的记号,,,,,,23()()()()()()()2!3!i i i i i i i s x s x s x y s x x x x x x x =+-+-+- 有数值逼近书上的推导，我们令：111i i i i i x x u x x -+--=-，111i i i i i x x x x λ++--=-，11111111f (,,)i i i i i i i i i i i i i y y y y x x x x x x x x x +-+--++------=- 由于未知数的数目多于方程的个数，我们需要增加两个条件才能唯一确定一个分段三次函数1）D1的三次样条插值a ．实验方案与原理：我们加上条件：,,,,11()(),()()n n s x f x s x f x ==我们建立三弯矩方程组：1211211111126(,,)26(,,),2,3.....126(,,)i i i i i i i i n n n n n M M f x x x u M M M f x x x i n M M f x x x λ-+-+--+=⎧⎪++==-⎨⎪+=⎩然后采用追赶法迭代求方程组，但是我们在程序中采用简单的方法（矩阵计算）直接求解降低编程难度,2）D2三次样条插值223123211i 11112111126(,,)26(,,),3,4.....226(,,)i i i i i i i n n n n n n n n M M f x x x u M u M M M f x x x i n u M M f x x x M λλλ-+-+----+-+=-⎧⎪++==-⎨⎪++=-⎩3）D3三次样条插值22321231i 111211126(,,)26(,,),3,4.....126(,,)n i i i i i i i n n n n n n n M M u M f x x x u M M M f x x x i n M u M M f x x x λλλ-+-+--+++=⎧⎪++==-⎨⎪++=⎩b ．编写程序：见附录调用test,每次显示一个图像，关闭后按回车继续显示下一幅。

三次样条实验报告范文三次样条插值的实验报告范文

三次样条实验报告范文三次样条插值的实验报告范文湘潭大学实验报告课程名称计算机图形学实验名称参数三次样条的绘制页数专业计算机科学与技术班级一班同组者姓名学号2022551208姓名刘兆臣实验日期2022.05.05实验目的使学生掌握三次参数样条的定义、画法和程序的编写。

实验内容和要求给定型值点，要求用VC++6.0画出通过给定型值点的参数三次样条曲线。

实验方案设计给定型值点，计算出参数三次样条曲线的每个区间段的代数式，由Hermit曲线定义画出每个区间的曲线。

开始开始给定型值点给定型值点计算出参数计算出参数ifn=1Yifn=1N给系数矩阵赋值给系数矩阵赋值Ifn=2YIfn=2给定系数值给定系数值N求解方程组，计算出各点导数求解方程组，计算出各点导数Ifi<ni=0Ifi<ni=0Y计算第i段的二次项，三次项系数计算第i段的二次项，三次项系数结束结束t=0t=0Ift<tt[i-1]Ift<tt[i-1]NYt=t+et=t+e计算出第i段中各点的某,y值计算出第i段中各点的某,y值连线段连线段i++i++程序运行和实验结果说明和分析。

使用VC++6.0运行程序后得到如下图形图为四段曲线组成的三次参数样条曲线，其中各段的绘制是通过给定了型值点的相关参数计算出型值点的导数，再由Hermit曲线知识，在各段上以直代曲绘制出每段的图形。

基本达到了实验目的，完成实验要求。

性能、扩展性等方面存在的不足和可能的改进之处。

不足：在源代码中的n个点采用了数据初始化的方法给出且给定了型值点的个数。

由于给定型值点较少，三次参数样条曲线看起来不够明显。

可改进：可将初始化的型值点数据去除，采用手动键盘输入或文件输入的方法导入多个型值点数据。

附件一源程序，执行程序，符号列表文件。

#include<graphic.h>#include<math.h>#include<conio.h>main(){intgdriver=DETECT,gmode;float某[100],y[100],a[100],b[100],c[100];floatp某[100],py[100],q某[100],qy[100],tt[100];floatd某[100],dy[100];inti,n=4,t,e=3;floatb某3,b某4=0.0,by3=0.0,by4,c某,cy;initgraph(&gdriver,&gmode,"");for(i=0;i<n;i++) {a[i]=0.0;b[i]=0.0;c[i]=0.0;p某[i]=0.0;py[i]=0.0;d某[i]=0.0;dy[i]=0.0;tt[i]=0.0;q某[i]=0.0;qy[i]=0.0;}p某[0]=1.0;py[0]=1.0;p某[4]=1.0;py[4]=1.0;某[0]=10.0;y[0]=110.0;某[1]=40.0;y[1]=100.0;某[2]=80.0;y[2]=90.0;某[3]=130.0;y[3]=95.0;某[4]=200.0;y[4]=105.0;moveto(某[0],y[0]);for(i=0;i<n;i++)putpi某el(某[i],y[i],15);putpi某el(某[0],y[0],15);for(i=0;i<n;i++)tt[i]=qrt((某[i]-某[i-1])某(某[i]-某[i-1])+(y[i]-y[i-1])某(y[i]-y[i-1]));if(n==1)gotopO;for(i=1;i<=n-1;i++){a[i]=2某(tt[i]+tt[i+1]);b[i]=tt[i+1];c[i]=tt[i];d某[i]=3某(tt[i]某(某[i+1]-某[i])/tt[i+1]+tt[i+1]某(y[i]-y[i+1])/tt[i]);}d某[i]=d某[1]-tt[2]某p某[0];d某[n-1]=d某[n-1]-tt[n-1]某p某[n];dy[1]=dy[1]-tt[2]某py[0];dy[n-1]=dy[n-1]-tt[n-1]某py[n];if(n==2){p某[1]=d某[1]/a[1];py[1]=dy[1]/a[1];gotopO;}c[1]=c[1]/a[1];for(i=2;i<=n-1;i++){a[i]=a[i]-b[i]某c[i-1];c[i]=c[i]/a[i];}q某[1]=d某[1]/a[1];qy[1]=dy[1]/a[1];for(i=2;i<=n-1;i++){q某[i]=(d某[i]-b[i]某q某[i-1])/a[i]; qy[i]=(dy[i]-b[i]某qy[i-1])/a[i];}p某[n-1]=q某[n-1];qy[n-1]=qy[n-1];for(i=n-2;i>=1;i--){p某[i]=q某[i]-c[i]某p某[i+1];py[i]=qy[i]-c[i]某py[i+1];}pO:for(i=0;i<=n-1;i++){b某3=(3某(某[i+1]-某[i])/tt[i+1]-2某p某[i]-p某[i+1])/tt[i+1];b某4=((2某(某[i]-某[i+1])/tt[i+1]+p某[i]+p某[i+1])/tt[i+1])/tt[i+1];by3=(3某(y[i+1]-y[i])/tt[i+1]-2某py[i]-py[i+1])/tt[i+1];by4=((2某(y[i]-y[i+1])/tt[i+1]+py[i]+py[i+1])/tt[i+1])/tt[i+1];t=0;do{t=t+e;c某=某[i]+(p某[i]+(b某3+b某4某t)某t)某t;cy=y[i]+(py[i]+(by3+by4某t)某t)某t;lineto(c某,cy);}while(t<tt[i+1]);}getch();cloegraph();}某[i]i型值点某坐标y[i]i型值点Y坐标a[i]初始赋值方程组系数矩阵mi，i的值b[i]初始赋值方程组系数矩阵mi，i+1的值c[i]初始赋值方程组系数矩阵mi，i-1的值p某[i]i型值点导函数某值py[i]i型值点导函数Y值q某[i],qy[i],d某[i],dy[i]均为解方程组中的各项系数tt[i]第i段参数范围(型值点i-1到型值点i的距离)b某,y3每段函数中二次项的系数b某,y4每段函数中三次项的系数c某每段函数中各点的某值cy每段函数中各点的Y值e作图时每段‘以直代曲’中的参数增量其中方程矩阵形式为：m1,1m1,2p1C1m2,1m2,2m2,3p2C2m3,2m3,3m3,4p3C3............=...mn-1,n-2mn-1,n-1mn-1,npn-1Cn-1mn,n-1mn,npnCn将系数矩阵改写为:其中方程矩阵形式为：(变量与代码变量不对应,如a2不等于a[2]) m1,1m1,2l11u1m2,1m2,2m2,3a2l21u2m3,2m3,3m3,4a3l31u3.........=.........mn-1,n-2mn-1,n-1mn-1,nan-1ln-11un-1mn,n-1mn,nanln1pi对应代码中p某[i],py[i]的值Ci对应代码中q某[i],qy[i]的值附件二运行结果文件。

关于三次样条插值函数的学习报告

关于三次样条插值函数的学习报告三次样条插值函数是一种广泛应用于数值分析领域的插值方法，用于逼近一组已知数据点构成的函数。

在这篇学习报告中，我将介绍三次样条插值函数的定义、原理、应用及其优缺点，并通过实际例子说明其如何在实际问题中使用。

一、三次样条插值函数的定义三次样条插值函数是指用分段三次多项式对一组已知数据点进行插值的方法。

具体来说，对于已知数据点$(x_0,y_0),(x_1,y_1),...,(x_n,y_n)$，三次样条插值函数会在每相邻两个数据点之间构造一个三次多项式，使得这些多项式在相应的数据点上满足插值条件，并且在相邻两个多项式之间满足一定的连续性条件。

二、三次样条插值函数的原理三次样条插值函数的原理是利用三次多项式在每个数据点上的取值和导数值来确定三次多项式的系数，从而构造出满足插值条件和连续性条件的插值函数。

具体来说，对于每个相邻的数据点$(x_i,y_i),(x_{i+1},y_{i+1})$，我们可以构造一个三次多项式$S_i(x)$，满足以下条件：1.$S_i(x_i)=y_i$，$S_i(x_{i+1})=y_{i+1}$，即在数据点上满足插值条件；2.$S_i'(x_{i+1})=S_{i+1}'(x_{i+1})$，$S_i''(x_{i+1})=S_{i+1}''(x_{i+1})$，即在数据点上满足连续性条件。

通过求解上述条件，可以得到每个相邻数据点之间的三次多项式$S_i(x)$，从而得到整个插值函数。

三、三次样条插值函数的应用三次样条插值函数在数值分析领域有广泛的应用，尤其在曲线拟合、数据逼近等问题中起到重要作用。

例如，当我们需要根据已知的离散数据点绘制平滑的曲线图形时，可以使用三次样条插值函数来进行插值，从而得到更加连续和光滑的曲线。

另外，在信号处理、图像处理等领域也常常会用到三次样条插值函数。

例如，在数字图像处理中，我们需要对像素点进行插值以得到更高分辨率的图像，三次样条插值函数可以很好地满足这个需求，使图像更加清晰和真实。

三次样条插值报告

三次样条插值多项式实验的目的及意义：为了取得理想结果：在不增加更多的插值条件下，能够求得一个插值多项式，既有良好的逼近效果，又有好的光滑性，引进三次样条插值多项式。

如果已知函数y=f(x)在节点a=x0<x1<…<xn=b 处的函数值和导数值：()i i x f y =，i=0,1,2,…,n如果S(x)满足条件：1. S(x)是一个分段的三次多项式且()i i y x S =；2. S(x)在[a,b]具有二阶连续导数。

则称S(x)是三次样条插值函数。

S(x)的具体形式为：()()()()⎪⎪⎩⎪⎪⎨⎧∈∈∈=-]12,121,01,[,...............][,][,n n n x x x x s x x x x s x x x x s x s其中()x S i 在[]i i x x ,1-上是三次多项式()iiiiid x c x b x a x S +++=23由插值条件()ii y x S =，i=0,1,2,…,n ，得n+1个条件。

边界条件一：()()nn y x S y x S '',''00== 边界条件二：()()nn y x S y x S '''',''''00==数学公式：()()2211133[2]()[2()]()i i i i i i i i i i ih x x x x h x x x x H x y y h h ---+-----=++2211122()()()()i i i i i i i ix x x x x x x x m m h h -------+ 算法描述：Step1:输入未知数X 及(xi,yi),i=0,1,…,n ； Step2:计算步长H[i]; Step3:计算[][]()⎪⎪⎪⎩⎪⎪⎪⎨⎧+=-=+=-+++i i i i i i i ii i i i i x x f x x f u g u hh h ,,311111λλλStep4:根据边界条件，求解相应的方程得到m0,m1,…, mn Step5:判断X 属于[]i i x x ,1-，i=1,2,…,n 中的哪一个Step6:计算()x s y i i ≈Step7:输出y. 程序原代码如下： #include "stdio.h" #define N 5 void main() { int i,k; float X,s,y0,yn;float a[N][N+1],h[N],u[N],v[N],g[N],m[N],p[N],q[N],w[N];printf("please input X:"); //X 为未知数的大小scanf("%f",&X);printf("please input x:"); //输入x的大小for(i=0;i<N;i++)scanf("%f",&a[i][0]);printf("please input y:"); //输入y的大小for(i=0;i<N;i++)scanf("%f",&a[i][1]);for(i=1;i<N;i++)h[i]=a[i][0]-a[i-1][0]; //计算步长for(i=1;i<N;i++){v[i]=h[i+1]/(h[i]+h[i+1]);u[i]=1-v[i];g[i]=3*u[i]*(a[i+1][1]-a[i][1])/h[i+1]+3*v[i]*(a[i][1]-a[i-1][1])/h[i]; }printf("\t(1)已知边界条件1\n");printf("\t(2)已知边界条件2\n");printf("请选择边界条件序号:");scanf("%d",&k);if(k==1){printf("请输入y0和yn的一阶导:"); //输入边界条件一scanf("%f%f",&m[0],&m[N-1]);p[0]=0; //用追赶法求解m[N]q[0]=0;g[1]=g[1]-v[1]*m[0];g[N-2]=g[N-2]-u[N-2]*m[N-1];for(i=1;i<N;i++){w[i]=2-u[i]*p[i-1];p[i]=v[i]/w[i];q[i]=(g[i]-u[i]*q[i-1])/w[i];}m[N-2]=q[N-2];for(i=N-3;i>0;i--)m[i]=q[i]-p[i]*m[i+1];printf("输出m[i]的值:\n");for(i=0;i<N;i++)printf("%f\n",m[i]);for(i=1;i<N;i++) //计算最终结果if(X>a[i-1][0]&&X<a[i][0])s=(h[i]+2*(X-a[i-1][0]))*(X-a[i][0])*(X-a[i][0])*a[i-1][1]/(h[i]*h[i]*h[i]) +(h[i]-2*(X-a[i][0]))*(X-a[i-1][0])*(X-a[i-1][0])*a[i][1]/(h[i]*h[i]*h[i])+ (X-a[i-1][0])*(X-a[i][0])*(X-a[i][0])*m[i-1]/(h[i]*h[i])+(X-a[i-1][0])*(X-a[i-1][0])*(X-a[i][0])*m[i]/(h[i]*h[i]);printf("s(%f)=%f\n",X,s);}if(k==2){printf("请输入y0和yn的二阶导:"); //输入边界条件二scanf("%f%f",&y0,&yn);g[0]=3*(a[1][1]-a[0][1])/h[1]-h[1]*y0/2;g[N-1]=3*(a[N-1][1]-a[N-2][1])/h[N-1]+h[N-1]*yn/2;q[0]=g[0];u[0]=1;v[N-1]=1;w[0]=2;for(i=1;i<N;i++){w[i]=2-v[i]*u[i-1]/w[i-1];q[i]=g[i]-v[i]*q[i-1]/w[i-1];}m[N-1]=q[N-1]/w[N-1];for(i=N-2;i>=0;i--)m[i]=(q[i]-u[i]*m[i+1])/w[i];printf("输出m[i]的值:\n");for(i=0;i<N;i++)printf("%f\n",m[i]);for(i=1;i<N;i++)if(X>=a[i-1][0]&&X<=a[i][0])s=(h[i]+2*(X-a[i-1][0]))*(X-a[i][0])*(X-a[i][0])*a[i-1][1]/(h[i]*h[i]*h[i]) +(h[i]-2*(X-a[i][0]))*(X-a[i-1][0])*(X-a[i-1][0])*a[i][1]/(h[i]*h[i]*h[i])+ (X-a[i-1][0])*(X-a[i][0])*(X-a[i][0])*m[i-1]/(h[i]*h[i])+(X-a[i-1][0])*(X-a[i-1][0])*(X-a[i][0])*m[i]/(h[i]*h[i]);printf("s(%f)=%f\n",X,s);}}数值计算：已知y=f(x)的如下数值求三次样条插值函数S（x），满足条件1.s’(0)=0,s’(4)=482.s’’(0)=0,s’’(4)=24Please input X:2.5Please input x:0 1 2 3 4 Please input y:-8 -7 0 19 56(1)已知边界条件1(2)已知边界条件2请选择边界条件的序号:1请输入y0和yn的一阶导:0 48 0.0000003.00000012.00000027.00000048.000000s(2.500000)=7.625000press any key tocontinue请选择边界条件的序号:2请输入y0和yn的二阶导:0 24 -0.0000003.00000011.99999927.00000248.0000007.625000press any key tocontinue s(2.500000)=7.625000 对计算结果进行评价分析：()()443845h M x S x f ≤-三次样条插值函数与三次Hermite 插值函数相比，不仅光滑度有提高，而且要求求解时还不需要增加内节点处的导数值，因此比较实用。

数值分析实验报告-插值、三次样条

实验目的：通过编程实现牛顿插值方法和三次样条方法，加深对多项式插值的理解。

应用所编程序解决实际算例。

实验要求：1．认真分析问题，深刻理解相关理论知识并能熟练应用； 2．编写相关程序并进行实验； 3．调试程序，得到最终结果； 4．分析解释实验结果； 5．按照要求完成实验报告。

实验原理：详见《数值分析第5版》第二章相关内容。

实验内容：（1）牛顿插值多项式1.1 当n=10时：在Matlab 下编写代码完成计算和画图。

结果如下：代码：clear allclcx1=-1:0.2:1;y1=1./(1+25.*x1.^2);n=length(x1);f=y1(:);for j=2:nfor i=n:-1:jf(i)=(f(i)-f(i-1))/(x1(i)-x1(i-j+1));endendsyms F x p;F(1)=1;p(1)=y1(1);for i=2:nF(i)=F(i-1)*(x-x1(i-1));p(i)=f(i)*F(i);endsyms PP=sum(p);P10=vpa(expand(P),5);x0=-1:0.001:1;y0=subs(P,x,x0);y2=subs(1/(1+25*x^2),x,x0);plot(x0,y0,x0,y2)grid onxlabel('x')ylabel('y')P10即我们所求的牛顿插值多项式，其结果为：P10（x）=-220.94*x^10+494.91*x^8-9.5065e-14*x^7-381.43*x^6-8.504e-14*x^5+123.36* x^4+2.0202e-14*x^3-16.855*x^2-6.6594e-16*x+1.0并且这里也能得到该牛顿插值多项式的在[-1，1]上的图形，并和原函数进行对比（见Fig.1）。

三次样条插值报告(三转角)

数值逼近实习报告题目：三次样条插值II班级：计算111班学号：3110811009姓名：刘艳平指导老师：秦新强2013-3题目：三次样条插值（三转角算法）姓名：刘艳平学号：3110811009 一、目的意义以上我们讨论的分段线性插值，逼近程度虽然好，但光滑性差，分段三次Hermite 插值，逼近程度好，光滑性也有所提高，但也增加了更多的条件，不太实用。

接下来我们研究的三次样条插值多项式，这就结合二者的优点，即逼近程度虽然好，光滑性强，不需要增加太多的条件，很实用。

当已知区间两个端点的导数值时，既满足边界条件1，用三弯矩算法来求解，相对简单，且很准确。

二、数学公式⎪⎪⎩⎪⎪⎨⎧∈∈∈-],[),(............],[),(],[),(1212101n n nx x x x s x x x x s x x x x s 具体：iii i i ii i iii i i i ii i i i m h x x x x m h x x x x y h x x x x h y h x x x x h x s 22112213211321))(())(())]((2[))]((2[)(--------+--+---+--+=其中：n i x x x x x h i i i i i ,...,2,1],,[,11=∈-=-- i i i i y x s y x s ==--)(,)(11需要求出：ii i i m x s m xs ='='--)(,)(11三、计算流程 Step 1:输入节点n x x x ,...,,10，函数值ny y y ,...,,10,边界条件及x.Step 2:计算⎪⎪⎪⎩⎪⎪⎪⎨⎧+=+=-=+=-++++]),[],[(3111111i i i i i i i ii i ii i i i i x x f x x f g h h h h h h λμλμλStep 3：根据边界条件，求解相应的方程得到nm m m ,...,,10.Step 4：判断x 属于],[1i i x x -，i=1,2,...,n中的那一个.Step 5：计算)(x s y i ≈.Step 6：输出y. 四、代码#include <iostream.h> //三次样条三转角边界条件1 #include "process.h" #define N 2 #define M 4double a[N],d[N],c[N],p[N],q[N],b[N],xx[N]; void p_q() { int k; if(d[0]==0) {cout<<"Method failed"<<endl;exit(0);}p[0]=d[0];q[0]=c[0]/p[0];for(k=1;k<N-1;k++){p[k]=d[k]-a[k]*q[k-1];if(p[k]==0){cout<<"Method failed"<<endl; exit(0);}q[k]=c[k]/p[k];}p[N-1]=d[N-1]-a[N-1]*q[N-2]; if(p[N-1]==0){cout<<"Method failed"<<endl; exit(0);}}void jisuan(){int k;double y[N];y[0]=b[0]/p[0];for(k=1;k<=N-1;k++)y[k]=(b[k]-a[k]*y[k-1])/p[k];xx[N-1]=y[N-1];for(k=N-2;k>=0;k--)xx[k]=y[k]-q[k]*xx[k+1];for(k=0;k<N;k++)cout<<"m的值为:"<<xx[k]<<'\t'<<endl;}void main(){double X[M],Y[M],h[M],v[M],u[M],m[M],g[M],mm[M]; int i;double x,y1,y2,w1,w2,w3,w4,w,hh;cout<<"输入已知点的X坐标:"<<endl;for(i=0;i<M;i++)cin>>X[i];cout<<"输入已知点的Y坐标:"<<endl;for(i=0;i<M;i++)cin>>Y[i];cout<<"输入边界条件:"<<endl;cin>>y1>>y2;cout<<"输入需要计算的x值:"<<endl;cin>>x;for(i=1;i<M;i++){h[i]=X[i]-X[i-1];}for(i=1;i<M-1;i++){v[i]=h[i+1]/(h[i]+h[i+1]);u[i]=h[i]/(h[i]+h[i+1]);g[i]=3*((v[i]*(Y[i+1]-Y[i])/(X[i+1]-X[i]))+(u[i]*(Y[i]-Y[i-1])/(X[i]-X[i-1])));}a[0]=0;b[0]=g[1]-v[1]*y1;for(i=1;i<N;i++)a[i]=v[i+1];for(i=0;i<N;i++)d[i]=2;for(i=0;i<N;i++)c[i]=u[i+1];for(i=1;i<N;i++)b[i]=g[i+1];p_q();jisuan();mm[0]=y1;mm[1]=xx[0];mm[2]=xx[1];mm[3]=y2;for(i=1;i<M;i++){hh=X[i]-X[i-1];if(x>=X[i-1]){w1=((hh+(x-X[i-1])*2)*(x-X[i])*(x-X[i])*Y[i-1])/(hh*hh*hh);w2=((hh-(x-X[i])*2)*(x-X[i-1])*(x-X[i-1])*Y[i])/(hh*hh*hh);w3=((x-X[i-1])*(x-X[i])*(x-X[i])*mm[i-1])/(hh*hh);w4=((x-X[i])*(x-X[i-1])*(x-X[i-1])*mm[i])/(hh*hh);w=w1+w2+w3+w4;}}cout<<"x 点的结果为"<<w<<endl; }五、数值实例已知函数y=f(x),的如下数据，试求其在区间]3,0[上的三次样条插值函数S(x). X 0 1 2 3 Y0 1 0 1 y '1解：在此题中，解得：151,15421=-=m m其中：0,130==m m ，一起带入上边公式即可得到所求的三次样条函数如下⎪⎪⎪⎩⎪⎪⎪⎨⎧∈--+---=∈--+------=∈---+--=]3,2[,)3)(2(151)2)](3(21[)(]2,0[),2()1(151)2)(1(154)2)](1(21[)(]1,0[),1(154)1()]1(21[)(22322222221x x x x x x s x x x x x x x x s x x x x x x x x s 假设，求出x=1.5处的函数值为：0.458333 程序运行结果截图如下：六、对计算机结果进行分析评价三次样条插值函数与三次Hermite插值函数相比，不仅光滑度有提高，而且要求求解时还不需要增加内节点处的导数值，因此比较实用。

第二类边界条件三次样条插值实验报告

数值计算实验—实验报告2一、实验项目：第二类边界条件三次样条插值二、实验目的和要求a.通过本实验深入地理解三次样条插值多项式的基本原理b.通过数值算例更好的领会三次样条插值多项式具有较高的准确性三、实验内容1.用调试好的程序解决如下问题：点中点处的函数值，并将计算结果与sinx在相应点的数值相比较。

n=8;p1=0.4794;pn=0.9463;u=[0.6,0.8,1.0,1.2,1.4,1.6,1.8];p=7;x=[0.5,0.7,0.9,1.1,1.3,1.5,1.7,1.9];y=[0.4794,0.6442,0.7833,0.8912,0.9636,0.9975,0.9917,0.9463];for i=1:n-1h(i)=x(i+1)-x(i);enda2(1)=1;g(1)=3*(y(2)-y(1))/h(1)-p1*h(1)/2;for k=2:n-1a1(k-1)=h(k)/(h(k)+h(k-1));a2(k)=h(k-1)/(h(k)+h(k-1));g(k)=3*a2(k)*(y(k+1)-y(k))/h(k)+3*a1(k-1)*(y(k)-y(k-1))/h(k-1); enda1(n-1)=1;g(n)=3*(y(n)-y(n-1))/h(n-1)+pn * h(n-1)/2;%追赶法求三转角方程b1(1)=2;m(1)=g(1)/2;b2(1)=a2(1)/b1(1);for i=2:nb1(i)=2-a1(i-1)*b2(i-1);if(i~=n)b2(i)=a2(i)/b1(i);endm(i)=(g(i)-a1(i-1)*m(i-1))/b1(i);endfor i=n-1:-1:1m(i)=m(i)-b2(i)*m(i+1);endfor j=1:pfor i=1:nif((u(j)>=x(i))&&(u(j)<x(i+1)))k=i;break;endends(j)=0;s(j)=s(j)+(h(k)+2*(u(j)-x(k)))*(u(j)-x(k+1))^2*y(k)/(h(k))^3;s(j)=s(j)+(h(k)-2*(u(j)-x(k+1)))*(u(j)-x(k))^2*y(k+1)/(h(k))^3;s(j)=s(j)+(u(j)-x(k))*(u(j)-x(k+1))^2*m(k)/(h(k))^2;s(j)=s(j)+(u(j)-x(k+1))*(u(j)-x(k))^2*m(k+1)/(h(k))^2;end（2）.运行结果3. 根据Lagrange插值多项式基本原理编制程序，并计算下面的数值算例：=-5+kh,其中h=10/n,n=10,20,40.给定函数f(x)=1/(1+x^2)（-5≤x≤5），取等距节点xk边界条件为S''(x0)=f''(x0),S''(x n)=f''(x n).用上述算法计算S10(x),S20(x), S40(x),并与函数f(x)以及10次Lagrange插值多项式L10(x)在给定点处的函数值进行比较。

实验一--三次样条插值的三弯矩法

实验一三次样条插值的三弯矩法
一、实验目的
用三次样条插值的三弯矩法,编制第一与第二种边界条件的程序。已知数据， , 及边界条件，编程计算的三次样条插值函数 .具体要求为：输出用追赶法解出的弯矩向量及的值;画出的图形,图形中描出插值点及分别用‘o’和‘＊’标记.
二、实验原理
（1）由，按公式
dS（2。841)=—0。068625
dS（3。67）=-0.034911
d2S(—4.166）=0。016506
d2S（—3.233)=0.04067
d2S(-2.495)=0.093623
d2S(-1.888）=0。20668
d2S(-1.469）=0.34754
d2S（-1.193)=0.46113
用追赶法求得的弯矩量为
要计算的20个节点出的值为
M0=0。0084
M1=0。019204
M2=0。047836
M3=0.10891
M4=0.2393
M5=0.38381
M6=0。48683
M7=0。50559
M8=0。081512
M9=—0.55789
M10=—2。2251
M11=-0.64145
M12=0。097766,， Nhomakorabea求出。
（2）由边界条件确定和。
（3)用追赶法解线性方程组
求出
（4）由
确定 ,得到S（x）。
三、实验结果
所用数据：
x=[—5,—3.8884，—3.0142,—2。3216，—1.7433，-1。3771,-1.1318,—0.86821,—0.58598，—0。4314，-0。0293，0。4314，0.58598，0.86821,1.1318，1。3771,1。7433,2.1242，2.3216,3.0142，3。8884,5]；

代数插值实验-三次样条曲线

实验7 代数插值实验-（2）实验类型：●验证性实验 ○综合性实验 ○设计性实验实验目的：进一步熟练掌握spline(样条)插值算法，提高编程能力和解决插值问题的实践技能。

实验内容：编写程序用计算机求解三次压紧样条曲线，经过点（0，0.0），（1，0.5），（2，2.0）和（3，1.5），而且一阶导数边界条件S ’(0)=0.2和S ’(3)=-1实验原理若已知N+1个点的{}N k k k y x 0),(=及其一阶导数的边界条件S ’(a)=0d 和S ’(b)=N d ,则存在唯一的三次样条曲线。

构造并求解下列线性方程组01112100111122211113(2)3('())22()(2,3,,2)3(2)3('())2k k k k k k k k N N N N N N N N h h m h m u d S x h m h h m h m u k N h m h h m u S x d ---+-------⎧++=--⎪⎪+++==-⎨⎪⎪++=--⎩ 构造分段函数S(x):23,0,1,2,31()()()()()[,]k k k k k k k k k k S x S x s s x x s x x s x x x x x +==+-+-+-∈ k k y s =0, ; 6)2(11,++-=k k k k k m m h d s ; 22,k k m s = ; kk k k h m m s 613,-=+ 实验说明需要建立两个数组X ，Y 分别存放（1x ，…，)n x ，（1y ，…，)n y ，建立两个变量dxo 、 dxn 分别存放左右端点的一阶导数 dxo = S'(x0) , dxn = S'(xn)，上述方程组中10,1,,1k k k h x x k N +=-=- 其中；1()/0,1,,1k k k k d y y h k N +=-=- 16(),1,2,,1k k k u d d k N -=-=- 实验步骤1 要求上机实验前先编写出程序代码2 编辑录入程序3 调试程序并记录调试过程中出现的问题及修改程序的过程4 经反复调试后，运行程序并验证程序运行是否正确。

LAB03_三次样条插值实验

LAB03_三次样条插值实验实验内容1.用三次样条插值的三弯矩法,编制第一与第二种边界条件的程序.已知数据如下：求)(x f 的三次样条插值函数)(x S 满足： (1)自然边界条件''(0.2)''(1.0)0;S S ==(2)第一种边界条件.55741.1)0.1(',20271.0)2.0('==S S 要求输出用追赶法解出的弯矩向量(0M ,1M ,2M ,3M ,4M )及)1.02.0(i S +(i=0,1,2,3,4,5,6,7,8)的值.并画出)(x S y =的图形.实验原理简述我们先利用一阶导数)(x s ''在节点i x （i=1,2,…n-1）上的连性以及边界条件，列出确定二阶导数)(i i x s M ''=（i=0,1,2,…n-1）的线性方程组（三弯矩方程），并由此解出i M 然后用i M 表达是)`(x s 即可i i i i i i g M M M =+++-112λμ（三弯矩方程） jj j i h h h +=--11μ，jj j i h h h +=-1λ， ,,,1,0n i =[][][]11111,,6,,6+---+=+-=j j j jj jj j j i x x xf h h x x f x x fg .其中有（1+n ）个未知数n M M M ,...,10，而方程只有（n-1）个，当满足第一种边界条件时，可的另两个方程[]()010010,62f x x f h M M '-=+， []()n n n n n n x x f f h M M ,62111----'=+如果令[]()[]()n n n n n n x x f f h g f x x f h d ,6,1,,6,111010000---'=='-==μλ，将上述方程综合后的一下矩阵形式：222211110nn n μλμλμλ--nn M M M M 110- =n n g g g g 110-可以证明此方程组满足追赶法的条件，我们用追赶法可得M 的)1(+n 值，将其带入公式即得)(x s对第二种边界条件，直接的端点方程n n f M f M ''=''=,00并且令n n n f d f d ''=''===2,2,0000μλ，则又得三弯矩方程同理即可求得解。

关于三次样条插值函数的学习报告

关于三次样条插值函数的学习报告三次样条插值函数是一种常用的插值方法，它利用多项式函数的特性来逼近一组数据点，并且具有较高的精度和平滑性。

本学习报告将对三次样条插值函数进行详细介绍，并探讨其在实际应用中的优缺点。

首先，三次样条插值函数的数学表达式可以表示为：S(x) = \sum_{i=0}^n {a_i(x-x_i)^3 + b_i(x-x_i)^2 + c_i(x-x_i) + d_i}其中，n是数据点的数量，a_i、b_i、c_i、d_i是通过求解一系列方程得到的系数。

这些方程的目标是使得插值函数在每个数据点之间的一阶和二阶导数连续。

对于每个区间[x_i,x_{i+1}]，我们可以得到一个关于未知系数a_i、b_i、c_i、d_i的线性方程组。

通过求解这些方程组，我们可以确定插值函数在每个区间的系数。

最终，我们得到一个全局的三次插值函数，它在整个插值区域内都具有较高的拟合精度。

三次样条插值函数的优点之一是它可以通过调整插值区间的数量来灵活控制插值的精度。

当插值区间数量增加时，插值函数与原始数据点之间的误差会减小，从而获得更精确的插值结果。

另外，三次样条插值函数还具有较好的平滑性，能够克服其他插值方法中可能出现的震荡现象，使得插值函数更加平滑。

然而，三次样条插值函数也存在一些缺点。

首先，它对于数据点分布较为密集的情况下，有时会出现振荡现象，导致插值函数不够平滑。

其次，三次样条插值函数在插值区间的两个端点附近可能无法很好地逼近原始数据。

这是因为在每个区间的端点处，插值函数需要满足特定的边界条件，通常是一阶或二阶导数为零。

这种约束可能导致插值函数在端点处的拟合程度较低。

为了解决上述问题，可以使用更高阶的样条插值函数，如五次样条插值函数。

五次样条插值函数通过增加插值函数的阶数，以获取更高的拟合精度和平滑性。

此外，还可以尝试使用其他插值方法，如非均匀插值、基函数插值等，以应对不同的插值需求。

总结来说，三次样条插值函数是一种精度较高且平滑的插值方法，具有较好的数学属性。

三次样条插值实验报告

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

画图：
x1=0:.01:1;y1=polyval(S1(1,:),x1-X(1)); x2=1:.01:2;y2=polyval(S1(2,:),x2-X(2)); x3=2:.01:3;y3=polyval(S1(3,:),x3-X(3)); x4=3:.01:4;y4=polyval(S1(4,:),x4-X(4)); x5=4:.01:5;y5=polyval(S1(5,:),x5-X(5)); x6=5:.01:6;y6=polyval(S1(6,:),x6-X(6)); >> plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5,x6,y6,X,Y,'.') >> hold on >> x1=0:.01:1;y1=polyval(S2(1,:),x1-X(1)); x2=1:.01:2;y2=polyval(S2(2,:),x2-X(2)); x3=2:.01:3;y3=polyval(S2(3,:),x3-X(3)); x4=3:.01:4;y4=polyval(S2(4,:),x4-X(4)); x5=4:.01:5;y5=polyval(S2(5,:),x5-X(5)); x6=5:.01:6;y6=polyval(S2(6,:),x6-X(6)); >> plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5,x6,y6,X,Y,'.') >> hold on >> x1=0:.01:1;y1=polyval(S3(1,:),x1-X(1)); x2=1:.01:2;y2=polyval(S3(2,:),x2-X(2)); x3=2:.01:3;y3=polyval(S3(3,:),x3-X(3)); x4=3:.01:4;y4=polyval(S3(4,:),x4-X(4)); x5=4:.01:5;y5=polyval(S3(5,:),x5-X(5)); x6=5:.01:6;y6=polyval(S3(6,:),x6-X(6)); >> plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5,x6,y6,X,Y,'.')
csfit2:第二类边界条件
function S=csfit2(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); C=H(2:N); U=6*diff(D); U(1)=U(1)-dx0; U(N-1)=U(N-1)-dxn;
三次样条插值实验报告
专业班级学号姓名
1、实验内容和要求
1、阅读上面的文字和程序，试运行，检验程序和上面叙述的正确性。
2、阅读上面的 MATLAB 程序；查资料，了解各 MATLAB 语句及命令。
3、画程序流程图，理解并描述算法。 4、修改上面的程序，能根据给定数据点，求（1）自然样条插值，边界 ,= = （2）第二类边界条件，是确定的。 5、使用上面的程序，根据数据点（0，1），（1，0），（2， 0），（3，1），（4，2），（5，2），和（6，1），求三种不同的三次样条插值，其中， =-1.8，在同一坐标系中，画出这 3 个三次样条插值和这些数据点。
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
（3）
将
代入上（3），并使用
，可分别得到包含
的方程：
（4）
求解这两个方程，求出
，而且将这些
值代入方程（3）中，可得到三次多项式方程：
（5）
表达式5可以简化成只包含未知系数
的
形式。
为求这些值，必须使用（5）式的导数，即
（6）
在处计算（6），并简化结果可得到
，其中
（7）
同理，在式（6）中用
并
计算在处的解，可得
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); End
6、写实验报告（实验内容+算法描述+程序+写成分段函数的结果描述+截图）。
2、算法说明
定义：设有N+1个点
，其中
。如果存在
N个三次多项式，系数为
，满足如下性质：（1）（2）（3）
则称函数为三次样条函数。因为是分段三次多项式，所以
区间上是分段线性的。
用上式，可得
在
（1）代入
（2）将（2）式积分两次，会引入两个积分常数，并得到
4、实验结果
第一类边界条件：
第二类边界条件：
自然边界条件：
图像：
5、说明与分析
实验结果显示，对同一组数据，不同边界条件，会导致形成曲线有一定差距，但总体趋势不变
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)-dx0)/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
（8）
ቤተ መጻሕፍቲ ባይዱ
利用节点处一阶导函数连续及方程（6）、
（7），可得到包含
的重要关系式
（9）
其中
方程组（9）中的未知数是要求的值
，其
他的项可通过数据点集
进行简单的计算
得到的常量。因此方程（9）是包含N+1个未知数，
具有N-1个线性方程的不定方程组。所以需要另外
两个方程组才能求解，即边界条件。
如果已知
，则
（10）
M(1)=dx0; M(N+1)=dxn; 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
（11）
根据（9）（10）求出公式计算的样条系数
，，，，，
后，可利用下面的。
3、源程序
csfit1:第一类边界条件
function S=csfit1(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); C=H(2:N); U=6*diff(D);