gmd分解，即几何均值分解

输入任意矩阵H，得到其gmd分解矩阵，具体算法参照论文An iterative geometric mean decomposition algorithm for MIMO communications systems。

function [QRS]=mygmd（H）
% gmd分解 H=Q*R*S‘
% K=rank（H） H \in NxM Q \in NxK S \in MxK R \in KxK
% Chen Chiao-En Yu-Cheng Tsai and Chia-Hsiang Yang. “An iterative geometric mean decomposition algorithm for MIMO communications systems.“ IEEE Transactions on Wireless Communications 14.1 （2014）: 343-352.
[usv]=svd（H）;
tm=diag（s）;
sigma=geomean（tm）;
K=rank（H）;
Q=u（:1:K）;R=diag（tm）;S=v（:1:K）;

for k=1:K-1
    % step 2
    diag_ele=diag（R）;
    if R（kk）>=sigma
        tm1=find（diag_ele        p=tm1（1）;
    end
    if R（kk）<=sigma
        tm2=fin