romberg积分法,gauss型积分法.docxVIP

数学软件实验任务书 课程名称 数学软件实验 班级 数 0901 Romberg积分法， Gauss 型积分法，高斯 - 勒让德积分 实验课题 法，高斯 - 切比雪夫积分法，高斯 - 拉盖尔积分法，高 斯- 埃尔米特积分法 熟悉 Romberg 积分法，Gauss型积分法，高斯 - 勒让德 实验目的 积分法，高斯 - 切比雪夫积分法， 高斯 - 拉盖尔积分法， 高斯 - 埃尔米特积分法 实验要求 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中 一种语言完成 Romberg积分法， Gauss 型积分法，高斯 - 勒让德积分 实验容 法，高斯 - 切比雪夫积分法，高斯 - 拉盖尔积分法，高 斯- 埃尔米特积分法 成绩 教师 实验 1Romberg积分法 1 实验原理 Romberg 方法是实用性很强的一种数值积分方法，其收敛速度是 很快的，这里给出 Romberg积分的计算方法。 （ 1）计算 R(0,0) 1 (b a)[ f ( a) f (b)] 2 （ 2）计算 1 hi 1 2 i 2 1 2 (i 1,0) 2 f (a (k 2 )hi 1) R(i,0) k 1 （ 3）计算 R( m, j ) 4j 1 R(m, j 1) R( m 1, j 1) 4 j 1 1 2 实验数据 用 Romberg 积分方法计算： 15 x 2 dx S 4 x 0 3 实验程序 程序 1 function s=rombg(a,b,TOL) n=1; h=b-a; delt=1; x=a; k=0; R=zeros(4,4); R(1,1)=h*(rombg_f(a)+rombg_f(b))/2; while deltTOL k=k+1; h=h/2; s=0; for j=1:n x=a+h*(2*j-1); s=s+rombg_f(x); end R(k+1,1)= R(k,1)/2+h*s; n=2*n; for i=1:k R(k+1,i+1)=((4^i)*R(k+1,i)-R(k,i))/(4^i-1); end delt=abs(R(k+1,k)-R(k+1,k+1)); end s=R(k+1,k+1); 程序 2 function f=rombg_f(x) f=x/(4+x^2); 程序 3 s=rombg(0,1.5,1.e-6) %作出图形 x=0:0.02:1.5; y=x./(4+x.^2); area(x,y) grid 实验结果 s = 0.2231 实验 2 高斯 - 勒让德积分法 实验原理 Gauss-Legendre 求积公式为 1 n f ( x)dx Ak f ( xk ) 1 1 k 其中 xk 为 Legendre 多项式在 1,1 区间上的零点。 n 阶 Legendre 多项式定义为： Pn (t ) 1 d n 2 n n n! dt n [( t 1) ] 2 Ak 为权系数， Ak 2 2(1 xk )2 2[ Pn ( xn )2 ] n2 [ Pn ( xk )]2 (1 xk ) 对于一般的积分区间为 a, b 问题，可以做变换 x a b b a 2 t 2 b b a n a b b a a f ( x)dx Ak f ( xk ) 2 k 1 2 2 实验数据 用 Gauss-Legendre 积分方法计算定积分 S 2 x2 cosxdx 0 实验程序 function s=gau_leg(a,b) %5阶 Legendre 多项式结点 node=[-0.9061798459,-0.5384693101,0,0.5384693101,0.90617 98459]; %结点对应的权 quan=[0.2369268851,0.4786286705,0.5688888889,0.478628670 5,0.2369268851]; %t 为（ 1， 5）的行向量，整个区间上的结点 t=(b+a)/2+(b-a)*node/2; s=((b-a)/2)*sum(quan.*gau_leg_f(t)); function f=gau_leg_f(x) f=(x.^2).*cos(x); disp( 计算结果为 :) s=gau_leg(0,pi/2) %画出图形 x=0:0.01:pi/2; y=(x.^2).*cos(x); bar(x,y) grid 实验结果计算结果为 : s = 0.4674 实验 3 高斯 - 拉盖尔积分法 实验原理 n 个结点 Gauss-Laguerre 求积公为： n S  Ak f ( xk ) k 1 其中 xk 为零点， Ak 为权系数 Ak xk 2