基于Matlab模拟独立瑞利衰落下双分支分集接收机QPSK的误码率

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

现代通信系统中,对于通信距离、容量和质量的要求,调制解调技术是其中的关键技术之一,一直是人们研究的重点和热点。

⛄ 完整代码

x = 1; % signal to transmit Eb = 1

TRIAL = 10000; %number of simulation runs per EbN0 %50000

for EbN0 = 0:1:20 %dB

    linear_EbN0 = 10^(EbN0/10); nvar = 1/(linear_EbN0); %calculation of N0, remember Eb = 1

    error1 = 0; %set error counter to 0

    error2 = 0; %set error counter to 0

    error3 = 0; %set error counter to 0

        for trial = 1:TRIAL % monte carlo trials.. count the errors

            n1 = sqrt(nvar/2)*randn; %noise for the first

            n2 = sqrt(nvar/2)*randn; %noise for the first

            h1 = sqrt(0.5)*abs(randn + j*randn); %rayleigh amplitude 1

            h2 = sqrt(0.5)*abs(randn + j*randn); %rayleigh amplitude 1

            %Equal Gain combining

            y1 = x*h1+n1; % Signal 1

            y2 = x*h2+n2; % Signal 2

            y_equal = 0.5*(y1+y2); 

            %Maximal Ratio combining

            a1 = (abs(h1))^2;

            a2 = (abs(h2))^2;

            y_maximal = x*(a1*h1+a2*h2)+a1*n1+a2*n2;

            %Selection combining

            P1 = chi2rnd(4);

            P2 = chi2rnd(4);

            as1 = P1*(abs(h1))^2;

            as2 = P2*(abs(h2))^2;

            if as1 >= as2

                y_selection = x*(as1*h1)+as1*n1;

            end

            if as1 < as2

                y_selection = x*(as2*h2)+as2*n2;

            end

            if y_equal < 0 %define decision region as 0 

                error1 = error1 + 1;

            end

            if y_maximal < 0 

                error2 = error2 + 1;

            end

            if y_selection < 0 

                error3 = error2 + 1;

            end

        end

    BER1(EbN0+1) = error1/(TRIAL);

    BER2(EbN0+1) = error2/(TRIAL);

    BER3(EbN0+1) = error3/(TRIAL);

end

% plot simulations

figure

EbNo=0:1:20; %changed from 10

mu = 10.^(EbNo./10);

ber_theory = (1/2)*(1 - sqrt(mu ./ (mu + 1))); 

semilogy(EbNo,BER1,'r*-',EbNo,BER2,'b--o',EbNo,BER3,'c-o',EbNo,ber_theory,'b'); % plot EG BER vs EbNo 

legend('EG','MR','SC','theory');

xlabel('EbNo(dB)') %Label for x-axis

ylabel('Bit error rate') %Label for y-axis

⛄ 运行结果

⛄ 参考文献

[1] 韩英娜. 基于FPGA的空时分组码的设计与实现[D]. 大连海事大学, 2010.

[2] 蔡莉莉, 王进华. 基于Matlab的独立瑞利衰落模型的仿真[J]. 江西通信科技, 2008(2):4.

[3] 刘芳. MATLAB多径衰落对误码率影响的仿真与分析[J]. 大陆桥视野, 2010(06X):2.

[4] 朱黄锁. 基于FPGA的QPSK中频接收机的研究与实现[D]. 西安电子科技大学, 2013.

[5] 马翘楚. 基于Matlab的QPSK系统的设计[D]. 吉林大学.

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料