ex5_13.m

源代码集合

A=[6,1,4,2,1; 3,0,1,4,2; -3,-2,-5,8,4]; rank(A)

iA = pinv(A)

B=iA*A*iA

norm(iA-B)   % Testing for iA*A*iA=iA

norm(A*iA*A-A) %Testing for A*iA*A=A

norm(iA*A-A'*iA')  % Testing for iA*A symmetric

norm(A*iA-iA'*A')  % Testing for A*iA symmetric

iiA = pinv(iA)

norm(iiA-A)