matlab 2psk调制与解调程序并画出误码率（参数码元宽度为39微秒，采样频率为1MHz）

以下是使用MATLAB进行2PSK调制解调并绘制误码率曲线的示例代码，其中参数码元宽度为39微秒，采样频率为1MHz：

% 2PSK调制解调误码率曲线

% 信号参数
fc = 1000; % 信号频率
fs = 1000000; % 采样频率
T = 1/fs; % 采样时间
A = 1; % 信号幅值
T_symbol = 39e-6; % 参数码元宽度

% 调制信号
t = 0:T:T_symbol-T; % 信号时间
m = A * sin(2 * pi * fc * t); % 基带信号
s = A * cos(2 * pi * fc * t); % 2PSK调制信号

% 绘制调制信号
figure;
plot(t, s);
title('2PSK调制信号');
xlabel('时间（s）');
ylabel('幅值');

% AWGN信道模拟
EbNo = 0:2:20; % 信噪比范围
ber = zeros(1, length(EbNo)); % 误码率
for i = 1:length(EbNo)
    % 计算噪声方差
    snr = EbNo(i) + 10*log10(fs*T_symbol);
    sigma = sqrt(A^2/(2*10^(snr/10)));
    % 添加噪声
    r = s + sigma*randn(size(s));
    % 2PSK解调
    r_demod = sign(r).*A;
    % 计算误码率
    [~, ber(i)] = biterr(m > 0, r_demod > 0);
end

% 绘制误码率曲线
figure;
semilogy(EbNo, ber);
title('2PSK调制解调误码率曲线');
xlabel('信噪比（dB）');
ylabel('误码率');
grid on;

运行代码后，将会得到绘制的2PSK调制解调误码率曲线图像。需要注意的是，由于参数码元宽度为39微秒，因此在调制时需要仅取每个码元的第一个采样点进行调制，而在解调时也需要将每个码元的采样点平均起来。