计算方法答案第三章

合集下载

计算方法第三章习题答案

计算方法第三章习题答案计算方法第三章习题答案计算方法是一门涵盖了数值计算和计算机编程的学科，它在现代科学和工程中扮演着重要的角色。

第三章是计算方法课程中的重要章节，主要涉及到数值计算中的误差分析和插值方法。

本文将为大家提供第三章习题的详细答案，帮助读者更好地理解和应用这些概念。

1. 误差分析误差分析是计算方法中非常重要的一部分，它帮助我们理解和评估数值计算中的误差来源。

以下是一些常见的误差类型：- 绝对误差：绝对误差是指数值计算结果与真实值之间的差异。

它可以通过计算两者之差来得到。

- 相对误差：相对误差是指绝对误差与真实值之间的比值。

通常以百分比的形式表示。

- 截断误差：截断误差是由于在计算过程中舍入或截断数字而引入的误差。

它通常是由于计算机的有限精度导致的。

- 舍入误差：舍入误差是由于将无限位数的小数截断为有限位数而引入的误差。

它通常是由于计算机的有限精度或计算方法的近似性质导致的。

2. 插值方法插值方法是一种用于通过已知数据点来估计未知数据点的技术。

以下是一些常见的插值方法：- 线性插值：线性插值是一种简单的插值方法，它假设两个已知数据点之间的未知数据点的取值在直线上。

通过已知数据点的斜率和截距，我们可以计算出未知数据点的值。

- 拉格朗日插值：拉格朗日插值是一种使用多项式来逼近已知数据点的方法。

它通过构造一个满足通过已知数据点的多项式来估计未知数据点的值。

- 牛顿插值：牛顿插值是一种使用差商来逼近已知数据点的方法。

它通过构造一个满足通过已知数据点的差商多项式来估计未知数据点的值。

3. 习题答案以下是一些第三章习题的答案，供大家参考：- 习题1：已知函数f(x)在区间[a, b]上连续，且在[a, b]上的导数存在且连续，证明存在一点c∈(a, b)，使得f(b) - f(a) = (b - a)f'(c)。

这是拉格朗日中值定理的一个特例，根据定理的条件，我们可以得到上述结论。

- 习题2：已知函数f(x)在区间[a, b]上连续，且在(a, b)内可导，证明存在一点c∈(a, b)，使得f'(c) = (f(b) - f(a))/(b - a)。

第三章物理学中定积分的数值计算方法

第三章物理学中定积分的数值计算方法一、填空题1、库仑常数k 等于 9×109mV/C ，真空中的介电常数ε0等于8.85×10-12F/m 。

2、对于电量为Q 的点电荷，在距离r 处产生的电场强度为21ˆˆ()4QrE r rrrπε==。

3、已知定积分()ba f x dx ⎰，被积分函数为()f x ，积分区间为[],ab 。

将该区间N 等分，步长()/x b a N ∆=-，用曲线下的虚矩形面积和近似替代积分值，该方法称为矩形法。

积分近似计算公式为1()()N bi ai I f x dx f x x -==≈∆∑⎰。

4、毕奥—萨伐尔定律所描述的公式为034Idl rdB r μπ⨯=。

5、玻尔兹曼常数是 k=1.38×1023 J/K 。

6、麦克斯韦速率分布律公式23/22/2()4()2v kTdN f v dv v e dv N kTμμππ-==。

7、在计算物理中求解定积分的方法有辛普森法、龙贝格法、高斯求积法等。

二、简答1、写出库仑常数、真空中的介电常数和玻尔兹曼常数的值。

答：库仑常数k= 9×109mV/C ，真空中的介电常数ε0= 8.85×10-12F/m ，玻尔兹曼常数是 k=1.38×1023 J/K 。

2、什么是矩形法？答：已知定积分()ba f x dx ⎰，被积分函数为()f x ，积分区间为[],ab 。

将该区间N 等分，步长()/x b a N ∆=-，用曲线下的虚矩形面积和近似替代积分值，该方法称为矩形法。

积分近似计算公式为1()()N bi ai I f x dx f x x -==≈∆∑⎰。

3、毕奥—萨伐尔定律和麦克斯韦速率分布律公式。

答：毕奥—萨伐尔定律所描述的公式为034Idl rdB rμπ⨯=。

麦克斯韦速率分布律公式23/22/2()4()2v kTdN f v dv v e dv N kTμμππ-==。

西北工业大学计算方法作业集答案及试题

答：（1）
2 则有 er ( S ) < er ( a * ) + er (b * ) + er (c * )
*
注意当 0 < c <
*
π
时， tgc * > c * > 0 ，即 (tgc * )
−1
< (c * ) 。
−1
7．设 y0 = 由
1 * * 2 ， y0 = 1.41 ， y0 − y0 ≤ × 10 − 2 = δ 2 * −1 * −1 y1 − y1 = 10 y0 − y0 ≤ 10 δ ，
η ∈ [ a, b]
1 f ′(η )(b − a ) 2 2
（2）右矩形公式将 f(x)在 b 处展开，并积分，得（3）中矩形公式将 f(x)在 a + b 处展开，得
2
∫
b
a
f ( x)dx = (b − a ) f (b) −
x * ( x > 0 )的相对误差约是 x * 的相对误差的 1/2 倍； * * n （2） ( x ) 的相对误差约是 x 的相对误差的 n 倍。 1 * * 1 * 1 * b sin c *e(a * ) a sin c *e(b* ) a b cos c *e(c * ) * 2 2 2 6．根据 er ( S ) ≤ + + 1 * * 1 * * 1 * * a b sin c * a b sin c * a b sin c * 2 2 2 * * * e(a ) e(b ) e(c ) = + * + a* b tgc *
I = 5.6308e −2.8882t
3.1781 4 3.1781 3.6092

计算方法-刘师少版课后习题答案

1.1 设3.14, 3.1415, 3.1416分别作为π的近似值时所具有的有效数字位数解近似值x =3.14＝0.314×101,即m =1，它的绝对误差是－0.001 592 6…，有31105.06592001.0-*⨯≤=- x x .即n =3，故x =3.14有3位有效数字. x =3.14准确到小数点后第2位.又近似值x =3.1416，它的绝对误差是0.0000074…，有5-1*10⨯50≤00000740=-.. x x即m =1,n ＝5，x =3.1416有5位有效数字.而近似值x =3.1415，它的绝对误差是0.0000926…，有4-1*10⨯50≤00009260=-.. x x即m =1,n ＝4，x =3.1415有4位有效数字.这就是说某数有s 位数，若末位数字是四舍五入得到的，那么该数有s 位有效数字 1.2 指出下列各数具有几位有效数字，及其绝对误差限和相对误差限： 2.0004 －0.00200 9000 9000.00解（1）∵ 2.0004＝0.20004×101, m=1绝对误差限：4105.0000049.020004.0-*⨯≤≤-=-x x x m -n =-4,m =1则n =5,故x =2.0004有5位有效数字1x =2，相对误差限000025.010221102151)1(1=⨯⨯=⨯⨯=---n r x ε（2）∵ －0.00200= -0.2×10-2, m =-25105.00000049.0)00200.0(-*⨯≤≤--=-x x xm -n =-5, m =-2则n =3,故x =-0.00200有3位有效数字1x =2，相对误差限3110221-⨯⨯=r ε=0.0025(3) ∵ 9000=0.9000×104, m =4，0105.049.09000⨯<≤-=-*x x xm -n =0, m =4则n =4,故x =9000有4位有效数字4110921-⨯⨯=r ε＝0.000056 (4) ∵9000.00=0.900000×104, m =4，2105.00049.000.9000-*⨯<≤-=-x x xm -n =-2, m =4则n =6,故x =9000.00有6位有效数字相对误差限为6110921-⨯⨯=rε＝0.000 00056由(3)与(4)可以看到小数点之后的0，不是可有可无的，它是有实际意义的.1.3 ln2=0.69314718…，精确到310-的近似值是多少？解精确到310-＝0.001，即绝对误差限是ε＝0.0005,故至少要保留小数点后三位才可以.ln2≈0.6932.1 用二分法求方程013=--x x在[1, 2]的近似根，要求误差不超过31021-⨯至少要二分多少？解：给定误差限ε＝0.5×10－3，使用二分法时，误差限为)(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211a b k 即可，亦即 96678.912lg 10lg 35.0lg 12lg lg )lg(=-+-=---≥εa b k只要取n ＝10.2.3 证明方程1 -x –sin x ＝0 在区间[0, 1]内有一个根，使用二分法求误差不超过0.5×10-4的根要二分多少次？证明令f (x )＝1－x －sin x ， ∵ f (0)=1>0，f (1)=－sin1<0∴ f (x )=1－x －sin x =0在[0，1]有根.又 f '(x )=-1－c os x<0 (x ∈[0.1]),故f (x ) 在[0，1]单调减少，所以f (x ) 在区间[0，1]内有唯一实根.给定误差限ε＝0.5×10－4，使用二分法时，误差限为)(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211a b k 即可，亦即7287.1312lg 10lg 45.0lg 12lg lg )lg(=-+-=---≥εa b k只要取n ＝14.2.4 方程0123=--x x 在x =1.5附近有根，把方程写成四种不同的等价形式，并建立相应的迭代公式：（1）211xx +=，迭代公式2111kk x x +=+ （2）231x x +=，迭代公式3211k k x x +=+ （3）112-=x x，迭代公式111-=+k k x x （4）13-=x x ，迭代公式131-=+k k x x试分析每种迭代公式的收敛性，并选取一种收敛迭代公式求出具有四位有效数字的近似根。

数值计算方法课后习题答案

第一章绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。

[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ，相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%，求n x 的相对误差。

[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。

3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。

[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。

4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。

（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。

计算方法各章习题及答案

第二章数值分析2.1 已知多项式432()1p x x x x x =-+-+通过下列点：试构造一多项式()q x 通过下列点：答案：54313()()()3122q x p x r x x x x x =-=-++-+． 2.2 观测得到二次多项式2()p x 的值：表中2()p x 的某一个函数值有错误，试找出并校正它．答案：函数值表中2(1)p -错误，应有2(1)0p -=．2.3 利用差分的性质证明22212(1)(21)/6n n n n +++=++．2.4 当用等距节点的分段二次插值多项式在区间[1,1]-近似函数xe 时，使用多少个节点能够保证误差不超过61102-⨯．答案：需要143个插值节点．2.5 设被插值函数4()[,]f x C a b ∈，()3()h H x 是()f x 关于等距节点01n a x x x b =<<<=的分段三次艾尔米特插值多项式，步长b ah n-=．试估计()3||()()||h f x H x ∞-．答案：()443||()()||384h M f x H x h ∞-≤．第三章函数逼近3.1 求()sin ,[0,0.1]f x x x =∈在空间2{1,,}span x x Φ=上最佳平方逼近多项式，并给出平方误差．答案：()sin f x x =的二次最佳平方逼近多项式为-522sin ()0.832 440 710 1.000 999 10.024 985 1x p x x x ≈=-⨯+-，二次最佳平方逼近的平方误差为0.122-1220(sin )())0.989 310 710x p x dx δ=-=⨯⎰．3.2 确定参数,a b c 和，使得积分2121(,,)[I a b c ax bx c -=++-⎰取最小值．答案：810, 0, 33a b c ππ=-== 3.3 求多项式432()251f x x x x =+++在[1,1]-上的３次最佳一致逼近多项式()p x ．答案：()f x 的最佳一致逼近多项式为323()74p x x x =++． 3.4 用幂级数缩合方法，求() (11)x f x e x =-≤≤上的３次近似多项式6,3()p x ，并估计6,3||()()||f x p x ∞-．答案：236,3()0.994 574 650.997 395 830.542 968 750.177 083 33p x x x x =+++， 6,3||()()||0.006 572 327 7f x p x ∞-≤3.5 求() (11)xf x e x =-≤≤上的关于权函数()x ρ=的三次最佳平方逼近多项式3()S x ，并估计误差32||()()||f x S x -和3||()()||f x S x ∞-．答案：233()0.994 5710.997 3080.542 9910.177 347S x x x x =+++，32||()()||0.006 894 83f x S x -=，3||()()||0.006 442 575f x S x ∞-≤．第四章数值积分与数值微分4.1 用梯形公式、辛浦生公式和柯特斯公式分别计算积分1(1,2,3,4)n x dx n =⎰，并与精确值比较．答案：计算结果如下表所示4.2 确定下列求积公式中的待定参数，使得求积公式的代数精度尽量高，并指明所确定的求积公式具有的代数精度．（１）101()()(0)()hh f x dx A f h A f A f h --≈-++⎰（２）11211()[(1)2()3()]3f x dx f f x f x -≈-++⎰ （３）20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-⎰答案：（１）具有三次代数精确度（２）具有二次代数精确度（３）具有三次代数精确度．4.3 设10h x x =-，确定求积公式12300101()()[()()][()()][]x x x x f x dx h Af x Bf x h Cf x Df x R f ''-=++++⎰中的待定参数,,,A B C D ，使得该求积公式的代数精确度尽量高，并给出余项表达式．答案：3711,,,20203020A B C D ====-，(4)6()[]1440f R f h η=，其中01(,)x x η∈．4.4 设2()P x 是以0,,2h h 为插值点的()f x 的二次插值多项式，用2()P x 导出计算积分30()hI f x dx =⎰的数值积分公式h I ，并用台劳展开法证明：453(0)()8h I I h f O h '''-=+．答案：3203()[(0)3(2)]4h h I p x dx h f f h ==+⎰．4.5 给定积分10sin xI dx x =⎰（１）运用复化梯形公式计算上述积分值，使其截断误差不超过31102-⨯．（２）取同样的求积节点，改用复化辛浦生公式计算时，截断误差是多少？（３）要求的截断误差不超过610-，若用复化辛浦生公式，应取多少个节点处的函数值？答案：（１）只需7.5n ≥，取９个节点，0.946I ≈（２）4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=⨯ （３）取７个节点处的函数值．4.6 用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分10sin xI dx x =⎰．要求用事后误差估计法时，截断误不超过31102-⨯和61102-⨯．答案：使用复化梯形公式时，80.946I T ≈=满足精度要求；使用复化辛浦生公式时，40.946 083I s ≈=满足精度要求．4.7（１）利用埃尔米特插值公式推导带有导数值的求积公式2()()[()()][()()][]212ba b a b a f x dx f a f b f b f a R f --''=+--+⎰，其中余项为 5(4)()[](), (,)4!30b a R f f a b ηη-=∈．（２）利用上述公式推导带修正项的复化梯形求积公式020()[()()]12Nx N N x h f x dx T f x f x ''≈--⎰，其中 0121[()2()2()2()()]2N N N hT f x f x f x f x f x -=+++++，而 00, (0,1,2,,), i N x x ih i N Nh x x =+==-．4.8 用龙贝格方法计算椭圆2214x y +=的周长，使结果具有五位有效数字．答案：49.6884l I =≈．4.9确定高斯型求积公式0011()()()x dx A f x A f x ≈+⎰的节点0x ，1x 及系数0A ，1A ．答案：00.289 949x =，10.821 162x =，00.277 556A =，10.389 111A =．4.10 验证高斯型求积公式00110()()()x e f x dx A f x A f x +∞-≈+⎰的系数及节点分别为0001 2 2A A x x ===-=+第五章解线性方程组的直接法5.1 用按列选主元的高斯-若当消去法求矩阵A 的逆矩阵，其中111210110A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭．答案： 1110331203321133A -⎛⎫ ⎪ ⎪⎪=- ⎪ ⎪⎪-- ⎪⎝⎭5.2 用矩阵的直接三角分解法解方程组1234102050101312431701037x x x x ⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案： 42x =，32x =，21x =，11x =．5.3 用平方根法(Cholesky 分解法)求解方程组12341161 4.25 2.750.51 2.75 3.5 1.25x x x -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪-=- ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭答案： 12x =，21x =，31x =-．5.4 用追赶法求解三对角方程组123421113121112210x x x x ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案：42x =，31x =-，21x =，10x =．第六章解线性代数方程组的迭代法６.１对方程1212123879897x x x x x x x -+=⎧⎪-+=⎨⎪--=⎩作简单调整，使得用高斯－赛得尔迭代法求解时对任意初始向量都收敛，并取初始向量(0)[0 0 0]T x =，用该方法求近似解(1)k x+，使(1)()3||||10k k x x +-∞-≤．答案：近似解为(4)[1.0000 1.0000 1.0000]Tx =．６.２讨论松弛因子 1.25ω=时，用SOR 方法求解方程组121232343163420412x x x x x x x +=⎧⎪+-=⎨⎪-+=-⎩ 的收敛性．若收敛，则取(0)[0 0 0]T x=迭代求解，使(1)()41||||102k k x x +-∞-<⨯．答案：方程组的近似解为*1 1.50001x =，*23.33333x =，*3 2.16667x =-．６.３给定线性方程组Ax b =，其中111221112211122A ⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭，证明用雅可比迭代法解此方程组发散，而高斯－赛得尔迭代法收敛．６.４设有方程组112233302021212x b x b x b -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭，讨论用雅可比方法和高斯－赛得尔方法解此方程组的收敛性．如果收敛，比较哪种方法收敛较快．答案：雅可比方法收敛，高斯－赛得尔方法收敛,且较快．６.５设矩阵A 非奇异．求证：方程组Ax b =的解总能通过高斯－赛得尔方法得到．６.６设()ij n nA a ⨯=为对称正定矩阵，对角阵1122(,,,)nn D diag a a a =．求证：高斯－赛得尔方法求解方程组1122D AD x b --=时对任意初始向量都收敛．第七章非线性方程求根例７.４对方程230xx e -=确定迭代函数()x ϕ及区间[,]a b ，使对0[,]x a b ∀∈，迭代过程1(), 0,1,2,k x x k ϕ+==均收敛，并求解．要求51||10k k x x -+-<．答案：若取2()x x ϕ=，则在[1,0]-中满足收敛性条件，因此迭代法121, 0,1,2,k x k x k +==在(1,0)-中有惟一解．取00.5x =-，*70.458960903x x ≈=-．取2()x x ϕ=，在[0,1上满足收敛性条件，迭代序列121, 0,1,2,k x k x k +==在[0,1]中有惟一解．取00.5x =，*140.910001967x x ≈=-在[3,4]上，将原方程改写为23xe x =，取对数得2ln(3)()x x x ϕ==．满足收敛性条件，则迭代序列21ln(3), 0,1,2,k k x x k +==在[3,4]中有惟一解．取0 3.5x =， *16 3.733067511x x ≈=．例７.６对于迭代函数2()(3)x x c x ϕ=+-，试讨论：（１）当c 为何值时，1()k k x x ϕ+=产生的序列{}k x（２）c 取何值时收敛最快？（３）取1,2c =-()x ϕ51||10k k x x -+-<．答案：（１）(c ∈时迭代收敛．（２）c =时收敛最快．（３）分别取1, 2c =--，并取0 1.5x =，计算结果如下表７.７所示表７.７例７.13 设不动点迭代1()k x x ϕ+=的迭代函数()x ϕ具有二阶连续导数，*x 是()x ϕ的不动点，且*()1x ϕ'≠，证明Steffensen 迭代式21(), (), 0,1,2,()2k k k k k k k k k k k y x z x k y x x x z y xϕϕ+==⎧⎪=-⎨=-⎪-+⎩二阶收敛于*x ．例7.15 设2()()()()()x x p x f x q x f x ϕ=--，试确定函数()p x 和()q x ，使求解()0f x =且以()x ϕ为迭代函数的迭代法至少三阶收敛．答案：1()()p x f x ='，31()()2[()]f x q x f x ''=' 例７.19 设()f x 在[,]a b 上有高阶导数，*(,)x a b ∈是()0f x =的(2)m m ≥重根，且牛顿法收敛，证明牛顿迭代序列{}k x 有下列极限关系：111lim2k kk k k k x x m x x x -→∞-+-=-+．第八章矩阵特征值8.1 用乘幂法求矩阵A 的按模最大的特征值与对应的特征向量，已知5500 5.51031A -⎛⎫ ⎪=- ⎪ ⎪-⎝⎭，要求(1)()611||10k k λλ+--<，这里()1k λ表示1λ的第k 次近似值．答案：15λ≈，对应的特征向量为[5,0,0]T-；25λ≈-，对应的特征向量为[5,10,5]T --． 8.2 用反幂法求矩阵110242012A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭的按模最小的特征值．知A 的按模较大的特征值的近似值为15λ=，用5p =的原点平移法计算1λ及其对应的特征向量．答案：(1) A 的按模最小的特征值为30.2384428λ≈(2) 1 5.1248854λ≈，对应的特征向量为(8)[0.242 4310, 1 ,0.320 011 7]T U =--．8.3 设方阵A 的特征值都是实数，且满足121, ||||n n λλλλλ>≥≥>，为求1λ而作原点平移，试证：当平移量21()2n p λλ=+时，幂法收敛最快． 8.4 用二分法求三对角对称方阵1221221221A ⎛⎫ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭的最小特征值，使它至少具有２位有效数字．答案：取5 2.234375λ≈-即有２位有效数字．8.5 用平面旋转变换和反射变换将向量[2 3 0 5]T x =变为与1[1 0 0 0]Te =平行的向量．答案：203/2/00001010/0T ⎛⎫⎪- ⎪=⎪--⎝0.324 442 8400.486 664 26200.811 107 1040.486 664 2620.812 176 04800.298 039 92200100.811 107 1040.298 039 92200.530 266 798H --⎛⎫⎪--⎪= ⎪ ⎪⎪--⎝⎭8.6 若532644445A -⎛⎫⎪=- ⎪ ⎪-⎝⎭，试把A 化为相似的上Hessenberg 阵，然后用QR 方法求A 的全部特征值．第九章微分方程初值问题的数值解法9.1 用反复迭代（反复校正）的欧拉预估－校正法求解初值问题0, 0<0.2(0)1y y x y '+=≤⎧⎨=⎩，要求取步长0.1h =，每步迭代误差不超过510-．答案： [4]11(0.1)0.904 762y y y ≈==，[4]22(0.2)0.818 594y y y ≈==9.2 用二阶中点格式和二阶休恩格式求初值问题2, 0<0.4(0)1dy x y x dx y ⎧=+≤⎪⎨⎪=⎩的数值解(取步长0.2h =，运算过程中保留五位小数)．答案：用二阶中点格式，取初值01y =计算得0n =时，1211.000 00, 1.200 00, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.298 72, (0.4)=1.699 74K K y y ==≈用二阶休恩格式，取初值01y =计算得0n =时，1211.000 00, 1.266 67, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.499 18, (0.4)=1.701 76K K y y ==≈9.3 用如下四步四阶阿达姆斯显格式1123(5559379)/24n n n n n n y y h f f f f +---=+-+-求初值问题, (0)1y x y y '=+=在[0,0.5]上的数值解．取步长0.1h =，小数点后保留８位．答案：4(0.4)0.583 640 216y y ≈=，5(0.5) 1.797 421 984y y ≈=． 9.4 为使二阶中点公式1(,(,))22n n n n n n h hy y hf x y f x y +=+++，求解初值问题 , (0)y y y aλλ'=-⎧⎨=⎩为实常数绝对稳定，试求步长h 的大小应受到的限制条件．答案：2h λ≤．9.5 用如下反复迭代的欧拉预估－校正格式(0)1(1)()111(,)[(,)(,)]2 0,1,2,; 0,1,2,nn n n k k n n n n n n y y hf x y h y y f x y f x y k n +++++⎧=+⎪⎪=++⎨⎪⎪==⎩，求解初值问题sin(), 01(0)1x y e xy x y '⎧=<≤⎨=⎩时，如何选择步长h ，使上述格式关于k 的迭代收敛．答案：2h e<时上述格式关于k 的迭代是收敛的．9.6 求系数,,,a b c d ，使求解初值问题0(,), ()y f x y y x a '==的如下隐式二步法221()n n n n n y ay h bf cf df +++=+++的误差阶尽可能高，并指出其阶数．答案：系数为142,,33a b d c ====，此时方法的局部截断误差阶最高，为五阶5()O h ．9.7 试用欧拉预估－校正法求解初值问题, (0)=1, 0<0.2()/, (0)2dyxy z y dxx dz x y z z dx⎧=-⎪⎪≤⎨⎪=+=⎪⎩，取步长0.1h =，小数点后至少保留六位．答案：由初值00(0)1, (0)2y y z z ====可计算得110.800 000z 2.050 000y =⎧⎨=⎩ ， 11(0.1)0.801 500(0.1) 2.046 951y y z z ≈=⎧⎨≈=⎩ 220.604 820z 2.090 992y =⎧⎨=⎩ ， 22(0.2)0.604 659(0.2) 2.088 216y y z z ≈=⎧⎨≈=⎩。

吉林大学工程数学计算方法第三章习题答案word资料9页

第三章习题答案1.分别用梯形公式、Simpson 公式、Cotes 公式计算积分1,I =⎰并估计误差。

解：1）用梯形公式有：事实上， 2）Simpson 公式事实上，()()()110.50.510.5410.000030462S E f f f f -⎡+⎤⎛⎫=-++= ⎪⎢⎥⎝⎭⎣⎦⎰3）由Cotes 公式有：事实上，()0.0000003C E f =2．证明Simpson 公式()2.8具有三次代数精度。

证明：而当()4f x x =时左侧：()()45515b b f x dx x dx b a a a ==-⎰⎰ 右侧：左侧不等于右侧。

所以Simpson 具有三次代数精度. 3.分别用复化梯形公式和复化公式Simpson 计算下列积分.(1)21,804x dx n x =+⎰，（3）,4n =⎰，6,sin 4602=-⎰n d ϕϕπ解：（1）用复化梯形公式有：10188b a h n --===，()()[]12345672128888888102(0.0311280.0615380.0905660.117650.142350.164380.18361)0.20.111416n h T f a f f f f f f f f ⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫=++++++++⎢⎥ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦=+⨯+++++++=由复化Simpson 公式有：解：删去解（3）：,4n =⎰由复化梯形公式有：由复化Simpson 公式有：（4）解：6,sin 4602=-⎰n d ϕϕπ由复化梯形公式：由复化Simpson 公式：4．给定求积节点012113,,,424x x x ===试推出计算积分()10f x dx ⎰的插值型求积公式，并写出它的截断误差。

解：考虑到对称性，有20A A =，于是有求积公式由于原式含有3个节点，故它至少有2阶精度。

数值计算方法(宋岱才版)课后答案

第一章绪论一本章的学习要求（1）会求有效数字。

（2）会求函数的误差及误差限。

（3）能根据要求进行误差分析。

二本章应掌握的重点公式（1）绝对误差：设x 为精确值，x *为x 的一个近似值，称e x x **=-为x *的绝对误差。

（2）相对误差：r e e x***=。

（3）绝对误差限：e x x ε***==-。

（4）相对误差限：r x x xxεε*****-==。

（5）一元函数的绝对误差限：设一元函数()()()0,df f x f x dx εε***⎛⎫==⋅ ⎪⎝⎭则。

（6）一元函数的相对误差限：()()1r df f x dx f εε****⎛⎫=⋅ ⎪⎝⎭。

（7）二元函数的绝对误差限：设一元函数()()(),0,f f x y f y y εε***⎛⎫∂==⋅ ⎪∂⎝⎭则。

（8）二元函数的相对误差限：()()()1r f f f x y x y f εεε******⎡⎤⎛⎫∂∂⎛⎫⎢⎥=⋅+⋅ ⎪ ⎪∂∂⎝⎭⎢⎥⎝⎭⎣⎦。

三本章习题解析1. 下列各数都是经过四舍五入得到的近似值，（1）试指出它们有几位有效数字，（2）分别估计1123A X X X ***=及224X A X **=的相对误差限。

12341.1021,0.031,385.6,56.430x x x x ****====解：（1）1x *有5位有效数字，2x *有2位有效数字，3x *有4位有效数字，4x *有5位有效数字。

（2）1111123231312123,,,,A A AA x x x x x x x x x x x x ∂∂∂====∂∂∂由题可知：1A *为1A 的近似值，123,,x x x ***分别为123,,x x x 近似值。

所以()()111rA A Aεε***=()()()12311111123A A A x x x A X X X εεε*******⎡⎤⎢⎥=++⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫⎛⎫∂∂∂ ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭43123131212311111010100.215222x x x x x x x x x **-**-**-***⎡⎤=⨯⨯+⨯⨯+⨯⨯=⎢⎥⎣⎦()222222424441,,,X A Ax A X x x x x ∂∂===-∂∂则有同理有2A *为2A 的近似值，2x *，4x *为2x ，4x 的近似值，代入相对误差限公式：()()222rA A Aεε***=()()24212224A A X X A X X εε*****⎡⎤⎢⎥=+⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫∂∂ ⎪ ⎪∂∂⎝⎭⎝⎭()33542224411*********X X X X X **--***⎡⎤⎢⎥=⨯⨯+⨯⨯=⎢⎥⎣⎦2. 正方形的边长大约为100cm ，怎样测量才能使其面积误差不超过21cm ? 解：设正方形的边长为x ，则面积为2S x =，2dsx dx=，在这里设x *为边长的近似值，S *为面积的近似值：由题可知：()()1ds s x dx εε***=≤⎛⎫ ⎪⎝⎭即：()21x x ε**⋅≤ 推出：()10.005200xcm ε*≤=。

计算方法的课后答案

《计算方法》习题答案第一章数值计算中的误差1．什么是计算方法？（狭义解释）答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。

2．一个实际问题利用计算机解决所采取的五个步骤是什么？答：一个实际问题当利用计算机来解决时，应采取以下五个步骤：实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值，并编程获得解。

解：400)(2345-+⋅+-⋅+=x x x x x x P ，从而 1 0 -1 0 1 -4 -3 -3 9 -24 72 -2191-38-2473-223所以，多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。

5．叙述误差的种类及来源。

答：误差的种类及来源有如下四个方面：（1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。

（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。

（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。

（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。

这样引起的误差称为舍入误差。

6．掌握绝对误差（限）和相对误差（限）的定义公式。

计算方法答案王能超

计算方法答案王能超【篇一：计算方法习题集及实验指导书】s=txt>计算机科学与技术系檀明2008-02-10课程性质及目的要求（一）课程性质自计算机问世以来，科学计算一直是计算机应用的一个重要领域，数值计算方法是解决各种复杂的科学计算问题的理论与技术的基础。

《计算方法》课程讨论用于科学计算中的一些最基本、最常用的算法，不但具有数学的抽象性与严密的科学性的特点，而且具有应用的高度技术性的特点。

它对于培养从事计算机应用的科技人才有着重要的作用，是计算机应用专业（本科段）的一门重要的技术基础课程。

（二）目的要求通过本课程的学习和上机实验，了解用计算机解决科学计算问题的方法特点，掌握计算方法中的一些基本概念、基本公式和相应的算法流程，提高根据算法描述设计高级语言程序并进行验证的技能。

在学习过程中，应注重理解和应用，在搞清基本原理和基本概念的基础上，通过习题、编程和上机等环节，巩固和加深已学的内容，掌握重要的算法及其应用。

注重理论与算法的学习和应用相结合，强调编程及上机计算的技能培养，是本课程不同于一般数学课程的重要特点。

（三）学习方法指导1．循序渐进逐章学习本课程从第二章开始，每章都讨论一个大类的算法。

虽然各算法是相对独立的，但是也存在相互联系与前后继承的关系。

前面的概念和算法学好了，后面的内容也就容易学，越学越感到容易。

前面的内容没有学好，后面就会感到难学，甚至会出现越来越感到困难、失去学习信心的情况。

2．稳扎稳打融会贯通学习要扎实、要讲求实效。

每一个重要的概念和公式，都会搞清楚，做到融会贯通。

只有这样，才能取得学习的学习效果。

3．多学练勤做习题教材及本习题集中的每一章都附有适量的习题，可以帮助考生巩固和加深理解所学的知识，提高解题能力。

因此，在学习过程中，应当适合习题进行思考，应当尽可能多做习题，遇到某些不会做的题，应三思之后再请老师给予提示。

4．抓住特点前后联系本课程只讲了五大类算法。

每类算法都是针对一类特定的计算问题，都有其自身的特点。

吉林大学工程数学计算方法第三章习题答案共7页word资料

第三章习题答案11. 分别用梯形公式、Simpson公式、Cotes公式计算积分I xdx,并估计误差。

0.5解：1)用梯形公式有：事实上，2) Simps on 公式事实上，E S(f) —[T Xdx —(0.5)+4f l f(1)]= 0.00003043) 由Cotes公式有：事实上，E C f = 0.00000032. 证明Simpson公式2.8具有三次代数精度。

证明: 而当f x =x4时b b 4 i 5 5左侧：f x dx x4dx b-a5'a 'a 5右侧：左侧不等于右侧。

所以Simps on具有三次代数精度.3•分别用复化梯形公式和复化公式Simps on计算下列积分.1 9 M r _____________________(1) dx, n=8 , (3) . xdx,n=4 , 仏4-sin2'd n = 604 x2 1 0解：(1)用复化梯形公式有:b - a 1_ 0 1h = ■n 8 8T n f a 210 2 (0.031128 0.061538 0.090566 0.11765 0.14235 0.16438 0.18361) 0.21 - 0.111416由复化Simpson公式有:解：删去解(3)：存dx,® 由复化梯形公式有: 由复化Simps on公式有:(4)解：\4_sin 2 d , n = 6 由复化梯形公式：由复化 Simps on 公式：113 14•给定求积节点x o ,x ),X 2,试推出计算积分.0f x dx 的插值型求积公式，并写出它的截断误差。

解:考虑到对称性，有氏=氏，于是有求积公式由于原式含有 3个节点，故它至少有 2阶精度。

考虑到其对称性，可以猜想到它可能有 3阶精度。

事实上，对 f =X 原式左右两端相等：此外，容易验证原式对 f =X 4不准确，故所构造出的求积公式有 3阶精度。

第三章计算题答案

1、假定经济中有C=100+0.75Yd(可支配收入),I=125-600r, G=50 ,T=20+0.2Y,TR=50.(1)推导IS曲线方程。

（2）IS曲线的斜率（3）当r=15%时，Y是多少？如果充分就业收入Yf=565，应增加多少政府支出才能实现充分就业？修改“TR=50”为“TR=0”答：（1）Y=C+I+G=[100+0.75（Y-20-0.2Y）]+ 125-600r+50=260+0.6Y-600r由此可知，r= -1/1500Y+650/1500 或Y=650-1500r（2）斜率为-1/1500（3）Y=650-225=425；与充分就业相比，Y的缺口是565-425=140，政府支出乘数是Kg=1/1-b（1-t）=1/1-0.75（1-0.2）=1/0.4=2.5，于是，政府应该增加的支出是140/2.5=56。

2、假定某经济中消费函数为C=0.8(1-t)Y,税率t=0.25,投资函数为：I=900-50r,政府购买G=800,货币需求为L=0.25Y+200-62.5r,实际货币供给位M/P=700,试求：(1) IS、LM曲线（2）两个市场同时均衡时的利率和收入。

答：（1）IS曲线：Y=C+I+G=0.8(1-t)Y+900-50r+800由此可知，Y=4250-125r。

LM曲线：由L=M可知，0.25Y+200-62.5r=700，Y=2000+250r。

（2）联立IS、LM曲线，可求出：Y=3500；r=6。

3、假定某经济中消费函数为C=100+0.9(1-t)Y,,投资函数为：I=200-50r,净出口为NX=100-0.12Y-500r,政府购买G=200,货币需求为L=0.8Y+200-2000r,税率t=0.2,名义货币供给位M=1000,价格水平不变，试求：(1) IS、LM曲线（2）两个市场同时均衡时的利率和收入。

（3）两个市场同时均衡时的消费、投资和净出口值。

北大版高等数学第三章积分的计算及应用答案第三章总练习题

第三章总练习题111121221.N ew to n -L eib n iz 1(1).[1,1],.tan (2).tan (0,2)2tan2.,x x xd de d x e e d x d x x d x d x u x xf F F ππ-⎛⎫⎛⎫⎛⎫=-- ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭=+⎰⎰为什么用公式于下列积分会得到不正确结果?无界从而不可积在的一些点不可导.证明奇连续函数的原函数为偶函数,而偶连续函数的原函数之一为奇函数.设奇连续函数的原函数为现在证明是偶证.()().(()())()()()()0,()(),(0)(0)0.()()0.,.()().(()())()()()()0,()().(0)0,(F x f x F x F x F x F x f x f x F x F x C C F F F x F x f F F F x f x F x F x F x F x f x f x F x F x C F C F ''''=--=---=---=--==--=--=''''=-+=--+=--+=--===函数设偶连续函数的原函数为现在证明是奇函数设则3003440010100)(0)0.()()0.sin ,0,3.()()()?0,0., 0,()()()sin co s |1co s .444.sin ().sin ()s b ab b b a aaabF F x F x x x f x f x f x d x a b x x f x d x f x d x f x d x x d x xd xxax b d x t d x d xd d x t d x d xd x--=-+=≥⎧==<>⎨<⎩=+=+=-=+-++=⎰⎰⎰⎰⎰⎰⎰⎰求定积分其中求微商解解()()110001201/2221210221in ()sin (1)sin ().5.lim()(),().1lim()()().6.lim(1).(2)!!(1)co s.(21)!!2x xh x h u xh u xnn nn n n u d u x x f x h t d x f x f x f x h t d u f t d tf x hx d x n x d x td t I n I π+→+→=→∞+++=+-+='+==--===+<⎰⎰⎰⎰⎰⎰⎰试证明其中是实轴上的连续函数求极限证解12210(2)!!(21)!!1,(21)!!(22)!!1100(),lim(1)0.sin co s 7..2sin 3co s sin co s (2sin 3co s )(2sin 3co s )(2sin 3co s )(2co s 3sin )(23)sin (32)co s nn n n n n n n I n x d x x x d x x xx x A x x B x x A x x B x x A B x A B +→∞+=+++<<→→∞-=+-'+=-+-=-++=++-+⎰⎰令解,x23115,,.3211313sin co s 2sin 3co s (2sin 3co s )(2sin 3co s )2sin 3co s ln |2sin 3co s |15ln |2sin 3co s |.1313A B A B A B x xd x x x A x x B x x d xx xA xB x xC x x x C +=⎧=-=⎨-+=⎩+=-'-+-=-=+-+=-+-+⎰⎰222228.:2(1),ln (2),.22222222.(2)(2)222xxu d u x u x u d x uu d u d u x u uu u C C xex x e xd==+=+⎛⎫==- ⎪++⎝⎭⎛⎛=-+=-+ ⎝⎝=-==⎰⎰⎰⎰⎰⎰⎰⎰通过适当的有理化或变量替换求下列积分()24.2).22(3)33.(4)(11ln .2x C x C x CC d x x C ⎛⎫=+⎝=-++==-⨯=-=+-==-+⎰⎰⎰⎰⎰224444224222244sec tan(1)9..sin co s1tan11112111.21111((arctan1)arctansin co sd x xd x u d ux x x uuuu ud xx x+==++++⎛==++⎝⎛⎫⎪⎪=+⎪⎛⎛⎪++-+⎪⎝⎝⎝⎭=+++⎰⎰⎰⎰)1).110.()(,),,()()(),:()()( 3.424)11.()[,],()0.:(,),()0.,()(,),,TxbaCf x Tg x f x f x d xTh x g t d t Tf x a b f x d x a b cf cf x a b ff-+-∞+∞=-===⎰⎰⎰设函数在上连续以为周期令证明函数也以为周期.此即习题第题设函数在区间上连续且证明在内至少存在一点使若不然在没有零点由的连续性和连续函数的中间值定理在证证(,).()0,(,).,,,[,]0.()()()()()0..b c d ba a c da b f x x a b c d a c d b f c dmf x d x f x d x f x d x f x d x m d c>∈<<<>=++≥->⎰⎰⎰⎰不变号不妨设取满足则在取最小值于是矛盾22222212.[,],()0,:()0,[,].[,][,],|()||()|,[,].2|()|()()()0.2.bab ea da b f x d x f x x a bd e a bf cf x x d ef cf x d x f x d x d e=≡∈∈≠⊆>∈≥>->⎰⎰⎰设函数f在区间上连续且证明若不然,存在c[a,b],f(c)0.由f在c的连续性,存在区间矛盾证00222/200013.()(-,),(1),()();sin(2);1co s4sin1(3)(1co s sin()()()()().a aa a aaf xa f x d x f a x d xx xd xxxd xx xf x d x x a t f a t d t f a t d t f a x d xπππ∞+∞=-=+=++=-=--=-=-⎰⎰⎰⎰⎰⎰⎰⎰设在上可积证明对于任意实数有证(1)22220211022222sin ()sin ()sin sin (2),1co s 1co s 1co s 1co s sin co s arctan |.21co s 21co s 14sin sin (/2)(3)co s sin co s(/2)s x x x x x x xI d x d x d x d x I x x xx xd xd x d uI u xxux x I d x x xx ππππππππππππππππ--==-==-++++==-===+++-==+-+⎰⎰⎰⎰⎰⎰⎰/2/20222/2/2/20/2/2/2in (/2)co s co s sin ,2co s sin co s sin co s sin 1|csc(/4)co t(/4)||co s sin 11ln1ln co s s 4d xx x xx d x I d x d xx xx xx xd x d xd x x x x xππππππππππππ==++++===+-++⎡⎤⎛⎫⎢⎥ ⎪=+-⎥ ⎪⎥ ⎪⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰⎰⎰11),in 411).I π⎛⎫⎪-= ⎪ ⎪⎝⎭=232221123222014.()(23)m /s .004m /s,(1);(2)3(1)()23,3,4,34,4,3230.()4.34(4)(1)0, 4.32(5)a t t t x v tx t t x t t C C x t t x t t C tC x t t t x t t t t t s x =-===-''''=-=-+-==--=--+'==--=--=-+===一质点作直线运动,其加速度若时且求质点改变动方向的时刻头5秒钟内质点所走的总路程.解3322543343(4)|(4)|424m.32322t t t t x x t t t t ==⎛⎫⎛⎫-+=-----=⎪⎪⎝⎭⎝⎭000200022002200215.100m ,10.2s,25m ,25m ., 0;(), 10.2., 0;()210.2./2/2253m /s .10010.2a t t t v t a t t t a tt t s t a t t C t t a t a t Ca t a a t C ≤≤⎧=⎨≤≤⎩⎧≤≤⎪=⎨⎪+≤≤⎩⎧=+⎪=≈⎨⎪=+⎩一运动员跑完共用了在跑头时以等加速度进行然后保持等速运动跑完了剩余路程.求跑头时的加速度解16.(1):利用积分的几何意义证明111111ln,1,2,111(2)1ln ,211111ln ,21,.111(3)lim 1ln E u ler 211ln |111ln (1)ln lnn n n n n n n n n nnn n n n nx n n y n n nx y n n n d x d x x n n xn n n n→∞++++<<=+=+++--=++++--⎛⎫++++- ⎪-⎝⎭=<=+++=+-=<⎰⎰令证明序列单调上升而序列单调下降证明极限存在(此极限称为常数).证 (1)1121.1111(2)1ln (1)1ln 22111ln 10((1)).111111ln (1)1ln 21211ln 10((1)).1(3)1ln 20(2)n n n n n n d x n n x x n n n n n n y y n n n n n n n y x x n ++=⎛⎫⎛⎫-=+++-+-+++- ⎪ ⎪-⎝⎭⎝⎭⎛⎫=-+> ⎪⎝⎭⎛⎫⎛⎫-=++++-+-+++- ⎪ ⎪+⎝⎭⎝⎭⎛⎫=-+< ⎪+⎝⎭>>=->>⎰由由n +1n11/22111/1/2222112.,lim .17.:0,11.111111(1/).111/118.()(,),,(2)().2,()0.(2)()n n n xxxxxay y x d t d t ttd t x u d x d t tuutf x x f x f x a f x d x f x f x →∞>=++==⨯=+++-∞+∞-=-≠=-=-⎰⎰⎰⎰⎰⎰单调下降有下界故有极限证明当时设在上连续(书上为可积，欠妥)且对一切实数均有求实数使条件证解(22220221(11)(11))()(2)(),()0.0.19.ln (1)arctan ,0 1.11,[0,1],[0,],1111ln (1)arctan ,0 1.20.(1)x x f x f x f x f x d x f u d u f u d u f x d x a x x x d t d t t x ttttx x x =+-=-+-=-=-==+≤≤≤≤∈≤+++++≤≤≤⎰⎰⎰⎰⎰⎰相当关于为奇函数取即可利用定积分的性质，证明不等式在上积分得设证()()[0,];()()2(2)a f x d x a f x a d x f x f a x =+-⎰在上可积，证明利用(1)中的公式求下列积分的值：22/2222sin ;22sin co s ()()(1)()()()()()()2()()()()()(-)1,.()()()()22a a a a a a a xx d x d xx x x xf x f a u I d x d uf x f a x f u f a u f x f a u I d x d uf x f a x f u f a u f x f a x a d x d x d x a I f x f a x f x f a x xx x π-++-==+-+--=++-+-=+===+-+--⎰⎰⎰⎰⎰⎰⎰⎰⎰证解(2)22222/20222 2.2(2)2sin /2.sin sin (/2)24xd x d x x x xd x x x ππππ==⨯=++-==+-⎰⎰⎰tan 2sin tan tan 22sin sin tan 22222sin tan 32222sin 2()21.()(1).()(1)tan sin ,()sec co s tan sec sin co s sec co s tan sec sin co s 3sec x xx x x xx xxxd f x f x xt d td x d xf x xt d t x x x t d t d f x x x x x x x x x t d td xt x x x x x x x x =+=+=-+=-+-+=-+-+==⎰⎰⎰⎰设求解()()22233222331co s tan sec sin co s tansin 31sec (1tan )co s (1sin )tansin .3x x x x x x x x x x x x x x x x x x -+-+-=+-++-()/22/2/22/220/20/2/2022.co s 3.11co s 3(1co s 6)sin 6|.2412423.|sin co s |.|sin co s ||sin co s ||sin co s ||sin co s ||sin ()2I d I d d I x x d x x x d xx x d x x x d x x x d x t ππππππππππθθππθθθθθπ===+=+==--=-+--++⎰⎰⎰⎰⎰⎰⎰⎰求定积分的值求定积分的值解解I =22=2()()()()/20/2/20/4/2/20/4/4/2/2/40co s()|2|sin co s ||co s sin |(co s sin )(sin co s )(co s sin )2sin co s |(co s sin )|(sin co s )|t d x x x d x t t d xx x d x x x d x t t d xx x x x x x ππππππππππππ⎛⎫-+ ⎪⎝⎭=-++=-+-++=++--+-=⎰⎰⎰⎰⎰⎰221010101010100110()/224.0,.2x x x x x x x x x x x x x x I x I xxx x x x u x --<<====+⎛⎫==- ⎪⎝⎭=⎰⎰⎰⎰⎰设求定积分的值解10()/210222102()2arcsin().28x x a aux x x a u a a x x aππ--==⎡⎤=⎢⎥⎣⎦-==⎰⎰4342(2)16468 4.y x x x y x x x =++-=++-与43323242224342442432164164684,680,6840,680,(2)(4)0,2,4.{(164)(684)][(684)(164)]y x x x x x x x x x x y x x x x x x x x x S x x x x x x d xx x x x x x d x⎧=++-⎪+-=+--+=⎨=++-⎪⎩=-+=--===++--++-+++--++-⎰⎰解222212552233532325.:(1)6827.682768,278150,(3)(5)0.3, 5.(27(68))(815)4415.33y x x y x y x x x x x y x x x x x x x S x x x d x x x d xx x x =-+=-⎧=-+-=-+⎨=-⎩-+=--====---+=-+-⎛⎫=-+-=⎪⎝⎭⎰⎰求下列曲线所围图形的面积与解/2/2/4/45/45/4/2/2(4)sin ,cos /2.(sin -cos )(cos sin )|1;(sin cos )(cos sin )| 1.y x y x x S x x dx x x S x x dx x x πππππππππ=====--==-=--=⎰⎰与解/21102/211226.co s ,1/2,,.,)2arcco s 2arcsin 2(1co s 2).242arcsin arcsin y x y x x V V V x d x y y d y y yd y V x d x V y yd y yd yππσπσππππππππ===⎛⎫=-==⎪⎝⎭=-====⎰⎰⎰⎰⎰⎰设区域由曲线及所围成将绕轴旋转一周得一旋转体试用两种不同的积分表示体积并且求的值. 2解V =(1-c o s 24323222444323202(68)(68)24248.44x x x d x x x x d xx x x x x x =-++-+-⎡⎤⎡⎤=-++-+-=⎢⎥⎢⎥⎣⎦⎣⎦⎰⎰2222212321(3)1 3.(3)1,7100,(2)(5)0,2,5.1,2.[(3)(1)]92.322y x y x x x x x x x x y S y y d x yy y --=-=--=--+=--===-=+-+⎛⎫=-++= ⎪⎝⎭⎰与解4π2π54π112212222arcsin ()1arcsin .22244y y y xyπππππππ=-⨯⎤=-=-=⎥⎦⎰22111/21/2200213397200257025027.:(1)arcsin |.4612(2)(918055801)400.28.()[0,7],()5,()6,() 3.(1)();(2)x xxx d x f x f x d x f x d x f x d x f x d x f πππ-===-=-++====⎰⎰⎰⎰⎰⎰求下列定积分的值设在上可积且一直已知求的值求7552500277557755().(3):(5,7),()0.(1)()()()5611.(2)()()()3118.(3),()0,(5,7),()0,()80,.x d x f x f x d x f x d x f x d x f x d x f x d x f x d x f x x f x d x f x d x <=+=+==-=-=-≥∈≥=-<⎰⎰⎰⎰⎰⎰⎰⎰⎰的值证明在内至少存在一点使若不然但是矛盾解证2/23/21/2/2/2/2/232111, 2,129.()sin ,(),()2, 2.(1)()();(2)()();(3)()().(1)()()sin 0.(2)()()()()()xx f x x h x g x x x f x g x d x g x h x d x f t g t d x f x g x d x xd x g x h x d x g x h x d x g x πππππππππ----≤≤⎧===⎨<≤⎩===+⎰⎰⎰⎰⎰⎰⎰设试求下列定积分的值或表达式:解322323221212/22/2/22()12125.6sin co s ,2(3)()()sin 2sin co s 22co s ,2.xxx h x d x d x d x xxxxtd t x t x f t g t d x td t td t x x πππππ=+=--=⎧=--≤≤⎪=⎨⎪+=-<≤⎩⎰⎰⎰⎰⎰⎰⎰()()30()[,](0),()()()()()().b f b af a f x a b ag y f x f x d x b f b a f a g y d x >=--⎰⎰设函数在区间上连续，严格单调递增是的反函数,利用定积分的几何意义证明下列公式并作图解释这一公式.解11()10031.(1)()[0,),()0,0()()()()()()()(*).00(0),()a B a aa x x x a B a B x d x x d xx x x d x x d x a a B B x x ϕϕϕϕϕϕϕϕϕϕϕϕϕ---+∞→∞→+∞≥≥≤++==>=→∞→+∞⎰⎰⎰⎰⎰设函数在上连续且严格单调递增又设当+时且(0)=0.证明:对于任意实数,下列不等式成立:其中是的反函数.由题时不等式显然成立.设由于+时证30,111101100()11()1()1,0,(),[0,],,0,().,*)()().,()()()()()()()()(())(a B aBa a B a Ba a a B a a a B a a a B x d x x d x a a x d x x d xx d x x d x x d xa a x d xa a a B ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ------'>''>>===+>+=++=++≥+⎰⎰⎰⎰⎰⎰⎰⎰存在在连续根据连续函数的中间值定理存在若则由(得若则1100()()11()11()).,()()()()()()()()(())(()).11(2)(1),0,,11,M in k o w sk i.1aB a a a Ba B pqa a B a a x d x x d xx d x x d x x d xa a x d xa a a a B a B ab p q p qaba b pqp ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ------=<+=+-=-≥--=≥≥+=≤+>⎰⎰⎰⎰⎰⎰若则利用中的不等式,对于任意实数证明下列不等式不妨设证111/(1)1/(1)1/(1)1/0.(),()..1/(1)1/(1)p p pp pp p pqa b ppx x x x abababa b x d x xd x pp pp p pqϕϕ----+-==≤+=+=+=+-+-⎰⎰在中取则(1)232.0,,1,112.xa a y y x a y π>=+===设求的值使由曲线及所围成的区域绕直线旋转所得之旋转体的体积等于b22222202222220002.)2,112,22,2,18,2ln 9,44aa x aa a x x u a d x d x xe d x e d x e d u e a a ππ=====-===⎰⎰⎰⎰⎰20解(y -1) 20033.1sin 21co s 43(1sin 2)(12sin 2).22r S d d ππθθπθθθθ=+-=+=++=⎰⎰作由极坐标方程所确定的函数的图形,并求它所围区域的面积.解。

计算方法各章作业答案

计算方法作业第二章插值1.（1（2）用二次Lagrange插值多项式求当X=0.15时Y的近似值。

（3）写出余项R(x)=f(x)-Pn(x)的表达式。

解：（1）Pn (x) =knknkjj jkj yxxxx)(00∑∏=≠=--n=3P 3(x)=321321))()(())()((yxxxxxxxxxxxx------+13121132))()(())()((yxxxxxxxxxxxx------+23212231))()(())()((yxxxxxxxxxxxx------+32313321))()(())()((yxxxxxxxxxxxx------x 0=0.0 x1=0.1 x2=0.2 x3=0.3y 0=0.0000 y1=0.0998 y2=0.1987 y3=0.2955P 3(x)=0000.0)3.00.0)(2.00.0)(1.00.0()3.0)(2.0)(1.0(⨯------xxx+0998.0)3.01.0)(2.01.0)(0.01.0()3.0)(2.0)(0.0(⨯------xxx+1987.0)3.02.0)(1.02.0)(0.02.0()3.0)(1.0)(0.0(⨯------xxx+2955.0)2.03.0)(1.03.0)(0.03.0()2.0)(1.0)(0.0(⨯------xxx(2) y(0.15) = P2(0.15) = 0.1494（3）R(x) = f(x)-Pn (x)=)!1()()1(++nf nξnk0=∏(x - x k)=!4)(4ξf(x – 0.0) (x – 0.1)(x – 0.2)(x – 0.3)第三章方程求根5．求解方程12-3x+2cosx=0的迭代法n n x x cos 3241+=+（1）证明对于任意的x 0€R 均有*lim x x n x =∞→ (x *为方程的根)（2）取x 0=4，用此迭代法求方程根的近似值，误差不超过10-３，列出各次的迭代值。

计算方法答案第三章

第三章插值法与最小二乘法1. 已知下列表值x 10 11 12 13 lnx 2.3026 2.3979 2.4849 2.5649用线形插值与二次Lagrange 插值计算ln11.75的近似值，并估计误差。

解：（1）线形插值说明：当插值点落在被插区间之内，这种方法称为内插法，此时插值精度较好。

x ],12,11[75.11∈=故选择x 0=11,x 1=12,求线形插值函数。

11001y x l y x l x P ⨯+⨯=∴）（）（）（=10100101y x x x x y x x x x ⨯--+⨯--=4849.21112113979.2121112⨯--+⨯--x x=2.4849（x-11）-2.3979（x-12）)1275.11(3979.2)1175.11(4849.2)75.11(75.11ln 1---=≈∴p =2.46315（2）二次拉格朗日插值选择插值结点：x 12,11,10210===x x P 2211002)()()()(y x l y x l y x l x ++= =212021012101200201021))(())(())(())(())(())((y x x x x x x x x y x x x x x x x x y x x x x x x x x ----+----+----=4849.2)1112)(1012()11)(10(3979.2)1211)(1011()12)(10(3026.2)1210)(1110()12)(11(----+----+----x x x x x x=1.1513（x-11）（x-12）-2.3979（x-10）（x-12）+1.24425（x-10）（x-11）)1175.11)(1011075(24245.1)1275.11)(1075.11(3979.2)1275.11)(1175.11(1513.1)75.11(75.11ln 2--+-----=≈∴P =1.15133125.124245.14375.03979.2)1875.0(⨯+⨯+-⨯ =2.4639282. 已知下列表值求f(x)在[0，2]之间零点近似值。

精选吉林大学工程数学计算方法第三章习题答案资料

第三章习题答案1.分别用梯形公式、Simpson公式、Cotes公式计算积分1,I=⎰并估计误差。

解：1）用梯形公式有：()()110.51[10.5]10.42678242f f⎛-≈+=+≈⎝⎭⎰()()()333333220.512.6042107.36571012124Tb aE f fηηη-----⎛⎫''=-=--=⨯≤⨯⎪⎝⎭事实上，()()()()()()110.430964410.50.510.4267767210.50.510.00418772Tf x II f fE f f f===-≈+=⎡⎤⎣⎦-∴=-+=⎡⎤⎣⎦⎰⎰2）Simpson公式()110.53111410.43093 64212f f f⎛-⎡⎤⎛⎫⎛⎫≈++=+=⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦⎝⎭⎰[]()()44744211111522 1.1837710180218028Sb a b aE f fηη--⎛⎫--⎪⎛⎫--⎛⎫=-=--≤⨯⎪ ⎪⎪⎝⎭⎝⎭⎪⎝⎭()312()''()48T f fb aE h=?--事实上，()()()110.50.510.5410.000030462SE f f f f-⎡+⎤⎛⎫=-++=⎪⎢⎥⎝⎭⎣⎦⎰3）由Cotes公式有：()() ()111537270.5321232719084814.9497525.2982210.3923029.9332670.43096180f f f f f-⎡⎤⎛⎫⎛⎫⎛⎫≈++++⎪ ⎪ ⎪⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦=++++=⎰17)180+++()6116211294522 2.697410945464C E f η--⎛⎫⨯ ⎪⎛⎫=-⨯-≤⨯ ⎪ ⎪⎝⎭⎪⎝⎭()7(6)945*4()()82Cf b a E f h =?-- 事实上，()0.0000003C E f =2．证明Simpson 公式()2.8具有三次代数精度。

数值计算引论(第二版)三四五章习题解答

close all clear all clc n=10; x=zeros(n+1,1); for k=1:n+1 x(k)=cos((2*k-1)*pi/2/(n+1)); end y=1./(1+25*x.^2); x0=-1:0.1:1; y0=interp1(x,y,x0,'spline'); plot(x0,y0,'r')
h2 1.5, h3 0.5, h4 1.5, h5 0.5
b [0
h2 h3 3 h3 A 6 0 h3 6 h3 h4 3 h4 6
2
0]
0
T
2 3 h4 1 12 6 h4 h5 0 3
（d）样条函数插值具有比较好的数值稳定性。 √
习题
3.以0.1，0.15，0.2为插值节点，计算 f ( x ) x 的二次Lagrange插值多项式 P2 ( x ) ，比较 P2 (0) 和 f (0) ，问定理4.1的结果是否适用于本问题。解答：首先构造二次Lagrange插值多项式
R=chol（A）
0 -0.8165 1.1547 0 0 0 -0.8660 1.1180
-0.7071 1.2247 0 0
方法2：利用Cholesky定义求解
6.矩阵
1 A1 1 2 2 1 2
2 2 1 , A2 2 1 1
2 2 0
(B) 0 1
2 1 2
Gauss-Seidel迭代
0 1 M (D L) U 0 0
(M ) 2 1

计算机组成与原理第三章答案

第三章3.1十进制数化成地进制数和八进制数（无法精确表示时，二进制数取3位小数，八进制取1位小数）。

7＋3/4,±3/64,73.5,725.9375,25.34解：（1）、（7+3/4）：（7）10=（111）2（3/4）10=（0.00）2∴（7+3/4）10=（111.11）=（7.6）8（2）、（±3/64）10=（±0.）2=（±0.03）8（3）、（73）10=64+8+1=（），（0.5）10=（0.1）2∴（73.5）10=（.1）2=（111.4）8（4）、（725）10=512+128+64+16+4+1=（）2=（1325）8（0.9375）10=（0.1111）2=（0.74）8∴（725.9375）10=（.1111）2=（1325.74）8（5）、（25）10=（11001）2=（31）8（0.34）10=（0.011）2=（0.3）8∴（25.34）10=（11001.011）2=（31.3）83.2 把下列各数化成十进制数：（101.10011）2,(22.2)8,(AD.4)16,解：（1）（101.10011）2＝２2＋２0＋２－１＋２-4＋２-5=(5.59375)10=(5+19/32)10（2）(123.123)4=42+2*41+3+4-1+2*4-2+3*4-3=(27.)10 =(27+27/64)10=(1+1/4)*(4+2*4+3)=27*(1/64+1)(3) (22.2)8=2*81+2+2*8-1=(10.25)10=(10+1/4)10(4)(AD.4) 16=10*16+13+4*16-1=(173.25)10=(173+1/4)10(5) (300.3)8=3*82+3*8-1=(192.375)10=(192+3/8)103.3 完成下列二进制运算：101.111+11.011,1001.10-110.01,101.11*11.01,÷11013.4写出下列各地进制数的原码、补码和反码：0.1010,0,-0,-0.1010,0.1111,-0.0100答：x [x]原[x]补[x]反0.1010 0.1010 0.1010 0.10100 0.0000 0.0000 0.0000-0 1.0000 0.0000 1.1111-0.1010 1.1010 1.0110 1.01010.1111 0.1111 0.1111 0.1111-0.0100 1.0100 1.1100 1.10113.5 已知[X]原为下述各值，求[X]补:0.10100,1.10111,1.10110答：[x]原0.10100 1.10111 1.10110[x]补0.10100 1.01001 1.010103.6 已知[X]补为下述各值，求X（真值）：0.1110,1.1100,0.0001,1.1111,1.0001答：[x]补0.1110 1.1100 0.0001 1.1111 1.0001x 0.1110 -0.0100 0.0001 -0.0001 -0.11113.7已知X＝0.1011,Y= -0.0101,试求：[X]补，[-X]补,[Y]补，[-Y]补,[X/2]补,[X/4]补,[2X]补,[Y/2]补,[Y/4]补,[2Y]补,[-2Y]补答：[x]补=0.1011; [-x]补=1.0101; [y]补=1.1011; [-y]补=0.0101;[x/2]补=0.0101(1); [x/4]补=0.0010(11); [2x]补=1.0110(溢出)；[y/2]补=1.1101(1); [y/4]补=1.1110(11); [2y]补=1.0110; [-2y]补=0.10103.8 设十进制数X＝（＋128.75）＊2－10（1）若（Y）2＝（X）10，用定点数表示Y值。

计算方法答案 第三章

计算方法第三章习题答案

第三章物理学中定积分的数值计算方法

西北工业大学计算方法作业集答案及试题

计算方法-刘师少版课后习题答案

数值计算方法课后习题答案

计算方法各章习题及答案

吉林大学工程数学计算方法第三章习题答案word资料9页

数值计算方法(宋岱才版)课后答案

计算方法的课后答案

计算方法答案王能超

吉林大学工程数学计算方法第三章习题答案共7页word资料

第三章计算题答案

北大版高等数学第三章 积分的计算及应用答案 第三章总练习题

计算方法各章作业答案

计算方法答案 第三章

精选吉林大学工程数学计算方法第三章习题答案资料

数值计算引论(第二版)三四五章习题解答

计算机组成与原理第三章答案

计算方法答案第三章

北大版高等数学第三章积分的计算及应用答案第三章总练习题

计算方法答案第三章