%% 参数设置 fs = 98.304e6; % 采样率 subcarriers = 16384; % 子载波数 subband_bw = 3e3; % 子带带宽 subband_fs = 6e3; % 子带采样率 n_symbols = 1024; % 符号数 rolloff = 0.75; % 滚降系数 span = 6; % 跨度 fc = 30e6; % 载波频率 max_doppler = 30; % 最大多普勒频偏 %% 生成随机的 BPSK 信号 data = randi([0 1], 1, subcarriers/4); modulated_data = 1 - 2 * data; %% 上采样 upsampled_data = upsample(modulated_data, subcarriers); %% 成型滤波 sps = subcarriers / n_symbols; h = rcosdesign(rolloff, span, sps, 'sqrt'); filtered_data = filter(h, 1, upsampled_data); %% 子带调制 tx_signal = zeros(1, subcarriers); for k = 1:subcarriers/4 % 计算中心频率 subband_fc = (k-1) / subcarriers * fs; % 考虑多普勒频偏 delta_f = rand * 2 * max_doppler - max_doppler; carrier = exp(1j * 2 * pi * (subband_fc + delta_f - fc) * (0:length(filtered_data)-1)/fs); % 调制 tx_signal(k) = real(filtered_data * carrier.'); end %% 作图 subplot(2,1,1); plot((0:length(modulated_data)-1)/subcarriers, modulated_data, 'LineWidth', 2); title('BPSK 调制后的信号'); xlabel('时间 (s)'); ylabel('幅度'); subplot(2,1,2); f = (-subcarriers/2:subcarriers/2-1) / subcarriers * fs; plot(f, abs(fftshift(fft(tx_signal))), 'LineWidth', 2); xlim([0 30e6]); title('调制后的信号频谱'); xlabel('频率 (Hz)'); ylabel('幅度');

这段代码主要是用来模拟 OFDM（正交频分复用）系统中的信号调制过程。

第一部分是设置参数，包括采样率、子载波数、子带带宽、符号数、滚降系数、跨度、载波频率和最大多普勒频偏等。

第二部分是生成随机的 BPSK 信号，即将信息进行调制。

第三部分是对信号进行上采样和成型滤波，即对信号进行调制前的处理。

第四部分是对信号进行子带调制，即将信号分成若干个子带，每个子带的信号分别进行调制。

最后一部分是对信号进行作图，分别绘制调制后的信号和调制后的信号频谱。