资源描述：
《matlab实现数值分析插值与积分》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、Matlab实现数值分析插值及积分摘要：数值分析(numericalanalysis)是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科，是数学的一个分支，它以数字计算机求解数学问题的理论和方法为研究对象。在实际生产实践中，常常将实际问题转化为数学模型来解决，这个过程就是数学建模。学习数值分析这门课程可以让我们学到很多的数学建模方法。分别运用matlab数学软件编程来解决插值问题和数值积分问题。题目中的要求是计算差值和积分，对于问题一，可以分别利用朗格朗日插值公式，牛顿插值公式，埃特金逐次线性插值公式来进行编程求解，具体matlab
2、代码见正文。编程求解出来的结果为：=+。其中Aitken插值计算的结果图如下：对于问题二，可以分别利用复化梯形公式，复化的辛卜生公式，复化的柯特斯公式编写程序来进行求解，具体matlab代码见正文。编程求解出来的结果为：0.6932其中复化梯形公式计算的结果图如下：问题重述问题一：已知列表函数表格101234121782257分别用拉格朗日，牛顿，埃特金插值方法计算。问题二：用复化的梯形公式，复化的辛卜生公式，复化的柯特斯公式计算积分，使精度小于5。问题解决问题一：插值方法对于问题一，用三种差值方法：拉格朗日，牛顿，埃特金差值方法来解决。一、拉
3、格朗日插值法：拉格朗日插值多项式如下： 首先构造个插值节点上的插值基函数，对任一点所对应的插值基函数，由于在所有取零值，因此有因子。又因是一个次数不超过的多项式，所以只可能相差一个常数因子，固可表示成：利用得：于是因此满足的插值多项式可表示为：从而次拉格朗日插值多项式为：matlab编程：编程思想：主要从上述朗格朗日公式入手：依靠循环，运用poly（）函数和conv（）函数表示拉格朗日公式，其中的poly（i）函数表示以i作为根的多项式的系数，例如poly（1）表示x-1的系数，输出为1-1，而poly（poly（1））表示（x-1）*（x-1
4、）=x^2-2*x+1的系数，输出为1-21；而conv（）表示多项式系数乘积的结果，例如conv（poly（1），poly（1））输出为1-21；所以程序最后结果为x^n+x^n-1+……+x^2+x+1（n的值据结果的长度为准）的对应系数。在命令窗口输入editlagran来建立lagran.m文件，文件中的程序如下：function[c,l]=lagran(x,y)w=length(x);n=w-1;l=zeros(w,w);fork=1:n+1v=1;forj=1:n+1ifk~=jv=conv(v,poly(x(j)))/(x(k)-
5、x(j));endendl(k,:)=v;endc=y*l;输入：>>x=[01234];>>y=[121782257];>>lagran(x,y)运行结果为ans=1.0000-0.0000-0.000001.0000结果为：=+。如图表1：图表1二．牛顿插值法newton插值多项式的表达式如下：其中每一项的系数ci的表达式如下：即为f(x)在点处的i阶差商，（，），由差商的性质可知：matlab编程：编程思想：主要从上述牛顿插值公式入手：依靠循环，运用poly（）函数和conv（）函数表示拉格朗日公式，其中的poly（i）函数表示以i作为根
6、的多项式的系数，例如poly（1）表示x-1的系数，输出为1-1，而poly（poly（1））表示（x-1）*（x-1）=x^2-2*x+1的系数，输出为1-21；而conv（）表示多项式系数乘积的结果，例如conv（poly（1），poly（1））输出为1-21；所以程序最后结果为x^n+x^n-1+……+x^2+x+1（n的值据结果的长度为准）的对应系数。在命令窗口输入editnowpoly来建立newpoly.m文件，文件中的程序如下：function[c,d]=newpoly(x,y)n=length(x);d=zeros(n,n);d
7、(:,1)=y';forj=2:nfork=j:nd(k,j)=(d(k,j-1)-d(k-1,j-1))/(x(k)-x(k-j+1));endendc=d(n,n);fork=(n-1):-1:1c=conv(c,poly(x(k)));m=length(c);c(m)=c(m)+d(k,k);end输入：>>x=[01234];>>y=[121782257];>>newpoly(x,y)运行结果为ans=10001所以=+。如图表2：图表2三．埃特金插值法：Aitken插值公式如下：递推表达式为：=+当n=1时，=+当n=2时，=+其中的
8、带入递推表达式求得。由此递推下去，最终得到的结果。matlab编程：编程思想：埃特金插值多项式又称作Aitken逐次线性插值多项式，根据公式的特点，可
