不同多普勒频移下的ofdm误码率代码

以下是一个MATLAB实现的OFDM系统的误码率仿真代码，考虑了不同多普勒频移的影响：

clear all; close all; clc;

% OFDM参数定义
N = 64; % 子载波数
cp_len = 16; % 循环前缀长度
frame_len = 100; % 帧长
mod_order = 4; % 调制阶数

% 信道参数定义
fd = [0 20 100]; % 多普勒频移
tau = [0 3 6] * 1e-6; % 时延
pdb = [0 -3 -6]; % 衰减

% 生成随机数据
tx_data = randi([0 mod_order-1], 1, N*frame_len);

% 调制
tx_mod = qammod(tx_data, mod_order);

% OFDM调制
tx_ofdm = reshape(tx_mod, N, frame_len);
tx_ofdm_ifft = ifft(tx_ofdm);

% 循环前缀
tx_ofdm_cp = [tx_ofdm_ifft(N-cp_len+1:N,:); tx_ofdm_ifft];

% 信道模型
for i = 1:length(fd)
    h(:,i) = rayleigh_chan(N, fd(i), tau(i), pdb(i), frame_len);
end

% 加噪声并接收
rx_ofdm_cp = zeros(size(tx_ofdm_cp));
for i = 1:length(fd)
    rx_ofdm_cp = rx_ofdm_cp + awgn(conv(h(:,i), tx_ofdm_cp), 10, 'measured');
end

% 去掉循环前缀
rx_ofdm_ifft = rx_ofdm_cp(cp_len+1:cp_len+N, :);

% OFDM解调
rx_ofdm = fft(rx_ofdm_ifft);

% 解调得到的符号
rx_mod = reshape(rx_ofdm, 1, N*frame_len);

% 误码率计算
num_err = sum(abs(tx_data-rx_mod) > 0);
ber = num_err / (N*frame_len);
fprintf('误码率为 %f\n', ber);

其中，`rayleigh_chan` 函数用于生成瑞利衰落信道。这里只是一个简单的实现，可以根据需要进行修改。