ea105.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 18 行

% 《工程线性代数（MATLAB版）》第一章例题1.5程序ea105
% 《工程线性代数（MATLAB版）》引用
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 消元法举例
% 
clear
A=[0,1,-2,6,7;1,1,2,-1,2;2,3,2,4,3;-3,-4,-4,-3,-5] ,pause
A1=A;A1([1,2],:)=A1([2,1],:),pause    % 交换第1，2行
%A1(2,:)=A1(2,:)-A1(1,:)*A1(2,1)/A1(1,1),pause			% 消去A(2,1)
A1(3,:)=A1(3,:)-A1(1,:)*A1(3,1)/A1(1,1),pause			% 消去A(3,1)
A1(4,:)=A1(4,:)-A1(1,:)*A1(4,1)/A1(1,1),pause			% 消去A(4,1)
A1(3,:)=A1(3,:)-A1(2,:)*A1(3,2)/A1(2,2),pause			% 消去A(3,1)
A1(4,:)=A1(4,:)-A1(2,:)*A1(4,2)/A1(2,2),pause			% 消去A(4,1)
A1(4,:)=A1(3,:)-A1(3,:)*A1(3,5)/A1(3,5),pause			% 消去A(4,5)
A1(2,:)=A1(2,:)-A1(3,:)*A1(2,5)/A1(3,5),pause			% 消去A(2,5)
A1(1,:)=A1(1,:)-A1(3,:)*A1(1,5)/A1(3,5),pause			% 消去A(1,5)
A1(1,:)=A1(1,:)-A1(2,:)*A1(1,2)/A1(2,2),pause			% 消去A(1,2)