128APSK代码 matlab

以下是一个简单的 128APSK 调制代码示例，使用 MATLAB 实现：

clear all;
close all;

M = 16; % 星座点数
k = log2(M); % 每个星座点所需比特数
N = 128; % 帧长
L = 32; % 每个符号的采样数
nsamp = 4; % 采样率
Fs = 10e6; % 采样频率
Fc = 2e6; % 载波频率
Es = 1; % 符号能量
Eb = Es/k; % 比特能量
gamma = 2; % 幅度比
constellation = gray2bin(0:M-1, 'qam', M); % 星座点

% 生成随机比特流
data = randi([0 1], N*k, 1);

% 将比特流转换为星座点
data = bi2de(reshape(data, k, N).', 'left-msb');
data = constellation(data+1).';

% 将星座点映射到 128APSK 符号
r = sqrt(Es./gamma).*data; % 半径
phi = (2*pi/N)*(0:N-1).'; % 角度
x = zeros(N*L, 1);
for ii = 1:M
    r_ii = r(constellation == constellation(ii));
    phi_ii = phi(constellation == constellation(ii));
    x_ii = repmat(r_ii.', L, 1).*exp(1i*repmat(phi_ii.', L, 1));
    x = x + x_ii;
end

% 调制
t = (0:length(x)-1).'/(L*nsamp*Fs);
carrier = exp(1i*2*pi*Fc*t);
x_mod = real(x.*carrier);

% 画图
subplot(2,1,1);
plot(real(x), imag(x), 'o');
title('128APSK 星座图');
xlabel('In-phase');
ylabel('Quadrature');
grid on;
axis square;

subplot(2,1,2);
plot(t, x_mod);
title('128APSK 调制波形');
xlabel('Time');
ylabel('Amplitude');
grid on;

其中，该示例使用 Gray 编码将星座点映射到比特序列中，并使用了均匀分布的星座点。您可以根据需要自行更改这些参数。