矩形运算方法
15页1、一败涂地、 解线性方程组(线性矩阵方程)解线性方程组是科学计算中最常见的问题。所说的“最常见”有两方面的含义:) 问题的本身是求解线性方程组;) 许多问题的求解需要或归结为线性方程组的求解。关于线性方程组()其求解方法有两类:) 直接法:高斯消去法(Gaussian Elimination);) 间接法:各种迭代法(Iteration)。、高斯消去法) 引例考虑如下(梯形)线性方程组:高斯消去法的求解思路:把一般的线性方程组()化成(上或下)梯形的形式。)高斯消去法示例考虑如下线性方程组:) 第一个方程的两端乘加到第二个方程的两端,第一个方程的两端乘加到第三个方程的两端,得)第二个方程的两端乘加到第三个方程的两端,得) 从上述方程组的第三个方程依此求解,得)高斯消去法的不足及其改进高斯(全、列)主元素消去法在上例中,由于建模、计算等原因,系数2.001而产生0.0005的误差,实际求解的方程组为注:数值稳定的算法高斯列主元素消去法就是在消元的每一步选取(列)主元素一列中绝对值最大的元取做主元素,高斯列主元素消去法是数值稳定的方法。列主元素消去法的基本思想:在每轮消元之前,选列主元素(绝
2、对值最大的元素),使乘数.列主元素消去法的步骤:设已经完成第1步到第步的按列选主元、交换两行、消元计算,得到矩阵.第步计算如下:对于,(1) 选列主元素,即确定使;(2) 如果,则方程组解不唯一,或者接近奇异矩阵,停止运算;(3) 如果,则交换第行与第行元素;(4) 消元计算: (5) 回代计算:完全主元素消去法即是每次选主元时,依次按行、列选取绝对值最大的元素作为主元素,然后交换两行、两列,再进行消元计算.完全主元素消去法的步骤:设已经完成第1步到第步的选主元、交换行和列、消元计算,得到矩阵.第步计算选主元素的范围为,即确定使.第步计算如下:对于,(1) 选主元素,即确定使;(2) 如果,则方程组解不唯一,或者接近奇异矩阵,停止运算;(3) 如果,则交换第行与第行元素;如果,则交换第列与第列元素;(4) 消元计算: (5) 回代求解.【注】 完全主元消去法是解低阶稠密矩阵方程组的有效方法,但完全主元消去法解方程组,在选主元素时要化费较多的计算机时间,行主元消去法与列主元消去法运算量大体相同,实际计算时,用列主元消去法即可满足一定的精度要求.对同一数值问题,用不同的计算方法,所得结果的
3、精度大不一样.对于一个算法来说,如果计算过程中舍入误差能得到控制,对计算结果影响较小,则称此算法是数值稳定的;否则,如果计算过程中舍入误差增长迅速,计算结果受舍入误差影响较大,则称此算法为数值不稳定的.因此,我们解数值问题时,应选择和使用数值稳定的算法,否则如果使用数值不稳定的算法,就可能导致计算失败.)高斯列主元素消去法的MATLAB实现:,意为例 LinearEquiation02.mopen LinearEquiation02 LinearEquiation02 一个典型的例子: Hilbert矩阵:注: 非奇异矩阵的条件数: )分解(Factorization)(高斯消去法、Doolittle分解)高斯消去法的消元过程,从代数运算的角度看就是用一个下三角矩阵左乘方程组的系数矩阵,且乘积的结果为上三角矩阵,即 ()可通过直接用A元素计算矩阵A的三角分解矩阵L和U.这种直接计算A的三角分解的方法有实用上的好处.下面利用矩阵乘法规则来确定三角矩阵L和U.第一步:利用A的第一行、第一列元素确定U的第一行、L的第一列元素.由矩阵乘法,得到,. (3.7)设已经计算出U的第1至r -1行元
4、素,L的第1至r -1列元素,现在要计算U的第r行元素及L的第r列元素.第r步:利用A的第r行、第r列剩下的元素确定U的第r行、L的第r列元素.由矩阵乘法,有,得U的第r行元素为. (3.8)由,得. (3.9)例5 用LU分解法求解方程组.解 对系数矩阵A进行LU分解,.由,有.,.因此.解方程组,得.解方程组,得.6) LU 分解的MATLAB实现:或例 A=rand(5);L,U,P=lu(A)A=rand(5);L,U,P=lu(A)L=PL 当是主对角占优的三对角矩阵时,基于Doolittle分解可得到解这类方程组的追赶法。、Cholesky分解 (Cholesky Factorization)对称正定矩阵的Cholesky分解和以为系数矩阵地的线性方程组的改进的平方根法:设阶方程组,是对称正定矩阵(Positive Definite Matrix),则有三角分解.再将分解为,则.(1) 对称正定矩阵有唯一的分解这是由于,且对称阵,则有再利用三角分解的唯一性,得.因此,对称正定矩阵有唯一的分解.(2) 是正定对角阵(即)由于对称正定的充要条件是对称正定,其中是阶可逆方阵.取,
《矩形运算方法》由会员n****分享,可在线阅读,更多相关《矩形运算方法》请在金锄头文库上搜索。
项目二财务管理价值观念
山东省安全生产风险分级管控与隐患排查治理信息化系统交流材料-2018.9.26
人教版高中地理必修3第一章地理环境与区域发展第二节《地理信息技术在区域地理环境研究中的应用》
第三章2房地产抵押贷款-固定利率抵押贷款
第八章工程质量法律制度
第25讲家庭电路与安全用电
餐厅点餐系统项目
项目7水箱水位控制
框架完整个人年度工作总结范文模板
科目名称-国土交通省
金融工程09课件
高校自主招生之结构化面试
房地产私募股权投资基金(PE)专题研究.
房地产基础知识培训2012
第一章食品检测技术基础知识
第10章网站设计与建设综合实例
第5章尝试迷人的机器人项目机器人灭火项目
自考英语二unit3
企业人力资源管理师第六章劳动法与劳动关系管理
第三章市场营销宏观环境分析
2022-08-31 17页
2023-04-16 12页
2022-12-18 4页
2022-08-24 4页
2023-08-09 12页
2023-04-07 15页
2023-12-24 14页
2024-01-28 14页
2023-03-31 14页
2023-03-28 17页