干扰场景下汽车雷达测距和测速分析【附MATLAB代码】

​微信公众号：EW Frontier QQ交流群：949444104

简介

该资源库包含一个综合分析项目，重点研究频率调制连续波（FMCW）雷达干扰对汽车应用中目标探测的影响。其目的是探索不同程度和类型的干扰如何影响基于雷达的测距和测速估计的准确性，而测距和测速是汽车安全和导航系统开发和性能的关键因素。

关键目标

影响分析：研究 FMCW 雷达干扰如何影响目标探测精度，特别是如何影响附近物体的测距和测速。

数据驱动的洞察力：利用收集到的数据，分析干扰如何干扰雷达精确探测和跟踪目标的能力，并将其可视化。

滤波器实施：测试并演示各种过滤技术的有效性，这些技术旨在减轻雷达干扰的负面影响，提高雷达的性能和可靠性。

数据描述

数据集包括受害雷达系统和干扰雷达系统的参数：

受害者雷达参数：

EIRP: -5 dBW

天线增益：10 dBi

起始频率：76.5 千兆赫

带宽：500 兆赫

扫描时间：50 µs

脉冲重复间隔：60 µs

相干处理间隔：10 毫秒

低通滤波器截止频率：3.5 MHz

信号调制：FMCW

天线配置：单节点（初始阶段），TDMA-MIMO（高级阶段）

速度：(0, 0) 米/秒

干扰雷达参数：

EIRP: -5 dBW

起始频率：76.45 千兆赫

带宽：900 兆赫

扫描时间：25 µs

脉冲重复间隔：30 µs

信号调制：FMCW

目标：

目标 1：尺寸（4 米 x 1.9 米 x 1.5 米），速度（10, 0）米/秒，位置（40, 3.6）米

目标 2：尺寸 (4.5 m x 1.9 m x 1.5 m)，速度 (-5, 0) m/s，位置 (20, 7.2) 米

MATLAB部分代码

clear all;
clc;
% Victim radar parameters
fc = 76.5e9; % Start frequency
bw_sv = 500e6; % Bandwidth
tm_sv = 50e-6; % Sweep time
fs_sv = bw_sv; % Sampling frequency
N_sv = 25000; % Number of samples in each chirp
Nd_sv = 135; % Number of chirps
% Initialize variables for later use
N=N_sv;
Nd=Nd_sv;
fs=fs_sv;
c=3e8;
slope=bw_sv/tm_sv; % Frequency slope
lambda = c/fc; % Wavelength (m)
tm=tm_sv;
% Interferer radar parameters
fc_int = 76.45e9; % Start frequency
bw_int = 900e6; % Bandwidth
tm_int = 25e-6; % Sweep time
fs_int = bw_int; % Sampling frequency (Hz)
% Time vectors
t_sv = 0:1/fs_sv:(tm_sv*Nd_sv)-1/fs_sv; % Victim radar
% Transmitted signal (Victim radar)
% --- Target Parameters ---
t1_vel = 10; % Velocity of target 1 (m/s)
t2_vel = -5; % Velocity of target 2 (m/s)
t1_pos = [40, 3.6]; % Initial position of target 1 (m)
t2_pos = [20, 7.2]; % Initial position of target 2 (m)
victim_pos = [0, 0]; % Position of victim radar (m)
% Calculate initial distances from radar to targets
r1 = pdist([victim_pos; t1_pos], 'euclidean');
r2 = pdist([victim_pos;t2_pos ], 'euclidean');
range1 = r1 + t1_vel * t_sv; %change of range with respect to time
delta_t1 = 2 * range1 / 3e8; %change in delay time
t_chirp = 0:1/fs_sv:(tm_sv)-1/fs_sv; %time vectors for single chirp
t_int1=0:1/fs_int:tm_int-1/fs_int;% time vectors for interference radar
t_int= repmat(t_int1,[1,150]);
% --- Power and Amplitude Calculations ---
L = 4; % Length of the longest side of target (m)
L_int = 4.5; % Length of the longest side of interferer (m)
EIRP_dBW = -5; % Effective Isotropic Radiated Power (dBW)
EIRP_dBW_int = -5; % EIRP for interferer (dBW)
G_dBi = 10; % Antenna gain (dBi)
lambda_val = 0.0039; % Wavelength (m)
% Calculate Radar Cross Section (RCS)
rcs = 4 * pi * (L^4) / (3 * (lambda^2)); % RCS for target (m^2)
rcs_int = 4 * pi * (L_int^4) / (3 * (lambda^2)); % RCS for interferer (m^2)
% Convert EIRP and Gr from dB to linear scale
EIRP_linear = 10^(EIRP_dBW/10); % in Watts
EIRP_linear_int = 10^(EIRP_dBW_int/10); % in Watts
G_r_linear = 10^(G_dBi/10); % Linear scale
% Radar equation to compute Pr
Pr = (EIRP_linear * G_r_linear * lambda_val^2 * rcs) / ((4*pi)^3 * r1^4);
A = sqrt(2*Pr); % Amplitude for victim radar
% Radar equation to compute Pr of interference
Pr_int = (EIRP_linear * G_r_linear * lambda_val^2*rcs_int ) / ((4*pi)^2 * r2^2);
A_int = sqrt (2*Pr_int); % Amplitude for recieved radar
% --- Signal Generation for Each Chirp ---
for chirp_idx = 1:Nd_sv-1
 % Update the range for target 1 for each chirp
 r1(chirp_idx+1) = r1(chirp_idx) + t1_vel * tm_sv;
end
for chirp_idx = 1:Nd_sv
 % Calculate the time offset for the current chirp
 time_offset = (chirp_idx - 1) * tm_sv;
 
 % Calculate the time delay for target 1
 delta_t1 = 2 * r1(chirp_idx) / c;
 % Generate the transmitted signal for one chirp
 st_chirp = exp(1j * 2 * pi * (fc * t_chirp + (bw_sv / tm_sv) * (t_chirp.^2)/2));
 % Generate the received signal for one chirp considering target 1
 sr_chirp = A * exp(1j * 2 * pi * (fc * (t_chirp - delta_t1) + (bw_sv / tm_sv)*(t_chirp - delta_t1).^2 / 2 - 2 * fc * t1_vel / c));
 % Add noise to the received signal
 noise_level = A * randn(size(sr_chirp));
 sr_chirp1 = sr_chirp + noise_level;
 % Store the transmitted and received signals for this chirp
 st_all_chirps(:, chirp_idx) = st_chirp;
 sr_all_chirps(:, chirp_idx) = sr_chirp1;
end
% --- Generate Interference Signal ---
slope_int = bw_int / tm_int; % Frequency slope for interferer
% Generate the interference signal
sr_interference = A_int * exp(1j * 2 * pi * (fc_int * t_int + (slope_int / 2) * (t_int.^2) - 2 * fc_int * t2_vel / c));
% Reshape the interference signal to match dimensions
sr_int = reshape(sr_interference, [N, Nd]);
% Add the interference to the received signal
sr = sr_all_chirps + sr_int;

MATLAB仿真结果

添加图片注释，不超过 140 字（可选）