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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

MPDR with DOA mismatch  

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

  下载LOFTER 我的照片书  |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.35
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
BWNN = 4/N;
signalRange = (0:0.25/100:0.25)*BWNN;

Vm = exp(j*n*pi*0);

k1 = 1;
for ASNR = 0:10:30
   M = 10^(ASNR/10);
   k2 = 1;
   for ua = signalRange
      Va = exp(j*n*pi*ua);
      Bc = 1/N*Vm'*Va;
      A(k1,k2) = N*abs(Bc)^2/( 1 + ( 2*M + M^2 )*( 1 - abs(Bc)^2 )  );
MPDR with DOA mismatch - 阿英 - Mr.Right
 
MPDR with DOA mismatch - 阿英 - Mr.Right
 
      k2 = k2 + 1;
   end
   k1  = k1 + 1;
end

A = 10*log10(A);
plot(signalRange/BWNN, A(1,:),'-',signalRange/BWNN, A(2,:),'--',signalRange/BWNN, A(3,:),'-.')
hold on
plot(signalRange/BWNN, A(4,:),':')


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Gain coresponding to CRB
hold on
 
h2 = 1;
for K = [2*N, 10*N]
   h1 = 1;
   for ASNR = 0:2:40
  M = 10^(ASNR/10);
  a = 3/(8*K);
    b = inv( 1 - 1/N^2);
  delta(h1) = sqrt(a*(1/M + 1/M^2)*b);
   
%   [x,I] = min( abs(signalRange/BWNN-delta(h1)) )
    Va = exp( j*n*pi*delta(h1)*BWNN );
  Bc = 1/N*Vm'*Va;
    CR_Gain(h2,h1) = N*abs(Bc)^2/( 1 + ( 2*M + M^2 )*( 1 - abs(Bc)^2 )  );
%   CR_Gain(h1) = A(h1,I);
   
  h1 = h1 + 1;
   end
   h2 = h2 + 1;
   
end

CR_Gain = 10*log10(CR_Gain);
plot(delta(4:length(delta)),CR_Gain(1,4:length(delta)), '-o',delta(6:length(delta)),CR_Gain(2,6:length(delta)), '-p');
axis([0 0.25 -60 10])
grid
h=legend('{\itASNR}=0 dB','{\itASNR}=10 dB','{\itASNR}=20 dB','{\itASNR}=30 dB','CRB,{\itK}=2{\itN}','CRB,{\itK}=10{\itN}',4);
set(h,'Fontsize',12)
xlabel('{\itu}_{\ite} /{\itBW}_{\itNN}','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
title('MPDR with DOA mismatch, N=10')
axis([0 0.1 -40 20])
set(gca, 'xtick',[0:0.01:0.25])
set(gca, 'ytick',[-40:5:20])
MPDR with DOA mismatch - 阿英 - Mr.Right
 
   
  评论这张
 
阅读(430)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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