计算方法-实验三列主元高斯消去法.doc

计算方法课程设计报告 实验三 高斯列主元消去法 姓名：黄仁化 学号：031010151551017 班级：计算机科学与技术2004班 日期 ：二○○六年六月十日 一、实验目的： 1、掌握高斯消去法的基本思路和迭代步骤。 培养编程与上机调试能力。 二、高斯列主元消去法的基本思路与计算步骤： 列主元高斯消去法计算步骤： 将方程组用增广矩阵表示。 步骤1：消元过程，对 选主元，找使得 如果，则矩阵奇异，程序结束；否则执行（3）。 如果，则交换第行与第行对应元素位置，，。 消元，对，计算对，计算 步骤 2：回代过程： 若则矩阵奇异，程序结束；否则执行（2）。 对，计算 三：程序流程图 四：程序清单： function X=uptrbk(A,b) % A是一个n阶矩阵。 % b是一个n维向量。 % X是线性方程组AX=b的解。 [N N]=size(A); X=zeros(1,N+1); Aug=[A b]; for p=1:N-1 [Y,j]=max(abs(Aug(p:N,p))); C=Aug(p,:); Aug(p,:)=Aug(j+p-1,:); Aug(j+p-1,:)=C; if Aug(p,p)==0 A是奇异阵，方程无惟一解 break end for k=p+1:N m=Aug(k,p)/Aug(p,p); Aug(k,p:N+1)=Aug(k,p:N+1)-m*Aug(p,p:N+1); end end % 这里用到程序函数backsub来进行回代。 X=backsub(Aug(1:N,1:N),Aug(1:N,N+1)); function X=backsub(A,b) % A是一个n阶上三角非奇异阵。 % b是一个n维向量。 % X是线性方程组AX=b的解。 n=length(b);%取b向量的个数。 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 五、测试数据与结果： 测试数据：（第8章习题三第2题） 求解线性方程组： A=[1,2,3;5,4,10;3,-0.1,1]; b=[1;0;2]; uptrbk(A,b) 计算结果如下: ans = 1.2000 2.0000 -1.4000 六、小结 本实验通过MATLAB程序编程实现了高斯列主元消去法的求解，能加深对高斯列主元消去法的基本思路与计算步骤的理解。 1 姓名：黄仁化　　　　　　　　　　　　学号：031010151551017　　　　　　　　　　实验三 高斯列主元消去法 开始 设置方程组的参数矩阵A，结果矩阵b 代入uptrbk(A,b)函数 调用backsub(A,b) 函数进行回代 程序自动进行计算 输出方程组的解 结束