ea104a.m

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

% 《工程线性代数（MATLAB版）》第一章例题1.4程序ea104a
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 用消元子程序消元法举例
% 
A=[2,-2,0,6;2,-1,2,4;3,-1,4,4;1,1,1,8] ,pause
b=[-2;-2;-3;2] ,pause			% 给系数矩阵A,b赋值
B=[A,b] ,pause					% 赋值B为增广矩阵
B1=gauss(B,1,2,1),pause			% 消去B(2,1)
B1=gauss(B1,1,3,1),pause		% 消去B1(3,1)
B1=gauss(B1,1,4,1),pause		% 消去B1(4,1)
B1=gauss(B1,2,3,2),pause		% 消去B1(3,2)
B1=gauss(B1,2,4,2),pause		% 消去B1(4,2)
B1([3,4],:)= B1([4,3],:) 		% 由于B1(3,3)=0,需交换B1中的第3,4行
B1=gauss(B1,3,4,3),pause		% 消去B1(4,3)
% 至此完成了行阶梯变换,最后的B1就是行阶梯形式的等价增广矩阵
B2= gauss(B1,4,3,4),pause		% 消去B1(3,4)
B2= gauss(B2,4,2,4),pause		% 消去B2(2,4)
B2= gauss(B2,4,1,4),pause		% 消去B2(1,4)
B2= gauss(B2,3,2,3),pause		% 消去B2(2,3)
B2= gauss(B2,3,1,3),pause		% 消去B2(1,3)
B2= gauss(B2,2,1,2),pause		% 消去B2(1,2)
% 至此完成了对角形的行阶梯变换 
B2(1,:)=B2(1,:)/B2(1,1),pause		% B2第一行对角元素归一
B2(2,:)=B2(2,:)/B2(2,2),pause		% B2第二行对角元素归一
B2(3,:)=B2(3,:)/B2(3,3),pause		% B2第三行对角元素归一
B2(4,:)=B2(4,:)/B2(4,4) 			% B2第四行对角元素归一
% 至此完成了最简形行阶梯变换，B2是所需的结果。