ea105.m

matlab 在线性代数中的应用 利用matlab实现可以提高对所学知识的理解 更有助于你考研

% 《工程线性代数（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)