计算方法 第四章 参考答案
- 格式:doc
- 大小:80.00 KB
- 文档页数:2
第四章作业参考答案4. 用推广的Euclid算法求67 mod 119的逆元解:初始化:(1,0,119), (0,1,67)1:Q=119/67=1,(0,1,67) , (1,-1,52)2:Q=67/52=1,(1,-1,52), (-1,2,15)3:Q=52/15=3,(-1,2,15), (4,-7,7)4:Q=15/7=2,(4,-7,7), (-9,16,1)所以67-1 mod 119=1610.设通信双方使用RSA加密体制,接收方的公开钥是(e,n)=(5,35),接收到的密文是C =10,求明文M。
解:由n=35,易知35=5×7,进而ϕ(n)=ϕ(35)=24,由RSA加密体制可知,ed≡1 mod ϕ(n),即5d≡1 mod 24,所以d=5∴M=C d mod n=105 mod 35=511. 已知c d mod n的运行时间是O(log3n),用中国剩余定理改进RSA的解密运算。
如果不考虑中国剩余定理的计算代价,证明改进后的解密运算速度是原解密运算速度的4倍。
证明:RSA的两个大素因子p,q的长度近似相等,约为模数n的比特长度log n的一半,即(log n)/2,而在中国剩余定理中要计算模p和模q两个模指数运算,与c d mod n的运行时间规律相似,每一个模指数运算的运行时间仍然是其模长的三次幂,即O[((log n)/2)3]= O(log3n)/8,这样在不考虑中国剩余定理计算代价的情况下,总的运行时间为两个模指数的运行时间之和,即O(log3n)/8+O(log3n)/8=O(log3n)/4,得证。
12. 设RSA加密体制的公开钥是(e,n)=(77,221)。
(1) 用重复平方法加密明文160,得中间结果为1602(mod 221)=185,1604(mod 221)=191,1608(mod 221)=16,16016(mod 221)=35,16032(mod 221)=120,16064(mod 221)=35,16072(mod 221)=118,16076(mod 221)=217,16077(mod 221)=23,若敌手得到以上中间结果就很容易分解n,问敌手如何分解n解:由以上中间结果得16016(mod 221)=35=16064(mod 221),此即16064-16016=0 (mod 221)即(16032-1608) (16032+1608)=0 (mod 221)(120-16)(120+16)=0 (mod 221)104×136=0 (mod 221)由gcd(104,221)=13及gcd(136,221)=17,可知221的分解为221=13×17(2) 求解密密钥dd=e-1mod ϕ(221)=77-1 mod 12×16由扩展Eucild算法可得d=5。
第四章网络层1.网络层向上提供的服务有哪两种?是比较其优缺点。
网络层向运输层提供“面向连接”虚电路(Virtual Circuit)服务或“无连接”数据报服务前者预约了双方通信所需的一切网络资源。
优点是能提供服务质量的承诺。
即所传送的分组不出错、丢失、重复和失序(不按序列到达终点),也保证分组传送的时限,缺点是路由器复杂,网络成本高;后者无网络资源障碍,尽力而为,优缺点与前者互易2.网络互连有何实际意义?进行网络互连时,有哪些共同的问题需要解决?网络互联可扩大用户共享资源范围和更大的通信区域进行网络互连时,需要解决共同的问题有:不同的寻址方案不同的最大分组长度不同的网络接入机制不同的超时控制不同的差错恢复方法不同的状态报告方法不同的路由选择技术不同的用户接入控制不同的服务(面向连接服务和无连接服务)不同的管理与控制方式3.作为中间设备,转发器、网桥、路由器和网关有何区别?中间设备又称为中间系统或中继(relay)系统。
物理层中继系统:转发器(repeater)。
数据链路层中继系统:网桥或桥接器(bridge)。
网络层中继系统:路由器(router)。
网桥和路由器的混合物:桥路器(brouter)。
网络层以上的中继系统:网关(gateway)。
4.试简单说明下列协议的作用:IP、ARP、RARP和ICMP。
IP协议:实现网络互连。
使参与互连的性能各异的网络从用户看起来好像是一个统一的网络。
网际协议IP是TCP/IP体系中两个最主要的协议之一,与IP协议配套使用的还有四个协议。
ARP协议:是解决同一个局域网上的主机或路由器的IP地址和硬件地址的映射问题。
RARP:是解决同一个局域网上的主机或路由器的硬件地址和IP地址的映射问题。
ICMP:提供差错报告和询问报文,以提高IP数据交付成功的机会因特网组管理协议IGMP:用于探寻、转发本局域网内的组成员关系。
5.IP地址分为几类?各如何表示?IP地址的主要特点是什么?分为ABCDE 5类;每一类地址都由两个固定长度的字段组成,其中一个字段是网络号net-id,它标志主机(或路由器)所连接到的网络,而另一个字段则是主机号host-id,它标志该主机(或路由器)。
第4章练习单项选择题:1.在Excel 2003中,关于图表中数据系列的长度和源数据表中的相应数据的动态关系,以下描述不正确的是()。
A、鼠标拖动改变数据系列的长度,不会影响源数据表中的数据变化B、图表中的系列长度增加,源数据表的数值大小也增加C、源数据表中的数值增大,图表中的系列长度将变长D、删除图表中某一个数据系列,源数据表中对应的行数据不会自动删除参考答案:A2.在Excel 2003的页面设置的“页面”选项卡中,不可以设置()。
A、打印的居中方式B、缩放比例C、打印质量D、纸张的方向参考答案:A3.在Excel 2003中,打开两个工作簿BOOK1和BOOK2,在工作簿Book1中引用工作簿Book2的Sheet1工作表中的第3行第5列单元格,可表示为()。
A、[Sheet2]$E$3B、Sheet1!E3C、Book1:Sheet1!E3D、[Book2]Sheet1!E3参考答案:D4.Excel工作簿存盘时默认的文件扩展名为()。
A、.SLXB、.XLSC、.DOCD、.GZB参考答案:B5.在Excel 2003中可以依据字母、数字或日期等数据类型按一定顺序进行排列,这种操作称为()A、筛选B、分类排序C、分类汇总D、排序参考答案:D6.在Excel中,A1单元格设定其数字格式为整数,当输入“33.51”时,显示为()。
A、33.51B、33C、34D、ERROR参考答案:C7.在Excel 2003的某个单元格中输入文字,若文字长度较长而列宽为固定值时利用“单元格格式”对话框中的()选项卡,实现文字自动换行。
A、数字B、对齐C、字体D、图案参考答案:B8.Excel是微软Office套装软件之一它属于()软件。
A、电子表格B、文字输入C、公式计算D、公式输入参考答案:A9.在Excel 2003中,运算符运算的优先级不同,如要修改计算顺序,则()。
A、应把公式中需要首先计算的部分括在( )内B、应把公式中需要首先计算的部分括在[ ]内C、应把公式中需要首先计算的部分括在< >内D、应把公式中需要首先计算的部分括在{ }内参考答案:A10.在Excel 2003中,当用户希望标题文字能够相对于表格居中时,以下操作正确的是()A、填充B、合并及居中C、分散对齐D、居中参考答案:B11.Excel 2003中数据清单最多可以同时使用()个关键字进行排序A、4B、5C、3D、2参考答案:C12.不能够对Excel 2003的工作表进行()操作。
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 。
1.*x 为精确值x 的近似值;()**x f y =为一元函数()x f y =1的近似值;()**,*y x f y =为二元函数()y x f y ,2=的近似值,请写出下面的公式:**e x x =-:***r x xe x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂ 2、 计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫 舍入误差 。
3、 分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有6 位和7 1.73≈(三位有效数字)-211.73 10 2≤⨯。
4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为 0.0055 。
5、 设121.216, 3.654x x ==均具有3位有效数字,则12x x +的误差限为 0.01 。
6、 已知近似值2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为 0.0000204 .7、 递推公式,⎧⎪⎨⎪⎩0n n-1y y =10y -1,n =1,2,如果取0 1.41y =≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3*=π,则近似值141.3*1=π和1415.3*2=π分别有 3 位和 4 位有效数字。
9、 若*2.71828x e x =≈=,则x 有 6 位有效数字,其绝对误差限为1/2*10-5。
10、 设x*的相对误差为2%,求(x*)n的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字;12、计算方法主要研究( 截断 )误差和( 舍入 )误差;13、为了使计算 ()()2334610111y x x x =++---- 的乘除法次数尽量地少,应将该表达式改写为11,))64(3(10-=-++=x t t t t y ,为了减少舍入误差,应将表达式19992001-改写为199920012+。
计算方法作业集及答案第一章数值计算基本常识一.填空题1.用四舍五入得到的近似数0.628,有_____位有效数字,其绝对误差限是____________。
2.用四舍五入得到的近似数0.586,有_____位有效数字,其绝对误差限是____________。
3.用四舍五入得到的近似数0.69,其绝对误差是__________,由此计算出的相对误差限是__________。
4.用四舍五入得到的近似数0.7960,其绝对误差是__________,由此计算出的相对误差限是__________。
5.设0.484是0.4900的近似值,那么0.484具有____位有效数字。
6.设某某=0.231是真值某=0.229的近似值,则某某有_____位有效数字。
7.设某某=0.23是真值某=0.229的近似值,则某某有_____位有效数字。
8.设某=2.3149541,取5位有效数字,则所得的近似值某某=_____。
9.设某=2.3149541,取4位有效数字,则所得的近似值某某=_____。
10.若近似数0.1100有4位有效数字,由有效数字计算出的相对误差是____________。
11.若近似数76.82有4位有效数字,由有效数字计算出的相对误差是____________。
12.若近似数576.00有5位有效数字,由有效数字计算出的相对误差是____________。
13.用3.15作为π的近似值有_____位有效数字。
14.用3.14作为π的近似值有_____位有效数字。
15.用3.1416作为π的近似值有_____位有效数字。
解答:1.3、0.5某10-32.3、0.5某10-33.0.5某10-2、0.725%4.0.5某10-4、0.00628%5.16.27.28.2.31509.2.31510.0.05%11.0.007%12.0.001%13.214.315.5二.选择题1.3.141580是π的近似值,有()位有效数字。
第二章 数值分析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]-上的3次最佳一致逼近多项式()p x .答案:()f x 的最佳一致逼近多项式为323()74p x x x =++. 3.4 用幂级数缩合方法,求() (11)x f x e x =-≤≤上的3次近似多项式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 确定下列求积公式中的待定参数,使得求积公式的代数精度尽量高,并指明所确定的求积公式具有的代数精度. (1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰(2)11211()[(1)2()3()]3f x dx f f x f x -≈-++⎰ (3)20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-⎰答案:(1)具有三次代数精确度(2)具有二次代数精确度(3)具有三次代数精确度.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 =⎰(1)运用复化梯形公式计算上述积分值,使其截断误差不超过31102-⨯. (2)取同样的求积节点,改用复化辛浦生公式计算时,截断误差是多少?(3)要求的截断误差不超过610-,若用复化辛浦生公式,应取多少个节点处的函数值? 答案:(1)只需7.5n ≥,取9个节点,0.946I ≈(2)4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=⨯ (3)取7个节点处的函数值.4.6 用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分10sin xI dx x =⎰.要求用事后误差估计法时,截断误不超过31102-⨯和61102-⨯. 答案:使用复化梯形公式时,80.946I T ≈=满足精度要求;使用复化辛浦生公式时,40.946 083I s ≈=满足精度要求.4.7(1)利用埃尔米特插值公式推导带有导数值的求积公式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 ηη-=∈. (2)利用上述公式推导带修正项的复化梯形求积公式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 =.第六章 解线性代数方程组的迭代法6.1 对方程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 =.6.2 讨论松弛因子 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 =-.6.3 给定线性方程组Ax b =,其中111221112211122A ⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭,证明用雅可比迭代法解此方程组发散,而高斯-赛得尔迭代法收敛.6.4 设有方程组112233302021212x b x b x b -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭,讨论用雅可比方法和高斯-赛得尔方法解此方程组的收敛性.如果收敛,比较哪种方法收敛较快.答案:雅可比方法收敛,高斯-赛得尔方法收敛,且较快.6.5 设矩阵A 非奇异.求证:方程组Ax b =的解总能通过高斯-赛得尔方法得到.6.6 设()ij n nA a ⨯=为对称正定矩阵,对角阵1122(,,,)nn D diag a a a =.求证:高斯-赛得尔方法求解方程组1122D AD x b --=时对任意初始向量都收敛.第七章 非线性方程求根例7.4 对方程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 ≈=.例7.6 对于迭代函数2()(3)x x c x ϕ=+-,试讨论:(1)当c 为何值时,1()k k x x ϕ+=产生的序列{}k x(2)c 取何值时收敛最快?(3)取1,2c =-()x ϕ51||10k k x x -+-<.答案:(1)(c ∈时迭代收敛.(2)c =时收敛最快.(3)分别取1, 2c =--,并取0 1.5x =,计算结果如下表7.7所示表7.7例7.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 ''=' 例7.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 ⎛⎫ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭的最小特征值,使它至少具有2位有效数字.答案:取5 2.234375λ≈-即有2位有效数字.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 =,小数点后保留8位.答案: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 ≈=⎧⎨≈=⎩。