資源簡介
輸入任意矩陣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
- 上一篇:matlab仿真車輛追逐跟蹤
- 下一篇:基于matlab的數字變聲器和濾波器
評論
共有 條評論