线性方程组的直接法.doc
《线性方程组的直接法.doc》由会员分享,可在线阅读,更多相关《线性方程组的直接法.doc(31页珍藏版)》请在咨信网上搜索。
1、第二章 线性方程组的直接法在近代数学数值计算和工程应用中,求解线性方程组是重要的课题。例如,样条插值中形成的关系式,曲线拟合形成的法方程等,都落实到解一个元线性方程组,尤其是大型方程组的求解,即求线性方程组(2.1)的未知量的数值. (2。1)其中ai j,bi为常数.上式可写成矩阵形式Ax = b,即 (2。2)其中,为系数矩阵,为解向量,为常数向量。当detA=D0时,由线性代数中的克莱姆法则,方程组的解存在且惟一,且有 为系数矩阵的第列元素以代替的矩阵的行列式的值。克莱姆法则在建立线性方程组解的理论基础中功不可没,但是在实际计算中,我们难以承受它的计算量。例如,解一个100阶的线性方程组
2、,乘除法次数约为(101100!99),即使以每秒的运算速度,也需要近年的时间。在石油勘探、天气预报等问题中常常出现成百上千阶的方程组,也就产生了各种形式方程组数值解法的需求。研究大型方程组的解是目前计算数学中的一个重要方向和课题.解方程组的方法可归纳为直接解法和迭代解法。从理论上来说,直接法经过有限次四则运算,假定每一步运算过程中没有舍入误差,那么,最后得到方程组的解就是精确解。但是,这只是理想化的假定,在计算过程中,完全杜绝舍入误差是不可能的,只能控制和约束由有限位算术运算带来的舍入误差的增长和危害,这样直接法得到的解也不一定是绝对精确的.迭代法是将方程组的解看作某种极限过程的向量极限的值
3、,像第2章中非线性方程求解一样,计算极限过程是用迭代过程完成的,只不过将迭代式中单变量换成向量而已。在用迭代算法时,我们不可能将极限过程算到底,只能将迭代进行有限多次,得到满足一定精度要求的方程组的近似解。在数值计算历史上,直接解法和迭代解法交替生辉。一种解法的兴旺与计算机的硬件环境和问题规模是密切相关的.一般说来,对同等规模的线性方程组,直接法对计算机的要求高于迭代法。对于中等规模的线性方程组,由于直接法的准确性和可靠性高,一般都用直接法求解。对于高阶方程组和稀疏方程组(非零元素较少),一般用迭代法求解.1 消元法一、 三角形方程组的解形如下面三种形式的线性方程组较容易求解。对角形方程组 (
4、2。3)设,对每一个方程,。显然,求解n阶对角方程的运算量为 。下三角方程组 (2。4)按照方程组的顺序,从第一个方程至第个方程,逐个解出。由方程,得。将的值代入到第二个方程得将的值代入到第个方程得计算需要次乘法或除法运算,。因此,求解过程中的运算量为上三角方程组 (2。5)与计算下三角方程组的次序相反,从第个方程至第一个方程,逐个解出。由第个方程。将的值代入到第个方程得将的值代入到第个方程得解的通式 计算需要次乘法或除法运算。因此求解过程中的运算量为消元法的基本思想就是通过对方程组做初等变换,把一般形式的方程组化为等价的具有上述形式的易解方程组.二、 高斯消元法与列主元消元法高斯消元法高斯消
5、元法是我们熟悉的古老、简单而有效的解方程组的方法。下面是中学阶段解二元方程组(高斯消元法)的步骤:(2.6)(2.7)方程(2。6)乘以3加到第(2。7)个方程中得代入(2。6)得。其方法相当于对方程组的增广矩阵做行的初等变换:已是上三角矩阵,而为原方程组的等价方程组,已化成易解的方程组形式.再用回代方法求解,得到:这就是高斯消元法解方程组的消元和回代过程.一般地,可对线性方程组(2。1)施行以下一系列变换;(1)对换某两个方程的次序;(2)对其中某个方程的两边同乘一个不为零的数;(3)把某一个方程两边同乘一个常数后加到另一个方程的两边。记变换后的方程组为: (2。8)显然方程组(2。1)与(
6、2。8)是等价方程组,或者说它们有相同的解.分别记方程组(2.1)与(2。8)的增广矩阵为:可以看出,实际上是由按一系列初等换后得到的(1)对换某两行元素;(2)中的某行乘一个不为零的数;(3)把的某一行乘一个常数后加到另一行。高斯消元法就是通过以上(3)的变换,把化为等价的上三角形式。下面我们以为例演示消元过程。设方程组: (2.9)其增广矩阵为:(1)若,则将第一行乘以加到第二行上;将第一乘以加到第三行上;将第一行乘以 加到第四行上得到 (2.10)即其中:(2)若则将第二行乘以加到第三行上;将第二行乘以加到第四行上,得到(2.11)其中:(3)若则将第三行乘以加到第四行上,得到 (2。1
7、2)其中:已是上三角矩阵,于是得到了与原方程等价的易解形式的方程组:(2.13)再对方程组(2。13)依次回代解出。从式(2.12)可以得到系数矩阵的行列式的值为的对角元素的乘积。即这也正是在计算机上计算方阵的行列式的常规方法。要将上述解方程步骤推广到阶方程组,只需将对控制量“4”的操作改成对控制量的操作即可。元方程组高斯消元法的步骤如下:对于约定有(2.14)经过以上消元,我们已得到与等价的方程组,其中已是一个上三角矩阵.为简单起见,仍记的元素为 (2.15)即已得到原方程组的解。高斯消元法算法在算法中略去了变量,矩阵和向量的定义部分.在消元过程中,将仍放在元素的位置上。1输入:方程组阶数n
8、,方程组系数矩阵A和常数向量b。2FOR k:=1 TO n1 /消元过程 FOR i:=k+1 TO n / 假定 FOR j:=k+1 TO n / ENDFORJ / ENDFORI / ENDFORK3FOR i:=n TO1/回代求解 s:=bi FOR j:=i+1 TO n DO 4输出方程组的解。高斯消元法的运算量整个消元过程即式(2.14)的乘法和除法的运算量为回代过程即式(2.15)的乘除运算量为故高斯消元法的运算量为 (2.16)高斯消元法的可行性在上面的消元法中,未知量是按照在方程组中的自然顺序消去的,也叫顺序消元法。在消元过程中假定对然元素,消元步骤才能顺利进行,由于
9、顺序消元不改变的主子式值,故高斯消元法可行的充分必要条件为的各阶主子式不为零。但是,实际上只要方程组就有解。故高斯消元法本身具有局限性。另一方面,即使高斯消元法可行,如果很小,在运算中用它作为除法的分母,会导致其它元素数量级的严重增长和舍入误差的扩散。这是高斯消元法的另一缺陷。例2。1 方程组(2。17)(2。18)的精确解为:。在高斯消元法计算中取5位有效数字。解:方程(2。17)(1)/0.0003+方程(2。18)得:,代入方程(2.17)得。由此得到的解完全失真,如果交换两个方程的顺序,得到等价方程组经高斯消元后有由此可看到,在有些情况下,调换方程组的次序对方程组的解是有影响的,对消元
10、法中抑制舍入误差的增长是有效的。如果不调换方程组的次序,取6位有效数字计算方程组的解,得到取9位有效数字计算方程组的解,得到由此可见有效数字在数值计算中的作用.列主元消元法如果在一列中选取按模最大的元素,将其调到主干方程位置再做消元,则称为列主元消元法.调换方程组的次序是为了使运算中做分母量的绝对值尽量地大,减少舍入误差的影响.用列主元方法可以克服高斯消元法的额外限制,只要方程组有解,列主元消元法就能畅通无阻地顺利求解,同时又提高了解的精确度。更具体地,第一步在第一列元素中选出绝对值最大的元素,交换第一行和第行的所有元素,再做化简为零的操作。对于每个在做消元前,选出中绝对值最大的元素,对行和行
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性方程组 直接
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。