消元法实验报告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=1:num
b(i)=input('');
end
b=b'
for k=1:num-1
%max(A(:,k))
for i=k+1:num
if A(k,k)~=0
l=A(i,k)/A(k,k);
A(i,:)=A(i,:)-A(k,:)*l;
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
b =
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=R