examp3_13.m

来自「大量matlab实例代码」· M 代码 · 共 19 行

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)   % 测试关系式 iA*A*iA=iA

norm(A*iA*A-A) % 测试关系式 A*iA*A=A

norm(iA*A-A'*iA')  % 测试 iA*A 的对称性

norm(A*iA-iA'*A')  % 测试 A*iA 的对称性

iiA = pinv(iA)  % 对广义逆结果再进行广义逆

norm(iiA-A)     % 和原矩阵进行比较