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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

MPDR波束形成器的阵列增益和SNR的关系  

2012-03-15 22:58:16|  分类: 学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.36a
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Ex 6.6.4 
%MPDR, with array perturbation
%Relation between signal mismatch, SNR and Array Gain is studied

clear all
close all

N = 10;
BWNN = 4/N;
n = (-(N-1)/2:(N-1)/2)';

% uncomment to get figure in text (takes a long time to run)
SNRrange = [-20:1:50];
%SNRrange = [-20:10:50];

totalTrialNumber = 500;

Ui1 = 0.30;
Us = 0;

INR1 = 10^(20/10);

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


k1 = 1;
for sigmaP = [0, 0.05, 0.1, 0.2]
   disp(['loop ' int2str(k1) ' of 4 ...'])
   k2 = 1;
   for SNR = 10.^(SNRrange/10)
      Gain(k1,k2) = 0;
      for trial = 1:totalTrialNumber
            ds = (n/2+sigmaP*randn(N,1))*Us  + sigmaP*randn(N,1)*sqrt(1-Us^2);
         d1 = (n/2+sigmaP*randn(N,1))*Ui1 + sigmaP*randn(N,1)*sqrt(1-Ui1^2);
        Vi1 = exp(j*2*pi*d1);
        Vs  = exp(j*2*pi*ds);
        Ss = SNR*Vs*Vs';
        Sn = eye(N) + INR1*Vi1*Vi1';
        Sx = Ss + Sn;    
        W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;           %consider White Noise only
        SINR0 = (W'*Ss*W)/(W'*Sn*W);
        SINRi = SNR/(1+INR1);
         Gain(k1,k2) = Gain(k1,k2) + SINR0/SINRi;
      end
      k2 = k2 + 1;  
   end
   k1 = k1 + 1;   
end
Gain = 10*log10(real(Gain)/totalTrialNumber);

plot(SNRrange,Gain(1,:),'-',SNRrange,Gain(2,:),'--',SNRrange,Gain(3,:),'-.',SNRrange,Gain(4,:),':')
grid
xlabel('{\itSNR} (dB)','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14);
h=legend('\sigma_{\itp} = 0','\sigma_{\itp} = 0.05\lambda','\sigma_{\itp} = 0.10\lambda','\sigma_{\itp} = 0.20\lambda',3);
set(h,'Fontsize',12)
axis([-20 50 -100 40])
title(['MPDR with perturbation., ui=',num2str(Ui1),...
      '(',num2str(10*log10(INR1)),'dB), ua=',num2str(Us),...
      ', N=',num2str(N)])


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.36b
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Ex 6.6.4 
%MPDR, with array perturbation
%Relation between signal mismatch, SNR and Array Gain is studied

clear all
close all

N = 10;
BWNN = 4/N;
n = (-(N-1)/2:(N-1)/2)';
% uncomment to get figure in text (takes a long time to run)
SNRrange = [-20:1:50];
%SNRrange = [-20:10:50];

totalTrialNumber = 500;

Ui1 = 0.30;
Us = 0;

INR1 = 10^(30/10);

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


k1 = 1;
for sigmaP = [0, 0.05, 0.1, 0.2]
   disp(['loop ' int2str(k1) ' of 4 ...'])
   k2 = 1;
   for SNR = 10.^(SNRrange/10)
      Gain(k1,k2) = 0;
      for trial = 1:totalTrialNumber
            ds = (n/2+sigmaP*randn(N,1))*Us  + sigmaP*randn(N,1)*sqrt(1-Us^2);
         d1 = (n/2+sigmaP*randn(N,1))*Ui1 + sigmaP*randn(N,1)*sqrt(1-Ui1^2);
        Vi1 = exp(j*2*pi*d1);
        Vs  = exp(j*2*pi*ds);
        Ss = SNR*Vs*Vs';
        Sn = eye(N) + INR1*Vi1*Vi1';
        Sx = Ss + Sn;    
        W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;           %consider White Noise only
        SINR0 = (W'*Ss*W)/(W'*Sn*W);
        SINRi = SNR/(1+INR1);
         Gain(k1,k2) = Gain(k1,k2) + SINR0/SINRi;
      end
      k2 = k2 + 1;  
   end
   k1 = k1 + 1;   
end
Gain = 10*log10(real(Gain)/totalTrialNumber);

plot(SNRrange,Gain(1,:),'-',SNRrange,Gain(2,:),'--',SNRrange,Gain(3,:),'-.',SNRrange,Gain(4,:),':')
grid
xlabel('{\itSNR} (dB)','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14);
h=legend('\sigma_{\itp} = 0','\sigma_{\itp} = 0.05\lambda','\sigma_{\itp} = 0.10\lambda','\sigma_{\itp} = 0.20\lambda',3);
set(h,'Fontsize',12)
axis([-20 50 -80 60])
title(['MPDR with perturbation., ui=',num2str(Ui1),...
      '(',num2str(10*log10(INR1)),'dB), ua=',num2str(Us),...
      ', N=',num2str(N)])

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

历史上的今天

评论

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

页脚

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