09级矩阵与数值分析上机作业.docVIP

《矩阵与数值分析》课程数值实验报告 矩阵05班马强 土木水利学院道路与铁道工程专业 任课教师：董波 一、解线性方程组 1．问题 (1)分别Jacobi迭代法和Gauss-Seidel迭; (2). 用Gauss列主元消去法、 2．程序及计算结果 (1) 源代码： clear all A=[3.4336 -0.52380 0.67105 -0.15270; -0.52380 3.28326 -0.73051 -0.26890; 0.67105 -0.73051 4.02612 -0.09835; -0.15272 -0.26890 0.01835 2.75702]; %A-方程组系数矩阵 b=[-1.0;1.5;2.5;-2.0]; %b-方程组的矩阵表达形式Ax=b的右端列向量表 c=0.000001; %终止条件 N=length(b); x0j=zeros(N,1);%Jacobi迭代法初始值 x0g=zeros(N,1);%Gauss-Seidel迭代法初始值 kj=0;kg=0;%迭代次数计数 D=zeros(N,N);L=zeros(N,N);U=zeros(N,N); %将A分解为D-L-U的形式 for i=1:N D(i,i)=A(i,i); end for i=2:N for j=1:i-1;