基于RMDCFT算法的天基雷达空间机动目标检测方法MATLAB仿真,对比FRFT和DFT算法

目录

✨1.课题概述

📊2.系统仿真结果

✅3.核心程序或模型

🚀4.系统原理简介

💢5.完整工程文件


✨1.课题概述

RMDCFT将Radon距离徙动校正与修正离散Chirp傅里叶变换耦合,构建包络对齐 - 相位匹配一体化相参积累框架,从数学层面同时抵消天基机动目标回波中的距离单元走动与二次多普勒调频相位,彻底解决DFT无法处理加速度、FRFT无法补偿距离跨徙的固有缺陷。整套算法从预处理、包络校正、多普勒匹配、恒虚警检测到参数精估形成完整闭环,每一步均具备严谨数学表达式支撑。虽然RMDCFT相比DFT、FRFT存在更高多维搜索计算开销,但通过分层搜索与快速变换优化可大幅降低运算负载,在天基雷达远距离、低信噪比、高速机动目标探测场景下具备不可替代的检测增益优势,是当前空间目标长时相参积累领域的主流改进变换类算法。

📊2.系统仿真结果

✅3.核心程序或模型

bash 复制代码
fc = 1e6;
Bs = 30e6;
v0 = [3050:-10:2850];
a0 = 400;
R0 = [0:10:600];
SNR= -10;
c  = 3e8;
L  = 800;
Ar = 0.9;
t  = 0e-06:80e-06/(length(v0)-1):80e-06;
Tr = 1/1500;
Ff = 0;
Sr = zeros(length(t),length(Ff));
M  = 64;
for i = 1:length(t)
    for j = 1:length(R0)
        tmps = 0;
        for k = 1:M
            tm(j)   = (((Ff)/2*(c/fc) - v0(i))/a0)/M;
            Rf(j)   = (R0(j) + v0(i)*tm(j) + 0.5*a0*tm(j)^2);
            tmps    = tmps + Ar*sinc(Bs*(t(i)-2*Rf(j)/c))*exp(sqrt(-1)*4*pi*fc/c*Rf(j));
        end
        Sr(i,j) = tmps;
    end
end
Sr2 = awgn(Sr,SNR,'measured');

[X,Y] = meshgrid(v0,R0);


figure;
subplot(121);
imagesc(abs(Sr1))
colorbar;
axis square
xlabel('脉冲数');
ylabel('距离单元');

subplot(122);
meshc(X,Y,abs(Sr2'))
axis square
axis([2850,3050,0,600,0,60]);
xlabel('速度');
ylabel('距离单元');
zlabel('幅度');
view([40,30]);
 
save R1.mat Sr1 X Y Sr2

🚀4.系统原理简介

天基雷达平台处于轨道高速运动状态,空间机动目标包含径向速度、径向加速度甚至加加速度,长时相参积累过程中会同时产生距离单元走动 RCM与多普勒单元徙动 DCM,回波慢时间维度呈现线性调频(Chirp)相位特征。设雷达脉冲重复周期为Tr​,积累脉冲总数M,快时间采样间隔Ts​,光速c,载波波长λ;目标初始距离R0​、径向速度v、径向加速度a,第m个脉冲时刻目标径向距离为:

RMDCFT全称Radon修正离散Chirp傅里叶变换(Radon Modified Discrete Chirp-Fourier Transform),融合Radon变换包络校正与MDCFT相位匹配滤波,同步完成距离走动包络对齐与线性调频多普勒相位相参积累,专门适配天基雷达高速机动目标二次相位耦合距离徙动回波。

MDCFT在传统DCFT基础上消除采样周期与调频率耦合约束,构造二维匹配核同时匹配多普勒中心频率f与调频率γ,单距离门离散MDCFT变换式:

匹配核共轭抵消回波慢时间一次、二次相位,当搜索参数(f,γ)与目标真实(fd​,γ0​)完全匹配时,求和项内相位全部归零,各脉冲回波同相叠加形成尖锐峰值,实现线性调频信号全能量积累。传统 DCFT存在周期模糊,MDCFT引入归一化修正因子消除调频率周期截断误差,旁瓣抑制能力提升 8~12dB。

Radon变换沿慢时间维度对距离 - 慢时间矩阵做直线积分,补偿由径向速度、加速度引发的距离单元跨徙。对距离慢时间二维矩阵s(k,m),给定速度搜索量vp​、加速度搜索量ap​,Radon包络映射函数为:

该步骤将分散在不同距离单元的同一目标包络对齐至同一距离门,彻底消除距离走动带来的能量分散,是RMDCFT区别于FRFT、MDCFT的核心创新。

恒虚警检测CFAR峰值提取

构造四维参数空间幅值矩阵,采用单元平均恒虚警CA-CFAR检测,检测阈值计算式:

遍历四维空间所有像素,若∣S(vp​,ap​,f,γ)∣>Tth​,判定该位置存在空间机动目标,同时输出对应搜索参数作为目标运动参数估计值。

对比

DFT仅适用于近程、短积累、无加速度平稳目标,无法用于天基预警雷达长时相参积累;FRFT适用于地面雷达近距离小幅机动目标,距离走动微弱可忽略场景;RMDCFT专门适配天基 / 高轨雷达高速、大加速度、长积累周期空间机动目标,同步处理距离弯曲与多普勒扩展,是三类算法中唯一适配天基复杂运动回波的处理框架。

相关推荐
noipp2 小时前
推荐题目:洛谷 B2099 矩阵交换行
线性代数·算法·矩阵
五条凪3 小时前
简单理解 BM25 与 TF-IDF
人工智能·算法·搜索引擎·全文检索·tf-idf
行者全栈架构师4 小时前
【码动四季】Spring Boot 可观测性体系:Micrometer + OpenTelemetry + Grafana 全链路搭建
java·算法·架构
TCW11214 小时前
AI底层系列:用C++实现线性代数的公式推导与算法设计-8.线性变化(3)
c++·人工智能·算法
皓月斯语4 小时前
B3849 [GESP样题 三级] 进制转换 题解
c++·算法·题解
中微极客4 小时前
剪枝与量化:让YOLO在边缘设备上高效部署
算法·yolo·剪枝
牢姐与蒯5 小时前
双指针算法
数据结构·算法
Hesionberger5 小时前
LeetCode406:重建身高队列精髓解析
开发语言·数据结构·python·算法·leetcode