2DPSK的matlab仿真

以下是一个简单的2DPSK Matlab仿真的示例代码：

% 数据生成
M = 2;  % 调制阶数
N = 1000;  % 数据长度
data = randi([0 M-1], N, 1);  % 生成随机数据

% 星座图
constellation = exp(1j*2*pi*(0:M-1)/M);

% 符号映射
sym = constellation(data+1);

% 差分编码
diff_sym = [sym(1); diff(sym)];

% 信道模型
snr_db = 10;
snr = 10^(snr_db/10);
noise = sqrt(1/2*1/snr)*(randn(size(sym))+1j*randn(size(sym)));
recv_sym = sym + noise;

% 差分解码
diff_recv_sym = [recv_sym(1); diff(recv_sym)];
diff_recv_sym(1) = [];

% 判决
recv_data = zeros(size(data));
for i = 1:length(data)
    [~, idx] = min(abs(constellation-diff_recv_sym(i)));
    recv_data(i) = idx-1;
end

% BER计算
ber = sum(recv_data~=data)/N;
disp(['SNR(dB)=' num2str(snr_db) ', BER=' num2str(ber)])

注：以上代码中的星座图为常规的QPSK星座图，即M=2时的星座图。若要绘制其他调制阶数的星座图，请自行修改constellation的取值。