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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

LCMV波束形成器  

2012-03-16 11:20:46|  分类: 学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.47
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Example 6.7.5 LCMV beamformer
%This routine is to test the Array Gain, when there is signal mismatch


clear all
close all

N = 10;
BWNN = 4/N;
n = (-(N-1)/2:(N-1)/2)';
SNR = 10^(20/10);
SINRi = SNR;

signalRange = 0.25*BWNN*(0:1/100:1);

detU = 0.0866;
C = exp(j*n*pi*[0 -detU detU]);
f1 = [1;1;1];
f2 = [1;sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];
f3= [1;1/2+1/2*sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));1/2+1/2*sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];
f4= [1;2/3+1/3*sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));2/3+1/3*sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];
f = [f1,f2,f3,f4];

for k1 = 1:4
    g = f(:,k1);
    k2 = 1;
    for ua = signalRange
        Vs = exp(j*n*pi*ua);
        Sn = eye(N); %LCMV, here only white noise is available
        W = inv(Sn)*C*inv(C'*inv(Sn)*C)*g;           %
        Ss = SNR*Vs*Vs';
        A(k1,k2) = ((W'*Ss*W)/(W'*Sn*W))/SINRi; % white noise
        k2 = k2+1;
    end
end


A = 10*log10(abs(A));
x = signalRange/BWNN;
plot(x,A(1,:),'-',x,A(2,:),'--',x,A(3,:),'-.',x,A(4,:),':');            %BWnn = 4/N
2012年03月16日 - 阿英 - Mr.Right
 
grid
xlabel('{\itu}_{\ita} /{\itBW}_{\itNN}','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
title('LCMV-DIR,WhiteNoise, SNR=20dB, uc=[0,+/-0.0866]')
axis([0 0.25 4 10])
h=legend('{\itg}_{{\itc}1}','{\itg}_{{\itc}2}','{\itg}_{{\itc}3}','{\itg}_{{\itc}4}',3);
set(h,'Fontsize',12)
line([0.0866/BWNN 0.0866/BWNN],[0 9])   
text(0.0866/BWNN+0.01,9,'{\itu}_{{\itc}1}','Fontsize',12)
set(gca,'xtick',0:0.02:0.25)
set(gca, 'ytick',4:0.5:10)
2012年03月16日 - 阿英 - Mr.Right
 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.48
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Example 6.7.6 LCMP beamformer, No Loading
%This routine is to test the Array Gain, when there is signal mismatch
%

clear all
close all

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


detU = 0.0866;
C = exp(j*n*pi*[0 -detU detU]);
f = [1;sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];

m = 1;
for SNR = 10.^([-20:10:30]/10)
   k = 1;
   for ua = signalRange
      Vs = exp(j*n*pi*ua);
      Ss = SNR*Vs*Vs';
      Sn = eye(N);
      Sx = Ss + Sn;
      W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;           %LCMP, white noise only
      A(m,k) = ((W'*Ss*W)/(W'*Sn*W) )/ SNR; % white noise
      k = k+1;
   end
   m = m + 1;
end

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

grid
xlabel('{\itu}_{\ita} /{\itBW}_{\itNN}','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
h=legend('{\itSNR}=-20 dB','{\itSNR}=-10 dB','{\itSNR}=0 dB','{\itSNR}=10 dB','{\itSNR}=20 dB','{\itSNR}=30 dB',3);
set(h,'Fontsize',12)
title('LCMP, White Noise, uc=[0 +/-0.0866], g=[1;Bc;Bc]')
axis([0 0.25 -50 20])     
line([0.0866/BWNN 0.0866/BWNN],[-80 15])  
text(0.0866/BWNN+0.01,14,'{\itu}_{{\itc}1}','Fontsize',12)
set(gca,'xtick',0:0.02:0.25)
  

LCMV波束形成器 - 阿英 - Mr.Right


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.49
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Example 6.7.7 LCMP beamformer, No Loading
%This routine is to test the Array Gain, when there is signal mismatch
% two interferers are added

clear all
close all

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


detU = 0.0866;
C = exp(j*n*pi*[0 -detU detU]);
f = [1;sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];

ui = [-0.3, 0.3];
Vi = exp(j*n*pi*ui);
INR = 10^(10/10);
Sn = INR*Vi*Vi' + eye(N);

m = 1;
for SNR = 10.^([0:10:30]/10)
   SINRi = SNR/(1+2*INR);
   k = 1;
   for ua = signalRange
      Vs = exp(j*n*pi*ua);
      Ss = SNR*Vs*Vs';
      Sx = Ss + Sn;
      W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;           %LCMP, white noise only
      A(m,k) = ((W'*Ss*W)/(W'*Sn*W) )/SINRi; % white noise
      k = k+1;
   end
   m = m + 1;
end

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

grid
xlabel('{\itu}_{\ita} /{\itBW}_{\itNN}','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
h=legend('{\itSNR}=0 dB','{\itSNR}=10 dB','{\itSNR}=20 dB','{\itSNR}=30 dB',3);
set(h,'Fontsize',12)
%title('LCMP-DIR, ui=+/-0.30(10dB each), uc=[0 +/-0.0866], g=[1;Bc;Bc]')
line([0.0866/BWNN 0.0866/BWNN],[-20 25])  
text(0.0866/BWNN+0.001,23,'{\itu}_{{\itc}1}','Fontsize',12)
set(gca,'xtick',0:0.02:0.25)
axis([0 0.25 -20 30])      
  

LCMV波束形成器 - 阿英 - Mr.Right

 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.50
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Example 6.7.7 LCMP beamformer, No Loading
%This routine is to test the Array Gain, when there is signal mismatch
% two interferers are added

clear all
close all

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


detU = 0.0866;
C = exp(j*n*pi*[0 -detU detU]);
f = [1;sin(-(N/2)*pi*detU)./(N*sin(-.5*pi*detU));sin((N/2)*pi*detU)./(N*sin(.5*pi*detU))];

ui = [-0.3, 0.3];
Vi = exp(j*n*pi*ui);
INR = 10^(30/10);
Sn = INR*Vi*Vi' + eye(N);

m = 1;
for SNR = 10.^([0:10:30]/10)
   SINRi = SNR/(1+2*INR);
   k = 1;
   for ua = signalRange
      Vs = exp(j*n*pi*ua);
      Ss = SNR*Vs*Vs';
      Sx = Ss + Sn;
      W = inv(Sx)*C*inv(C'*inv(Sx)*C)*f;           %LCMP, white noise only
      A(m,k) = ((W'*Ss*W)/(W'*Sn*W) )/SINRi; % white noise
      k = k+1;
   end
   m = m + 1;
end

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

grid
xlabel('{\itu}_{\ita} /{\itBW}_{\itNN}','Fontsize',14)
ylabel('Array gain (dB)','Fontsize',14)
h=legend('{\itSNR}=0 dB','{\itSNR}=10 dB','{\itSNR}=20 dB','{\itSNR}=30 dB',3);
set(h,'Fontsize',12)
title('LCMP-DIR, ui=+/-0.30(30dB each), uc=[0 +/-0.0866], g=[1;Bc;Bc]')
line([0.0866/BWNN 0.0866/BWNN],[0 48])  
text(0.0866/BWNN+0.01,45,'{\itu}_{{\itc}1}','Fontsize',12)
set(gca,'xtick',0:0.02:0.25)
axis([0 0.25 0 50])     
  

  评论这张
 
阅读(565)| 评论(1)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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