解线性方程组的直接法.ppt

《解线性方程组的直接法.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第二章解线性方程组的直接法主要内容：高斯消去法；高斯主元消去法；直接三角分解法；解对称正定方程组的平方根法；行列式和矩阵的求逆运算；方程组的状态和条件数.§2.1高斯消去法线性代数方程组的求解问题：各种各样的科学与工程问题往往最终都要归结为一个线性方程组的求解问题，例如结构分析、网络分析、大地测量、数据分析、最优化及非线性方程组和微分方程组数值解等，常遇到线性方程组的求解问题.求解线性方程组的数值方法大体上可分为直接法和迭代法两大类，其中直接法是指在没有舍入误差的情况下经过有限次运算可求得精确解；而迭代法则是从一个初始向量出发按照

2、一定的计算格式逐次逼近精确解.解第1步：将方程（2.1）乘上（-3/2）加到方程（2.2）将方程（2.1）乘上（-1/2）加到方程（2.3）则得到与原方程等价的方程组（2.1）（2.3）（2.2）高斯消去法是一个古老的求解线性方程组的方法，但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法．例1用加减消元法解方程组一、高斯消去法高斯消去法的思路：系数矩阵的消元和回代通过系数矩阵A的LU三角分解（高斯变换）转化为较容易求解的三角形矩阵高斯消去法就是加减消元法,将给定方程组通过加减消元化为较容易求解

3、的上三角形方程组.n阶线性代数方程组的一般形式为到此原方程组化为(1)到此原方程组化为(n-1)原方程组化为以上为消元过程.(上三角方程组)(n)回代求解公式是回代过程.二.高斯-约当消去法高斯--约当消去法就是将给定方程组通过加减消元化为对角形方程组的方法.是无回代的消去法,所进行的乘除运算次数比高斯消去法(有回代的消去法)多.解第1步：将方程（2.1）乘上（-3/2）加到方程（2.2）将方程（2.1）乘上（-1/2）加到方程（2.3）则得到与原方程等价的方程组（2.1）（2.3）（2.2）用高斯--约当消去法解方程组（2.4）

4、（2.5）（2.6）例2第2步：消去（2.4）和（2.6）式中未知数x2，得到等价的方程组不用回代，即可求得原方程组的解为：第3步：消去（2.7）和（2.8）式中未知数x3，得等价的对角形方程组（2.7）（2.8）三.高斯消去法的计算工作量(略)教材34页高斯—约当消去法消元过程中总的乘除法次数约为算法的时间复杂度:高斯消去法消元和回代过程中总的乘除法次数约为总的加减法次数约为高斯消去法解20阶方程组大约需做3060次乘除法.克莱姆（Cramer）法则大约需做5×1019次乘法．若用矩阵来描述消去法的约化过程，即为高斯消去法解方程

5、组的基本思想是用矩阵的初等变换将系数矩阵约化为具有简单形式的矩阵（上三角矩阵，单位矩阵等），从而容易求解.四.消去法进行到底的条件【定理2.1】当系数矩阵的各阶顺序主子式不为零时,高斯消去法能进行到底.(略)教材35页算法的稳定性:当

6、A(k,k)

7、=0时，运算就会中断；当

8、A(k,k)

9、0但很小时，损失精度和计算溢出.§2.2高斯主元消去法主元：消元过程中的若,则可利用高斯消元法.例3设有方程组解精确解为[方法1]得到回代求解．对于用具有舍入的3位浮点数进行运算，这是一个很好的计算结果．[方法2]用具有行交换的高斯消去法(避免

10、小主元)．用绝对值很小的aii作除数，会导致矩阵元素数量级严重增长和舍入误差的扩散,使得最后的计算结果不可靠．选主元法的思路：为减小方程组求解误差可通过置换A来选主元方法使A的LU三角分解没有太大的元素.列选主元方法：（行交换的技巧）设A∈Rnxn，通过行交换在其原主元以下部分所在列中选主元但其计算过程可能不稳定，不适于求解较大规模的线性代数方程组.用aii作除数，会导致结果矩阵元素数量级严重增长和舍入误差的扩散,使得最后的计算结果不可靠．高斯列主元素消去法计算结果可靠原方程组化为到此原方程组化为是回代过程.(上三角方程组)(n)

11、回代求解公式(n-1)原方程组化为以上为消元过程.解先交换方程1与方程2用列主元消去法解方程组例4全选主元法全选主元消去法：优点:计算结果更可靠；缺点:挑主元花机时更多，次序有变动，程序复杂.（行列交换的技巧）设A∈Rnxn，通过行交换与列交换在其原主元以下部分所在行列（子阵）中选主元.作业P70习题1(1),8OK!Let’shaveabreak!重要的matlab函数[L,U,P,Q]=lu(A)(LU分解，可返回四个参量)；cond(A,p)(条件数)，rcond(A)(条件数倒数)；Ab(Ax=b的解，=inv(A)*b

12、)；hilb(n)(n阶Hilbert矩阵)。