线性方程组的直接法和迭代法文档格式.docx
- 文档编号:13904692
- 上传时间:2022-10-14
- 格式:DOCX
- 页数:13
- 大小:377.90KB
线性方程组的直接法和迭代法文档格式.docx
《线性方程组的直接法和迭代法文档格式.docx》由会员分享,可在线阅读,更多相关《线性方程组的直接法和迭代法文档格式.docx(13页珍藏版)》请在冰豆网上搜索。
所以当时,方程组有唯一解。
定理2当齐次线性方程组,时该方程组有唯一的零解。
定理3齐次线性方程组有非零解。
1.2Gauss消元法
Gauss消元法是线性代数中的一个算法,可用来为线性方程组求解,求出矩阵的秩,以及求出可逆方阵的逆矩阵。
当用于一个矩阵时,高斯消元法会产生出一个“行梯阵式”。
1.2.1用Gauss消元法为线性方程组求解
eg:
Gauss消元法可用来找出下列方程组的解或其解的限制:
这个算法的原理是:
首先,要将以下的等式中的消除,然后再将以下的等式中的消除。
这样可使整个方程组变成一个三角形似的格式。
之后再将已得出的答案一个个地代入已被简化的等式中的未知数中,就可求出其余的答案了。
在刚才的例子中,我们将和相加,就可以将中的消除了。
然后再将和相加,就可以将中的消除。
方程组则变为:
现在将和相加,就可将中的消除,方程组变为:
这样就完成了整个算法的初步,一个三角形的格式(指:
变量的格式而言,上例中的变量各为3,2,1个)出现了。
第二步,就是由尾至头地将已知的答案代入其他等式中的未知数。
第一个答案就是。
然后直接带入,立即就可得出第二个答案:
和最后一个答案。
这样,我们利用高斯消元法解决了这个方程组。
2.线性方程组迭代法
2.1Jacobi迭代法
对于线性方程组则,即将A分解为一个严格下三角矩阵、一个对角阵和一个严格上三角矩阵之和,从而可写出Jacobi迭代格式的矩阵表示形式为:
,其迭代矩阵)称为雅可比迭代矩阵.
将线性方程组变为一个通解方程组,对其进行迭代式改写,矩阵B为迭代矩阵
由方程组(I)的第i个方程解出,得到一个同解方程组:
构造相应的迭代公式
取初始向量,利用(III)反复迭代可以得到一个向量序列,利用此迭代格式求解方程组的解法称为Jacobi迭代法。
用Jacobi迭代求解下列方程组
输入
A=[430;
33-1;
0-14];
b=[24;
30;
-24];
[x,k,index]=Jacobi(A,b,1e-5,100)
输出:
x=
-2.9998
11.9987
-3.0001
k=
100
index=
0
所以解为:
=-2.9998,=11.9987,=-3.0001
2.2Gauss-Seide迭代
若L、U、D为上述的L、U、D。
则Gauss—Seidel迭代法的矩阵表示为:
,现将显示化由得:
,令,,则得:
,此即为Gauss—Seidel迭代法的矩阵表示形式,G称为迭代阵。
由Jacobi迭代法中,每一次的迭代只用到前一次的迭代值,若每一次迭代充分利用当前最新的迭代值,即在计算第个分量时,用最新分量,代替旧分量,,就得到所谓解方程组的Gauss-Seidel迭代法。
其迭代格式为(初始向量),
或者写为
用Gauss-Seide迭代求解下列方程组
x0=[0;
0;
0];
[v,sN,vChain]=gaussSeidel(A,b,x0,0.00001,11)
v=
0.6169
11.1962
-4.2056
sN=
11
vChain=
6.000010.0000-6.0000
-1.50002.0000-3.5000
4.500010.3333-5.5000
-1.75003.6667-3.4167
3.250010.6111-5.0833
-1.95835.0556-3.3472
2.208310.8426-4.7361
-2.13196.2130-3.2894
1.340311.0355-4.4468
-2.27667.1775-3.2411
0.616911.1962-4.2056
000
所以结果为:
=0.6169,=11.1962,=-4.2056。
2.3SOR迭代
在很多情况下,Jacobi和Gauss—Seidel法收敛速度较慢,SOR法是Gauss—Seidel法的一种加速方法,需要施加合适的松弛因子。
若L、U、D为上述的L、U、D。
SOR迭代公式为,其中
,。
用SOR迭代求解下列方程组。
取初始点松弛因子,精度要求。
A=[0.76-0.01-0.14-0.16;
-0.010.88-0.030.06;
-0.14-0.031.01-0.12;
-0.160.06-0.120.72];
B=[0.681.180.120.72];
X0=[0;
0];
W=1.05
[x,n]=SOR(A,b,x0,w)
输出
x=
1.2715
1.2844
0.4858
1.2843
n=
7
有上述结果得出:
经过7次迭代后,该方程组的解为
x1=1.2715,x2=1.2844,x3=0.4858,x4=1.2843
2.4迭代法收敛
引理:
设A为n阶方阵,则的充要条件为(为普半径)。
证明:
必要性若
由矩阵收敛的定义知有因为所以
由夹逼定理可得出,又因为所以由可得出
充分性:
若,取,存在矩阵范数,使得则有:
,
由算子范数相容性可得:
由夹逼定理可得出。
定理1:
迭代公式收敛的充分必要条件是迭代矩阵的谱半径
必要性
设存在n维向量,使得,则满足
由迭代公式得出
所以
因为为任意n维向量,因此上式成立必须
由引理可得出。
若,则不是B的特征值,因而有,于是对任意的n维向量f,方程组有唯一解,记为。
即:
且,又因为
所以,对任意初始向量都有
即迭代公式收敛。
注解:
代矩阵的谱半径这里的是矩阵的特征值。
定理2:
若迭代矩阵的一种范数,则对任意的初始向量和任意。
迭代格式均收敛。
迭代法收敛与否只决定于迭代矩阵的普半径,于初始向量及右端项无关。
对于同一方程组,由于不同的迭代法迭代矩阵不同,可能出现有的方法收敛,有的发散的情形。
三种常用矩阵范数为:
,,(是的最大特征值)。
由上述两个定理可得:
Jacobi迭代收敛的充分必要条件是,收敛的充分条件是任一种范数,Gauss—Seidel迭代收敛的充分必要条件是,收敛的充分条件是任一种范数。
定理3:
对角占优线性方程组的Jacobi迭代格式和Gauss—Seidel迭代格式均收敛。
n阶方阵A,如果其主对角线元素的绝对值大于同行其他元素绝对值之和,则称A是对角占优的.
虽然定理1是充分必要条件,可是需要计算迭代矩阵的特征值,我们在线性代数上看到一个大型矩阵的特征值是非常难求的,计算量是很大的。
定理2和定理3的计算量相对定理1来说是相当小的,因为它们只是简单的加减乘除;
但是他们是充分条件,不满足时需要再改用定理1来判断。
所以收敛的判断可先采用定理2和定理3,不行再选择定理1。
预处理是用来加速迭代法收敛的一个重要手段。
值得注意的是,经典的求解线性方程组的迭代法也都能够看作是求解采用不同的因子预处理之后得到的线性方程组的迭代法。
换句话来说,原线性方程组的松弛迭代法等价于预处理之后的方程组的定点迭代法。
谱条件数是反映预处理因子性态是否良好的一个有效指标。
在估计特征值和条件数的界的研究领域,虽然已经有了很多的研究成果,但是支持理论还是一个全新的概念。
支撑理论是一个用于分析预处理方程组的最大(或最小)特征值和条件数的代数架构,它最初产生于对称正定的线性方程组(7)
2.5迭代法收敛的应用
用Jacobi迭代,Gauss—Seidel迭代两种方法求解下列方程组是否收敛。
因为迭代法收敛与否只决定于迭代矩阵的普半径。
所以求普半径是否小于1.
Jacobi迭代矩阵
其特征方程
所以,所以所以Jacobi迭代法收敛。
Gauss—Seidel迭代矩阵
所以特征值为所以所以Gauss—Seidel迭代法发散。
上述例子说明对于同一方程组,由于不同的迭代法迭代矩阵不同,可能出现有的方法收敛,有的发散的情形。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性方程组 直接 迭代法