动目标显示与动目标检测

《动目标显示与动目标检测》由会员分享，可在线阅读，更多相关《动目标显示与动目标检测（18页珍藏版）》请在装配图网上搜索。

1、引言1.目标回波频谱6.5.1目标回波和杂波的频谱j 2.杂波频谱1原理递归传统非递归零点分配算法滤波器设计-优化预测误差算法 结语原理MTI+FFT6.5.3 MTD滤波器滤波器设计Y点最佳、等间隔最佳结语6.5.4改善因子分析J MTII MTD6.5动目标显示与动目标检测雷达探测的运动目标如飞机，导弹，舰艇，车辆等周围存在各种背景，包括不动的地物和运动着的云雨，海浪或金属丝干扰等。动目标 显示(Moving Target Indicator :MTD与动目标检测(Moving Target Detection: MTD)就是使用各种滤波器，滤去这些背景产生的杂波而取出 运动目标的回波。此

2、外也可以通过把雷达安装在山上、增加雷达天线的倾角、安装防杂波网来阻止杂波进入天线；或通过调整雷达天线的波 束形式、采用极化技术、降低雷达的分辨单元、在时域采用CFAR检测、自适应门限、杂波图来抑制杂波。在频域上应用MTI与MTD技术 可以提高信杂比，改善杂波背景下检测运动目标的能力。本节首先分析目标回波和杂波的频谱特性；然后分别讨论MTI与MTD原理及滤波 器设计方法；最后分析MTI与MTD对改善因子的提高。6.5.1目标回波和杂波的频谱运动目标回波和杂波在频谱结构上有所差别，运动目标检测就是利用这种差别，从频率上将它们区分，以达到抑制杂波而显示目标回波 的目的。为此，应首先弄清楚目标和杂波的

3、回波的特性。(1)目标回波的频谱雷达发射相参脉冲串，其脉冲宽度为T，脉冲重复频率为fr。当天线不扫描而对准目标时所得脉冲为无限脉冲串。调制信号)及其频谱Ui(f)分别为u (t) = A 党 rect(t nT )件Tn=-3V 七 J(6.5.1)AT Sin( fT )尹 eU 1( f )=兀 fT 5 ( f FDre n=3(6.5.2)A为信号振幅。而高频载波(t)及其频谱(f)为u (t) = cosW tU (f)=邹(f - f ) + 6 (f + f )2200(6.5.3)(6.5.4)发射的相参脉冲串u(t) - ui(t) u2(t)故其频谱U (f)为(6.5.5

4、)脉冲发射A发射信号频谱(a)回波K固定)f(b)回波(运动)U ( f) - U i( f) U 2( f)-(f七(c)f0+fdW天线扫描动目标回波频谱U现f)1(a)T固定回波频谱祖 f0 fo+1 f气-M (b)f0-1 f0动目标回波频谱(c)3一飞飞飞_._：.：!t一_(d)-(f+1(d)1 4+fdf-(f0 + fd)+,(d)f0+ fd-T图6.5.1发射接收信号时域与频域特性雷达发射信号通常是窄带信号(如图6.5.1(a)，因而运动目标回波频谱的特征是将发射信号的频谱位置在频率轴上平移一个多普勒频率 厂2v。fd =了(如图6.5.1(c)，fd的符号由目标运动的

5、方向决定，靠近为正，远离为负。固定回波的频谱与发射信号一致，只是幅度有衰减(如多普勒频率f可以直观地解释为：振荡源发射的电磁波以恒速c传播，当接收者相对振荡源不动，则它的接收频率等于发射频率。当 d接收者与振荡源之间有相对接近的运动时，则它接收电磁波的频率大于振荡源发射频率，当两者背向运动时，结果相反。运动目标与雷达有”-2 R相对运动，所以存在多普勒频率fd。定义固定回波的波程为2R0，R0为雷达到目标的距离，则回波到达雷达的时间为宣；由于运动目标c/ A 、2(R -AR)与雷达有相对运动，其波程为2(R0-AR)，时间为一。雷达工作时，天线以各种方式进行扫描，这时收到的回波脉冲为有限数，

6、且其振幅受天线方向图调制。设天线方向图可用高斯函数来 表示，则收到的回波脉冲串的包络函数可写为(6.5.6)m (t) =、. 2兀 b exp ( 2 兀 2 b 212)b是和天线波瓣宽度及扫描速度有关的参数。减小，表示观察时间增加。天线扫描时收到的回波信号，可以用m(t)和无限脉冲串七(t)的乘积表示。七(t)为天线不扫描时的回波脉冲串，即(6.5.7)u (t) = m(t)u (t)其包络函数m(t)的频谱为图6.5.2天线扫描条件下回波频谱(6.5.8)天线扫描时回波信号的频谱U (f)为 mUm (f) = M (f) Ur (f)(6.5.9)即无限回波脉冲串频谱气(f)的每一

7、根谱线均按M(f)的形状展宽(如图6.5.1(d)。谱线展宽的程度反比于天线波束照射目标的时间。已 求出当天线方向图为高斯形时谱线展宽的均方值为0.265 0.265f(6.5.10)=rt9式中：fr为雷达重复频率，n为在单程天线方向图3dB宽度内收到的脉冲数。中频回波信号经过相位检波器后，相当于把中频信号的频谱搬移到零频率附近，根据目标多普勒频移fd的不同，相位检波后谱线nfr 土 fd的具体位置也有差异，每根谱线均按脉冲串包络的频谱形状展宽。杂波频谱对于固定点杂波，当天线不扫描时，固定杂波的频谱是位于nf上的谱线，可以用对消器全部滤去。当天线扫描时，由于回波数目有r限，谱线将展宽。由于天

8、线扫描引起双程天线方向图对回波信号调幅，杂波谱展宽可用高斯函数表示为(6.5.11)G( f) = G e -f 2 2b20其中b广0.265/n，n为在单程天线方向图3dB宽度内的脉冲数。设为天线照射目标的等效时间，则n = Tq f，即b广0.265,气，即b，与目标照射时间成正比。杂波信号的功率谱的实验公式可近似为(6.5.12)W( / ) = |g (/ )|2 = |g J2expi |_ V f 0 J _其中：W（ f ）作为频率函数的杂波功率谱，g（ f ）杂波的傅立叶变换，f雷达载波，a和杂波相关的参数。杂波频谱可以用杂波频率分布的均方根值。（Hz）或速度分布的均方根值。

9、（ms）来表示，式（2-12）可写为CV2b其中：Wo = go2，W，rW (f) = WexpV2b 2 JCr=W exp 一0 V(6.5.13)C ,C 2f 可得a =航J0Vb 2为杂波功率谱方差。Cb V为杂波内部起伏运动速度的均方根值，和工作波长无关。相同的b,值，对不同的雷达工作波长产生的杂波谱线的宽度也是不同的。工作波长越短，杂波谱的展宽越严重。图5.4.1中的载波频率fo 为 1GHz。(pI10ff(c)%fff(b)0ff f滤波器频率特性0fdf图6.5.3杂波的功率谱在接收机前端引入发射信号作为基准电压，可得到收发频率的差频电压，即多普勒电压。图6.5.4给出各

10、主要点的频谱图。(%(b)f0f + fdf(c)0fdf0ff0 - fd0f0+ f ff0(d)图6.5.4主要点频谱图及滤波器特性f0f点，图中给出的是一条谱线。实际中，由 度的展宽，由于系统不稳等原因，杂波谱 设计的滤波器特性如图6.5.5地杂波的功率谱滤波器要满足的要求：（1）凹口适当扩宽，相位检波后的频谱如图（d），固定杂波的频谱在零 于天线扫描收到的回波脉冲数有限，谱线会有一定程 中还会出现杂乱分量，把它近似看成均匀谱。本节要图6.5.4（d）所示，取出动目标频谱，滤除杂波频谱。 滤波器凹口和通带的平坦程度是关注的特性。动目标 与杂波梳状谱宽度相当。（2）杂波有多普勒频移，即不

11、在零频时，滤波器凹口要对准杂波谱平均多普勒位置。（3）凹口深度能使杂波尽量多的 衰减，目标回波能在尽可能大的速度范围内有较大输出。改善因子是综合评定滤波器性能的参数。6.5.2 MTI滤波器（1）MTI滤波器原理当杂波和运动目标回波在雷达显示器上同时显示时，会使目标的观察变得很困难。如果目标处在强杂波背景内，弱的目标淹没在强杂 波中，特别是当强杂波使接收机发生过载时，将很难发现目标。目标回波和杂波在时间域上难以区分，但由于目标的速度远大于背景的速度， 目标回波的多普勒频移远大于背景的多普勒频移，从而可在频域上区分目标与杂波。动目标显示滤波器MTI）利用运动目标回波和杂波在 频谱上的区别，有效地

12、抑制杂波而提取信号。在雷达上加装MTI滤波器，大大的改善了雷达在强杂波背景中检测运动目标的能力MTI有 多种实现方法，包括传统的相消器和各种优化的FIR滤波器。采用重复参差和时变加权的MTI体制可以克服盲速。MTI和MTD是对多个回波数据进行处理，因此数据需要进行存储。其中每一行的数据是沿距离单元采样值，反映了某一距离单元的信 号特征。每一列的数据为从相同的距离单元，依次间隔一个脉冲重复周期的采样值，这些数据的变化反映了在同一距离单元目标的变化情况。 MTI和MTD都是对同一距离单元的数据，即同一列的数据进行处理。通常MTI雷达滤波器如图6.5.6所示图，图中Tr为雷达重复周期，在 这里作为延

13、迟线的延迟时间，巧为滤波器权系数值。i图6.5.6 MTI滤波器的组成传统的MTI相消器可以滤除零频杂波，性能不高，改善因子在20dB左右。优化的MTI滤波器可以满足上一节提出的动目标滤波器要 求，现在就滤波器的设计讨论如下。(2) MTI滤波器设计1.传统MTI滤波器设计：在相位检波器输出端，固定目标的回波是一串振幅不变的脉冲，而运动目标的回波是一串振幅调制的脉冲。在把回波信号送到终端显示 器前，必须先消除固定目标回波。最直接的方法是将相邻重复周期的回波信号相减，则固定目标回波由于振幅不变而互相抵消，运动目标回 波相减后剩下相邻重复周期振幅变化的部分。传统的MTI滤波器有两种形式：非递归形和

14、递归形。(a)非递归滤波器不带反馈的滤波器称为非递归型滤波器。下面以一次对消器为例进行说明。一次对消器，即二脉冲对消。其结构图如图6.5.7(a)，对消器的输入X(z)相位检波器的输出信号。它是一个单零点系统，零点位置在Z = 1， 令s =加，即z = ej回在Z平面上是单位圆。(a)框图(b)零点图图6.5.7 一次对消滤波器由相位检波器输出的脉冲包络为u(t) = U0 cos中(t)(6.5.14)中为回波与基准电压之间的相位差中(t) = -o t =-o 2(0_tl) =w t 中(6.5.15)0 r 0 cd 0回波信号按重复周期Tr出现，将回波信号延迟一周期后，其包络是u(

15、t) = U0 coswd (t T )甲0 (6.5.16)u =u - u = 2U sinT )fT)I 2 Jsin气 一平20 J(6.5.17)输出包络为一多普勒频率的正弦信号，其频率为(6.5.18)为多普勒频率的函数。当JJ2 = n兀(n=1,2,3,.)时，输出振幅为零。这时的目标速度正相当于盲速，盲速是运动目标回波在相位检波 器的输出端与固定目标回波相同，因而经对消设备后输出为零。下面从频率域来说明对消器的工作原理。对消器的输出为u0 =气(1 e-件)(6.5.19)对消器的频率响应特性为. uJH (爬)= 1 e -沁 Tr = (1 cos T) + j sin

16、T = 2sin 兀 fTej(2仃)(6.5.20)i图6.5.8对消器框图及其输出响应其频率响应特性如图6.5.8。对消器等效于一个梳状滤波器，其频率特性在f = f 各点处均为零。固定目标频谱的特点是它的谱线位于nfr点上，因而在理想情况下，通过对消器后输出为零。当目标的多普勒频率为重复频率的整数倍时，其频谱结构也有相同的特点，故通过上述梳状滤波器后无输出。表6.5.1二项式权系数和增益表对消器阶数增益G = 1 w.2 (线性)j=0二项式系数121, -1261, -2, 13201, 3, 3,-14701, -4, 6, -4, 152521, -5,10, -10,5, -1采

17、用级联二项式形的N次对消器(N +1脉冲对消)的权系数为G X展开式的系数，其计算式为(6.5.21)* =(T)-1 N iN G 1)!，i = 1,2,.“N +1非递归滤波器中的权系数变化它的输出也会变化。MTI滤波器中的权系数可选二项式系数，也可选其它的系数。在不同的准则下，系 数也不同。在没有先验知识的情况下，二项式系数式一个很好的选择。如果已知信号和杂波的协方差阵，在最优MTI滤波器的准则是改善因子最大化时，则可以求其最优化的权系数。如果杂波协方差阵为：A = E(c *cT)，信号协方差阵为：S = E(s * st )。其中：权系数向量为w = (N x1)，信号向量为s =

18、 (N x1)，杂波向量为c = (N x 1)，*代表共轭，T代表转置。则FIR滤波器的最优化权系数为w = kA-1 s *(6.5.22)i式中k为任意常数，A -1为杂波协方差阵的逆，为目标回波信号向量。准则不同结果也不同。如果没有先验知识或先验知识不准确，则求得的最优权系数的实现效果可。递归型一次对消器的传递函数为最优化权系数是相对于一定准则的， 能会小于二项式权系数实现的结果。(b)递归滤波器带有反馈的滤波器称为递归型滤波器(6.5.23)Z )_ 1 - Z -1H (z)X( z)1 - Kz -11递归滤波器的框图和频响如图6.4.4所示。虽然递归滤波器有较好的频响，可以减小

19、滤波器的阶数，但它的暂态响应差，限制了它在雷达信号处理中的应用。2.5I.D10.5频率0.D1(b)递归滤波器频率响应图6.5.9波器框图及其输出响应2.优化MTI滤波器设计：在自适应杂波抑制的应用环境下，我们需要的理想滤波器是在杂波分布的频率点处，能使杂波得到最大限度的抑制，而在其他频率点具 有最大平坦幅度。(a)零点分配算法零点分配算法是在设计带阻滤波器时，在凹口处设置频率响应零点的一种方法，他不需要计算杂波的相关系数。的延迟时间是匚FIR滤波器在t时刻的输出对于FIR滤波器，如图6.5.10，设第i阶 是前几个时刻输入信号的加权和。N阶滤波器 的传输函数为：2，N，滤波器H(f)=W

20、w e - j 2fii=0(6.5.24)式中，T0=0且T，=L k ； i=0,1,2,N。 H (f )在f=f0处展开成泰勒级数: k=1H (f )=H 侦)+H (f0)(f-f+H(”/2顷-状+L (皿25)式中:H (k)(f )= ( j2 兀 刃 Tkwe - j 2fT0z(6.5.26)i=0要在f=f0处设计带阻滤波器，为使在孔处附近幅度最大平坦，则必须使泰勒级数展开式中(f - f) k(k=0,1,2,N-1)的系数为零，即HMf ) = 0，k=0，1，2，N-1(6.5.27)J 0这样，就产生了 N个关于绪，的齐次线形方程。刃 Tkwe-j2f. = 0

21、k=0,1,2,N-1(6.5.28)i=0其中，To = 1， w0为一个常数。将式(6.5.28)写成矩阵形式:A W =- wU(6.5.29)其中：W = (* %L , WN ，U =(L0, L ,0 )te-j2fTe-j2兀 fT1Me-2fLe-fT e-fLT e-f2NMMTN*-j2兀暮TN-1e-j2兀值 L12-150-1012归一化频率f/fr3-150-10123归一化频率f/fr取t . =Tr解矩阵方程(6.5.29)得权系数向量W带入式(6.5.24)即得滤波器的频率响应如下:(6.5.30)(a)凹口对准0频点杂波(b)凹口对准动杂波图6.5.11零点分

22、配算法MTI滤波器频率响应目标速度为对应f的整数时倍称为盲速。山图6.5.11，MTI滤波器会在整数倍fr的多普勒频率处形成凹口。因此，当运动目标的多普勒频率等于整数倍f时，这些运动目标也会被MTI滤波器滤掉。使用参差变T可以将第一个盲速扩大到可以接受的范围，即使其大于 r可能出现的目标最大速度。如图G50，设计一个3阶3参差变T滤波器，L： T 2： T 3=11：12：13，得MTI滤波器频率响应如下:o O 心 。心 5 0 5 0 -112- - -:应响率频.)0应响率频 _ rrrrr _ _ _ rrrrr-250-25002468100246810归一化频率f/fr归一化频率f

23、/fr(a)凹口对准0频点杂波(b)凹口对准动杂波图6.5.11参差变T零点分配算法MTI滤波器频率响应(b)预测误差算法利用预测误差算法设计MTI滤波器时利用到杂波谱的相关系数:f 7 * 2-f X(以厂旧 T )-f 2r = R(i, j) = exp=i=i(6.5.31)ij2b 2f其中T =irrLrw r 2223L2 n112rrrwr32333n2=13MMOMMMr_rKrLwLr L n 2n 2nnn1n(6.5.32)1，W2吃通过解下面方程得 T i，匕是杂波的多普勒频率，改变匕的取值可以改变滤波器凹口的位置。N阶滤波器的权值w i=1040123归一化频率f/

24、fro o o O 2 4 6 8 - - - -0-80一1c1110246810归一化频率f/fr000246-,)0-应响率频采用预测误差算法得到的MTI滤波器频率响应特性如图6.5.12所示:-100 c-1(a)等 T(b)参差变T图6.5.12预测误差算法MTI滤波器频率响应6.5.3 MTD滤波器(1)MTD滤波器原理早期的动目标显示雷达性能不高，其改善因子一般在20dB左右。随着在系统设计与实现技术的改进、数字技术的提高，主要依靠信号 处理的潜在能力，MTI雷达的性能还将进一步的改善和提高： 增大信号处理的线性动态范围； 增加一组多普勒滤波器，使之更接近于最佳滤波，提高改善因子

25、； 能抑制地杂波(其平均多普勒频移通常为零)且能同时抑制运动杂波(如云雨、鸟群、箔条等)； 增加一个或多个杂波图，可有帮助检测切向飞行大目标等作用； 增加CFAR检测电路。做了上述改进的系统称之为动目标检测(MTD)系统。本节介绍的多普勒滤波器组就是一种MTD滤波器根据最佳滤波理论，当杂波功率谱C( f )和信号频谱S( f )已知时，最佳滤波器的频响是:式中，七是使滤波器能够实现而附加的延迟时间，(6.5.33)式(6.5.34)的滤波器可分成两个级联的滤波器H1(f)和H2(f)其传递函数形式为：(6.5.34)H 2( f) = S *( f )e - f(6.5.35)粗略的认为，其中

26、气(f)用来抑制杂波，相当于MTI滤波器，如图6.5.13(a)；H2( f)用来对脉冲信号匹配。对于相参脉冲串，H (f)又可以表示为： 2H 2( f) = H 21( f) H 22( f)(6.5.36)式中，H21(f)为单个脉冲的匹配滤波器，通常在接收机中放实现，H22(f)是对相参脉冲串进行匹配，它利用了回波脉冲串的相参性进行相参积累。H22(f)是梳齿形滤波器，齿的间隔为脉冲重复频率f，如图6.5.13(b)中0号滤波器；齿的位置取决于回波信号的多普勒频移，9 / 15而齿的宽度应和回波谱线的宽度一致。由于实际中fd不能预知，因此要用一组乡里且部分重叠的滤波器组(如图6.5.1

27、3(b)中0-7号滤波器)， 覆盖整个多普勒频率范围，这就是多普勒滤波器组所要完成的功效。w)固定杂波 气象杂波 动目标回波 Jr fHff杂波剩余r(a)动目标显示滤波(b)多普勒滤波器组特性输入N个脉冲的横向滤波器组有N-1根延 冲不同的加权并求和后，可以做成N个相邻的 同时得到N个滤波器的响应，则在图6.5.14中 输出并相应的加权。图6.5.14 MTD滤波设加在第k个滤波器的第i个输出端头的加权值图6.5.13 MTI与MTD滤波器特性w = e-j2顶i-1)k n,i -1,2,L N ik迟线，每根延迟时间孔=1 f，经过各脉 窄带滤波器组，频率覆盖0到fr。如果要 横向滤波器

28、的每一抽头应该有N个分开的 器组成框图为：(6.5.37)k表示标号从0到N-1的滤波器,每一个k值对应一组不同的加权值，相应地对应一个不同的多普勒滤波器响应。图6.5.13(b)中所示滤波器响 应是N=8时按式(6.5.38)加权所得各标记k的滤波器频率响应，k取0-7。第2个滤波器，即当k=1时，峰值响应产生在f/8,以及f + fJ8， 2f + fk取其他值，可以此类推。可写出横向滤波器按式(6.5.38)加权时的脉冲响应及其频响函数:(6.5.38)(6.5.39)h (t) = IL5 t - (i - 1)Te-j冲-1)k ni=1H (f)= e- j2兀 ft ej源(i-

29、1) fT-k n ki=1滤波器振幅特性;H (f )|= 丈 ej2(/-1)fT-kN ki=1sinK N(fT - k N)sinK (fT - k N)(6.5.40)滤波器峰值产生于sin兀(fT -kN) = 0。0号滤波器的中心位置在零频以及重复频率的整数倍nf处，通过没有多普勒频移的杂波，因此对地杂波没有抑制能力。所以MTD滤波器常常和MTI滤波器配合使用。用横向滤波器实现窄带滤波器可不采用式(6.5.38)所示的权值，而是根据特定的需要灵活地选用不同的加权矢量，可以根据不同频率处 设置特性相异的滤波器。每个滤波器都有一定的副瓣，幅瓣的大小决定着杂波抑制能力的大小。为了压低

30、副瓣，一般都需要对数据做加窗处理。目前常用的窗函 数有矩形窗、三角形窗、Hanning窗、Hamming窗、Blackman窗、Gaussian窗、Kaiser窗、Bartlett窗。MTD系统的实现通常采用MTI+FFT 的形式。MTD滤波器组对于动杂波的抑制效果较好。由于MTD滤波器组中每个窄带滤波器的主瓣宽度只占对消器通频带的1/N宽度，因而输 出端的信噪比有相应提高：白噪声提高近N倍；有色杂波各滤波器输出的信噪比提高程度因距离杂波中心及杂波谱宽的不同而不同。MTD 滤波器组能较MTI对消器提高改善因子是因为细分频带后，各滤波器的杂波输出功率只有各自通带范围内的杂波谱部分，而不是整个多普

31、 勒频带内的杂波功率。副瓣的影响会降低改善因子，解决办法是MTD滤波之前先采用MTI对消处理；或者采用加权法降低各个滤波器的副 瓣，但代价是滤波器主瓣有所加宽。(改善因子的分析可以放在后一节)(2)MTD滤波器设计(a)基于MTI+FFT的MTD滤波器组基于MTI+FFT的MTD处理分两步进行：首先进行MTI滤波，主要用于抑制地物杂波，一般用三脉冲对消器；然后再按式(6.5.38)进行 FFT处理。FFT的每点输出，相当于N点数据在这个频率上的积累，也可以说是以这个频率为中心的一个带通滤波器的输出。对N个数据 的FFT的运算量中复乘法次数为N- 10g2N，比DFT的复乘法次数N 2要少，特别

32、是N很大时，FFT的运算量要少得多。每个滤波器都有一定的副瓣，副瓣的大小决定着杂波抑制能力的大小。为了降低副瓣，一般都需要加窗。目前常用的窗函数为海明窗 (Hamming)，加窗可降低副瓣电平，如图6.5.15，副瓣电平从-20dB左右下降到-40dB左右，但各滤波器的主瓣有一定展宽。应响率频-100 00.20.40.60.81归一化频率f/fr0,rno-)E应响率频00002468 00.20.40.60.81归一化频率f/fr(b)加hamming窗8点FFT滤波器组图6.5.15 FFT滤波器组频率特性0o o o O2 4 6 8 - - - -100(a)8点FFT滤波器组(b)

33、点最佳MTD滤波器组点多普勒横向滤波器只在所需的多普勒处理波段中的某一点上达到最佳，而在所有其它点都是不匹配的。点最佳MTD滤波器组就是用 许多点横向滤波器填满感兴趣的多普勒区域。尽管实际应用中常用n个点滤波器等间隔填满多普勒范围，使n等于处理的相干脉冲数，但我 们也可以设置任意个横向滤波器个数。不模糊多普勒覆盖波段0工2旦4宣&二旦210 prf101010101010101010零多普勒x=n fd/fr图6.5.16点最佳MTD滤波器组多普勒覆盖波段示意图点多普勒横向滤波器复数输入信号表示为:(6.5.41)s = Aejd ZW (t -iT)式中A是幅度i =0是多普勒角频率，n是处

34、理脉冲数，T是雷达脉冲间的间距。信号矢量可以表示为 d(6.5.42)这里:s广 EK-1)T； kT，2.n(6.5.43)协方差矩阵:Ms = E (ss *t )(6.5.44)对它归一化(处以A2)1O e- jwdT L e 0 ML L L Oe- jwd (n-1)TMejwd (n-1)T L归一化学方差矩阵的元素m为:s(6.5.46)m = ej%( k-i)Te -加 d( j-i)t = ed( k - j)tskj式中k为行数，j为列数。k=1,2,n; j=1,2n对应于R -1M的最大本征值的矢量即为最佳权值。但不一定非要求解本征矢量，现给出了用统计检测理论求出的

35、最佳权值: n sW0 = R tS*(6.5.47)其中Rn为杂波为干扰矩阵，可以表示为杂波协方差矩阵Rc和接收机噪声矩阵之和R = R +Q 21其中R =A】+R 2式中R 1地物杂波，R 2为气象杂波Pn-11R =b 2c cPn-1(6.5.48)(6.5.49)P1L(6.5.50)其中P是杂波相关系数，如式(6.5.31)。200-80-1001-120 00.20.40.60.8归一化频率f/fr-20B -40H -60200-20-40-60-80-100-12000.20.40.60.81归一化频率f/fr(a)N=10,10脉冲对消滤波器组频响(b)k=5时单个滤波器

36、频响图6.5.16点最佳MTD滤波器组频率响应(c)等间隔最佳MTD滤波器组 等间隔多普勒横向滤波器在某个间隔内检测性能最佳，这个间隔是多普勒处理段的一个子集。把多普勒覆盖段分成相等的间隔，每个间隔内都设计一个最佳横向滤波器，这就构成了一个多普勒处理机。 I I I I I I I I I I II1 I2 I3 I4 I5 :6 I7 I8 I9 I1 I多普勒偏移/雷达PRF1.0图6.5.17等间隔最佳多普勒覆盖波段示意图等间隔多普勒处理机反映出一个同输入信号多普勒偏移的分布有关的补偿性的假设。这种处理机在比点多普勒处理机在更宽的区域上得到最佳，并且具有MTI处理机所没有的S/N改善。但

37、是它的最佳权值和改善因子的计算比其它类型的处理机要困难。 等间隔多普勒处理机的目标多普勒概率分布为：JnT,(2k-1)/2nT f wtM w*=#fwtrw*，可得:图6.4.5二项式系数的MTI改善因子(N=2,3,4)点最佳MTD滤波器组设计算法的改善因子是将式(6.5.48)代入I|w|2 + 2云云 Re(w*w )cos( j - i) T - Im(w*w )sin( j - i)w TiL i jdi jd iL=j=1i=0 j=i+1 弓.i j=12 + 2习研(i- j)Re(w*w)i=0 j=i+1(6.5.65)上式就是计算任意一个多普勒频率。处改善因子的一般表

38、达式。但必须记住上式中的权值要在横向滤波器的最佳点多普勒频率处计算。最 佳点的改善因子：I =& +2云 a cos( j 一 i)s Ti=1i=1 j=i+1(6.5.66)等间隔最佳MTD滤波器组的改善因子，是用(6.5.55)的平均加权求解改善因子: i c 一n一2 *-乙 w 2 + 2乙乙Re(w w)cos( / - i)T - Im(w w)sin( / - i)TI i|_i j J d i j J 】dI = 4=1i=0 j=i+1町 + 2 习 p(i-j)Re(W*W7) j=1i=0 j=i+1(6.5.67)对于相等数目的分布滤波器，等间隔最佳横向滤波器具有比点最佳滤波器更均匀的响应且在零多普勒频率附近时有更好的响应，但点最 佳设计的改善因子峰值超过等间隔设计的峰值。