Matlab求解非方阵的广义逆矩阵程序_zhangyin0913

个人资料

微博

正文字体大小：大中小

Matlab求解非方阵的广义逆矩阵程序

(2016-11-14 21:20:44)

标签：

clear all

close all

clc

G = [0 1 0 0 1 0 0 0;0 0 1 0 0 0 -1 0;1 0 0 0 0 1 0 0;

-1 0 0 0 0 1 0 0;0 0 0 1 0 0 0 -1;0 1 0 0 -1 0 0 0];

% G = [2 4 1 1;1 2 -1 2;-1 -2 -2 1]; %非行满秩非列满秩示例

k = rank(G) %求解矩阵秩

[m,n] = size(G); %提取矩阵行数和列数

if k == m %判断矩阵是否行满秩

G0 = G'*inv(G*G'); %行满秩时广义逆求解

elseif k == n %判断矩阵是否列满秩

G0 = inv(G'*G)*G'; %列满秩时广义逆求解

else %既非行满秩也非列满秩时广义逆求解

A = orth(G); %行满秩分解

B = A'*G; %行满秩分解

G0 = B'*inv(B*B')*inv(A'*A)*A';

end

GG = vpa(G0) %化简结果

%G1 = pinv(G); %调用matlab函数pinv（）求解

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report