(1)(2)(3)
??h?h2hf(x)dx?A?1f(?h)?A0f(0)?A1f(h); ;
??2h1f(x)dx?A?1f(?h)?A0f(0)?A1f(h)?1hf(x)dx??f(?1)?2f(x1)?3f(x2)?/3;
.
(4)?0f(x)dx?h?f(0)?f(h)?/1?ah2?f?(0)?f?(h)?1?x22. 分别用梯形公式和辛普森公式计算下列积分:
1(1?e)xdx,n?8dx,n?10??04?x20x(1); (2);
1(3)1; (4)
3. 直接验证柯特斯公式(2.4)具有5次代数精度. 4. 用辛普森公式求积分0并计算误差. 5. 推导下列三种矩形求积公式:
(1)(2)(3)
?9?xdx,n?4?60?sin2?dx,n?6.
?1e?xdx??baba?baf?(?)(b?a)22; f?(?)f(x)dx?(b?a)f(b)?(b?a)22;
a?bf?(?)f(x)dx?(b?a)f()?(b?a)3224. f(x)dx?(b?a)f(a)?6. 证明梯形公式(2.9)和辛普森公式(2.11)当n??时收敛到积分7. 用复化梯形公式求积分a超过?(设不计舍入误差)?
?baf(x)dx.
?bf(x)dx1,问要将积分区间?a,b?分成多少等分,才能保证误差不
28. 用龙贝格方法计算积分??0e?xdx,要求误差不超过10.
??5cS?a?21?()2sin2?d?0a9. 卫星轨道是一个椭圆,椭圆周长的计算公式是,这里a是椭圆
的半长轴,c是地球中心与轨道中心(椭圆中心)的距离,记h为近地点距离,H为远地点距离,R?6371公里为地球半径,则a?(2R?H?h)/2,c?(H?h)/2.我国第一颗人造
卫星近地点距离h?439公里,远地点距离H?2384公里,试求卫星轨道的周长. 10. 证明等式
法求?的近似值.
nsin?n????33!n2??55!n4??试依据nsin(?/n)(n?3,6,12)的值,用外推算
11. 用下列方法计算积分
(1) 龙贝格方法;
(2) 三点及五点高斯公式;
(3) 将积分区间分为四等分,用复化两点高斯公式.
?31dyy并比较结果.
f(x)?12. 用三点公式和五点公式分别求
差.f(x)的值由下表给出: 1.0 1.1 x 1(1?x)2在x?1.0,1.1和1.2处的导数值,并估计误
1.2 0.2066 1.3 0.1890 1.4 0.1736 f(x) 0.2500 0.2268 第五章 常微分方程数值解法
1. 就初值问题y??ax?b,y(0)?0分别导出尤拉方法和改进的尤拉方法的近似解的表达式,并与准确解
2. 用改进的尤拉方法解初值问题
y?12ax?bx2相比较。
?y??x?y,0?x?1;??y(0)?1,
x取步长h=0.1计算,并与准确解y??x?1?2e相比较。
3. 用改进的尤拉方法解
?y??x2?x?y;??y(0)?0,
?x2y??e?x?x?1相比较。 y(0.5)取步长h=0.1计算,并与准确解
4. 用梯形方法解初值问题
证明其近似解为
?y??y?0;??y(0)?1,
n?2?h?yn???,2?h??
?x并证明当h?0时,它原初值问题的准确解y?e。
5. 利用尤拉方法计算积分
??y??x?y,0?x?1;? 1)?y(0)?1,
?y??3y/(1?x),0?x?1;? 2)?y(0)?1.
x0etdt2
在点x?0.5,1,1.5,2的近似值。
6. 取h=0.2,用四阶经典的龙格-库塔方法求解下列初值问题:
7. 证明对任意参数t,下列龙格-库塔公式是二阶的:
8. 证明下列两种龙格-库塔方法是三阶的:
h?y?y?(K2?K3);n?n?12??K?f(x,y);nn?1?K2?f(xn?th,yn?thK1);???K3?f(xn?(1?t)h,yn?(1?t)hK1).
h?y?y?(K1?3K3);n?n?14??K1?f(xn,yn);??hhK?f(x?,y?K1);nn?233??K?f(x?2h,y?2hK);nn2?3331) ? h?y?y?(2K1?3K2?4K3);n?n?19?K?f(xn,yn);??1?hhK?f(x?,y?K1);nn?222??K?f(x?3h,y?3hK).nn2?3442) ?
9. 分别用二阶显式亚当姆斯方法和二阶隐式亚当姆斯方法解下列初值问题:
y??1?y,y(0)?0,
?xh?0.2,y?0,y?0.181,y?1?ey(1.0)01取计算并与准确解相比较。
10. 证明解y??f(x,y)的下列差分公式
yn?1?是二阶的,并求出截断误差的首项。 11. 导出具有下列形式的三阶方法: 12. 将下列方程化为一阶方程组:
1h??1?yn??3yn??1)(yn?yn?1)?(4yn24
??b1yn??1?b2yn??2). yn?1?a0yn?a1yn?1?a2yn?2?h(b0yny???3y??2y?0,1)y(0)?1,y?(0)?1;
y???0.1(1?y2)y??y?0,2)y(0)?1,y?(0)?0;
xy??,y(t)??,r?x2?y2,33rr3)
x(0)?0.4,x?(0)?0,y(0)?0,y?(0)?2.
x??(t)??13. 取h=0.25,用差分方法解边值问题
14. 对方程y???f(x,y)可建立差分公式
?y???y?0;??y(0)?0,y(1)?1.68.
试用这一公式求解初值问题
yn?1?2yn?yn?1?h2f(xn,yn),
验证计算解恒等于准确解
?y???1;??y(0)?y(1)?0,
15. 取h=0.2用差分方法解边值问题
x2?xy(x)?.2
?(1?x2)y???xy??3y?6x?3;??y(0)?y?(0)?1,y(1)?2.
第六章 方程求根
1. 用二分法求方程x?x?1?0的正根,要求误差<0.05。
2. 用比例求根法求f(x)?1?xsinx?0在区间[0,1]内的一个根,直到近似根xk满足精度
2|f(xk)|?0.005时终止计算。
323. 为求方程x?x?1?0在x0?1.5附近的一个根,设将方程改写成下列等价形式,并建立相应的迭代公式。
22x?1?1/xk?1kx?1?1/x1),迭代公式;
2332x?1?xk?1kx?1?x2),迭代公式;
1x2?x?1,迭代公式xk?1?1/xk?1。 3)
试分析每种迭代公式的收敛性,并选取一种公式求出具有四位有效数字的近似根。 4. 比较求e?10x?2?0的根到三位小数所需的计算量;
1)在区间[0,1]内用二分法;
xkx?(2?e)/10,取初值x0?0。 k?12) 用迭代法
5. 给定函数f(x),设对一切x,f?(x)存在且0?m?f?(x)?M,证明对于范围内0???2/M的任意定数λ,迭代过程xk?1?xk??f(xk)均收敛于f(x)的根x?。
x6. 已知x??(x)在区间[a,b]内只有一根,而当a |??(x)|?k?1, 试问如何将x??(x)化为适于迭代的形式? 将x?tgx化为适于迭代的形式,并求x=4.5(弧度)附近的根。 7. 用下列方法求f(x)?x?3x?1?0在x0?2附近的根。根的准确值x=1.87938524…,要求计算结果准确到四位有效数字。 1) 用牛顿法; ?32)用弦截法,取x0?1,x1?1.9; 3)用抛物线法,取x0?1,x1?3,x2?2。 8. 用二分法和牛顿法求x?tgx?0的最小正根。 9. 研究求a的牛顿公式 xk?1?证明对一切k?1,2,?,xk?1a(xk?),x0?0,2xk a且序列x1,x2,?是递减的。 10. 对于f(x)?0的牛顿公式xk?1?xk?f(xk)/f?(xk),证明 Rk?(xk?xk?1)/(xk?1?xk?2)2 收敛到?f??(x)/(2f?(x)),这里x为f(x)?0的根。 11. 试就下列函数讨论牛顿法的收敛性和收敛速度: ?????x,x?0;f(x)??????x,x?0; 1) 23??x,x?0;f(x)??23??x,x?0. ?2) 3212. 应用牛顿法于方程x?a?0,导出求立方根a的迭代公式,并讨论其收敛性。 13. 应用牛顿法于方程值。 f(x)?1?a?0x2,导出求a的迭代公式,并用此公式求115的 f(x)?1?a?0nxn,分别导出求a的迭代公 14. 应用牛顿法于方程f(x)?x?a?0和 式,并求 k??nlim(na?xk?1)/(na?xk)2.15. 证明迭代公式 xk?1x(x?3a)?kk23xk?a 2?是计算a的三阶方法。假定初值x0充分靠近根x,求 lim(a?xk?1)/(a?xk)3.k?? 第七章 解线性方程组的直接方法 1. 考虑方程组: ?0.4096x1?0.1234x2?0.2246x?0.3872x?12??0.3645x1?0.1920x2??0.1784x1?0.4002x2?0.3678x3?0.2943x4?0.4043;?0.4015x3?0.1129x4?0.1550;?0.3781x3?0.0643x4?0.4240;?0.2786x3?0.3927x4??0.2557; (a) 用高斯消去法解此方程组(用四位小数计算), (b) 用列主元消去法解上述方程组并且与(a)比较结果。 2. (a) 设A是对称阵且a11?0,经过高斯消去法一步后,A约化为