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

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

MPDR特征空间波束形成器  

2012-03-16 14:48:39|  分类: 学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.73 Example 6.8.1
% EigenSpace MPDR beamformer
% Array Gain vs. ua/BWnn with correct (Dr=3) estimation
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all

N = 10;
BWNN = 4/N;
n = (-(N-1)/2:(N-1)/2)';
signalRange = [0:BWNN/100:BWNN/2];
Dr = 3;         %Number of eigenvectors selected

ui = [-0.30, 0.30];
Vi = exp(j*n*pi*ui);

INR = 10^(30/10);      

C = exp(j*n*pi*[0]);      %MPDR
f = [1];

k1 = 1;
for SNR = 10.^([-10 0 10 20 30]/10)
    k2 = 1;
    for ua = signalRange
        Va = exp(j*n*pi*ua);
        Ss = SNR*Va*Va';
        Sn = eye(N) + INR*Vi*Vi';
        Sx = Ss + Sn;

        [U,G,V] = svd(Sx);
        Ur = V(:,1:Dr);
        Ces = Ur'*C;
        Lr = G(1:Dr,1:Dr);
        Wes = inv(Lr)*Ces*inv(Ces'*inv(Lr)*Ces)*f;
        W = Ur*Wes;

        SINRo = real((W'*Ss*W)/(W'*Sn*W));
        SINRi = SNR/(1 + 2*INR);
        Gain(k1,k2) = SINRo/SINRi;
        k2 = k2 + 1;
    end
    k1 = k1 + 1;
end

Gain = 10*log10(abs(Gain));
x = signalRange/BWNN;
plot(x,Gain(1,:),'-',x,Gain(2,:),'--',x,Gain(3,:),'-.',x,Gain(4,:),':',x,Gain(5,:),'-*')
grid
xlabel('{\itu}_{\ita}/{\itBW}_{\itNN}')
ylabel('Gain (dB)')
title('Eigenspace MPDR, N=10, Dr=3, ui=+/-0.3, INR=30dB(each), No Loading')
legend('{\itSNR} = -10 dB','{\itSNR} = 0 dB','{\itSNR} = 10 dB','{\itSNR} = 20 dB','{\itSNR} = 30 dB',3)
axis([0 0.25 0 50])

MPDR特征空间波束形成器 - 阿英 - Mr.Right

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.74 Example 6.8.2
% EigenSpace MPDR beamformer
% Array Gain vs. ua/BWnn with under (Dr=2) estimation
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

N = 10;
BWNN = 4/N;
n = (-(N-1)/2:(N-1)/2)';
signalRange =[0:BWNN/200:BWNN/4];
Dr = 2;         %Number of eigenvectors selected

ui = [-0.3, 0.3];
Vi = exp(j*n*pi*ui);

INR1 = 10^(30/10);      %
INR2 = 10^(30/10);      %

C = exp(j*n*pi*[0]);      %MPDR
f = [1];

k1 = 1;
for SNR = 10.^([-10 0 10 20 30]/10)
    k2 = 1;
    SINRi = SNR/(1 + INR1+INR2);
    for ua = signalRange
        Va = exp(j*n*pi*ua);
        Ss = SNR*Va*Va';
        Sn = eye(N) + INR1*Vi(:,1)*Vi(:,1)'+INR2*Vi(:,2)*Vi(:,2)';
        Sx = Ss + Sn;

        [U,G,V] = svd(Sx);
        Ur = V(:,1:Dr);
        Ces = Ur'*C;
        Lr = G(1:Dr,1:Dr);
        Wes = inv(Lr)*Ces*inv(Ces'*inv(Lr)*Ces)*f;
        W = Ur*Wes;

        SINRo = real((W'*Ss*W)/(W'*Sn*W));
        Gain(k1,k2) = SINRo/SINRi;
        k2 = k2 + 1;
    end
    k1 = k1 + 1;
end

Gain = 10*log10(abs(Gain));
x = signalRange/BWNN;
plot(x,Gain(1,:),'-',x,Gain(2,:),'--',x,Gain(3,:),'-.',x,Gain(4,:),':',x,Gain(5,:),'-*')
grid
xlabel('{\itu}_{\ita}/{\itBW}_{\itNN}')
ylabel('Array Gain (dB)')
legend('{\itSNR} = -10 dB','{\itSNR} = 0 dB','{\itSNR} = 10 dB','{\itSNR} = 20 dB','{\itSNR} = 30 dB',3)
axis([0 0.25 -20 5])

MPDR特征空间波束形成器 - 阿英 - Mr.Right

 

  评论这张
 
阅读(540)| 评论(0)
推荐 转载

历史上的今天

评论

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

页脚

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