使用matlab编写2psk的调制解调，并绘制调制解调波形图和蒙特卡罗仿真

以下是使用MATLAB编写2PSK调制解调的代码，包括调制函数、解调函数以及蒙特卡罗仿真代码。在调制解调波形图中，我们绘制了原始信号、调制信号和解调信号的波形。在蒙特卡罗仿真中，我们生成了1000个随机的二进制信息比特流，并对其进行2PSK调制和解调，然后计算误码率并绘制BER曲线。

% 2PSK调制解调

% 生成随机的二进制信息比特流
data = randi([0 1], 1, 1000);

% 调制函数
function modulated_signal = modulate_2psk(data)
    T = 1;               % 信号持续时间
    fc = 2;              % 载波频率
    Ts = 0.01;           % 采样时间间隔
    t = 0:Ts:T-Ts;       % 时间序列
    w = 2*pi*fc;         % 角频率
    modulated_signal = sqrt(2/T) * cos(w*t + pi*(data==0));
end

% 解调函数
function demodulated_data = demodulate_2psk(received_signal)
    T = 1;               % 信号持续时间
    fc = 2;              % 载波频率
    Ts = 0.01;           % 采样时间间隔
    t = 0:Ts:T-Ts;       % 时间序列
    w = 2*pi*fc;         % 角频率
    demodulated_data = (received_signal * sqrt(2/T) .* cos(w*t)) > 0;
end

% 生成调制信号和解调信号
data = randi([0 1], 1, 1000);
modulated_signal = modulate_2psk(data);
received_signal = awgn(modulated_signal, 10);

% 绘制波形图
subplot(3,1,1);
plot(modulated_signal);
title('Modulated signal');
subplot(3,1,2);
plot(received_signal);
title('Received signal');
subplot(3,1,3);
plot(demodulate_2psk(received_signal));
title('Demodulated signal');

% 蒙特卡罗仿真
nBits = 1000;               % 比特流长度
SNR = 0:2:20;               % 信噪比范围
nSimulations = 1000;        % 仿真次数
ber = zeros(size(SNR));     % 误码率

for i = 1:length(SNR)
    for j = 1:nSimulations
        % 生成随机的二进制信息比特流
        data = randi([0 1], 1, nBits);
        
        % 调制信号
        modulated_signal = modulate_2psk(data);
        
        % 添加高斯白噪声
        noise = sqrt(0.5/(10^(SNR(i)/10))) * randn(size(modulated_signal));
        received_signal = modulated_signal + noise;
        
        % 解调信号
        demodulated_data = demodulate_2psk(received_signal);
        
        % 统计误码率
        ber(i) = ber(i) + sum(data ~= demodulated_data) / nBits;
    end
    ber(i) = ber(i) / nSimulations;
end

% 绘制BER曲线
semilogy(SNR, ber);
xlabel('SNR (dB)');
ylabel('BER');
title('BER vs SNR for 2PSK modulation');