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

w360366

已于 2023-11-01 10:15:42 修改

阅读量1.9k

点赞数 6

文章标签： matlab 开发语言

于 2022-11-13 10:07:00 首次发布

本文链接：https://blog.csdn.net/wangwttwtt/article/details/127828905

版权

等增益合并

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:20
            error=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 channel
                                
            H=abs(h1)+abs(h2);
            for i=1:N
                if(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))
                    s(i)=-1;
                else
                    s(i)=1;
                end

                if(x(i)~=s(i))
                    error=error+1;
                else 
                    error=error+0;
                end
            end
                    
            ber_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:20
            error=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:N
                if(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))
                    s(i)=-1;
                else
                    s(i)=1;
                end

                if(x(i)~=s(i))
                    error=error+1;
                else 
                    error=error+0;
                end
            end
                    
            ber_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:20
            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');
         
           hPower1 = sum(h1.*conj(h1));
           hPower2 = sum(h2.*conj(h2));
           if(hPower1>hPower2) 
               y_rec=y1;
               H=h1;
           else
               y_rec=y2;
               H=h2;
           end
                
           error=0;          
           for i=1:N
                if(abs(y_rec(i)-H(i))>abs(y_rec(i)+H(i)))
                    s(i)=-1;
                else
                    s(i)=1;
                end

                if(x(i)~=s(i))
                    error=error+1;
                else 
                    error=error+0;
                end
            end
                    
            ber_SC(snr_dB)=error/N;
end

figure;
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');