matlab编写2psk、qpsk、8psk调制相干解调的蒙特卡罗仿真，并将误码曲线绘制在同一个图中

以下是一个简单的 Matlab 代码，可以用于模拟2PSK、QPSK、8PSK的相干调制和解调，并绘制误码率曲线。

clear all;
close all;
clc;

% 设置模拟参数
SNR_dB = 0:2:20; % 信噪比范围
num_bits = 1e6; % 传输比特数
M = [2 4 8]; % 调制阶数

% 初始化误码率向量
BER = zeros(length(M), length(SNR_dB));

for m = 1:length(M)
    % 生成调制映射表
    switch M(m)
        case 2
            hMod = comm.PSKModulator(M(m), 'BitInput', true);
            hDemod = comm.PSKDemodulator(M(m), 'BitOutput', true);
        case 4
            hMod = comm.QPSKModulator('BitInput', true);
            hDemod = comm.QPSKDemodulator('BitOutput', true);
        case 8
            hMod = comm.PSKModulator(M(m), 'BitInput', true, 'SymbolMapping', 'gray');
            hDemod = comm.PSKDemodulator(M(m), 'BitOutput', true, 'SymbolMapping', 'gray');
    end
    
    for i = 1:length(SNR_dB)
        % 生成随机比特流
        data = randi([0 1], num_bits, 1);
        
        % 调制
        modSignal = step(hMod, data);
        
        % 加噪声
        snr = 10^(SNR_dB(i)/10);
        noisePower = 1/sqrt(2*snr);
        noisySignal = modSignal + noisePower*randn(size(modSignal));
        
        % 解调
        demodData = step(hDemod, noisySignal);
        
        % 计算误码率
        [~, BER(m, i)] = biterr(data, demodData);
    end
end

% 绘制误码率曲线
figure;
semilogy(SNR_dB, BER(1,:), 'o-', 'LineWidth', 2);
hold on;
semilogy(SNR_dB, BER(2,:), 's-', 'LineWidth', 2);
semilogy(SNR_dB, BER(3,:), 'd-', 'LineWidth', 2);
grid on;
title('误码率曲线');
xlabel('信噪比 (dB)');
ylabel('误码率');
legend('2PSK', 'QPSK', '8PSK');

运行上述代码，将生成一个包含三条误码率曲线的图像，对应于2PSK、QPSK和8PSK。

需要注意的是，在实际应用中，误码率曲线的形状可能会受到许多因素的影响，包括信道特性、调制方式、码率、编码方式等。因此，需要根据具体的应用场景进行调整和优化。