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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

MVDR波束方向图--对称干扰  

2012-03-15 00:12:11|  分类: 默认分类 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.21
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
vs = ones(N,1);

u = -1:0.01:1;
v = exp(j*n*pi*u);   % exp(j*2*pi/lambda * d * u), 若 d = 0.5*lambda
SNR = 1;
INR = 10.^([20]/10);

ui = [0.3 0.18 0.02];  % 干扰的位置,传波方向的方向余弦是-au是传播方向的反方向的方向余弦,是个单位向量

for i=1:3
    VI = [exp(j*n*pi*ui(i)) exp(-j*n*pi*ui(i))];   % 对称干扰是这样的加的,我以前不知道啊
    R = INR*VI*VI'+eye(N);
   
    w = inv(vs'*inv(R)*vs)*inv(R)*vs;
    B = w'*v;
   
    figure
    plot(u,10*log10(abs(B).^2),'-')
    hold on
    if i==3
        plot(ui(i)*[1 1],[-10 20],'--')
        plot(-ui(i)*[1 1],[-10 20],'--')
        axis([-1 1 -10 30])
    else
        plot(ui(i)*[1 1],[-30 -10],'--')
        plot(-ui(i)*[1 1],[-30 -10],'--')
        axis([-1 1 -30 0])
    end
    xlabel('{\itu}','Fontsize',14)
    ylabel('Beam pattern (dB)','Fontsize',14)
    grid on
end

2012年03月14日 - 阿英 - Mr.Right

 

2012年03月14日 - 阿英 - Mr.Right

 

2012年03月14日 - 阿英 - Mr.Right

 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.22
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

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

ui = BWNN*(0.01:0.01:2.5);
ni = length(ui);

SNR = 1;
INR = 10.^([-20 -10 0 10 20]/10);
nn = length(INR);

Ao = zeros(nn,ni);
for m=1:nn   % length(INR) = 5
    for i=1:ni
        VI = [exp(j*n*pi*ui(i)) exp(-j*n*pi*ui(i))];   % 注意对称干扰怎样加入
        R = INR(m)*VI*VI'+eye(N);
        SINR = SNR/(2*INR(m)+1);
       
        w = inv(vs'*inv(R)*vs)*inv(R)*vs;   %MVDR 波束形成器
        Ao(m,i) = SNR/(N*SINR*real(w'*R*w));  % 最优阵列增益为什么这么算的,我没有看懂
    end
end

plot(ui/BWNN,10*log10(Ao),'-');
xlabel('{\itu}_{\itI} /{\itBW}_{\itNN}','FontSize',14)
ylabel('{\itA}_{\ito} /{\itN} (dB)','FontSize',14)
axis([0 2.5 -10 25])

for m=1:nn
    if m==1
        text(1.5, 10*log10(Ao(m,ni))-1, ['{\itINR}=' num2str(10*log10(INR(m))) ' dB'],'FontSize',12)
    else
        text(1.5, 10*log10(Ao(m,ni))+1, ['{\itINR}=' num2str(10*log10(INR(m))) ' dB'],'FontSize',12)
    end
end

hold off

2012年03月14日 - 阿英 - Mr.Right
  评论这张
 
阅读(542)| 评论(1)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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