2904203015.李亚韦.实验三.docVIP

 MACROBUTTON MTEditEquationSection2 Equation Chapter 1 Section 1 SEQ MTEqn \r \h \* MERGEFORMAT  SEQ MTSec \r 1 \h \* MERGEFORMAT  SEQ MTChap \r 1 \h \* MERGEFORMAT 三角线性方程组的数值计算 2904203015 李亚韦 实验描述 求解三角线性方程组，可以通过已知函数从只有一个未知函数的式子开始，依次求解出每个未知量即采用消元法和LU分解法求解；也可以通过迭代近似的方式求解出，而迭代的方法大致有雅克比迭代、高斯-赛德尔迭代等方法。 实验内容 假设d1=d2=…=dn=12，b1=b2=…=bn=5;c1=c2=..=cn=a1=a2=…=an=-2;的情况下分别用高斯消元法，LU分解法，雅可比迭代法和高斯-赛德尔迭代法求解  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1)  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) A矩阵为： B矩阵为：B=[5;5;5;...;5;5;5]; 高斯消元法: 根据题意，构造上述的A矩阵和B矩阵，然后通过高斯消元的过程求解出最后的变量值。 流程图： 开始 i=1,j=1 消去i行j列的元 i≤9 得消去i行后的矩阵 j≤9 输出消元后的矩阵 结束 Y Y N N 流程图 1 实验程序： n=input(input a n=: ); A=zeros(n,n); A(1,1)=12; A(1,2)=-2; A(n,n-1)=-2; A(n,n)=12; for i=2:n-1 A(i,i)=12; A(i,i-1)=-2; A(i,i+1)=-2; end %构造A矩阵 B=ones(n,1)*5; %构造B矩阵 R=[A B]; %构造增广矩阵R for i=1:n-1 if R(i,i)==0 break end for j=i+1:n m=R(j,i)/R(i,i); R(j,i:n+1)=R(j,i:n+1)-m*R(i,i:n+1); end end %构造三角矩阵 x=zeros(n,1); x(n)=B(n)/A(n,n); for k=n-1:-1:1 x(k)=(B(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); end x(1:n) %求解X向量 用LU分解求解: LU分解为非直接求解，设A为n×n矩阵，由一般线性方程组如果存在置换矩阵P，使得可得以下关系式  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 我们有替换矩阵对方程组求得Y然后通过会带求得X。 通过以上分析，先通过上三角求解法求得L，U矩阵，然后求得Y ,最后即可得到X。 实验流程图： 开始 I=1,j=1 消去i行，j列的元 i≤9 得到消去i行的矩阵 j≤9 得到L,U矩阵 得到Y向量 输出X向量 结束 Y Y N N 流程图 2 实验程序： n=input(input a n=: ); A=zeros(n,n); A(1,1)=12; A(1,2)=-2; A(n,n-1)=-2; A(n,n)=12; for i=2:n-1 A(i,i)=12; A(i,i-1)=-2; A(i,i+1)=-2; end