MATLAB中初等變換與初等矩陣

已知方陣A與置換向量P

方法一：

對A進行行置換ans=A(p,:)

對A進行列置換ans=A(:,p)

對A進行元素置換ans=A(p,p)

方法二：

列置換：

A = gallery('wathen',1,1);
P = [1 3 5 7 2 4 6 8];
% % % 列置換
Index1= eye(size(A,1));
Q1 = Index1(:,P);
eps=A*Q1-A(:,P);
max(max(eps))

 

行置換：

A = gallery('wathen',1,1);
P = [1 3 5 7 2 4 6 8];
% % % 行置換
Index= eye(size(A,1));
Q2 = Index(P,:);
eps2=Q2*A-A(P,:);
max(max(eps2))

元素置換

A = gallery('wathen',1,1);
P = [1 3 5 7 2 4 6 8];
% % % 元素置換
Index= eye(size(A,1));
Q1 = Index(:,P);
Q2 = Index(P,:);
eps3=Q2*A*Q1-A(P,P);
max(max(eps3))