重庆大学数值计算实验.doc

WORD完美整理版 范文范例 参考指导 重 庆 大 学 学 生 实 验 报 告 实验课程名称 数值计算 开课实验室 数统学院楼104实验室 学 院 年级 专业班 班 学 生 姓 名 学 号 开 课 时 间 至 学年第 学期 总 成 绩 教师签名 课程 名称 数值计算实验 实验项目 名 称 线性方程组的数值解法 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 √ 教师签名 年 月 日 实验目的： （1）高斯列主元消去法求解线性方程组的过程 （2）熟悉用迭代法求解线性方程组的过程 （3）设计出相应的算法，编制相应的函数子程序 实验内容 分别用高斯列主元消去法 ,Jacobi迭代法,Gauss--Saidel迭代法,超松弛迭代法求解线性方程组 实验步骤 高斯列主元消去法 代码： format long;%设置为长格式显示,显示15位小数 A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]%方程系数矩阵 b=[b1,b2,b3] %方程右侧系数矩阵 [m,n]=size(A); %先检查系数正确性 if m~=n error(矩阵A的行数和列数必须相同); return; end if m~=size(b) error(b的大小必须和A的行数或A的列数相同); return; end %再检查方程是否存在唯一解 if rank(A)~=rank([A,b]) error(A矩阵的秩和增广矩阵的秩不相同,方程不存在唯一解); return; end c=n+1; A(:,c)=b; %(增广) for k=1:n-1 [r,m]=max(abs(A(k:n,k))); %选主元 m=m+k-1; %修正操作行的值 if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去 end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end disp(X=); disp(x); format short;%设置为默认格式显示,显示5位 Jacobi迭代法 代码： clc; A=[ a11,a12,a13,b1;a21, a22,a23,b2; a31,a32,a33,b3] %增广矩阵 MAXTIME=50;%最多进行50次迭代 eps=1e-5;%迭代误差 [n,m]=size(A); x=zeros(n,1);%迭代初值 y=zeros(n,1); k=0; %进入迭代计算 disp(迭代过程X的值情况如下:) disp(X=); while 1 disp(x); for i=1:1:n s=0.0; for j=1:1:n if j~=i s=s+A(i,j)*x(j); end y(i)=(A(i,n+1)-s)/A(i,i); end end for i=1:1:n maxeps=max(0,abs(x(i)-y(i))); %检查是否满足迭代精度要求 end if maxeps=eps%小于迭代精度退出迭代 for i=1:1:n x(i)=y(i);%将结果赋给x end return; end for i=1:1:n%若不满足迭代精度要求继续进行迭代 x(i)=y(i); y(i)=0.0; end k=k+1; if kMAXTIME%超过最大迭代次数退出 error(超过最大迭代次数,退出);