基于先验信息的变凹口宽度mtd滤波器设计方法
【专利摘要】本发明公开了一种基于先验信息的变凹口宽度MTD滤波器设计方法，其步骤包括：1.查寻参数；2.获得地杂波的协方差矩阵；3.获得初始化干扰信号功率；4.设计数字综合法的FIR滤波器；5.搜索地杂波功率谱方差；6.获得最优杂波改善因子的FIR滤波器；7.设计变凹口宽度的MTD滤波器。本发明利用雷达周围场景的先验信息，使滤波器输出端杂波改善因子保持最优的方法。克服了现有技术中的MTD滤波器零频处凹口宽度都固定的问题。使得本发明设计的MTD滤波器对不同的通带中心频率滤波器有不同的期望凹口宽度，而且增强了对于低速运动目标检测的优点。
【专利说明】
基于先验信息的变凹口宽度MTD滤波器设计方法
技术领域
[0001] 本发明属于信号处理技术领域，更进一步涉及杂波抑制技术领域中的一种基于先 验信息的变凹口宽度动目标检测MTD(Moving Target Detection)滤波器设计方法。本发明 可用于雷达对杂波的抑制，实现对目标的检测。
【背景技术】
[0002] 杂波信号的有效抑制是雷达系统进行正常目标检测的前提和基础。在自适应杂波 抑制方法中，动目标检测MTD(Moving Target Detection)是一种常用的技术。动目标检测 MTD滤波器是利用多个带通滤波器组成的滤波器组对雷达回波进行滤波处理，再对滤波器 组输出进行检测来发现目标的。动目标检测MTD滤波器组最简单的实现方法是采用离散傅 里叶变换DFT滤波器组。在离散傅里叶变换DFT滤波器组之前加上动目标显示MTI处理，这样 可以先用动目标显示MTI消掉地杂波，再用离散傅里叶变换DFT滤波器组进行滤波处理，离 散傅里叶变换DFT滤波器组可以用快速傅里叶变换FFT实现。因为快速傅里叶变换FFT可以 节省很大的运算量，所以这种方法至今仍在大量使用，尤其是动目标检测MTD滤波器组阶数 较大时。但是由于快速傅里叶变换FFT的阶数一定要是2的整数幂，使得动目标显示MTI加快 速傅里叶变换FFT实现法的应用受到了一定的限制，并且离散傅里叶变换DFT滤波器组位于 动目标显示MTI之后，滤波器组的各滤波器增益受到动目标显示MTI滤波器频率响应的调 制，所以目前常用的动目标检测MTD滤波器组实现方法是采用限长脉冲响应FIR滤波器组。 对于有限长脉冲响应FIR滤波器组中的每个滤波器来说其幅度频率响应都要在零频附近有 较深的零陷，用于抑制杂波其中主要是地杂波，这就是常规自适应动目标检测MTD滤波器。 但常规自适应动目标检测MTD滤波器通常具有较高的副瓣电平，这样容易造成滤波器组中 不同滤波器的目标之间的互相影响，会带来虚警。
[0003] 谷泓，赵永波等人在其发表的论文"一种基于数字综合算法的MTD滤波器组设计方 法"中提出了一种数字综合算法的动目标检测MTD滤波器设计方法。该方法首先基于自适应 阵列理论，然后在滤波器的副瓣区放置大量的干扰信号，改变干扰信号的功率强度，从而自 适应的控制天线方向图的副瓣电平。该方法所设计的动目标检测MTD组滤波器，其频率响应 在零频处在其零频附近有较宽而深凹口，并且在其他副瓣区域可以为任意形状。该方法存 在的不足之处是，其没有利用地基雷达周围场景杂波的先验信息，其杂波的相关矩阵是通 过查不同场景的地杂波功率谱方差和地杂波功率表得到的，导致所有滤波器零陷的凹口宽 度都是固定的，在实际应用中对于不同的雷达杂波环境，滤波器的凹口宽度不能同杂波信 号的谱宽相匹配，并且当滤波器通带中心频率较低时滤波器零频的凹口会使滤波器的主瓣 发生畸变，滤波器的主瓣将会偏离中心频率，使主瓣增益下降，输出的目标信号有损失，而 不能有效的检测到目标。
【发明内容】
[0004] 本发明的目的在于克服上述已有技术的不足，提出了一种基于先验信息的变凹口 宽度动目标检测MTD滤波器设计方法，以解决了上述方法中的对于动目标检测MTD滤波器中 的各个有限长脉冲响应FIR滤波器在其零频附近的凹口宽度都是固定的问题。该方法在数 字综合算法的动目标检测MTD滤波器的基础上，通过搜索最优的杂波功率谱方差使得输出 的杂波改善因子最大，从而得到最优的有限长脉冲响应FIR滤波器权矢量。
[0005]实现本发明的基本思路是，首先查寻雷达工作场景的杂波功率谱方差和杂波功 率，计算出地杂波的协方差矩阵。然后设计数字综合法的有限长脉冲响应FIR滤波器，在[0， 0。2]范围内搜索杂波功率谱方差，来改变数字综合法的有限长脉冲响应FIR滤波器的权矢 量，计算出对应的地杂波改善因子，从所有地杂波改善因子的值中搜索最大的地杂波改善 因子值，将最大地杂波改善因子值所对应的有限长脉冲响应FIR滤波器的权矢量作为期望 的有限长脉冲响应FIR滤波器的权矢量。重复执行上述步骤，直至得到满足设计者要求的所 有有限长脉冲响应FIR滤波器，最后将所有有限长脉冲响应FIR滤波器组成变凹口宽度的动 目标检测MTD滤波器。
[0006] 为实现上述目的，本发明包括如下步骤：
[0007] (1)查寻参数：
[0008] 通过查寻不同场景的地杂波功率谱方差和地杂波功率表，得到雷达工作场景地杂 波的功率谱方差和地杂波功率；
[0009] (2)获得地杂波的协方差矩阵：
[0010] (2a)利用高斯功率谱密度公式，计算雷达工作场景地杂波功率谱密度：
[0011] (2b)利用维纳辛钦公式，计算地杂波协方差矩阵的元素值；
[0012] (3)获得初始化干扰信号功率：
[0013] (3a)将L个干扰信号添加在有限长脉冲响应FIR滤波器副瓣区，覆盖在有限长脉冲 响应FIR滤波器频率响应的整个脉冲周期，每个干扰信号所对应的频率分别为心，f 2，…，fL， L的取值范围为有限长脉冲响应FIR滤波器的阶数的整数倍以上；
[0014] (3b)将有限长脉冲响应FIR滤波器副瓣区的L个干扰信号功率初始化为0;
[0015] (4)设计数字综合法的有限长脉冲响应FIR滤波器：
[0016] (4a)将当前搜索迭代次数k初始化为1;
[0017] (4b)按照下式，计算当前搜索迭代次数的干扰、杂波、加噪声的协方差矩阵：
[0019] 其中，Rb(k)表示当前搜索迭代次数的干扰、杂波、加噪声的协方差矩阵，k表示当 前搜索迭代次数，R。表示地杂波的协方差矩阵，心 2表示雷达内部接收单元的噪声功率，I表 示单位矩阵，Σ表示求和操作，L表示添加在有限长脉冲响应FIR滤波器副瓣区干扰信号的 个数，1表示第1个干扰，Si 2表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号功 率，fi表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率，a(fi)表 示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的导频矢量，Η表示共辄转置操 作；
[0020] (4c)利用自适应权公式，计算有限长脉冲响应FIR滤波器在当前搜索迭代次数时 的权矢量；
[0021] (4d)按照下式，计算有限长脉冲响应FIR滤波器在当前搜索迭代次数时的频率响 应：
[0022] H(k,f)=ffH(k)a(f)
[0023] 其中，H(k，f)表示在当前搜索迭代次数时有限长脉冲响应FIR滤波器的频率响应， k表示当前搜索迭代次数，f表示多普勒频率，WH(k)表示在当前搜索迭代次数时有限长脉冲 响应FIR滤波器权矢量的共辄转置，a(f)为多普勒频率f处的导频矢量；
[0024] (4e)按照下式，计算在当前搜索迭代次数时添加到有限长脉冲响应FIR滤波器副 瓣区干扰期望的电平值：
[0025] d(k,fi)=h(k)/10D(fl)/2°
[0026] 其中，d(k，f〇表示在当前搜索迭代次数时添加到有限长脉冲响应FIR滤波器副瓣 区第1个干扰期望的电平值，k表示当前搜索迭代次数，h表示添加到有限长脉冲响应FIR滤 波器副瓣区第1个干扰的多普勒频率，h(k)表示在当前搜索迭代次数时有限长脉冲响应FIR 滤波器的响应，D (fi)表示有限长脉冲响应FIR滤波器频率响应的副瓣电平在第1个干扰信 号处低于主瓣峰值的大小；
[0027] (4f)按照下式，计算下次搜索迭代时，有限长脉冲响应FIR滤波器副瓣区L个干扰 信号的功率：
[0029] 其中，δ^α+Ι)表示下次搜索迭代时，有限长脉冲响应FIR滤波器副瓣区第1个干扰 信号的功率，k表示当前搜索迭代次数，h表示添加在有限长脉冲响应FIR滤波器副瓣区第1 个干扰信号的多普勒频率，firft(k)，f right(k)分别表示在当前搜索迭代次数时有限长脉冲 响应FIR滤波器主瓣宽度的左右边界值，fzl(k)，f z2(k)分别表示在当前搜索迭代次数时有 限长脉冲响应FIR滤波器凹口宽度的左右边界值，U表示或操作，F r表示脉冲重复周期，max 表示求最大值操作，δ^α)表示在当前搜索迭代次数时，添加到有限长脉冲响应FIR滤波器 副瓣区第1个干扰信号的功率，K表示一常数标量，称为迭代增益，H(k，f〇表示在当前搜索 迭代次数时有限长脉冲响应FIR滤波器在频率fi时的频率响应，d(k，fi)表示在当前搜索迭 代次数时添加到有限长脉冲响应FIR滤波器副瓣区第1个干扰信号期望的电平值，h表示添 加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率；
[0030] (4g)判断W(k)-W(k-1) = 0是否成立，若是，则得到有限长脉冲响应FIR滤波器副瓣 区干扰信号功率后执行步骤(5)，否则，将当前搜索迭代次数k加1后执行步骤(4b);其中，W (k)表示在当前搜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，W(k-l)表示上一次搜 索迭代次数时有限长脉冲响应FIR滤波器的权矢量，k表示当前搜索迭代次数；
[0031] (5)搜索地杂波功率谱方差：
[0032] (5a)将当前的地杂波功率谱方差〇2初始化为0;
[0033] (5b)利用自相关公式，计算当前搜索的地杂波协方差矩阵元素值；
[0034] (5c)由当前搜索的地杂波协方差矩阵和步骤(4)的干扰信号功率，计算出当前搜 索的干扰、杂波、加噪声的协方差矩阵；
[0035] (5d)利用自适应权公式，计算当前搜索的有限长脉冲响应FIR滤波器权矢量；
[0036] (5e)将目标信号和地杂波输入到由新权矢量构成的有限长脉冲响应FIR滤波器 中，获得当前搜索的输出端的目标信号功率和地杂波功率；
[0037] (5f)利用杂波改善因子公式，计算当前搜索的地杂波改善因子；
[0038] (5g)判断|〇2_~2| <0.1是否成立，若是，则得到所有地杂波改善因子的值执行步 骤(6)，否则，将地杂波的功率谱方差〇2加0.1后执行步骤(5b);其中，| · |表示取绝对值操 作，σ2表示当前搜索的地杂波的功率谱方差，σ。2表示地杂波功率谱方差；
[0039] (6)获得最优杂波改善因子的有限长脉冲响应FIR滤波器：
[0040] (6a)从所有地杂波改善因子的值中搜索最大的地杂波改善因子值，将最大地杂波 改善因子值所对应的有限长脉冲响应FIR滤波器的权矢量作为期望的有限长脉冲响应FIR 滤波器的权矢量;其中，σ。2表示地杂波功率谱方差；
[0041] (6b)将期望的有限长脉冲响应FIR滤波器的权矢量，加入有限长脉冲响应FIR滤波 器中，获得最优杂波改善因子的一个有限长脉冲响应FIR滤波器；
[0042] (7)设计变凹口宽度的动目标检测MTD滤波器：
[0043] 重复执行上述6个步骤，直至得到满足设计者要求的所有有限长脉冲响应FIR滤波 器，将所有有限长脉冲响应FIR滤波器组成变凹口宽度的动目标检测MTD滤波器。
[0044] 本发明与现有技术相比具有以下优点：
[0045] 第一，本发明方法通过搜索地杂波功率谱方差，从所有地杂波改善因子的值中搜 索最大的地杂波改善因子值，将最大地杂波改善因子值所对应的有限长脉冲响应FIR滤波 器的权矢量作为期望的有限长脉冲响应FIR滤波器的权矢量，克服了现有技术中没有利用 雷达周围场景的杂波先验信息，其杂波的相关矩阵是通过查不同场景的地杂波功率谱方差 和地杂波功率表得到的，导致所有滤波器零陷的凹口宽度都是固定的，对于不同的雷达杂 波环境，滤波器的凹口宽度不能同杂波信号的谱宽相匹配的问题，使得本发明具有对于设 计的MTD滤波器中不同的通带中心频率有限长脉冲响应FIR滤波器有不同的期望凹口宽度， 能够在保证好的杂波抑制性能的同时，降低滤波器对信号的滤波损失的优点。
[0046] 第二，由于本发明使有限长脉冲响应FIR滤波器输出端杂波改善因子保持最优的 方法，克服了现有技术中的当滤波器通带中心频率较低时滤波器零频的凹口宽度会使滤波 器的主瓣发生畸变，滤波器的主瓣将会偏离中心频率，使得输出的目标信号有损失，而不能 有效的检测到目标的问题。使得本发明具有增强了对于低速运动目标检测的优点。
【附图说明】
[0047]图1为本发明的流程图；
[0048]图2为使用本发明对传统MTD滤波器的信杂比改善值的仿真效果图；
[0049] 图3为在与图2不同的杂波背景下，使用本明对传统MTD滤波器的信杂比改善值的 仿真效果图。
【具体实施方式】
[0050] 下面结合附图对本发明的步骤做进一步详细描述。
[0051]参照图1，本发明的具体实施步骤如下。
[0052] 步骤1，查寻参数。
[0053]通过查寻不同场景的地杂波功率谱方差和地杂波功率表，得到雷达工作场景地杂 波的功率谱方差和地杂波功率。
[0054]本发明的实例采用对空情报雷达，收发共用等距线阵。查寻得到的雷达工作场景 的地杂波功率谱方差用σ。2表示，查寻得到的雷达工作场景的地杂波功率用δ。2表示。
[0055]步骤2,获得地杂波的协方差矩阵。
[0056]按照下式高斯功率谱密度公式，计算雷达工作场景地杂波功率谱密度：
[0058] 其中，Sjf)表示雷达工作场景地杂波功率谱密度，c表示雷达工作场景地杂波，f 表示雷达工作场景地杂波功率谱密度的频率，·/:表示取根号操作，31表示圆周率，exe表示 以e为底数的指数操作，σ。 2表示雷达工作场景地杂波功率谱方差。
[0059] 按照下式维纳辛钦公式，计算地杂波协方差矩阵的元素值：
[0060] Rc (X, y) = IFT [ Sc (f) ] = 5〇2exp ( -2jt2〇c2 (y-χ) 2Tr2)
[0061] 其中，Κ(χ，5〇表示地杂波协方差矩阵第xRy列的元素值，c表示雷达工作场景地 杂波，X和y分别表示协方差矩阵的行和列，χ = 1，2···η，γ = 1，2···η，η表示有限长脉冲响应 FIR滤波器的阶数，IFT表示逆傅里叶变换操作，Sdf)表示雷达工作场景地杂波功率谱密 度，f表示雷达工作场景地杂波功率谱密度的频率，δ。 2表示雷达工作场景地杂波功率，exp表 示以e为底数的指数操作，π表示圆周率，〇。 2表示雷达工作场景地杂波功率谱方差，Tr表示脉 冲重复周期。
[0062]步骤3,获得初始化干扰信号功率。
[0063] 将L个干扰信号添加在有限长脉冲响应FIR滤波器副瓣区，覆盖在有限长脉冲响应 FIR滤波器频率响应的整个脉冲周期，每个干扰信号所对应的频率分别为fi，f2，···，&，L的 取值范围为有限长脉冲响应FIR滤波器的阶数的整数倍以上。
[0064] 将有限长脉冲响应FIR滤波器副瓣区的L个干扰信号功率初始化为0。
[0065]步骤4，设计数字综合法的有限长脉冲响应FIR滤波器。
[0066] 第一步，将当前搜索迭代次数k初始化为1。
[0067] 第二步，按照下式，计算当前搜索迭代次数的干扰、杂波、加噪声的协方差矩阵：
[0069]其中，Rb(k)表示当前搜索迭代次数的干扰、杂波、加噪声的协方差矩阵，k表示当 前搜索迭代次数，R。表示地杂波的协方差矩阵，心2表示雷达内部接收单元的噪声功率，I表 示单位矩阵，Σ表示求和操作，L表示添加在有限长脉冲响应FIR滤波器副瓣区干扰信号的 个数，1表示第1个干扰，Si 2表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号功 率，fi表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率，a(fi)表 示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的导频矢量，Η表示共辄转置操 作。
[0070] 第三步，按照下式自适应权公式，计算有限长脉冲响应FIR滤波器在当前搜索迭代 次数时的权矢量：
[0071] ff(k)=Rb_1(k)a(fi)
[0072] 其中，W(k)表示在当前搜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，k表 示当前搜索迭代次数，R^a)表示在当前搜索迭代次数时干扰、杂波、加噪声的协方差矩阵 的逆，aaO表示通带中心频率为h的有限长脉冲响应FIR滤波器的导频矢量。
[0073] 第四步，按照下式，计算有限长脉冲响应FIR滤波器在当前搜索迭代次数时的频率 响应：
[0074] H(k,f)=ffH(k)a(f)
[0075] 其中，H(k，f)表示在当前搜索迭代次数时有限长脉冲响应FIR滤波器的频率响应， k表示当前搜索迭代次数，f表示多普勒频率，WH(k)表示在当前搜索迭代次数时有限长脉冲 响应FIR滤波器权矢量的共辄转置，a( f)为多普勒频率f处的导频矢量。
[0076] 第五步，按照下式，计算在当前搜索迭代次数时添加到有限长脉冲响应FIR滤波器 副瓣区干扰期望的电平值：
[0077] d(kJ^ = h[k)![Qnu：]：ii
[0078] 其中，d(k，f〇表示在当前搜索迭代次数时添加到有限长脉冲响应FIR滤波器副瓣 区第1个干扰期望的电平值，k表示当前搜索迭代次数，h表示添加到有限长脉冲响应FIR滤 波器副瓣区第1个干扰的多普勒频率，h(k)表示在当前搜索迭代次数时有限长脉冲响应FIR 滤波器的响应，D (fi)表示有限长脉冲响应FIR滤波器频率响应的副瓣电平在第1个干扰信 号处低于主瓣峰值的大小。
[0079] 第六步，按照下式，计算下次搜索迭代时，有限长脉冲响应FIR滤波器副瓣区L个干 扰信号的功率：
[0081] 其中，δ^α+Ι)表示下次搜索迭代时，有限长脉冲响应FIR滤波器副瓣区第1个干扰 信号的功率，k表示当前搜索迭代次数，h表示添加在有限长脉冲响应FIR滤波器副瓣区第1 个干扰信号的多普勒频率，firft(k)，f right(k)分别表示在当前搜索迭代次数时有限长脉冲 响应FIR滤波器主瓣宽度的左右边界值，fzl(k)，f z2(k)分别表示在当前搜索迭代次数时有 限长脉冲响应FIR滤波器凹口宽度的左右边界值，U表示或操作，F r表示脉冲重复周期，max 表示求最大值操作，δ^α)表示在当前搜索迭代次数时，添加到有限长脉冲响应FIR滤波器 副瓣区第1个干扰信号的功率，K表示一常数标量，称为迭代增益，H(k，f〇表示在当前搜索 迭代次数时有限长脉冲响应FIR滤波器在频率fi时的频率响应，d(k，fi)表示在当前搜索迭 代次数时添加到有限长脉冲响应FIR滤波器副瓣区第1个干扰信号期望的电平值，h表示添 加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率。
[0082] 第七步，判断W(k)-W(k-1)=0是否成立，若是，则得到有限长脉冲响应FIR滤波器 副瓣区干扰信号功率后执行步骤5,否则，将当前搜索迭代次数k加1后执行本步骤的第二 步;其中，W(k)表示在当前搜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，W(k-l)表 示上一次搜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，k表示当前搜索迭代次数。 [0083]步骤5,搜索地杂波功率谱方差。
[0084] 第一步，将当前的地杂波功率谱方差σ2初始化为0。
[0085] 第二步，按照下式自相关公式，计算当前搜索的地杂波协方差矩阵元素值：
[0086] R(x，y) =5c2exp(-2JT2o2(y-x) 2Tr2)
[0087] 其中，R(x，y)表示当前搜素的地杂波协方差矩阵第xRy列的元素值，c表示雷达工 作场景地杂波，X和y分别表示协方差对角矩阵的行和列，χ=1，2···η，γ = 1，2···η，η表示有限 长脉冲响应FIR滤波器的阶数，δ2。表示地杂波功率，exp表示以 e为底数的指数操作，π表示圆 周率，〇2表示当前搜索的地杂波功率谱方差，T r表示脉冲重复周期。
[0088] 第三步，由当前搜索的地杂波协方差矩阵和步骤4的干扰信号功率，计算出当前搜 索的干扰、杂波、加噪声的协方差矩阵。
[0089] 第四步，按照下式自适应权公式，计算当前搜索的有限长脉冲响应FIR滤波器权矢 量。
[0090] ff=R-Wfi)
[0091] 其中，W表示当前搜索的有限长脉冲响应FIR滤波器的权矢量，IT1表示当前搜索的 干扰、杂波、加噪声的协方差矩阵的逆，aaO表示通带中心频率为h的有限长脉冲响应FIR 滤波器的导频矢量。
[0092] 第五步，将目标信号和地杂波输入到由新权矢量构成的有限长脉冲响应FIR滤波 器中，获得当前搜索的输出端的目标信号功率和地杂波功率。
[0093] 其中目标信号采用的为仿真数据，地杂波采用的为雷达工作场景的地杂波数据。
[0094] 第六步，按照下式杂波改善因子公式，计算当前搜索的地杂波改善因子：
[0096]其中，I表示当前搜素的地杂波改善因子，So表示当前搜素的有限长脉冲响应FIR 滤波器输出端的目标信号功率，Co表示当前搜素的有限长脉冲响应FIR滤波器输出端的地 杂波功率，Si表示有限长脉冲响应FIR滤波器当前搜素的输入端的目标信号功率，(^表示当 前搜素的有限长脉冲响应FIR滤波器输入端的地杂波功率。
[0097]第七步，判断|〇2-〇。2| <0.1是否成立，若是，则得到所有地杂波改善因子的值执行 步骤6,否则，则将地杂波的功率谱方差σ2加0.1后执行本步骤的第二步;其中，| · |表示取 绝对值操作，σ2表示当前搜索的地杂波的功率谱方差，σ。2表示地杂波功率谱方差。
[0098]步骤6,获得最优杂波改善因子的有限长脉冲响应FIR滤波器。
[0099] 从所有地杂波改善因子的值中搜索最大的地杂波改善因子值，将最大地杂波改善 因子值所对应的有限长脉冲响应FIR滤波器的权矢量作为期望的有限长脉冲响应FIR滤波 器的权矢量。
[0100] 将期望的有限长脉冲响应FIR滤波器的权矢量，加入有限长脉冲响应FIR滤波器 中，获得最优杂波改善因子的一个有限长脉冲响应FIR滤波器。
[0101] 步骤7，设计变凹口宽度的动目标检测MTD滤波器。
[0102] 重复执行上述6个步骤，直至得到满足设计者要求的所有有限长脉冲响应FIR滤波 器，将所有有限长脉冲响应FIR滤波器组成变凹口宽度的动目标检测MTD滤波器。
[0103] 下面结合附图对本发明的效果做进一步描述。
[0104] 1.仿真条件：
[0105] 本发明仿真实验的软件运行系统为64位windows操作系统，仿真软件为MATLAB (R2007b)。
[0106] 仿真本发明的MTD滤波器对传统MTD滤波器的信杂比改善值所采用的仿真参数为： 有限长脉冲响应FTR滤波器的阶数为32,脉冲重复周期PRT = 2370us，接收单元的噪声功率 为1，滤波器的迭代增益K=l，目标信号为根据雷达参数模拟产生的多普勒频率由35Hz变化 到220HZ的仿真数据，杂波数据采用某型号对空情报雷达采集的某段杂波信号作为杂波背 景，得到仿真图2。
[0107] 在与图2不同的杂波背景下，仿真本发明的MTD滤波器对传统MTD滤波器的信杂比 改善值所采用的仿真参数为:有限长脉冲响应FTR滤波器的阶数为32,脉冲重复周期PRT = 2370us，接收单元的噪声功率为1，滤波器的迭代增益K= 1，目标信号为根据雷达参数模拟 产生的多普勒频率由35Hz变化到220Hz的仿真数据，杂波数据采用与图2相同的某型号对空 情报雷达采集的另一段杂波信号作为杂波背景，得到仿真图3。
[0108] 2.仿真内容及结果分析：
[0109] 图2为本发明方法对传统MTD滤波器的信杂比改善值的仿真效果图。其中横坐标表 示多普勒频率，纵坐标表示本发明方法对传统MTD滤波器的信杂比改善值，该条信杂比改善 值曲线是通过让杂波信号和目标信号分别通过基于先验信息的变凹口宽度的动目标检测 MTD滤波器和传统的数字综合法MTD滤波器，比较两个滤波器输出端的信杂比改善值得到 的，由图2可以得到信杂比改善平均值在2.8dB以上，所以本发明方法具有较传统MTD滤波器 设计方法有更好的信杂比改善性能。
[0110] 图3为在与图2不同的杂波背景下，本发明方法对传统MTD滤波器的信杂比改善值 的仿真效果图。其中横坐标表示多普勒频率，纵坐标表示本发明方法对传统MTD滤波器的信 杂比改善值，该条信杂比改善值曲线是通过让杂波信号和目标信号分别通过基于先验信息 的变凹口宽度的动目标检测MTD滤波器和传统的数字综合法MTD滤波器，比较两个滤波器输 出端的信杂比改善值得到的，由图2可以得到信杂比改善平均值在3.8dB以上，所以本发明 方法具有较传统MTD滤波器设计方法有更好的信杂比改善性能。
【主权项】
1. 一种基于先验信息的变凹口宽度MTD滤波器设计方法，包括如下步骤： (1) 查寻参数： 通过查寻不同场景的地杂波功率谱方差和地杂波功率表，得到雷达工作场景地杂波的 功率谱方差和地杂波功率； (2) 获得地杂波的协方差矩阵： (2a)利用高斯功率谱密度公式，计算雷达工作场景地杂波功率谱密度； (2b)利用维纳辛钦公式，计算地杂波协方差矩阵的元素值； (3) 获得初始化干扰信号功率： (3a)将L个干扰信号添加在有限长脉冲响应FIR滤波器副瓣区，覆盖在有限长脉冲响应 FIR滤波器频率响应的整个脉冲周期，每个干扰信号所对应的频率分别为fi，f2，-'，fL，L的 取值范围为有限长脉冲响应FIR滤波器的阶数的整数倍，W上； (3b)将有限长脉冲响应FIR滤波器副瓣区的L个干扰信号功率初始化为0; (4) 设计数字综合法的有限长脉冲响应FIR滤波器： (4a)将当前捜索迭代次数k初始化为1; (4b)按照下式，计算当前捜索迭代次数的干扰、杂波、加噪声的协方差矩阵：其中，化化)表示当前捜索迭代次数的干扰、杂波、加噪声的协方差矩阵，k表示当前捜索 迭代次数，R。表示地杂波的协方差矩阵，Sn2表示雷达内部接收单元的噪声功率，I表示单位 矩阵，Σ表示求和操作，L表示添加在有限长脉冲响应FIR滤波器副瓣区干扰信号的个数，1 表示第1个干扰，Si2表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号功率，fi表 示添加在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率，a(fi)表示添加 在有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的导频矢量，Η表示共辆转置操作； (4c)利用自适应权公式，计算有限长脉冲响应FIR滤波器在当前捜索迭代次数时的权 矢量； (4d)按照下式，计算有限长脉冲响应FIR滤波器在当前捜索迭代次数时的频率响应： Η 化，f)=wH 化)a(f) 其中，Η化，f)表示在当前捜索迭代次数时有限长脉冲响应FIR滤波器的频率响应，k表 示当前捜索迭代次数，f表示多普勒频率，wH(k)表示在当前捜索迭代次数时有限长脉冲响 应FIR滤波器权矢量的共辆转置，a(f)为多普勒频率f处的导频矢量； (4e)按照下式，计算在当前捜索迭代次数时添加到有限长脉冲响应FIR滤波器副瓣区 干扰期望的电平值：其中，cKk，fi)表示在当前捜索迭代次数时添加到有限长脉冲响应FIR滤波器副瓣区第1 个干扰期望的电平值，k表示当前捜索迭代次数，fi表示添加到有限长脉冲响应FIR滤波器 副瓣区第1个干扰的多普勒频率，h化)表示在当前捜索迭代次数时有限长脉冲响应FIR滤波 器的响应，D(fi)表示有限长脉冲响应FIR滤波器频率响应的副瓣电平在第1个干扰信号处 低于主瓣峰值的大小； (4f)按照下式，计算下次捜索迭代时，有限长脉冲响应FIR滤波器副瓣区L个干扰信号 的功率：其中，δι2化+1)表示下次捜索迭代时，有限长脉冲响应FIR滤波器副瓣区第1个干扰信号 的功率，k表示当前捜索迭代次数，fi表示添加在有限长脉冲响应FIR滤波器副瓣区第1个干 扰信号的多普勒频率，fleft化），fright化)分别表示在当前捜索迭代次数时有限长脉冲响应 FIR滤波器主瓣宽度的左右边界值，fzl(k)，fz2化)分别表示在当前捜索迭代次数时有限长 脉冲响应FIR滤波器凹口宽度的左右边界值，U表示或操作，Fr表示脉冲重复周期，max表示 求最大值操作，Si2化)表示在当前捜索迭代次数时，添加到有限长脉冲响应FIR滤波器副瓣 区第1个干扰信号的功率，K表示一常数标量，称为迭代增益，Η化，fi)表示在当前捜索迭代 次数时有限长脉冲响应FIR滤波器在频率fi时的频率响应，cKk，fi)表示在当前捜索迭代次 数时添加到有限长脉冲响应FIR滤波器副瓣区第1个干扰信号期望的电平值，fi表示添加在 有限长脉冲响应FIR滤波器副瓣区第1个干扰信号的多普勒频率； (4g)判断W化)-W化-1)=0是否成立，若是，则得到有限长脉冲响应FIR滤波器副瓣区干 扰信号功率后执行步骤(5)，否则，将当前捜索迭代次数k加1后执行步骤(4b);其中，W化)表 示在当前捜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，W化-1)表示上一次捜索迭 代次数时有限长脉冲响应FIR滤波器的权矢量，k表示当前捜索迭代次数； (5) 捜索地杂波功率谱方差： 巧a)将当前的地杂波功率谱方差σ2初始化为0; 巧b)利用自相关公式，计算当前捜索的地杂波协方差矩阵元素值； (5c)由当前捜索的地杂波协方差矩阵和步骤(4)的干扰信号功率，计算出当前捜索的 干扰、杂波、加噪声的协方差矩阵； 巧d)利用自适应权公式，计算当前捜索的有限长脉冲响应FIR滤波器权矢量； (5e)将目标信号和地杂波输入到由新权矢量构成的有限长脉冲响应FIR滤波器中，获 得当前捜索的输出端的目标信号功率和地杂波功率； 巧f)利用杂波改善因子公式，计算当前捜索的地杂波改善因子； (5g)判断I 1 <0.1是否成立，若是，则得到所有地杂波改善因子的值执行步骤(6)， 否则，则将地杂波的功率谱方差σ2加0.1后执行步骤(5b);其中，I · I表示取绝对值操作，曰2 表示当前捜索的地杂波的功率谱方差，σ。2表示地杂波功率谱方差； (6) 获得最优杂波改善因子的有限长脉冲响应FIR滤波器： (6a)从所有地杂波改善因子的值中捜索最大的地杂波改善因子值，将最大地杂波改善 因子值所对应的有限长脉冲响应FIR滤波器的权矢量作为期望的有限长脉冲响应FIR滤波 器的权矢量； (6b)将期望的有限长脉冲响应FIR滤波器的权矢量，加入有限长脉冲响应FIR滤波器 中，获得最优杂波改善因子的一个有限长脉冲响应FIR滤波器； (7)设计变凹口宽度的动目标检测MTD滤波器： 重复执行上述6个步骤，直至得到满足设计者要求的所有有限长脉冲响应FIR滤波器， 将所有有限长脉冲响应FIR滤波器组成变凹口宽度的动目标检测MTD滤波器。2. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(2a)所述的高斯功率谱密度公式如下：其中，ScXf)表示雷达工作场景地杂波功率谱密度，C表示雷达工作场景地杂波，f表示雷 达工作场景地杂波功率谱密度的频率，λΓ表示取根号操作，η表示圆周率，exe表示We为底 数的指数操作，σ。2表示雷达工作场景地杂波功率谱方差。3. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(2b)所述的维纳辛钦公式如下：其中，Rt(x，y)表示地杂波协方差矩阵第X行y列的元素值，C表示雷达工作场景地杂波，X 和y分别表示协方差矩阵的行和列，χ=1，2'''n，y = l，2···η，η表示有限长脉冲响应FIR滤波 器的阶数，IFT表示逆傅里叶变换操作，Sc(f)表示雷达工作场景地杂波功率谱密度，f表示 雷达工作场景地杂波功率谱密度的频率，δ。2表示雷达工作场景地杂波功率，exp表示We为 底数的指数操作，η表示圆周率，σ。2表示雷达工作场景地杂波功率谱方差，Tr表示脉冲重复 周期。4. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(4c)所述的自适应权公式如下：其中，W化)表示在当前捜索迭代次数时有限长脉冲响应FIR滤波器的权矢量，k表示当 前捜索迭代次数，Rb^i(k)表示在当前捜索迭代次数时干扰、杂波、加噪声的协方差矩阵的 逆，a(fi)表示通带中屯、频率为fi的有限长脉冲响应FIR滤波器的导频矢量。5. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(5b)所述的自相关公式如下：其中，R(x，y)表示当前捜素的地杂波协方差矩阵第X行y列的元素值，C表示雷达工作场 景地杂波，X和y分别表示协方差对角矩阵的行和列，又=1，2-'11，7=1，2-'11，11表示有限长脉 冲响应FIR滤波器的阶数，δ2。表示地杂波功率，exp表示We为底数的指数操作，31表示圆周 率，表示当前捜索的地杂波功率谱方差，Tr表示脉冲重复周期。6. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(5d)所述的自适应权公式如下：其中，W表示当前捜索的有限长脉冲响应FIR滤波器的权矢量，扩1表示当前捜索的干扰、 杂波、加噪声的协方差矩阵的逆，a(fi)表示通带中屯、频率为fi的有限长脉冲响应FIR滤波器 的导频矢量。7. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(5e)所述的目标信号采用的为仿真数据，地杂波采用的为雷达工作场景的地杂波 数据。8. 根据权利要求1所述的基于先验信息的变凹口宽度MTD滤波器设计方法，其特征在 于，步骤(5f)所述的杂波改善因子公式如下：其中，I表示当前捜素的地杂波改善因子，So表示当前捜素的有限长脉冲响应FIR滤波器 输出端的目标信号功率，Co表示当前捜素的有限长脉冲响应FIR滤波器输出端的地杂波功 率，Si表示有限长脉冲响应FIR滤波器当前捜素的输入端的目标信号功率，Cl表示当前捜素 的有限长脉冲响应FIR滤波器输入端的地杂波功率。
【文档编号】G01S7/36GK106093888SQ201610511506
【公开日】2016年11月9日
【申请日】2016年7月1日
【发明人】赵永波, 孙林, 程增飞, 刘宏伟, 何学辉, 苏洪涛
【申请人】西安电子科技大学