基于Alamouti 编码的 M-PSK 信号通过莱斯平坦衰落信道传输附matlab代码

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

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

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

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

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

⛄ 内容介绍

基于Alamouti 编码的 M-PSK 信号通过莱斯平坦衰落信道传输附matlab代码

⛄ 完整代码

clear all; close all; clc

snapshots = 100000;

EbNo = 0:10;

K = [4.0; 0.6];

M = [4; 8]; %Positions of modulation (M-PSK)

Mt = 2;

Mr = [1; 2];

ostbcEnc = comm.OSTBCEncoder('NumTransmitAntennas', Mt); % Alamouti

ric_ber = zeros(length(EbNo), length(M), length(K), length(Mr));

sum_BER = zeros(length(EbNo), length(M), length(K), length(Mr));

for mr = 1:length(Mr)

    ostbcComb = comm.OSTBCCombiner('NumTransmitAntennas', Mt, 'NumReceiveAntennas', Mr(mr));

    H = zeros(Mr(mr), Mt, snapshots);

    ric_msg = zeros(snapshots, Mr(mr));

    for k = 1:length(K)

        mu = sqrt( K(k)/(K(k)+1));

        s = sqrt(1/(K(k)+1));

        for m = 1:length(M)

            hModulator = comm.PSKModulator('ModulationOrder', M(m), 'BitInput', false); 

            hDemod = comm.PSKDemodulator('ModulationOrder', M(m), 'BitOutput', false);

            ric_ber(:,m,k,mr) = berfading(EbNo, 'psk', M(m), Mr(mr)*Mt, K(k));

            snr = EbNo+10*log10(log2(M(m)));

            message = randi([0,M(m)-1],100000,1);

            mod_msg = step(hModulator,message);

            Es = mean(abs(mod_msg).^2);

            alam_msg = step(ostbcEnc, mod_msg);

            % Channel

            h = mu + s*(1/sqrt(2))*(randn(Mr(mr),Mt,snapshots/Mt)...

            + 1j*randn(Mr(mr),Mt, snapshots/Mt));

            H(:,:,1:2:end-1) = h;

            H(:,:,2:2:end) = h;

            pathGainself = permute(H,[3,2,1]);

            % Transmit through the channel

            for q = 1:snapshots;  

                ric_msg(q,:) = (sqrt(Es/Mt)*H(:,:,q)*alam_msg(q,:).').';

            end

            for c = 1:100

                for jj = 1:length(EbNo)

                    noisy_mod = awgn(ric_msg,snr(jj),'measured','dB');

                    decodeData = step(ostbcComb,noisy_mod,pathGainself);

                    demod_msg = step(hDemod,decodeData);

                    [number,BER(c,jj)] = biterr(message,demod_msg);

                end

            end

            sum_BER(:,m, k, mr) = sum(BER)./c;

        end

    end

end

figure(1)

semilogy(EbNo,sum_BER(:,1,1,1),'r-o',EbNo,sum_BER(:,2,1,1),'g-o',...

         EbNo,ric_ber(:,1,1,1),'r-',EbNo,ric_ber(:,2,1,1),'g-',...

         EbNo,sum_BER(:,1,1,2),'b-o',EbNo,sum_BER(:,2,1,2),'y-o',...

         EbNo,ric_ber(:,1,1,2),'b-',EbNo,ric_ber(:,2,1,2),'y-',...

         'LineWidth', 1.5) 

title('Rician model (K = 4.0)') 

legend('QPSK(simulated) 2x1', '8-PSK(simulated) 2x1',...

    'QPSK(theory) 2x1','8-PSK(theory) 2x1',...

    'QPSK(simulated) 2x2', '8-PSK(simulated) 2x2',...

    'QPSK(theory) 2x2','8-PSK(theory) 2x2') 

xlabel('EbNo (dB)') 

ylabel('BER')

grid on

figure(2) 

semilogy(EbNo,sum_BER(:,1,2,1),'r-o',EbNo,sum_BER(:,2,2,1),'g-o',...

         EbNo,ric_ber(:,1,2,1),'r-',EbNo,ric_ber(:,2,2,1),'g-',...

         EbNo,sum_BER(:,1,2,2),'b-o',EbNo,sum_BER(:,2,2,2),'y-o',...

         EbNo,ric_ber(:,1,2,2),'b-',EbNo,ric_ber(:,2,2,2),'y-',...

         'LineWidth', 1.5) 

title('Rician model (K = 0.6)') 

legend('QPSK(simulated) 2x1', '8-PSK(simulated) 2x1',...

    'QPSK(theory) 2x1','8-PSK(theory) 2x1',...

    'QPSK(simulated) 2x2', '8-PSK(simulated) 2x2',...

    'QPSK(theory) 2x2','8-PSK(theory) 2x2')  

xlabel('EbNo (dB)') 

ylabel('BER')

grid on

⛄ 运行结果

⛄ 参考文献

[1]于子甲, 董全德, 徐旭. 基于Alamouti的MIMO系统仿真设计[J]. 贵阳学院学报：自然科学版, 2020, 15(1):5.

[21]李万臣, 鄂颖. 基于Alamouti空时码的网络编码协作技术的研究[J].  2013.

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

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