信息计算科学实验报告

1、信息计算:岸实验报土1=1实验报告专业：信息与计算科学年级：大三 班级：ap08102 学号：ap0810227姓名：庞锦芬、实验目的1、了解lagrange插值法的基本原理和方法；2、了解多项式拟合的基本原理和方法；3、了解数值积分的基本原理和方法；实验题目：实验三插值法与拟合实验1、插值效果比较：将区间【-5,5】10等分，对下列函数分别计算插值节点Xk的值，进行不同类型的插值，作出插值函数的图形并与y f(x)的图形进行比较:f(x) 11x2 ；f (x) arctanx ；f(x)2x1x4 .Xiyi分别对上述数据作三次多项式和五次多项式拟合，并求平方误差，作出离散函数(xi, %

2、)和拟合函数的(1)做拉格朗日插值；2、拟合多项实验： 给定数据如下表所示:图形。实验四数值微积分实验1、复化求积公式计算定积分 ：用复化梯形公式、复化辛普森公式、龙贝格公式求下列定积分，要求绝对误 差为 0.5 108，并将计算结果与精确解进行比较：(1) e42x3ex2dx，3(2) ln62x dx2 x2 3三、实验原理(将实验所涉及的基础理论、算法原理详尽列出。拉个朗日插值原理：经过n 1个点(X。，y°),(X1, %) ,(xn,yn)，构造一个n次多项式，形pn(x) 使得 Pn(xQyk (k 0,1,2 ,n)成立。ykik(x)k 0lk(x)(x xj (x

3、 Xk1)(X xk 1) (x xn)(Xk X0)(Xk Xk 1)(Xk Xk 1) (Xk Xn)0 x Xk1 X Xk其中为插值基函数。拟合多项式原理：假设给定数据点(Xi,yi)(i=0,1，,m)，nPn(X)akX现求一k 0，使得mPn(Xi)i 0yi为所有次数不超过n(nm)的多项式构成的函数类,nkakXik 0yimin当拟合函数为多项式时，称为多项式拟合， 特别地，当n=1时，称为线性拟合或直线拟合。 显然满足式(1)(1)Pn (X)称为最小二乘拟合多项式。为a0,a1， an的多元函数, 极值的必要条件，得Im n kI(akXii 0 k 0因此上述问题即为

4、求ajmJ kXi)aki 0mXiJyi,i 0(3)是关于a。®, an的线性方程组，用矩阵表示为mXii 0mXii 0m2Xii 0mnXii 0mn 1Xii 0aoa1yi)2I(ao,a1an)的极值问题。由多元函数求0,1,n0,1,nmnXii 0m n 1Xii 0anmyii 0mXiYii 0mnXi yii 0m 2nXii 0式(3)或式(4)称为正规方程组或法方程组。可以证明，方程组(4)的系数矩阵是一个对称正定矩阵，故存在唯一解。出ak(k=0,1，n)，从而可得多项式从式(4)中解nPn(X)akXkk 0可以证明，式(5)中的Pn(X)满足式(1)

5、，即Pn(X)为所求的拟合多项式。我们把i 0 称为最小二乘拟合多项式Pn(X)的平方误差，记作llrlPn(X)m2 22Pn(Xi) yii 0mPn(Xi)yi由式可得2 mnm|r2y：ak(Xikyi 0k 0 i 0(硬件和软件)。)四、实验内容(列出实验的实施方案、步骤、数据准备、算法流程图以及可能用到的实验设备实验步骤：1、先编写好matlabM文件，然后在命令窗口编辑程序并运行；2、运行，观察结果；3、根据运行结果进行结果分析。实验三各个实验在 matlab窗口输进的主要程序如下：拉格朗日插值：x=-5:1:5;y1=1./(1+x.A2);y2=ata n( x);y3=x

6、.A2丿(1+x44);L1=malagr(x,y1,x);L2=malagr(x,y2,x);L3=malagr(x,y3,x);plot(x,y1,'r',x,y2,'g',x,y3,'b',x ,L 1,'rp',x ,L2,'gd',x ,L 3,'b*'); xlabel('x');ylabel('y');lege nd('y1','y2','y3','L1','L2','

7、;L3')拟合多项式：作三次多项式拟合的程序：x=0;y=;y1=mafit(x,y,3)作五次次多项式拟合的程序：x= 0 ;y=;y2=mafit(x,y,5)求平方误差，作出离散函数(冷)和拟合函数的图形，程序为：% san ci ni he duo xia ng shi de xi shux= 0 ;y2=;y=mafit(x,y,3);P1=.*X.*X.A.*X+% san ci ni he duo xia ng shi de xi shux= 0 ；y=;y2=mafit(x,y,5);p2=.*x.A5+.*x.A4+.*x.A.*x.A.*x+norm(p1-y)n

8、orm(p2_y)plot(x,y, 's',x,p1,'d',x,p2,'p')xlabel('x');ylabel('y');legend('y', 'p1' , 'p2')% s fangxing; d lingxing; p wujiaox ing实验四复化求积公式计算定积分用复化梯形公式相关的程序和相关的注释如下：(1)% numerical integrate formulation 1(1) f2=diff( '2/3*xA3*exp(xA2)&

9、#39;, 'x',2)% 求对变量x的二阶偏导数%f2=4*x*exp(xA2)+28/3*xA3*exp(xA2)+8/3*xA5*exp(xA2)=exp(xA2)*8/3*x(2*x+1)(x+3)x=2;a=4*x*exp(xA2)+28/3*xA3*exp(xA2)+8/3*xA5*exp(xA2)%求2 的函数值f=i nline('1/12*hA2*+003)', 'h':);%复化梯形余项减去误差h=fzero(f,0)%求满足精度的h值n=abs(1/h);%求满足精度的n值fun=in li ne('2/3.*x.A

10、3*exp(x.A2)')T=matrap(fu n,1,2 ,n)% T=b=exp % exp (4)=(2)f2=diff( '2*x/(xA2-3)', 'x',2)% 求对变量 x 的二阶偏导数% f2 =-12/(xA2-3)A2*x+16*xA3/(xA2-3)A3 x=2;a=-12/(xA2-3)A2*x+16*xA3/(xA2-3)A3%求f2 的函数值f=i nline('1/12*hA2*', 'h')；%复化梯形余项减去误差h=fzero(f,0)%求满足精度的h值n=abs(1/h);%求满足精

11、度的n值fun=in li ne('2*x/(xA2-3)')T=matrap(fu n,2,3 ,n)% T=b=log (6)% log( 6)=用复化辛普公式相关的程序和相关的注释如下：%用复化辛普公式f2=diff( '2/3*xA3*exp(xA2)', 'x' ,4)%求对变量 x 的 4 阶偏导数%f2=80*x*exp(xA2)+200*xA3*exp(xA2)+96*xA5*exp(xA2)+32/3*xA7*exp(xA2) x=4;a=80*x*exp(xA2)+200*xA3*exp(xA2)+96*xA5*exp(xA2

12、)+32/3*xA7*exp(xA2)求f2的函数值f=inline('1/2880牛人4*+012)', 'h' ); %复化辛普余项减去误差h=fzero(f,0)%求满足精度的h值n=abs(1/h);%求满足精度的n值fun=in li ne( '2/3.*x.A3*exp(x.A2)')T=masimp(fu n,1,2, n)% T=b=exp % exp (4)=(2)%用复化辛普公式(2)f2=diff( '2*x/(xA2-3)'，乂 ,4)%求对变量 x 的 4 阶偏导数% f2 =-960/(xA2-3)A4

13、*xA3+240/(xA2-3)A3*x+768*xA5/(xA2-3F5 x=4;a=-960/(xA2-3)A4*xA3+240/(xA2-3F3*x+768*xA5/(xA2-3)A5%求f2 的函数值f=inline('1/12牛人2* ', 'h' ); %复化辛普余项减去误差h=fzero(f,0)%求满足精度的h值n=abs(1/h);%求满足精度的n值fun=inlin e('2*x/(xA2-3)')T=matrap(fu n,2,3 ,n)% T=b=log (6)% log( 6)=用龙贝格公式相关的程序和相关的注释如下：(

14、1)T仁maromb(inline('2./3*x.A3.*exp(x.A2)'),1,2,b=exp(4)(2)T仁 maromb( in li ne('2*x./(x.A2-3)'),2,3,b=log (6)五、实验结果(实验结果应包括试验的原始数据、中间结果及最终结果，复杂的结果可以用表格或图形形式实现，较为简单的结果可以与实验结果分析合并出现。)实验三 拉格朗日插值:y1 y2 y3L1L2L31.510.5-0.5-1-1.5x作三次多项式拟合结果为:y1 =因此三次拟合函数为：P1=.*X.*X.A.*X+作五次多项式拟合结果为：y2 =因此五次次

15、拟合函数为：p2=.*x.A5+.*x.A4+.*x.A.*x.A.*x+求平方误差，作出离散函数（冷）和拟合函数的图形，结果如下:拟合函数图形：相关数据输出:p1 =p2 = ans = ans =由此可知道：三次拟合离散函数为：(,),(,),(,),(,), (,),(,),(,)三次拟合离散函数为：(,),(,),(,),(,),(,),(,),(,)实验四求积公式计算定积分相关的结果：利用复化梯形公式：(1)f2=4*x*exp(x2)+28 /3*xA3*exp(xA2)+8 /3*xA5*exp(xA2)a =+003h = fun =Inline fun ctio n:fun(

16、x) = 2/3.*x.A3*exp(x.A2)T = b =f2 =-12/(xA2-3)A2*x+16*xA3 /仪人2-3)人3a =104 fun =Inline fun ctio n: fun(x) = 2*x/(xA2-3)T =用复化辛普公式相关的结果如下：(1)f2 =80*x*exp(xA2)+200*xA3*exp(xA2)+96*xA5*exp(xA2)+32/3*xA7*exp(xA2)a =+012h = fun =Inline fun ctio n:fun(x) = 2/3.*x.A3*exp(x.A2)T =(2)f2 =-960/(xA2-3)A4*xA3+240/(xA2-3F3*x+768*xA5心人2-3)人5a =fun =Inline fun ctio n:fun(x) = 2*x/(xA2-3)用龙贝格公式相关的结果如下（1）T =T