河南科技大学数值分析(计算方法)期末复习画题资料.docxVIP

1.数值积分公式形如(15) °f(x)dx 二 Aof(O) Af(1) B°f (0) (1)试确定求积公式中的参数 Ao,A,B°，使其代数精度尽可能高?并求出其代数精度。 ⑵ 已知该求积公式余项 R[f ]二kf「)，1三(0,1),试求出余项中的参数 k o (1)解：1f(x (1)解： 1 f(x^ —1 时，左=(f (x)dx =1，右=代+A，左=右得： A)+ A=1 B0 Al - _ 2 f (x) =x 时，左 2 f(x) =x 时，左 联立上述三个方程, 1 1 =I f (x)dx =3，右=Bo +A，左=右得： 1 1 =(0 f (x)dx =—，右=A，左=右得： A 3 解得： TOC \o 1-5 \h \z A 2 1 1 AO , Bo , A1 -— 3 6 3 3f(x)二 3 f(x)二 X 时，左 f (x)dx ，右二A ，左尸右 0 4 3 所以，该求积公式的代数精度是 2 (2)解：过点0, 1构造f(x)的Hermite插值H2(x)，因为该求积公式代数精度为 2,所 以有： 1 °H2(x)dx 二人出(0) 的2(0) B°H2(0) (0) Af (1) f (0) 其求积余项为： 1 R(f^f0f (x)dx —［人 f (0)十 A f ⑴十 B。f (0)］ 1 1 1 f ( ) 2 =0 f(x)dx - ° H2(x)dx = ° 3! x (x-1)dx f ( ) 1 2 )J°x2(x-1)dx f () 1 所以，k「72 72 y = 3x + 2y 2?设初值问题［y(°)=1叶涂1 写出用改进的Euler法解上述初值问题数值解的公式， 若h = 0.2 ,求解 浙，丫2 ,保留两位小 数。(10分) .解：改进的 Euler公式是： Yn 1 二 yn hf (Xn, %) j h % 1 = yn 2【f (Xn’yn) f(Xn1，Vn 1)] 具体到本题中，求解的公式是: % 1 二 yn 0.2(3Xn 2yn) =1.4yn 0.6Xn yn 1 二 yn 0.1[3Xn 2yn 3Xn 1 2% 1] y(0) =1 代入求解得：y1 =1.4, w =1.54 y2 =2.276,y2 =2.4832 3.分别用梯形公式，复化梯形公式计算积分： 1 1 dx 01 x 其中在用复化梯形公式求积分时，步长 h = 0.2 o（ 10分） 梯形公式为： TOC \o 1-5 \h \z b 1』 b — a q T (f(a) f (b)) 2 =0.75 2 2 复化梯形公式为： — n」 T”S(f(a) 2：f(xi) f(b)) 具体到本题中，可知h =：0.2,a =0,b =1 h 4 T6r —(f(0) 2、f (xi) f (1))=0.1 (1.5 5.456) =0.6956 2 y 4?用改进的欧拉方法求解初值问题: f 2 y =x +x — y (0Ex^0.4) y(0) =0 取步长h =0.2，计算过程中保留到小数点后四位。 （10 分） .改进的Euler公式为: yn 1 二 yn hf(Xn,yn) TOC \o 1-5 \h \z j h yn 1 =yn ；(f(Xn, yn) f(Xn 1, ?n 1)) L 2 具体到本题中，则为 2 yn 1 =yn h(Xn X. - y.) i h 2 2 2 yn1=yn ?[Xn 人 -y. (x. h) (x. h) - y. - h(x Xn-yn)] 经化简为： yn1 =0.82yn 0.18xn2 0.22xn 0.024 所以： y(0.2) 0.024 0 y(0.4) 0.0949 5?证明：设f 5?证明：设f(X)=1,左二b -a,右 f (x) = x 左=2 b2 (1 1) = b -a,左=右 2 右,=b?* b (a = j b 々a 2左=)右 .3 32 t .3 3 2 t. b — a f (x)二 X ,左二一3— 右二空迴(a2 b2) b3 b—,… 所以，该公式具有一次代数精度6. 所以，该公式具有一次代数精度 6.用两点Gauss-Legendre求积公式求积分 1 /dx 解：两点 Gauss-legrende求积公式为: 1」 1 」f (x)dx 二 f (- .3f(5) 所以7.用欧拉法求解常微分方程组初值问题：M =3% +2y?W y2 所以 7.用欧拉法求解常微分方程组初值问题： M =3% +2y? W y2 %(0) =0必(0)=1 在［0,0.4］上的数值解，取步长h =0.2，计算过程中保留两位小数。 Euler公式为: (10 分