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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

协方差矩阵锐化  

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

  下载LOFTER 我的照片书  |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figures 6.68 and 6.69
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Analysis of covariance matrix taper

close all
clear all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
SNR = 10^(0/10);
INR = 10^(30/10);

ui = 0.3;
Vi = exp(j*n*pi*ui);

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

BWNN = 4/N;

Ss = SNR*Va*Va';
Sx = Ss + INR*Vi*Vi' + eye(N);

delta_uiRange = -0.1:1/1000:0.1;
u = delta_uiRange;

gammaRange = [0:0.01:0.03];

k1 = 1;
for gamma = gammaRange 
   
   for h1 = 1:N
      for h2 = 1:N
         Sx(h1,h2)  = Sx(h1,h2)*sinc((h1-h2)*gamma/pi);   %here, sinc(x) = sin(x)/x;
      end
   end %Tapered covariance matrix
   
   W = inv(Sx)*Va/(Va'*inv(Sx)*Va);
   b(k1,:) = W'*exp(j*n*pi*(ui+u));  %波束方向图
   
   k2 = 1;
   for uia = ui + delta_uiRange;
   
      Via = exp(j*n*pi*uia);
      Sn = INR*Via*Via' + eye(N);
      SINRo = real(W'*Ss*W)/real(W'*Sn*W);
      SINRi = SNR/(1+INR);
      Gain(k1,k2) = SINRo/SINRi;
      
      k2 = k2 + 1;
   end
   
   k1 = k1 + 1;
end

   
x = delta_uiRange/BWNN;
   
Gain = 10*log10(Gain);
b = 20*log10(abs(b));
协方差矩阵锐化 - 阿英 - Mr.Right
 
plot(x, Gain(1,:),'-',x, Gain(2,:),'--',x, Gain(3,:),'-.',x, Gain(4,:),':')
grid
axis([-0.25 0.25 0 45])
h=legend('\gamma=0','\gamma=0.01','\gamma=0.02','\gamma=0.03',4);
set(h,'Fontsize',12)
xlabel('\Delta{\itu}_{\iti} /{\itBW}_{\itNN}','Fontsize',14)   
ylabel('Gain (dB)','Fontsize',14)
title('CMT, ua=0(0dB), ui(normial)=0.3(30dB), N=10')


协方差矩阵锐化 - 阿英 - Mr.Right
 
figure
plot(x, b(1,:),'-',x, b(2,:),'--',x, b(3,:),'-.',x, b(4,:),':')
grid
axis([-0.25 0.25 -100 0])
h=legend('\gamma=0','\gamma=0.01','\gamma=0.02','\gamma=0.03',4);
set(h,'Fontsize',12)
xlabel('\Delta{\itu}_{\iti} /{\itBW}_{\itNN}','Fontsize',14)   
ylabel('Beam pattern (dB)','Fontsize',14)
title('CMT, ua=0(0dB), ui(normial)=0.3(30dB), N=10')

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

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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