matlab 4PSK调制代码

以下是一个简单的 MATLAB 4PSK 调制代码：

% 定义调制参数
M = 4; % 调制阶数
k = log2(M); % 每个符号的比特数
EbNo = 10; % 信噪比
numBits = 10000; % 发送的比特数

% 生成随机比特序列
dataIn = randi([0 1], numBits, 1);

% 平均能量 per symbol
Es = 1;

% 计算每个符号的能量
Eb = Es/k;

% 生成一个 4PSK 映射表
grayCode = [0 1 3 2];
phaseOffset = pi/4;
constellation = sqrt(Es) * exp(1i*(grayCode*2*pi/M + phaseOffset));

% 将比特序列转换为符号序列
dataSymbolsIn = bi2de(reshape(dataIn, k, length(dataIn)/k).', 'left-msb');

% 将符号映射到 4PSK 星座图
txSignal = constellation(dataSymbolsIn+1).';

% 添加高斯白噪声
snr = EbNo + 10*log10(k) - 10*log10(2);
rxSignal = awgn(txSignal, snr, 'measured');

% 将接收信号解调为符号序列
[~, dataSymbolsOut] = min(abs(rxSignal - constellation), [], 2);
dataOut = de2bi(dataSymbolsOut-1, k, 'left-msb').';
dataOut = dataOut(:);

% 计算误比特率
[numErrors, ber] = biterr(dataIn, dataOut);

% 输出 BER 和调制星座图
disp(['误比特率 = ', num2str(ber)]);
scatterplot(rxSignal);
title('4PSK 星座图');