% OFDM参数N = 200; % 子载波数目cp_len = 64; % 循环前缀长度fs = 1e6; % 采样率fc = 800e6; % 载波频率T = 1/fs; % 采样时间Ts = N*T; % OFDM符号时间T_cp = cp_len*T; % 循环前缀时间f = (0:N-1)*fs/N; % 子载波频率t = (0:N-1)*T; % 子载波时间t_cp = (N-cp_len:N-1)*T; % 循环前缀时间% 生成随机数据data = randi([0 15], 1, N);% 将数据进行16QAM调制data_mod = qammod(data, 16);% 将数据进行IFFT变换data_ifft = ifft(data_mod, 512);% 添加循环前缀data_cp = [data_ifft(end-cp_len+1:end) data_ifft];% 将OFDM符号变为基带信号baseband = real(exp(1j*2*pi*fc*t).*data_cp);% 绘制OFDM信号频谱figure;plot(f/1e6, 20*log10(abs(fftshift(fft(baseband)))));xlabel('频率（MHz）');ylabel('幅度（dB）');title('OFDM信号频谱');对于此运算，数组的大小不兼容。 出错 ofdm (第 21 行) baseband = real(exp(1j*2*pi*fc*t).*data_cp);怎么改

您可以尝试将 `t` 和 `data_cp` 的大小改为 `512`，即：

t = (0:511)*T;
data_cp = [data_ifft(end-cp_len+1:end) data_ifft];
data_cp = repmat(data_cp, 1, 2);

这样可以使 `t` 和 `data_cp` 的大小与 `exp(1j*2*pi*fc*t)` 相匹配，从而避免数组大小不兼容的错误。同时，由于循环前缀的长度为 `64`，因此我们需要将 `data_cp` 重复两次，以便 OFDM 符号和循环前缀的长度都为 `512`。

完整的程序如下：

% OFDM参数
N = 200; % 子载波数目
cp_len = 64; % 循环前缀长度
fs = 1e6; % 采样率
fc = 800e6; % 载波频率
T = 1/fs; % 采样时间
Ts = N*T; % OFDM符号时间
T_cp = cp_len*T; % 循环前缀时间
f = (0:N-1)*fs/N; % 子载波频率
t = (0:511)*T; % 子载波时间
t_cp = (N-cp_len:N-1)*T; % 循环前缀时间

% 生成随机数据
data = randi([0 15], 1, N);

% 将数据进行16QAM调制
data_mod = qammod(data, 16);

% 将数据进行IFFT变换
data_ifft = ifft(data_mod, 512);

% 添加循环前缀
data_cp = [data_ifft(end-cp_len+1:end) data_ifft];
data_cp = repmat(data_cp, 1, 2);

% 将OFDM符号变为基带信号
baseband = real(exp(1j*2*pi*fc*t).*data_cp);

% 绘制OFDM信号频谱
figure;
plot(f/1e6, 20*log10(abs(fftshift(fft(baseband)))));
xlabel('频率（MHz）');
ylabel('幅度（dB）');
title('OFDM信号频谱');