3 线性方程组解法2

  • 格式:doc
  • 大小:623.50 KB
  • 文档页数:34

35 第3章 线性方程组的解法 本章讨论线性方程组 11112211211222221122nnnn

nnnnnn

axaxaxbaxaxaxbaxaxaxb





LLML

的求解问题. 线性方程组的矩阵表示 Axb 式中A称为系数矩阵,b称为右端项。 36

数值分析中,线性方程组的数值解法主要分为直接法和迭代法两大类。 直接法是用有限次计算就能求出线性方程组“准确解”的方法(不考虑舍入误差);迭代法是由线性方程组构造出迭代计算公式,然后以一个猜测的向量作为迭代计算的初始向量逐步迭代计算,来获得满足精度要求的近似解。 迭代法是一种逐次逼近的方法。 37

1 线性方程组的迭代解法 线性方程组迭代解法有Jocobi迭代法、Gauss-Seidel迭代法及Sor法等 基本思想(与简单迭代法类比) 将线性方程组Axb等价变形为 xBxg 以构造向量迭代格式 1kkxBxg



用算出的向量迭代序列12,,xxL去逼近解。 38

1. 构造原理 (1) Jacobi迭代法 将线性方程组的第i个变元ix用其他n-1个变元表出,可得

11122133111222112332221122111()1()1()nnnn

nnnnnnnnnxbaxaxaxaxbaxaxaxaxbaxaxaxa







LL

L…… 39

Jacobi迭代格式: (1)()()()11122133111(1)()()()22211233222(1)()()()1122111()1()1()kkkknnkkkknnkkkknnnnnnnnnxbaxaxaxaxbaxaxaxaxbaxaxaxa



…-

…-………- (3.6)

(3)取定初始向量000012,,,TnxxxxL,代入,可逐次算出向量序列12,,,kxxxLL,这里12,,,TkkkknxxxxL。 40

(2)Gauss-Seidel迭代法 Seidel迭代格式: (1)()()()11122133111(1)(1)()()22211233222(1)(1)(1)(1)1122111()1()1()kkkknnkkkknnkkkknnnnnnnnnxbaxaxaxaxbaxaxaxaxbaxaxaxa



…-

…-………-

例1对线性方程组 123123123

+22=1+=22+2=3xxxxxxxxx



写出Jacobi迭代格式和Gauss-Seidel迭代格式. 41

3)SOR法 SOR法的迭代格式

1111,1,2,,inkkkkiiiijjijjjjiiixxbaxaxina









L

式中参数称为松弛因子,当 =1时,SOR法就是Seidel迭代法. 42

2.迭代分析及向量收敛 1) 三种迭代法的向量迭格式 对 Ax=b,将系数矩阵A作如下分解 ADLU

1122

12121212,00000000,0000nnnn

nn

aaDaaaaaLUaa



















OLL

LL

MMOMMMOM

LL

则Ax=b可以写成 DLUxb

Jacobi迭代的向量迭代格式 1kkJJxBxg



1()JBDLU

,1JgDb. JB为Jacobi迭代法的

迭代矩阵. Seidel向量迭代格式 1kkSSxBxg



1SBDLU,1SgDLb

.sB为Seidel迭代法的

迭代矩阵. 43

SOR法的向量迭代格式 1kkxBxg



11BDLDU



,1gDLb.B为

超松弛迭代法的迭代矩阵。

三种迭代格式可写成迭代格式 1kkxBxg

 44

2)向量收敛定义 定义1 设向量序列12,,,TkkkknxxxxL及向量****12,,,TnxxxxL都是nR中的向量,如果有 *lim,1,2,,kiikxxinL

成立,则称()kx收敛于*x.简记为 *limkkxx

。 45

3)范数定义与科学计算中的常用范数 定义2 设L是数域K上的一个线性空间,如果定义在L上的实值函数()Px满足

1) xL,有()0Px, 且()00Pxx; 2) ,xLK,有()PxPx; 3) ,xyL,有()PxyPxPy, 则称()P是L上的一个范数,称()Px为x的一个范数。 46

范数的定义很象绝对值函数,故常用P

或表示范数,而范数()Px常记为Px或x。

这样,上面范数定义中的3个条件常写为

1)xL,有0x, 且00xx; 2),xLK,有xx; 3),xyL,有xyxy

将其与绝对值比较,是否很象? 实际上,很多有关绝对值的运算和结论可以平行引进到有关范数的运算和证明问题中。 47

数值分析中常用的线性空间有  n维向量空间 12|,,,,nnkRaaaaaaRL

 矩阵空间 |,mnmnmnijijmnRAAaaR

连续函数空间 



,|[,]Cabfxfxab在上连续

函数空间,Cab是由闭区间,ab上所有连续函数组成的集合,其线性运算定义为 加法:fgfgxfxgx 数乘 :ffxfx,为数 48

在这些空间上,数值分析中常用的范数有 (1)nR的向量范数

1) 11nkkxx 2) 221nkkxx 3) 1maxkknxx 式中向量12,,,TnxxxxL.

例2 计算向量1,2,3Tx的各种范数. 49

(2) nnR的矩阵范数 矩阵范数要满足如下四条 1)nnAR,有0A,且00AA; 2),nnARK,有AA; 3),nnABR,有ABAB; 4),nnABR,有ABAB. 50

由于线性方程组求解问题中,系数矩阵总是与向量联系在一起的,为描述这种联系,引入如下的算子范数概念. 定义3 设矩阵nnAR,称

0maxnPPxRpxAxAx 为矩阵A的算子范数。 容易证明,矩阵A的算子范数也是矩阵范数,且满足不等式关系

PPPAxAx. 51

例3设为矩阵的算子范数,证明若1B,则IB为非奇异矩阵,且 111IBB

证:用反证法。 若IB为奇异矩阵,则其对应的方程组 0IBx

有非零解x,即有0x,使0IBx,得出 Bxx 两边取范数并作范数运算

BxBxxx