实验八、解线性方程组的迭代法.doc

《实验八、解线性方程组的迭代法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验八、解线性方程组的迭代法解线性方程组的迭代法是用某种极限过程去逐步逼近线性方程组精确解的方法，即是从一个初始向量出发，按照一定的迭代格式产生一个向量序列，使其收敛到方程组的解。迭代法的优点是所需计算机存储单元少，程序设计简单，原始系数矩阵在计算过程中始终不变等。但迭代法存在收敛性及收敛速度问题。迭代法是解大型稀疏矩阵方程组的重要方法。一、实验目的1、熟悉迭代法的有关理论和方法；2、会编制雅可比迭代法、高斯-塞德尔迭代法的程序；3、注意所用方法的收敛性及其收敛速度问题。二、实验任务1、用雅可比迭代法解方程组.注意：若用高斯-塞德尔迭代法则发散。解：输

2、入主程序：functionX=jacdd(A,b,X0,P,wucha,max1)[nm]=size(A);forj=1:m%a(j)=sum(abs(A(:,j)))-2*(abs(A(j,j)));%end%fori=1:n%ifa(i)>=0%disp('请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛')%return%end%end%ifa(i)<0disp('请注意：系数矩阵A是严格对角占优的，此方程组有唯一解，且雅可比迭代收敛')%endfork=1:max1kforj=1:mX(j)=(b(j)-A(j,[1:j-1,j+1

3、:m])*X0([1:j-1,j+1:m]))/A(j,j);endX,djwcX=norm(X'-X0,P);xdwcX=djwcX/(norm(X',P)+eps);X0=X';X1=Ab;if(djwcXwucha)&(xdwcX>wucha)disp('请注意：雅可比迭代次数已经超过最大迭代次数max1')enda,X=X;jX=X1',在MATLAB中输入：A=[12-2；111；2

4、21]；b=[725]’X=jacdd(A,b,[0;0;0],inf,0.01,20)输入结果得：k=1X=725k=2X=13-10-13k=3X=12-1k=4X=12-1请注意：雅可比迭代收敛,此方程组的精确解jX和近似解X如下：X=12-12、用高斯-塞德尔迭代法解方程组.注意：若用雅可比迭代法则发散。解：输入主程序：functionX=gsdddy(A,b,X0,P,wucha,max1)D=diag(diag(A));U=-triu(A,1);L=-tril(A,-1);dD=det(D);ifdD==0disp('请注意：因为对角矩阵D

5、奇异，所以此方程组无解.')elsedisp('请注意：因为对角矩阵D非奇异，所以此方程组有解.')iD=inv(D-L);B2=iD*U;f2=iD*b;jX=Ab;X=X0;[nm]=size(A);fork=1:max1X1=B2*X+f2;djwcX=norm(X1-X,P);xdwcX=djwcX/(norm(X,P)+eps);if(djwcX

6、(xdwcX

7、(xdwcX

8、代收敛,此A的分解矩阵D,U,L和方程组的精确解jX和近似解X如下：')elsedisp('请注意：高斯-塞德尔迭代的结果没有达到给定的精度，并且迭代次数已经超过最大迭代次数max1,方程组的精确解jX和迭代向量X如下：')X=X';jX=jX'endendX=X';D,U,L,jX=jX'然后再MATLAB中输入：A=[10.90.9;0.910.9;0.90.91];b=[1.9;2.0;1.7];X=gsdddy(A,b,[0;0;0],inf,0.01,20)请注意：因为对角矩阵D非奇异，所以此方程组有解.k=1ans=1.90000.2900

9、-0.2710k=2ans=1.88290.5493-0.4890k=3ans=1.84570.7789-0.6622k=4ans=1.79490.9805-0.7979k=5ans=1.73561.1560-0.9025k=6ans=1.67181.3076-0.9815k=7ans=1.60651.4375-1.0396k=8ans=1.54191.5479-1.0808k=9ans=1.47961.6411-1.1086k=10ans=1.42081.7191-1.1259k=11ans=1.36611.7838-1.1349k=12ans=1.

10、31601.8370-1.1377k=13ans=1.27061.8804-1.1359k=1