本发明专利技术涉及利用稀疏傅里叶变换计算外辐射源雷达互模糊函数的方法。属于雷达目标探测处理领域。首先，将外辐射源雷达天线接收到的经滤波、下抽处理后的直达波信号与回波时延信号点乘，构造新向量；然后，将该新向量进行稀疏傅里叶变换，从而得到互模糊函数在该延时点的多普勒切面结果；通过稀疏傅里叶变换结果可估计目标多普勒频移等参数。与传统的外辐射源雷达借助傅里叶变换进行互模糊函数运算的方法相比，本发明专利技术根据实际情况目标在空中出现个数有限从而呈现稀疏性的特点，利用稀疏傅里叶变换进行互模糊函数求解，可大为降低长时间积累情况下互模糊函数的运算量。
【技术实现步骤摘要】
【技术保护点】
一种利用稀疏傅里叶变换计算外辐射源雷达互模糊函数方法，其特征在于具体步骤为：（1）参数设置：构造信号：er(n)=echo(n+m)·refr*(n)，n∈[1,L]；er(n)表示回波信号echo(n)的m点延时与直达波信号refr(n)共轭点乘；其中，echo(n)为n时刻采样回波信号，refr(n)为n时刻采样直达波信号，*为共轭运算符，m为延时单元数；定义信号：x(i)是对信号er经滤波、下抽处理后数据的第i采样点，i∈[1,N]；其中N为信号长度，N是2的整数次幂；信号x的频域表示为X，假设X包含k个大值点数，大值对应可能的目标；k取整数；（2）频谱排列：定义信号：s(i)=x(σ·i)，i∈[1,N]；其中σ∈[1,N]且是对N存在模逆的奇数。若σ·i>N，则s(i)=x((σ·i)modN)，(σ·i)modN表示σ·i对模N求余数；定义σ?1为σ模N的逆，两者满足关系(σ×σ?1)modN=1；（3）信号滤波定义平坦窗函数g(i)（i∈[1,N]），令g为对称向量，设其有效长度为w，w1且为整数；J中坐标看作以Z为观测域得到的大值“像”位置集合；将J映射回以X为观测域的大值“原像”；“原像”位置集合I={i∈[1,N]|hashσ(i)∈J}；其中，（i∈I）定义为“哈希函数”，表示“原像”与“像”的映射关系；每次映射得到“原像”集合新元素个数d·k·N/B；定义偏移函数：oσ(i)=σ·i?hashσ(i)·N/B，i∈I；执行步骤（5）loc次；（6）估值循环改变σ，执行（2）～（4）；记录步骤（2）中σ，将步骤（4）所得Z归入集合Z_SUM；执行步骤（6）est次；（7）频谱大值估计Z_SUM包含所有步骤（4）得到的子采样FFT中d·k个大值，I与Z_SUM对应以X为观测域的“原像”；对i∈I，X(i)估计值为hashσ(i)、oσ(i)中使用的参数σ来自步骤（5）、（6）的每次循环；执行步骤（7）loops次，其中，loops=est+loc为定位与估计的总循环次数；步骤（7）得到loops个值，取第个大值作为最终对的估计结果。FDA00003359538600011.jpg,FDA00003359538600021.jpg,FDA00003359538600023.jpg,FDA00003359538600024.jpg,FDA00003359538600022.jpg,FDA00003359538600031.jpg,FDA00003359538600033.jpg,FDA00003359538600032.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员：陶然，刘升恒，张果，单涛，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。