消元法实验报告14.doc

消元法实验报告14.doc西京学院数学软件实验任务书
课程名称
数学软件实验
班级
***
学号
***
姓名
***
实验课题
线性方程组高斯消去法，高斯列主元消去法，高斯全
主元消去法
实验目的
熟悉线性代数方程组高斯消去法,高斯列主元消去法，
高斯全主元消去法
实验要求
运用 Matlab/C/C++/Java/Maple/Mathematica 等其中
一种语言完成
实验内容
线性方程组高斯消去法
线性方程组高斯列主元消去法
线性方程组高斯全主元消去法
成绩
教师
一•线性方程组高斯消去法
1 •线性方程组高斯消去法程序
function maintest
clc
clear all
num=input('please input the order n=')
A=zeros(num,num)；
for i = 1：num
for j =1:num
A(i,j)=input('');
end
end
A
if det(A)~=0
for i=l:num
b (i)=input('')；
end
b=b'
for k=l:num-1
%max (A (:, k))
for i=k+l:num
if A(k,k)~=0
l=A(i,k)/A(k,k);
A(i,:)=A(i,:)-A(k,:)*1;
b(i) =b(i) -b(k) *1；
end
end
end
A
b
x(num)=b(num)/A(num,num)；
for i=num-1：-1：1
sum=0；
for j = i +1：num
sum=sum+A(i,j)*x(j)；
end
x (i) = (b (i) - sum) /A(i, i)；
end
x
end
end
2-线性方程组高斯消去法程序运行结果示例
please input the order n=3
num =
3
1
2
3
1
3
5
1
3
6
A =
1 2 3
1 3 5
1 3 6
2
3
4
A =
1 2 3
0 1 2
0 0 1
b =
2
1
1
x =
1 -1 1
»
二•线性方程组高斯列主元消去 法
1 -线性方程组高斯列主元消去法程序
function [Ra,Rb,n,X]=GaussXQLineMain(A,b) B=[A b];
n=length(b)；
Ra=rank(A)；
Rb=rank(B)； if Ra==Rb if Ra==n
disp ('因为Ra=Rb=n,所以此方程组有唯一•解.')
X=zeros(nz1)；
C=zeros(1,n+1)； for p= 1：n-1
[Y,j]=max(abs(