Guass-Legendre求积公式近似求定积分报告..docVIP

中国矿业大学（北京）理学院 数值分析实验报告 实验名称 Guass-Legendre求积公式近似求定积分 实验时间 5月15日 组长签名 龙纯鹏 班级 信息与计算科学（1）班 学号 11107200110 成绩 组员签名 11107200101 11107200102 11107200103 11107200119 11107200120 一、实验目的,内容 二、相关背景知识介绍 三、代码 四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会 一、实验目的、内容： 1.掌握Guass-Legendre求积公式 ； 2.学会编程用Guass-Legendre求积公式近似计算定积分； 3.利用两点、三点、四点Guass-Legendre求积公式计算下列积分； (1) ，并与准确值比较精度并估计误差； (2) 求的近似值。 相关背景知识介绍： Guass-Legendre求积公式是在高斯求积公式中，取权函数ρ（x）=1， 区间为[-1,1],则得公式 我们知道勒让德多项式是区间上的正交多项式，因此，勒让德多项式的零点就是求积公式(上式)的高斯点．形如(上式)的高斯公式特别地称为高斯－勒让德求积公式． 若取的零点做节点构造求积公式 令它对准确成立，即可定出． 这样构造出的一点高斯－勒让德求积公式是中矩形公式．再取的两个零点构造求积公式 令它对都准确成立，有 ． 由此解出，从而得到两点高斯－勒让德求积公式 ． 三点高斯－勒让德求积公式的形式是 ． 如表列出高斯－勒让德求积公式的节点和系数． 0 0.0000000 2.0000000 1 0.5773503 1.0000000 2 0.7745967 0.0000000 0.5555556 0.8888889 3 0.8611363 0.3399810 0.3478548 0.6521452 4 0.9061798 0.5384693 0.0000000 0.2369269 0.4786287 0.5688889 代码（Matlab）： %先定义积函数 function y=f(t) a=2; if(a==1) y=1.5*(1.5*t+1.5)*sqrt(1+(1.5*t+1.5)^2) else y=0.5*exp(-1*(0.5+0.5*t)^2); end end %当a=1时表示（1）的运算结果，当a=2时表示（2）的运算结果； q=4; if(q==2) ans=f(-1/(3^0.5))+f(1/(3^0.5)) end if(q==3) ans=(5/9)*f(-3/(15)^0.5)+(8/9)*f(0)+(5/9)*f(3/(15)^0.5) end if(q==4) ans=0.2369269*f(-0.9061798)+0.2369269*f(0.9061798)+0.4786287*f(0.5384693)+0.4786287*f(-0.5384693)+0.5688889*f(0) end %当a=1,q分别取2,3,4表示（1）的两点、三点、四点Guass-Legendre求积运算结果，当a=2时表示（2）的两点、三点、四点Guass-Legendre求积运算结果； 四、数值结果及误差分析 （1） 的运算结果 先计算真实结果： =10.2075922 两点Guass-Legendre求积 三点Guass-Legendre求积 四点Guass-Legendre求积 运算值 10.242281168751251 10.210631347268638 10.207552151929656 误差 0 0 0 （2）的运算结果如下表： 方法 两点Guass-Legendre求积 三点Guass-Legendre求积 四点Guass-Legendre求积 结果 0.746210796131749 0.746824183875915 0.746821079461317 五、 计算中出现的问题，解决方法及体会 此次上机实验，再次熟悉了matlab，通过编程利用Guass-Legendre求积公式进行积分计算的估计，