注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

being23

写给未来的自己

 
 
 

日志

 
 
关于我

真正的坚定,就是找到力量去做自己喜欢的事情,并为之努力,这样才会觉得生活是幸福的。

网易考拉推荐

MI Algorithm  

2010-07-30 15:40:46|  分类: 研一下 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%iterative algorithm for interference alignment
%communication scenario:3 users 2 antennas per user
%author 柒.Smilence
%2010/07/26 at room 904
clear all
close all
%degree of freedom equals 1
DoF = 1;
%transmit antennas
M = 2;
%receive antennas
N = 2;
%assuming the transmit power P per user equals 1
P = 1;
%generate channel
H11=1/sqrt(2)*randn(N,M);
H12=1/sqrt(2)*randn(N,M);
H13=1/sqrt(2)*randn(N,M);

H21=1/sqrt(2)*randn(N,M);
H22=1/sqrt(2)*randn(N,M);
H23=1/sqrt(2)*randn(N,M);

H31=1/sqrt(2)*randn(N,M);
H32=1/sqrt(2)*randn(N,M);
H33=1/sqrt(2)*randn(N,M);

%reverse channel
Hr11 = conj(H11');
Hr12 = conj(H21');
Hr13 = conj(H31');
Hr21 = conj(H12');
Hr22 = conj(H22');
Hr23 = conj(H32');
Hr31 = conj(H13');
Hr32 = conj(H23');
Hr33 = conj(H33');


%V transmit precoding matrix
%U receive suppression matrix
[U S V] = svd(H11);
V1 = V(:,2);
[U S V] = svd(H22);
V2 = V(:,2);
[U S V] = svd(H33);
V3 = V(:,2);

%V1 = [1;0];
%V2 = [1;0];
%V3 = [1;0];

counter = 0;
threshold = 1e-6;
figure
hold on
%iterative process
while(1)
 
    %the suppression matrix
    [U1] = Calc_Min_EigVec(M,P,DoF,H12,H13,V2,V3);
    [U2] = Calc_Min_EigVec(M,P,DoF,H21,H23,V1,V3);
    [U3] = Calc_Min_EigVec(M,P,DoF,H31,H32,V1,V2);

 
    %reverse precoding matrix
    Vr1 = U1;
    Vr2 = U2;
    Vr3 = U3;

    [Ur1] = Calc_Min_EigVec(M,P,DoF,Hr12,Hr13,Vr2,Vr3);
    [Ur2] = Calc_Min_EigVec(M,P,DoF,Hr21,Hr23,Vr1,Vr3);
    [Ur3] = Calc_Min_EigVec(M,P,DoF,Hr31,Hr32,Vr1,Vr2);
   
    V1 = Ur1;
    V2 = Ur2;
    V3 = Ur3;
   
    %equivalent channel
    h11 = U1'*H11*V1;
    h12 = U1'*H12*V2;
    h13 = U1'*H13*V3;
    h21 = U2'*H21*V1;
    h22 = U2'*H22*V2;
    h23 = U2'*H23*V3;
    h31 = U3'*H31*V1;
    h32 = U3'*H32*V2;
    h33 = U3'*H33*V3;
   
    counter = counter+1;
    WLI = h12 + h13 + h21 + h23 + h31 +h32 ; 
    plot(counter,abs(WLI),'r*');
    if abs(WLI) < threshold
        break;
    end
end

ylabel('abs(WLI)')
xlabel('lterative times ')


MI Algorithm - 柒.smilence - 此地花開成海
  评论这张
 
阅读(228)| 评论(1)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017