基于MATLAB的三种分集合并技术的实现

目录

等增益合并

 最大比合并

选择合并

 三种合并方式对比分析

等增益合并

N=10^6;
data=randi([0,1],1,N);
x=2*data-1;s=zeros(1,N);
ber_EGC=zeros(1,N);for snr_dB=1:20error=0;h1=randn(1,N)+(randn(1,N)*1i);h2=randn(1,N)+(randn(1,N)*1i);c1=h1.*x;c2=h2.*x;y1=awgn(c1,snr_dB,'measured');y2=awgn(c2,snr_dB,'measured');y_rec = y1.*exp(-1i*angle(h1)) + y2 .*exp(-1i*angle(h2)); % removing the phase of the channelH=abs(h1)+abs(h2);for i=1:Nif(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))s(i)=-1;elses(i)=1;endif(x(i)~=s(i))error=error+1;else error=error+0;endendber_EGC(snr_dB)=error/N;end         
figure;
semilogy(ber_EGC,'go-');
legend('EGC');
xlabel('SNR');
ylabel('BER');
title('BER for 1tx vs 2rx');

 最大比合并

N=10^6;
data=randi([0,1],1,N);
x=2*data-1;s=zeros(1,N);
ber_MRC=zeros(1,N);for snr_dB=1:1:20error=0;h1=randn(1,N)+(randn(1,N)*1i);h2=randn(1,N)+(randn(1,N)*1i);c1=h1.*x;c2=h2.*x;y1=awgn(c1,snr_dB,'measured');y2=awgn(c2,snr_dB,'measured');y_rec =conj(h1).*y1+conj(h2).*y2;H=abs(h1)+abs(h2);for i=1:1:Nif(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))s(i)=-1;elses(i)=1;endif(x(i)~=s(i))error=error+1;else error=error+0;endendber_MRC(snr_dB)=error/N;end         figure;
semilogy(ber_MRC,'k*-');
legend('MRC');
xlabel('SNR');
ylabel('BER');
title('BER for 1tx vs 2rx'); 

选择合并

N=10^6;
data=randi([0,1],1,N);
x=2*data-1;s=zeros(1,N);
ber_SC=zeros(1,N);
for snr_dB=1:20h1=randn(1,N)+(randn(1,N)*1i);h2=randn(1,N)+(randn(1,N)*1i);c1=h1.*x;c2=h2.*x;y1=awgn(c1,snr_dB,'measured');y2=awgn(c2,snr_dB,'measured');hPower1 = sum(h1.*conj(h1));hPower2 = sum(h2.*conj(h2));if(hPower1>hPower2) y_rec=y1;H=h1;elsey_rec=y2;H=h2;enderror=0;          for i=1:Nif(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))s(i)=-1;elses(i)=1;endif(x(i)~=s(i))error=error+1;else error=error+0;endendber_SC(snr_dB)=error/N;
endfigure;
semilogy(ber_SC,'m*-');hold on;
legend('SC');
xlabel('SNR');
ylabel('BER');
title('BER for 1tx vs 2rx');     

 三种合并方式对比分析

figure(2);
semilogy(ber_EGC,'go-');hold on;
semilogy(ber_MRC,'k*-');hold on;
semilogy(ber_SC,'m*-');hold on;
legend('EGC','MRC','SC');
ylabel('BER');
xlabel('SNR');
title('BER for 1tx vs 2rx');

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！