ofdm系统中多普勒扩展引起的ICI MATLAB代码

OFDM系统中多普勒扩展会引起干扰，称为多普勒扩展引起的干扰（Inter-Carrier Interference，ICI）。下面是OFDM系统中多普勒扩展引起的ICI的MATLAB代码实现。

% OFDM系统中多普勒扩展引起的ICI的MATLAB代码实现
clc; clear all; close all;

% 设置参数
N = 64;                 % 子载波数量
M = 16;                 % 星座数
SNR = 15;               % 信噪比
delay = [0 1 3 5];      % 多径时延
gain = [0 -2 -3 -4];    % 多径增益
fd = 1.5;               % 多普勒频移

% 生成OFDM信号
data = randi([0 M-1], N, 1);     % 生成随机数据
pilot = repmat([1;-1;1;-1],N/4,1);  % 生成导频
x = ifft(data);                  % IFFT变换
x = [pilot;x];                   % 插入导频
x = x(:);                        % 转为列向量

% 生成多普勒频移信号
Fs = 100;                        % 采样频率
t = (0:length(x)-1)/Fs;          % 时域采样点
fshift = exp(1j*2*pi*fd*t.');   % 生成频移信号

% 信号传输
h = rayleighchan(1/Fs, 100, delay, gain);   % 生成多径信道
h.StoreHistory = true;                      % 保存信道历史
y = filter(h,x.*fshift);                    % 多普勒频移和信道叠加

% 接收端处理
r = y.*conj(fshift);                        % 多普勒频移补偿
r = r(length(pilot)+1:end);                 % 去除导频

% 解调
rx = fft(r);                                % FFT变换
data_rx = qamdemod(rx, M);                  % 解调

% 处理ICI
f = -N/2:N/2-1;                             % 频域采样点
H_f = freqz(h.PathGains, 1, f, Fs);         % 频率响应
H_f_shift = freqz(h.PathGains, 1, f-fd, Fs);% 频率响应（多普勒频移后）
H_i = H_f_shift ./ H_f;                     % 计算ICI影响
I = fftshift(ifft(H_i));                    % ICI影响时域响应
ICI = filter(I, 1, r);                      % ICI处理

% 处理结果
err = sum(data_rx ~= data);                 % 比较接收数据和原始数据的误码率
BER = err/length(data);                     % 比特误率

% 输出结果
fprintf('误码率为: %.4f\n', BER);

上述代码中首先生成了一个包含随机数据和导频的OFDM信号，然后生成了一个多普勒频移信号。接下来将OFDM信号和多普勒频移信号进行叠加，并通过多径信道进行传输。接收端接收到信号后，首先进行多普勒频移补偿，然后进行FFT解调。接下来计算ICI影响，并进行ICI处理。最后计算误码率并输出结果。