当前位置：首页 > news >正文

基于OFDM数据符号的多普勒频率估计与移动台速度提取方法

news 2026/3/26 22:20:52

一、核心原理与数学模型

1.1 OFDM信号的多普勒效应

OFDM信号的每个子载波在传输过程中会受到目标运动引起的多普勒频移影响。假设目标径向速度为v，载波频率为fc，则多普勒频移为：

c为光速，λ为波长。多普勒频移会导致接收信号相位随时间线性变化。

1.2 相位演化与速度关系

对于第k个子载波，接收信号在N个OFDM符号周期内的相位累积为：

其中Ts为符号周期。通过分析相邻符号间的相位差，可提取多普勒频移：

二、算法实现步骤

2.1 信号预处理

% 参数设置
fc = 30e9;      % 载波频率 (Hz)
c = 3e8;        % 光速 (m/s)
lambda = c/fc;  % 波长 (m)
N = 64;         % 子载波数
Ts = 1e-6;      % 符号周期 (s)
N_sym = 100;    % 累积符号数% 接收信号处理
rx_signal = rx_data; % 去除CP后的时域信号
Y = fft(rx_signal, N); % FFT变换到频域

2.2 多普勒频移估计

方法1：FFT频域分析法

% 频域相位差计算
phase_diff = angle(Y(:,2:end) .* conj(Y(:,1:end-1)));
f_d_est = (phase_diff * fc) ./ (2*pi*N*Ts); % 估计多普勒频移

方法2：自相关时域分析法

% 自相关计算
R = xcorr(rx_signal, rx_signal, 'biased');
[~, lag] = max(R(N_sym:end)); % 延迟估计
f_d_est = (lag * fc) / (2*N_sym*Ts); % 多普勒频移计算

2.3 速度计算

v_{est} = \frac{f_d_est \cdot \lambda}{2} \quad \text{(单位：m/s)}

三、MATLAB仿真示例

% 参数设置
c = 3e8; fc = 2.4e9; lambda = c/fc;
v_true = 20; % 真实速度 (m/s)
f_d = 2*v_true*fc/lambda; % 理论多普勒频移% 生成带多普勒的OFDM信号
tx_data = randi([0 1], N, N_sym);
tx_ofdm = ifft(tx_data, N);
tx_signal = zeros(N + 100, N_sym);
for i = 1:N_symtx_signal(:,i) = filter([1 0.5], 1, [tx_ofdm(:,i); zeros(100,1)]);
end% 添加高斯噪声
SNR = 20; noise_power = var(tx_signal(:))/(10^(SNR/10));
rx_signal = awgn(tx_signal(:), SNR, 'measured');% 速度估计
Y = fft(rx_signal(1:N,:), N);
phase_diff = angle(Y(:,2:end) .* conj(Y(:,1:end-1)));
f_d_est = (phase_diff * fc) ./ (2*pi*N*Ts);
v_est = (f_d_est * lambda)/2;% 结果可视化
figure;
subplot(2,1,1);
plot(0:N_sym-1, cumsum(real(tx_signal)), 'b', 0:N_sym-1, cumsum(real(rx_signal)), 'r--');
legend('发射信号', '接收信号');
title('时域信号对比');subplot(2,1,2);
stem(f_d, 1, 'r', 'LineWidth', 2);
hold on;
stem(f_d_est, 0.8, 'b', 'LineWidth', 2);
xlabel('多普勒频移 (Hz)');
legend('理论值', '估计值');