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

Mr.Right

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

 
 
 

日志

 
 
关于我

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

网易考拉推荐

DMR波束形成器  

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

  下载LOFTER 我的照片书  |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Example 6.8.3  Dm=3
% Dominant Mode Rejection: Array Gain vs. us/BWnn
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all

N=10;
BWnn=4/N;
INR1=10.^([30]/10);
INR2=10.^([30]/10);
SNR=10.^([-10:10:30]/10);

n=(-(N-1)/2:(N-1)/2)';
us=[0:0.002:0.25];
ui1=0.3;
ui2=-0.3;
vi1=exp(i*n*pi*ui1);
vi2=exp(i*n*pi*ui2);

vm=ones(N,1);
Dm=3;
e=1;
w1=ones(N,1)/N;
%............................ MPDR and DMR ..............................
Ri=INR1*vi1*vi1'+INR2*vi2*vi2'+eye(N);
for numsnr=1:size(SNR,2)
    for num=1:size(us,2)
        vs=exp(i*n*pi*us(num));
        Rs=SNR(numsnr)*vs*vs';
        Rx=Rs+Ri;
        %....................DMR.............
        [E,eigv]=eig(Rx);
        [Y,index]=sort(diag(real(eigv)));
        alpha=sum(Y(1:N-Dm))/(N-Dm);
        Rd=0;
        for numd=1:Dm
            Rd=Rd+e*(Y(N-numd+1)-alpha)*E(:,index(N-numd+1))*E(:,index(N-numd+1))';
        end
        Rd=Rd+alpha*eye(N);

        w3=inv(Rd)*vm/(vm'*inv(Rd)*vm);
        ArrayGain(numsnr,num)=real(w3'*Rs*w3/(w3'*Ri*w3))/(SNR(numsnr)/(INR1+INR2+1));
    end;
end
figure
x=us/BWnn;
y1=10*log10(ArrayGain(1,:));
y2=10*log10(ArrayGain(2,:));
y3=10*log10(ArrayGain(3,:));
y4=10*log10(ArrayGain(4,:));
y5=10*log10(ArrayGain(5,:));
plot(x,y1,x,y2,'--',x,y3,'-.',x,y4,':',x,y5,'-*');
grid on
axis([0 0.25 -30 50])
legend('{\itSNR} = -10 dB','{\itSNR} = 0 dB','{\itSNR} = 10 dB','{\itSNR} = 20 dB','{\itSNR} = 30 dB')
xlabel('{\itu}_{\ita}/{\itBW}_{\itNN}')
ylabel('Array Gain (dB)')
DMR波束形成器 - 阿英 - Mr.Right
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Example 6.8.4 Dm=2
% Dominant Mode Rejection: Array Gain vs. us/BWnn
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all

N=10;
BWnn=4/N;
INR1=10.^([30]/10);
INR2=10.^([30]/10);
SNR=10.^([-10:10:30]/10);

n=(-(N-1)/2:(N-1)/2)';
us=[0:0.002:0.25];
ui1=0.3;
ui2=-0.3;
vi1=exp(i*n*pi*ui1);
vi2=exp(i*n*pi*ui2);

vm=ones(N,1);
Dm=2;
e=1;
w1=ones(N,1)/N;
%............................ MPDR and DMR ..............................
Ri=INR1*vi1*vi1'+INR2*vi2*vi2'+eye(N);
for numsnr=1:size(SNR,2)
    for num=1:size(us,2)
        vs=exp(i*n*pi*us(num));
        Rs=SNR(numsnr)*vs*vs';
        Rx=Rs+Ri;
        %....................DMR.............
        [E,eigv]=eig(Rx);
        [Y,index]=sort(diag(real(eigv)));
        alpha=sum(Y(1:N-Dm))/(N-Dm);
        Rd=0;
        for numd=1:Dm
            Rd=Rd+e*(Y(N-numd+1)-alpha)*E(:,index(N-numd+1))*E(:,index(N-numd+1))';
        end
        Rd=Rd+alpha*eye(N);

        w3=inv(Rd)*vm/(vm'*inv(Rd)*vm);
        ArrayGain(numsnr,num)=real(w3'*Rs*w3/(w3'*Ri*w3))/(SNR(numsnr)/(INR1+INR2+1));
    end;
end
figure
x=us/BWnn;
y1=10*log10(ArrayGain(1,:));
y2=10*log10(ArrayGain(2,:));
y3=10*log10(ArrayGain(3,:));
y4=10*log10(ArrayGain(4,:));
y5=10*log10(ArrayGain(5,:));
plot(x,y1,x,y2,'--',x,y3,'-.',x,y4,':',x,y5,'-*');
grid on
axis([0 0.25 -30 50])
legend('{\itSNR} = -10 dB','{\itSNR} = 0 dB','{\itSNR} = 10 dB','{\itSNR} = 20 dB','{\itSNR} = 30 dB','location','southeast')
xlabel('{\itu}_{\ita}/{\itBW}_{\itNN}')
ylabel('Array Gain (dB)')

DMR波束形成器 - 阿英 - Mr.Right
 
  评论这张
 
阅读(492)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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