《2022年数值计算方法期末复习答案终结版》由会员分享,可在线阅读,更多相关《2022年数值计算方法期末复习答案终结版(10页珍藏版)》请在金锄头文库上搜索。
1、一、名词解释1误差:设*x为准确值x的一个近似值, 称*()e xxx为近似值*x的绝对误差, 简称误差。2有效数字: 有效数字是近似值的一种表示方法,它既能表示近似值的大小,又能表示其精确程度。如果近似值*x的误差限是1102n,则称*x准确到小数点后 n 位,并从第一个不是零的数字到这一位的所有数字均称为有效数字。3. 算法: 是指解题方案的准确而完整的描述,是一系列解决问题的清晰指令,算法代表着用系统的方法描述解决问题的策略机制。计算一个数学问题,需要预先设计好由已知数据计算问题结果的运算顺序,这就是算法。4. 向量范数: 设对任意向量nxR,按一定的规则有一实数与之对应,记为|x,若|
2、x满足(1)| 0x,且|0x当且仅当0x;(2)对任意实数,都有| |x|x;(3)对任意,nx yR,都有| |xyxy则称|x为向量 x的范数。5. 插值法: 给出函数( )f x的一些样点值,选定一个便于计算的函数形式,如多项式、分段线性函数及三角多项式等,要求它通过已知样点,由此确定函数( )x作为( )f x的近似的方法。6 相对误差: 设*x为准确值x的一个近似值,称绝对误差与准确值之比为近似值*x的相对误差,记为*()re x,即*()()re xexx7. 矩阵范数: 对任意 n 阶方阵 A,按一定的规则有一实数与之对应,记为|A。若|A满足(1)| 0A,且| 0A当且仅当
3、0A;(2)对任意实数,都有| |A|A;(3)对任意两个 n 阶方阵 A,B,都有| |ABAB;(4)| |ABA|B称|A为矩阵 A 的范数。8 算子范数 :设 A 为 n 阶方阵,| |是nR中的向量范数,则0|maxxAxAx是一种矩阵范数,称其为由向量范数| |诱导出的矩阵范数,也称算子范数。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 10 页9. 矩阵范数与向量范数的相容性:对任意 n 维向量 x ,都有| |AxA|x这一性质称为矩阵范数与向量范数的相容性。10. 1范数,范数和2范数:(1)1范数11| | |ni
4、ixx(2)范数1| |m ax | ii nxx(3)2范数222212|nxxxx二、简答题1高斯消元法的思想是:先逐次消去变量,将方程组化成同解的上三角形方程组,此过程称为消元过程。然后按方程相反顺序求解上三角形方程组,得到原方程组的解,此过程称为回代过程。2. 迭代法的基本思想是:构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解得规则,由不同的计算规则得到不同的迭代法。3. 雅可比( Jacobi )迭代法的计算过程(算法) :(1)输入()ijAa,1(,)nbbb,维数 n,(0)(0)(0)(0)12(,)nxxxx,最大容许迭代次数 N。(2)置1k(3)对1,2,
5、in( 0 )1() /niii jji ijjixba xa(4)若(0)xx,输出 x 停机;否则转 5。(5) kN ,置(0)1,(1,2, )iikk xxin,转 3,否则,输出失败信息,停机。4. 插值多项式的误差估计: (P102)由(1)(1)101( )( )( )( )()()()(1)!(1)!nnnnnffRxxxxxxxxnn当(0,1, )ixx in时,上式自然成立,因此,上式对 , a b上的任意点都成立,这就叫插值多项式的误差估计。5. 反幂法的基本思想: 设 A 为阶非奇异矩阵,u 为 A 的特征值和相应的特征向量,则1A的特征值是 A 的特征值的倒数,而
6、相应的特征向量不变,即精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 10 页11Auu因此,若对矩阵1A用幂法, ,即可计算出1A 的按模最大的特征值,其倒数恰为 A 的按模最小的特征值。6. 雅可比( Jacobi )迭代法是: 选取初始向量(0)x代入迭代公式(1 )()kkixBxg(0 , 1 , 2 ,k产生向量序列( )kx, 由上述计算过程所给出的迭代法。7. 数值计算中应注意的问题是:(1)避免两个相近的数相减(2)避免大数“吃”小数的现象(3)避免除数的绝对值远小于被除数的绝对值(4)要简化计算,减少运算次数,提高效
7、率(5)选用数值稳定性好的算法8. 高斯消去法的计算量: 由消去法步骤知,在进行第k 次消元时,需作除法nk 次,乘法()nk (1)nk次,故消元过程中乘除运算总量为乘法次数121()(1)(1)3nknnknkn除法次数11()(1)2nknnkn在回代过程中,计算kx需要(1)nk次乘除法,整个回代过程需要乘除运算的总量为1(1)(1)2nknnkn,所以,高斯消去法的乘除总运算量为322(1)(1)(1)32233nnnnnNnnnn9. 迭代法的收敛条件: 对任意初始向量(0)x和右端项g,由迭代格式(1)()kkxMxg(0 , 1 , 2 ,k产生的向量序列()kx收敛的充要条件
8、是()1M。10. 迭代法的误差估计: 设有迭代格式(1)( )kkxMxg,若| 1M,( )kx收敛于*x,则有误差估计式()*(1)(0)|1 |KkMxxxxM。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 10 页二、计算题1. 假定运算中数据都精确到两位小数,试求*1.21 3.659.81x的绝对误差限和相对误差限,计算结果有几位有效数字?解:由式12121212121212()()()()()()rrre xxe xe xxxe xxe xe xxxxx和1221121212()()()()()()rrre x xx
9、e xx e xe x xe xe x得*()3.65(1.21)1.21(3.65)(9.81)e xeee因为式中数据都精确到两位小数,即其误差限均为21102,故有*| ()|3.65 | (1.21)| 1.21 | (3.65) | (9.81) |e xeee*| ()|0.0293|() |0.0054|5.3935re xexx所以,*x的绝对误差限为 0.0293,相对误差限为 0.0054,计算结果有两位有效数字。2. 求矩阵223477245A的三角分解。解:由式111111(1,2, )(2, , )() /(1,2,1,1, )jjiijijikkjkjijijikk
10、jjjkuajnual uin jinlal uujnijn,12122ua,13133ua2121114/22lau,3131112/12lau222221 127223ual u,2323211372 31ual u3232311222()/4( 1)2 / 32lal uu333331 133223()5(1)32 16ual ul u所以100223210031121006A21(3.651.211)100.02932精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 10 页3用幂法(2k)求矩阵210021012A的按模最大的特
11、征值和相应的特征向量。取(0)(0,0,0)Tx. (P77)解:(0)(0)(0,0,1)Tyx(1)(0)(0,1,2)TxAy, 2(1)(1)(0,0.5,1)Txy(2)(1)(0.5,2,2.5)TxAy,2. 54. 已知函数lnyx,x的值是 10,11,12,13,14 对应的lnyx的值分别是2.3026,2.3979, 2.4849, 2.5649, 2.6391 。用 Lagrange线性插值求 ln11.5的近似值。解:取两个节点011x,112x,插值基函数为1001( )(12)xxlxxxx0110( )11xxlxxxx由式011010110( )xxxxxy
12、yxxxx得1( )2.3979(12)2.4849(11)Lxxx将 x=11.5 代入,即得1ln11.5(11.5)2.39790.52.48490.52.4414L按式(1)1( )( )( )(1)!nnnfRxxn(,)a b得1(ln)( )(11)(12)2!xR xxx因为21(ln)xx,在 11 和 12 之间,故2211|(ln)|0.008264511x于是311|(11.5) |0.00826450.50.51.03306 102R精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 10 页5. 用 Jacobi
13、迭代法(1k)求解线性方程组1231231231027210283542xxxxxxxxx. 解:由 Jacobi迭代法得计算公式(1)( )11nkkiiijjjiiiij ibxa xaa得(1)( )( )123(1)()( )213(1)( )( )3120.10.27.20.10.28.30.20.28.4kkkkkkkkkxxxxxxxxx取(0)(0,0,0)Tx,代入上式得(1)17.2x( 1 )28 . 3x( 1 )38 . 4x(2)10.1 8.30.2 8.47.29.71x(2)20.1 7.20.28.48.310.70x(2)30.27.20.28.38.41
14、1.50x6. 设有方程组Axb,其中111221112211122A,讨论用 Jacobi迭代法求解的收敛性。解:因为 A 为对称矩阵,且其各阶主子式皆大于零,故A 为对称正定矩阵, A 不是弱对角占优阵,故不能判别Jacobi迭代的收敛性。易算出Jacobi迭代法的迭代矩阵为1110221102211022BIDA其特征方程311221113|22441122IB精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 10 页21() (1)02有根1212,31,因而( )1B。由向量序列( )kx收敛的充要条件是()1B,故Jacobi
15、迭代法不收敛。7用反幂法(1k)求矩阵210021012A接近 2.93 的特征值,并求相应的特征向量,取(0)(0,0,0)Tx. 解:对2.93AI 作三角分解得0.93102.9300.931010.93AI1000 . 9 31001000 . 9 311101000.930.930.938. 已知函数lnyx,x的值是10,11,12,13,14 对应的lnyx的值分别是2.3026,2.3979, 2.4849, 2.5649, 2.6391 。用 Lagrange抛物线插值求 ln11.5的近似值。解:取011x,112x,213x,插值多项式为2(12)(13)(11)(13)
16、(11)(12)( )2.39792.48492.5649(11 12)(1113)(1211)(1213)(1311)(13 12)xxxxxxLx1.19895(12)(13)2.4849(11)(13)1.28245(11)(12)xxxxxx所以2ln11.5(11.5)L1.19895( 0.5)( 1.5)2.48490.5( 1.5)1.282450.5( 0.5)2.442275因为32(ln)xx,于是2311132max | (ln) |0.15031011xx因此用抛物线插值法计算的误差为2|(ln)|(11.5) |(11.511)(11.512)(11.513) |3
17、!xR精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 10 页2510.1503 100.5 0.51.59.3938106查表可得 ln11.52.442347三、证明题1. 若x的近似值x1210.10 (0)mna aaa有n位有效数字,则111102na为其相对误差限。反之,若x的相对误差限r满足111102(1)nra,则x至少具有n位有效数字。证明:由式*1|102m nxx得*1| () | |102m ne xxx从而有*1*121110()12|()| |100.102m nnrmne xe xxa aaa所以1111
18、02na是*x的相对误差限。若111102(1)nra,由式*()|() | |rre xe xx得*12| () | |() | 0.10mrnre xx e xa aa111111(1) 1010102(1)2mnm naa由式*1|102m nxx,*x至少有 n 位有效数字。2. 设01,nxxx为1n个互异节点,( ),(0,1,)ilxi,n为这组点上的 Lagrange 插值基函数,试证明0( )1niilx。证明:上式的左端为插值基函数的线性组合,其组合系数均为1。显然,函数( )1f x在这n+1 个节点处取值均为1,即()1iiyf x(0 , 1 ,in,由式0( )(
19、)nni iiLxy lx知,它的 n次 Lagrange插值多项式为0( )( )nniiLxlx对任意 x,插值余项为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 10 页(1)1( )( )( )( )( )0(1)!nnnnfRxf xLxxn所以0()()()1nniiLxlxf x3 设A为任意n阶方阵,为任意由向量范数诱导出的矩阵范数,则( )AA证明:对 A 的任一特征值i及相应的特征向量iu,都有|i| | | |iiiiuuAuA|iu因为iu为非零向量,于是有| | | |iA由i的任意性即得()| | |AA4
20、. 设A为n阶方阵,则 lim0kkA的充分必要条件为()1A。证明:必要性。若lim0kkA由相关定义得l i m |kkA而0()() | |KKKAAA于是由极限存在准则,有l i m () kkA所以( )1A。充分性。若()1A,取1( )02A,由|( )AA,存在一种矩阵范数,使得1()|()12AAA而| |kkAA,于是l i m | | |l i m |kkkkAA所以l i m0kkA五、应用题1. 平面桁架是由刚性元件通过结点互相联结而组成的力学结构,它通常出现在桥梁结构和其他需要力学支撑的结构中。如图是一个简单的静力桁架结构,其中刚性元件(5m)通过精选学习资料 -
21、- - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 10 页结 点,A B C D相 连 。 求 各 个 结 点 的 合 力 方 程 , 并 求 出 当,36外 部 负 荷12250,1500gN gN时,求各个节点内力。解:设五个刚性元件的内力为125,fff,它们都处理为压力,如果解是负的,表明该力是张力。桁架的左边由固定结点A 支撑,右边由滑轮 D 支撑,678,fff是外部支撑力,12,gg是外部负荷。由于在静力平衡时,每个结点处的水平方向合力与垂直方向的合力为零,那么有结点 A12617cos0sin0fffff结点 B141134coscos0s
22、insin0ffgfff结点 C253200fffg结点 D4548cos0sin0ffff设f表示未知力向量,上述方程组可用矩阵表示为12cos10001000sin00000100cos00cos0000sin01sin0000001001000000100000000cos10000000sin00010gfg若取,36,外部负荷12250,1500gN gN。采用列主元素法,得各结点的内力如下:( 1174,837,1500,966.5,837, 250,1017, 483.3)Tf精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 10 页