数值计算实验.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数值计算 开课实验室 数学与统计学院 学 院 数学与统计学院 年级 X级 专业班 数学与应用数学一班 学 生 姓 名 X 学 号 X 开 课 时 间 2012 至 2013 学年第 一 学期 总 成 绩 教师签名 课程 名称 数值计算 实验项目 名 称 Gauss消元法 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 √ 一、实验目的 （1）高斯列主元消去法求解线性方程组的过程 （2）熟悉用迭代法求解线性方程组的过程 （3）设计出相应的算法，编制相应的函数子程序 二、实验内容 分别用高斯列主元消去法 ,Jacobi迭代法,Gauss--Saidel迭代法,超松弛迭代法求解线性方程组 三、实验原理 （1）高斯列主元消去法：设有方程组AX=b，将矩阵的初等行变换作用于方程组大的增广矩阵B=[A:b]，将A变换成一个上三角矩阵，然后对方程组进行求解。 （2）Jacobi迭代法：设方程组Ax=b 满足aii ≠0, 将方程组变形为: x=Bx+f, 则雅可比(Jacobi)迭代法是指x(k+1)=Bx(k)+f，即由初始解逐步迭代即可得到方程组的解。 （3）Gauss--Saidel迭代法: 在雅可比迭代法计算过程中,在迭代的每一步计算过程中是用的全部分量来计算的所有分量,显然在计算第i个分量时,已经计算出的最新分量没有被利用，对这些最新计算出来的第次近似的分量加以利用,就得到所谓解方程组的Gauss-Seidel迭代法。 （4）超松弛迭代法：在Gauss--Saidel迭代法的基础上，对每次迭代是的改变量乘以一个因子w，是收敛速度加快。 四、实验步骤 （1）高斯列主元消去法：根据实验原理编写matlab程序，首先判断解的个数，再分情况进行列主元消去求解。 Matlab程序如下： clc A=input(A); b=input(b); B=[A b]; n=length(b); RA=rank(A); RB=rank(B); c=RB-RA; if c0, disp(方程组无解) return end if RA==RB if RA==n fprintf(方程组有唯一解,n) X=zeros(n,1); for p=1:n-1 t=find(abs(B(p:end,p))==max(abs(B(p:end,p))))+p-1; if abs(B(t,p))~=abs(B(p,p)) l=B(t,:); B(t,:)=B(p,:); B(p,:)=l; end for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); end end b=B(1:n,n+1); A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); end else disp(方程组有无穷多解) return end end disp(X);X (2) Jacobi迭代法: 根据实验原理编写matlab程序，调用Jacobi_iterative函数，首先给定误差限，确定最大迭代次数，确认矩阵A的阶，然后进行迭代过程，最后输出结果。 Matlab程序如下： A=[2 10 0 -3;-3 -4 -12 13;1 2 3 -4;4 14 9 -13]; b=[10 5 -2 7]; e=0.000001; %控制误差 n=max(size(A)); %测定维数 for i=1:n if A(i,i