《数值分析》实验报告
一、问题的提出
由于计算机的实质计算是在一个有限的浮点数集合上进行的,在大量运算
中的误差积累会使得计算过程出现不稳定,通过实验深入了解高斯消元法从理
论算法到数值算法,能稳定的关键是选主元。
二、实验名称
运用 MATLAB 编程实现高斯消去法和高斯列主元消去法。
三、实验目的
1、 熟悉了解高斯消去法和高斯列主元消去法的算法。
2、 学习 MATLAB 软件的功能。
四、基本原理
五、实验环境
操作环境:Windows XP
实验平台:Matlab7.1 软件
六、实验设计
1 高斯顺序消去法
(1)算例:课本 p10 例 1
(2)程序清单
高斯消去法的 MATLAB 函数文件 gauss.m 如下:
clc
clear all
A=[1,1,1;0,4,-1;2,-2,1];
r=[6;5;1];
[row,col]=size(A);
n=row;
%-------------------------------------------------------------------
B=A;
b=r;
C=B;
c=b;
for i=1:n
[maxi,row]=max(abs(C(:,i)));
if C(row,i)<0
maxi=-maxi;
end
B(i,:)=C(row,:);
B(row,:)=C(i,:);
b(i)=c(row);
b(row)=c(i);
B(i,:)=B(i,:)/maxi;
b(i)=b(i)/maxi;
C=B;