【雷达干扰】基于CFAR的FMCW雷达干扰抑制附Matlab代码

最新推荐文章于 2024-07-19 15:52:23 发布

天天Matlab代码科研顾问

最新推荐文章于 2024-07-19 15:52:23 发布

阅读量1.3k

点赞数 27

分类专栏：雷达文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/m0_57702748/article/details/136097448

版权

雷达专栏收录该内容

140 篇文章 22 订阅

订阅专栏

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，代码获取、论文复现及科研仿真合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab完整代码及仿真定制内容点击👇

智能优化算法神经网络预测雷达通信无线传感器电力系统

信号处理图像处理路径规划元胞自动机无人机

🔥 内容介绍

1. 前言

随着雷达技术的发展，雷达干扰也变得越来越普遍。雷达干扰可以分为主动干扰和被动干扰。主动干扰是指干扰机主动向雷达发射干扰信号，以降低雷达的探测性能。被动干扰是指干扰机利用雷达发射的信号进行反射或散射，以产生虚假目标，从而干扰雷达的正常工作。

FMCW 雷达是一种新型的雷达，具有抗干扰能力强、分辨率高、测距精度高等优点。然而，FMCW 雷达也存在着一些弱点，例如容易受到连续波干扰的影响。连续波干扰是一种常见的主动干扰，它可以使 FMCW 雷达的接收机饱和，从而降低雷达的探测性能。

2. 基于 CFAR 的 FMCW 雷达干扰抑制方法

为了抑制连续波干扰，可以采用基于 CFAR（恒虚警率）的干扰抑制方法。CFAR 是一种自适应滤波技术，它可以根据雷达接收信号的统计特性，自动调整滤波器的参数，以抑制干扰信号，同时保持对目标信号的检测性能。

基于 CFAR 的 FMCW 雷达干扰抑制方法主要包括以下几个步骤：

信号预处理：对雷达接收信号进行预处理，包括滤波、去噪等。
背景噪声估计：估计雷达接收信号中的背景噪声功率。
门限计算：根据背景噪声功率和虚警率，计算干扰抑制的门限。
干扰信号检测：将雷达接收信号与门限进行比较，检测出干扰信号。
干扰信号抑制：对检测出的干扰信号进行抑制，可以采用滤波、抵消等方法。

📣 部分代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CFAR-based interference suppression%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%close all; clear;  clc;oldpath = path;addpath('./Util_CFAR_IM/','./Util_FMCW/','./Util_plot/')flag_plot = 0;  % boolean flag that turns on/off the figure plotting:                  %        flag_plot=1, plot figure;                %        flag_plot=0, do not plot.Ftsz = 12;fig_dir = './figs/';%% FMCW radarii = sqrt(-1);c      = 3e8;f_c    = 77e9; %3e9;             % center frequency 3GHzT_sw   = 100e-6;%500e-6;          % Sweep time 200 usBW     = 600e6;%40e6;Num    = 1;6, real(sig_CFAR_Burg(ind_inset_cfar)),'g-.')grid onylim([-3.2 3.2])if flag_plot == 1    set(gcf, 'PaperUnits', 'inches', 'PaperSize', [4 3], 'PaperPosition', [0 0 4 3])    print(hf, '-dpdf', '-r300', [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_Recv.pdf'])    print(hf, '-dpng', '-r300', [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_Recv.png'])    saveas(hf, [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_Recv.fig'])end%% Range profilen_ft = 2^(nextpow2(length(sig_wd)) + 2);r_axis_Aft = f_s/n_ft*(0:n_ft-1)/abs(sweep_slope)*c/2;%================ANC=================[rp_ANC, sig_fft_pos,sig_fft_neg] ...    = AdaptiveNoiseCanceller(sig_full, 10, 80, 1.8, n_ft);%rp_CFAR_Z  = fft(sig_CFAR_Z, n_ft);rp_CFAR_AC = fft(sig_CFAR_AC, n_ft);rp_CFAR_Burg = fft(sig_CFAR_Burg, n_ft);rp_wd      = fft(sig_wd, n_ft);max_rp = max(abs( [ rp_CFAR_Z; rp_CFAR_AC; rp_wd ; rp_ANC] ));rp_CFAR_Z_norm  = db( abs(rp_CFAR_Z) / max_rp );rp_CFAR_AC_norm = db( abs(rp_CFAR_AC) / max_rp );rp_wd_norm      = db( abs(rp_wd) / max_rp ); rp_ANC_norm     = db( abs(rp_ANC) / max_rp);rp_CFAR_Burg_norm = db( abs(rp_CFAR_Burg) / max_rp);hf = figureplot(r_axis_Aft(1:n_ft/2), rp_ANC_norm, 'g--',...     r_axis_Aft, rp_wd_norm,'k-',...        r_axis_Aft, rp_CFAR_Z_norm,'r--',...     r_axis_Aft, rp_CFAR_AC_norm,'b-.',...     r_axis_Aft, rp_CFAR_Burg_norm,'m:')axis([0 250 -50 10]);grid onlegend('ANC','WD','CFAR-Z','CFAR-AC','CFAR-Burg','Location','south','NumColumns',2)xlabel('Range [km]', 'fontsize', Ftsz)ylabel('Normalized amplitude [dB]', 'fontsize', Ftsz)title('RP after IM', 'fontsize', Ftsz)ind_inset = r_axis_Aft>=29 & r_axis_Aft<=31;r_axis_inset = r_axis_Aft(ind_inset);axes('position', [0.3 0.65 0.25 0.25])plot(r_axis_inset, rp_ANC_norm(ind_inset),'g--',...     r_axis_inset, rp_wd_norm(ind_inset),'k-',...     r_axis_inset, rp_CFAR_Z_norm(ind_inset),'r--',...     r_axis_inset, rp_CFAR_AC_norm(ind_inset),'b-.',...     r_axis_inset, rp_CFAR_Burg_norm(ind_inset),'m:')grid onylim([-20 2])ind_inset = r_axis_Aft>=148.5 & r_axis_Aft<=154.5;r_axis_inset = r_axis_Aft(ind_inset);axes('position', [0.66 0.65 0.23 0.25])plot(r_axis_inset, rp_ANC_norm(ind_inset),'g--',...     r_axis_inset, rp_wd_norm(ind_inset),'k-',...     r_axis_inset, rp_CFAR_Z_norm(ind_inset),'r--',...     r_axis_inset, rp_CFAR_AC_norm(ind_inset),'b-.',...     r_axis_inset, rp_CFAR_Burg_norm(ind_inset),'m:')grid onylim([-25 0])if flag_plot==1    set(gcf, 'PaperUnits', 'inches', 'PaperSize', [4 3], 'PaperPosition', [0 0 4 3])    print(hf, '-dpdf', '-r300', [fig_dir 'pt_RP_AftSup_' num2str(SNR) 'dB.pdf'])    print(hf, '-dpng', '-r300', [fig_dir 'pt_RP_AftSup_' num2str(SNR) 'dB.png'])    saveas(hf, [fig_dir 'pt_RP_AftSup_' num2str(SNR) 'dB.fig'])end%% Performance evaluationlen_cfar_z = length(sig_CFAR_Z);% SINRSINR_cfar_z = db(norm(sig_Rx_trc(1:len_cfar_z)) ...                 / norm(sig_CFAR_Z - sig_Rx_trc(1:len_cfar_z).' ) );SINR_cfar_ac = db( norm(sig_Rx_trc(1:len_cfar_z) ...                 / norm(sig_CFAR_AC - sig_Rx_trc(1:len_cfar_z).' ) ) );SINR_cfar_burg = db( norm(sig_Rx_trc(1:len_cfar_z) ...                 / norm(sig_CFAR_Burg - sig_Rx_trc(1:len_cfar_z).' ) ) );             SINR_wd = db( norm(sig_Rx_trc) / norm(sig_wd - sig_Rx_trc.') );sig_ANC = ifft(rp_ANC, n_ft);sig_ANC = sig_ANC(1:N_sig);SINR_anc = db( norm(sig_Rx_trc) / norm(sig_ANC - sig_Rx_trc.' ) );% correlation coefficientrho_cfar_z = sig_CFAR_Z' * sig_Rx_trc(1:len_cfar_z).'/norm(sig_CFAR_Z)/norm(sig_Rx_trc(1:len_cfar_z));rho_cfar_ac = sig_CFAR_AC' * sig_Rx_trc(1:len_cfar_z).'/norm(sig_CFAR_AC) /norm(sig_Rx_trc(1:len_cfar_z));rho_cfar_burg = sig_CFAR_Burg' * sig_Rx_trc(1:len_cfar_z).'/norm(sig_CFAR_Burg) /norm(sig_Rx_trc(1:len_cfar_z));rho_wd      = sig_wd' * sig_Rx_trc.'/norm(sig_wd) / norm(sig_Rx_trc);rho_anc     = sig_ANC' * sig_Rx_trc.' / norm(sig_ANC) / norm(sig_Rx_trc);[SINR_cfar_z, SINR_cfar_ac,SINR_cfar_burg, SINR_wd, SINR_anc]rho =[rho_cfar_z, rho_cfar_ac, rho_cfar_burg, rho_wd, rho_anc];rho_amp = abs(rho)rho_phase = angle(rho)% rho = rho_amp .* exp(ii*rho_phase)%close allhf = figureplot(t*1e6, real(sig_ANC));grid onaxis tight% ylim([-2.5 2.5])xlabel('Time [\mus]', 'fontsize', Ftsz)ylabel('Amplitude', 'fontsize', Ftsz)title('Beat signal recovered with ANC', 'fontsize', Ftsz)if flag_plot == 1    set(gcf, 'PaperUnits', 'inches', 'PaperSize', [4 3], 'PaperPosition', [0 0 4 3])    print(hf, '-dpdf', '-r300', [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_ANC.pdf'])    print(hf, '-dpng', '-r300', [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_ANC.png'])    saveas(hf, [fig_dir 'pt_RawSig_RealPart_AftSup_' num2str(SNR) 'dB_ANC.fig'])endif flag_plot==1save([fig_dir 'simu_pt_data_' num2str(SNR) '.mat'],...    'sig_Rx_trc','sig_CFAR_AC','sig_CFAR_Z', 'sig_ANC','sig_wd',...    'SINR_cfar_ac','SINR_cfar_z','SINR_anc','SINR_wd',...    'rho_cfar_ac', 'rho_cfar_z', 'rho_anc', 'rho_wd',...    'rp_CFAR_AC', 'rp_CFAR_Z', 'rp_ANC', 'rp_wd',...    'SNR', 'SINR_0', 'r_axis_Aft', 'n_ft', 't');end%% restore the search directorypath = oldpath;