计算方法-3线性方程组数值解法.ppt

《计算方法-3线性方程组数值解法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2学习要点消去法：包括高斯顺序消去法和列主元高斯消去法；三对角方程组的追赶法；各种向量范数和矩阵范数的概念；解线性方程组的迭代法：Jacobi迭代法；Gauss-Seidel迭代法；迭代法的收敛性判断；33.1问题的提出在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。例如:建筑工程中的结构力学问题;用最小二乘法求实验数据的曲线拟合问题;用差分法或者有限元法解常微分方程，偏微分方程边值问题等都导致求解线性方程组而且后面几种情况常常归结为求解大型线性方程组。4线性代数方面的计算方法就是研究求解线性方程组的一些数值解法与研究

2、计算矩阵的特征值及特征向量的数值方法。5式中，aij，bi为已知常数，xi为待求的未知量。记设有线性方程组6则式（3.1）可写出矩阵形式Ax=b也可以把式（3.1）写成增广矩阵形式7关于线性方程组的数值解法一般有两类：直接法：经过有限步算术运算，可求得方程组的精确解的方法（若在计算过程中没有舍入误差）迭代法：用某种极限过程去逐步逼近线性方程组精确解的方法，迭代法具有占存储单元少，程序设计简单，原始系数矩阵在迭代过程中不变等优点，但存在收敛性及收敛速度等问题83.2消去法一上三角方程组的解法设9写出矩阵形式为：Ux=y其中U称为上三角矩

3、阵10若uii≠0（i=1,2,……n）,则由下至上依次回代得11二高斯消去法其中aij(1)=aij,ai,n+1(1)=bi1≤i,j≤n1、第1步消元。设a11(1)≠0,将的第一行乘以-li1(li1=ai1(1)/a11(1))，加到i行记原线性方程组为12其中aij(2)=aij(1)-li1a1j(1)，2≤i≤n；2≤j≤n+1得到同解方程组13从运算过程来看，上面相当于用矩阵左乘即14若已经进行了k-1步次消元，同解方程组为：15其中，aij(k+1)=aij(k)-likakj(k)，k+1≤i≤n；k+1≤j≥n

4、+1第k步消元：16即第k步消元，也相当于用矩阵左乘17一直作到n-1步，得到同解方程组记则线性方程组为：Ux=y18例1用高斯消去法解线性方程组19解：20因此得到与原方程同解的三角方程组为通过回代，可以得到解为：X3=-3,x2=2,x1=121三列主元高斯消去法≠0（k=1,2,……n-1）;若很小，也不适用；高斯消去法的问题：列主元消去法：与高斯消去法相似；防止舍入误差的增长；选列主元，进行行交换。22具体过程：设k-1步消元得到：23在进行第k步消元前，选出第k列中位于对角线及其以下元素绝对值中的最大值者，即确定t使得将第t

5、行和第k行互相交换；再按高斯消去法进行第k步消元；若系数矩阵为对称正定或严格对角占优的方程组，不必选主元24例2用列主元高斯消去法解线性方程组253.3追赶法其系数矩阵为三对角形，元素满足以下条件：

6、b1

7、>

8、c1

9、>0

10、bi

11、≥

12、ai

13、+

14、ci

15、，且aici≠0i=2，3，……n-1;

16、bn

17、≥

18、an

19、>0。可以采用追赶法求解26消元过程:追赶273.4向量范数和矩阵范数为了研究线性方程组近似解的误差估计和迭代法的收敛性，我们需要对Rn(n维向量空间)中的向量或Rnxn中矩阵的“大小”引入一种度量——向量和矩阵的范数。28向量范数对

20、空间直角坐标系任意向量其长度为：29它满足如下3个条件对任意x属于R3，

21、x

22、≥0，

23、x

24、=0当且仅当x=0;对任意常数c属于R和任意x属于R3，有

25、cx

26、=

27、c

28、·

29、x

30、;对任意x属于R3，y属于R3，有

31、x+y

32、≤

33、x

34、+

35、y

36、;这三个条件分别称为向量长度的非负性、齐次性和三角不等式，推广此概念即可得到向量范数30定义设f(x)=

37、

38、x

39、

40、是定义在Rn上的实函数，如果它满足如下条件：1对对任意x属于Rn，

41、

42、x

43、

44、≥0，

45、

46、x

47、

48、=0当且仅当x=0;2对任意实常数c和任意x属于Rn，有

49、

50、cx

51、

52、=

53、c

54、·

55、

56、x

57、

58、;3对任意x属

59、于Rn，y属于Rn，有

60、

61、x+y

62、

63、≤

64、

65、x

66、

67、+

68、

69、y

70、

71、;则称

72、

73、·

74、

75、为Rn上的向量范数31最常用的是三种范数：向量的1-范数：向量的2-范数：向量的∞-范数：示例(-5，-3，4)T的三种范数分别为：32定义设向量x,y属于Rn上，则称

76、

77、x-y

78、

79、为x和y之间的距离三种范数之间的关系33矩阵范数定义:设A属于Rn×n，

80、

81、·

82、

83、是Rn上的任一向量范数，称为A的矩阵范数，记作Rn×n表示所有n×n阶矩阵组成的实线性空间。34最常用的三种矩阵范数：35矩阵范数的五个性质：36定理:设A∈R，

84、

85、·

86、

87、为一种矩阵范数，则ρ(A)≤

88、

89、

90、A

91、

92、证明：设λ为A的按模最大的特征值，x为相对应的特征向量，则37例5设A=383.5迭代法迭代法及其收敛性思路：与方程求根的迭代法思想类似，对于线性方程组：39通过变形，找到同解方程组：x=Bx+f建立迭代公式x