数值分析编程及运行结果（高斯顺序消元法） 

高斯消元法

1.程序:

clear format rat

A=input('输入增广矩阵A=') [m,n]=size(A); for i=1:(m-1) numb=int2str(i);

disp(['第',numb,'次消元后的增广矩阵']) for j=(i+1):m

A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end %回代过程 disp('回代求解') x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1

x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m)')/A(i,i); end x

.

2.运行结果：

.

.

高斯选列主元消元法

1. 程序:

clear format rat

A=input('输入增广矩阵A=') [m,n]=size(A); for i=1:(m-1) numb=int2str(i);

disp(['第',numb,'次选列主元后的增广矩阵']) temp=max(abs(A(i:m,i))); [a,b]=find(abs(A(i:m,i))==temp); tempo=A(a(1)+i-1,:); A(a(1)+i-1,:)=A(i,:); A(i,:)=tempo

disp(['第',numb,'次消元后的增广矩阵']) for j=(i+1):m

A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end %回代过程

disp('回代求解')

.

.

x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1

x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m)')/A(i,i); end x

2.运行结果：

.

.

.