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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

LCMP-DL波束形成器的阵列增益  

2012-03-16 13:32:17|  分类: 学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.56
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Example 6.7.8
% LCMV with Directional Constraints.
% Relations between Array Gain, INR, SNR and Optimal LNR are studied

clear all
close all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
ud = 0.1;                       %constraint region

INR = 10^(30/10);
SNRrange = [0 10 20 30];

%LNRrange = -20:50;
LNRrange = -20:2:50;
%signalRange = -ud:ud/200:ud;
signalRange = -ud:ud/50:ud;

ui1 = -0.30;       %interferer configuration
ui2 = 0.30;
Vi1 = exp(j*n*pi*ui1);
Vi2 = exp(j*n*pi*ui2);

uc = 0.0866;
C = exp(j*n*pi*[0 -uc uc]);   %LCMP with directional constraints
f = [1;sin(-(N/2)*pi*uc)./(N*sin(-.5*pi*uc));sin((N/2)*pi*uc)./(N*sin(.5*pi*uc))];

k1 = 1;
for SNR = 10.^(SNRrange/10)
   disp(['loop ' int2str(k1) ' of 4 ...'])
   k2 = 1;
   for LNR = 10.^(LNRrange/10)
      Gain(k1,k2) = 0;
      for ua = signalRange
         Va = exp(j*n*pi*ua);
         Ss = SNR*Va*Va';
         Sn = eye(N) + INR*Vi1*Vi1' + INR*Vi2*Vi2';
         Sx = Ss + Sn + LNR*eye(N);
        
         W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;            %LCMP
         SINR0 = (W'*Ss*W)/(W'*Sn*W);
         SINRi = SNR/(1+INR+INR);
         Gain(k1,k2) = Gain(k1,k2) + SINR0/SINRi;    % 我以前不知道阵列增益还可以这样算
      end     %end of ua
      k2 = k2 + 1;
   end      %end of LNR
   k1 = k1 + 1;
end       %end of SNR

A = 10*log10( abs(Gain)/length(signalRange) );

plot(LNRrange,A(1,:),'-',LNRrange, A(2,:),'--',LNRrange, A(3,:),'-.',LNRrange, A(4,:),':');
xlabel('{\itLNR} (dB)','Fontsize',14)
ylabel('{\itE}(gain) (dB)','Fontsize',14)
title('LCMP\_DL(DIR), u_i=+/-0.30, u_a~U(-0.1, 0.1), u_c=[0, -0.0866, 0.0866], g=[1, Bc, Bc], INR=20dB')
grid
axis([-20 50 0 50])
h=legend('{\itSNR}=0 dB','{\itSNR}=10 dB','{\itSNR}=20 dB','{\itSNR}=30 dB');
set(h,'Fontsize',12)

LCMP-DL波束形成器的阵列增益 - 阿英 - Mr.Right

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.57
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Example 6.7.8
% LCMV with Directional Constraints.
% Relations between Array Gain, INR, SNR and Optimal LNR are studied

clear all
close all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
ud = 0.1;                       %constraint region

INR = 10^(20/10);
SNRrange = [0 10 20 30];

%LNRrange = -20:50;
LNRrange = -20:2:50;
%signalRange = -ud:ud/200:ud;
signalRange = -ud:ud/50:ud;

ui1 = -0.30;       %interferer configuration
ui2 = 0.30;
Vi1 = exp(j*n*pi*ui1);
Vi2 = exp(j*n*pi*ui2);

uc = 0.0866;
C = exp(j*n*pi*[0 -uc uc]);   %LCMP with directional constraints
f = [1;sin(-(N/2)*pi*uc)./(N*sin(-.5*pi*uc));sin((N/2)*pi*uc)./(N*sin(.5*pi*uc))];

k1 = 1;
for SNR = 10.^(SNRrange/10)
   disp(['loop ' int2str(k1) ' of 4 ...'])
   k2 = 1;
   for LNR = 10.^(LNRrange/10)
      Gain(k1,k2) = 0;
      for ua = signalRange
         Va = exp(j*n*pi*ua);
         Ss = SNR*Va*Va';
         Sn = eye(N) + INR*Vi1*Vi1' + INR*Vi2*Vi2';
         Sx = Ss + Sn + LNR*eye(N);
        
         W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;            %LCMP
         SINR0 = (W'*Ss*W)/(W'*Sn*W);
         SINRi = SNR/(1+INR+INR);
         Gain(k1,k2) = Gain(k1,k2) + SINR0/SINRi;
      end     %end of ua
      k2 = k2 + 1;
   end      %end of LNR
   k1 = k1 + 1;
end       %end of SNR

A = 10*log10( abs(Gain)/length(signalRange) );

plot(LNRrange,A(1,:),'-',LNRrange, A(2,:),'--',LNRrange, A(3,:),'-.',LNRrange, A(4,:),':');
xlabel('{\itLNR} (dB)','Fontsize',14)
ylabel('{\itE}(gain) (dB)','Fontsize',14)
title('LCMP\_DL(DIR), u_i=+/-0.30, u_a~U(-0.1, 0.1), u_c=[0, -0.0866, 0.0866], g=[1, Bc, Bc], INR=20dB')
grid
axis([-20 50 0 50])
h=legend('{\itSNR}=0 dB','{\itSNR}=10 dB','{\itSNR}=20 dB','{\itSNR}=30 dB');
set(h,'Fontsize',12)

LCMP-DL波束形成器的阵列增益 - 阿英 - Mr.Right
  评论这张
 
阅读(424)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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