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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

噪声环境下阵列增益与方位角的关系  

2012-03-15 20:47:01|  分类: 学习 |  标签: |举报 |字号 订阅

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

clear all
close all

%************************
% Array
%************************
N = 10;                                 % Elements in array
d = 0.5;                                % sensor spacing half wavelength wrt wc
D = [-(N-1)/2:1:(N-1)/2].';
BWNN = 2/(N*d);
u=[-1:0.001:1];
nu=length(u);
vv = exp(j*2*pi*d*D*u);

%************************
% Source
%************************

INR = 10^(20/10);
alpha = 1;

theta_s = [0:0.01:1]*pi;
us = cos(theta_s);
AS = exp(j*2*pi*d*D*us);
ns = length(us);
A = zeros(1,ns);
Ac = zeros(1,ns);
for n=1:ns
   p = [0:1:N-1];
   pI = [1:1:N-1];
   r = sinc(p*2*d)+[0 ((j*alpha)./(pi*pI*2*d)).*(sinc(pI*2*d)-cos(pi*2*d*pI))];
噪声环境下阵列增益与方位角的关系 - 阿英 - Mr.Right
    Sn = INR*toeplitz(r,conj(r))+eye(N); %以前不知道共轭对称的谱矩阵Sn是用toeplitz函数生成的!
   Sninv = inv(Sn);
   Ac(n) =N*N*Sn(1,1)/ real(AS(:,n)'*Sn*AS(:,n)); 
   A(n) = real(AS(:,n)'*Sninv*AS(:,n))*Sn(1,1); 
end
figure
h1=plot(theta_s*180/pi,10*log10(A),'-');
hold on
h2=plot(theta_s*180/pi,10*log10(Ac),'--');
h=legend('MVDR', 'Conventional',4) ;
set(h,'Fontsize',12)
xlabel('\theta_{\its} (degrees)','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
grid on
hold off
%axis([0 180 5 20])
噪声环境下阵列增益与方位角的关系 - 阿英 - Mr.Right
噪声环境下阵列增益与方位角的关系 - 阿英 - Mr.Right
 

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.25
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

%************************
% Array
%************************
N = 10;                                 % Elements in array
d = 0.5;                                % sensor spacing half wavelength wrt wc
D = [-(N-1)/2:1:(N-1)/2].';
BWNN = 2/(N*d);
u=[-1:0.001:1];
nu=length(u);
vv = exp(j*2*pi*d*D*u);

%************************
% Source
%************************
INR = 10^(20/10);
alpha = 1;

theta_s = [0:0.01:1]*pi;
us = cos(theta_s);
AS = exp(j*2*pi*d*D*us);
ns = length(us);
A = zeros(1,ns);
Ac = zeros(1,ns);
for n=1:ns
   p = [0:1:N-1];
   pI = [1:1:N-1];
   r = sinc(p*2*d)+[0 alpha*((3./(pi*pI*2*d).^2)-1).*(sinc(pI*2*d))-alpha*cos(pi*2*d*pI)./(pi*pI*2*d).^2];
噪声环境下阵列增益与方位角的关系 - 阿英 - Mr.Right
 
   Sn = INR*toeplitz(r,conj(r))+eye(N); %以前不知道共轭对称的谱矩阵Sn是用toeplitz函数生成的!
   Sninv = inv(Sn);
   Ac(n) =N*N*Sn(1,1)/ real(AS(:,n)'*Sn*AS(:,n)); 
   A(n) = real(AS(:,n)'*Sninv*AS(:,n))*Sn(1,1); 
end
figure
h1=plot(theta_s*180/pi,10*log10(A),'-');
hold on
h2=plot(theta_s*180/pi,10*log10(Ac),'--');
h=legend('MVDR', 'Conventional',4) ;
set(h,'Fontsize',12)
xlabel('\theta_{\its} (degrees)','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
grid on
hold off
噪声环境下阵列增益与方位角的关系 - 阿英 - Mr.Right
 

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

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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