巴克码多普勒容限带宽的扩展杨文华 宋力平 王其扬 (航天工业总公司八院804所 上海200082)文 摘 介绍和分析了多普勒容限带宽扩展的几种方法,提出用设计包络约束(EC)滤波器的方法设计巴克码的多普勒频移补偿滤波器,可以同时进行旁瓣抑制和多普勒容限带宽的扩展,并且使信噪比损失较小。
经计算机仿真证明有较好的结果。
主题词 巴克码 多普勒容限带宽 包络约束滤波器 多普勒频移补偿滤波器收稿日期:1997 05 29系上海市宇航学会会员。
0 引言相位编码信号有许多优点,但它有多普勒频移敏感性。
当雷达信号回波存在多普勒频移时,脉压后主瓣加宽,旁瓣增大,甚至起不到脉冲压缩的作用,这就限制了它的应用,因此多普勒容限带宽的扩展日益受到重视。
一方面,人们通过寻找新的码型来解决这个问题,如NLFM 码(Non Linear FM Sig nals)就是通过取样和量化非线性调频信号而来的,它具有低旁瓣特性和多普勒频移的不敏感性;另一方面,人们在雷达系统和旁瓣抑制滤波器上进行了研究。
文献[2]提出了采用非相干处理结构来降低多普勒频移的敏感性,文献[3]在非相干处理结构的基础上,在某一多普勒频率上采用不同的最佳准则设计多普勒频移补偿滤波器,得到了比较好的效果。
在实际应用中,也采用了多路匹配滤波器等各种各样的处理方法。
1 几种多普勒补偿滤波器的设计在信号检测链路中,为了降低旁瓣电平,采用了旁瓣抑制滤波器。
旁瓣抑制滤波处理在求模运算以前进行,称之为相干处理;反之,旁瓣抑制滤波处理在求模运算以后进行,称之为非相干处理[2],如图1所示。
图1 信号检测链路框图分析表明:采用非相干处理结构,可以降低多普勒频移的敏感性和扩展多普勒容限带宽[3]。
但是这样扩展是有限的:当目标高速飞行或旁瓣抑制滤波器应用于多普勒(PD)雷达时,非相干处理的多普勒容限带宽扩展已不能满足需要。
由于非相干处理结构的旁瓣抑制滤波器所考虑的只是多普勒频率的绝对值,其性能是关于零多普勒频率对称的,所以,就能以某一多普勒频率为中心频率,设计非相干处理结构的旁瓣抑制滤波器,得到较大的多普勒容限带宽。
根据文献[2]的分析,不论旁瓣抑制滤波器是非相干处理型的,还是相干处理型的,其主瓣附近的第一旁瓣随多普勒频移增大而增加得非常快。
若考虑第一旁瓣,则旁瓣抑制滤波器的多普勒容限带宽是有限的。
因此,放宽对第一旁瓣的约束,把第一旁瓣的上升看作是主瓣的加宽是很有意义的。
当然,这会使雷达的距离分辨力下降。
曾有人运用最小平方近似逆滤波(LS)的方法设计出了补偿式非相干处理结构的旁瓣抑制滤波器。
下面介绍用LS法设计补偿滤波器的原理。
对给定的输入序列,通过滤波器后得到的输出序列与希望得到的输出序列之间的均方差最小。
假设实际输入序列为(a0,a1, !,a M)T,滤波器的加权序列为(b0,b1,!, b N)T,实际输出序列为(c0,c1,!,c M+N)T:c=a b(1)对于旁瓣抑制问题而言,我们希望输出序列为:在L处为1,而在别的时刻为0,则最小均方逆滤波就是使下式为最小V=c20+c21+!+(c L-1)2+!+c M+N2(2)式中V∀∀∀实际响应和理想输出序列之间的能量差。
以巴克码为例,当码元时宽T=0.7 s,多普勒频移f d=20kH z时,同时放宽第一旁瓣的约束,并对理想输出的约束加以适当处理,则可以得到长度为25的LS补偿式旁瓣抑制滤波器的权系数。
表1为对最大值归一化的权系数。
表1 LS补偿式旁瓣抑制滤波器的权系数-0.04964 0.00010 -0.04879 -0.00180 -0.04833 -0.03657-0.05363-0.04003-0.04757-0.03791 -0.057090.21673 1.000000.21673-0.05709 -0.03791-0.04757-0.04003-0.05363-0.03657 -0.04833-0.00180-0.048790.00010-0.04964同样,文献[3]用线性规划(LP)的方法设计了13位巴克码的补偿滤波器。
设子码宽度T=0.7 s,选用15kHz为多普勒频移补偿点,运用线性规划。
a目标函数:使主瓣峰极大化,即max#M-1i=0h i A l i-(M-2N+1)/2(3) b约束条件:限制最大旁瓣电平;满足所要求多普勒容限带宽下尽可能限制主瓣峰的展宽#M-1i=0|h i A l i-k|∃11-2N∃k∃M-1k%(M-2N+1)/2k%(M-2N+1)/2l(4) #M-1i=0|h i A li-(M-2N+1)/2|∃s(5)式中 N∀∀∀巴克码的位数;{h i}∀∀∀旁瓣抑制滤波器的权序列,i=0,1,!,M -1;{A l }∀∀∀巴克码经匹配滤波后的模序列,其多普勒频移f d =15kHz , l =1,2,!,2N -1。
在满足多普勒容限带宽下,M 应尽可能小,同时选择s 使主瓣峰作最小展宽。
经过计算、比较,以及经适当的合并、简化,得到M =29,s =20。
此时,可获得对理想滤波器逼近的最大多普勒容限带宽。
表2为LP 补偿式旁瓣抑制滤波器的权系数。
表2 LP 补偿式旁瓣抑制滤波器的权系数0.033208 0.019531 -0.0527340.011719 -0.0488280.011719-0.042969-0.042969-0.042968-0.042969-0.033208-0.042969-0.0527340.148438 1.0000000.148438-0.052734-0.042969-0.033208-0.042969-0.042969-0.042969-0.0429690.011719-0.0488280.011719-0.0527340.0195310.033208在以上两种方法中,考虑的是最佳的主旁瓣比和最佳的均方旁瓣。
如果信号在通过旁瓣抑制滤波器以后,有很大的信噪比损失,那末就会严重影响后面的检测性能。
因此,信噪比损失也是旁瓣抑制滤波器的一项重要指标。
设计包络约束滤波器的出发点就是使信噪比损失最小。
由此可见,用设计EC 滤波器的方法设计非相干处理结构的补偿滤波器是有实际意义的。
2 包络约束(EC)滤波器包络约束滤波器是Evans J 等人提出的。
其原理如下:假设一输入信号为s ,通过一滤波器后,其输出( )为一介于上包络e +和下包络e -之间的脉冲形信号,如图2所示。
设U 为滤波器有限长时不变冲激响应,则=sU(6)在实际中,噪声是普遍存在的。
以最普遍的白噪声为例,定义最优的包络抑制滤波器为在满足包络约束的条件下,噪声能量最小的滤波器。
因为噪声能量正比于范数&U &的平方,所以可以用数学公式表示minimize &U &2(7)条件是e -(t)∃ (t)∃e +(t)(8)这里 (t)=∋U ( ) s (t - )d(9)图2 EC 滤波器(a)框图;(b)脉冲包络示意图;(c)用于雷达脉冲压缩的脉冲包络示意图。
由于规定了上下包络,在物理意义上可以认为限制了主旁瓣比。
在此条件下,EC 滤波器使滤波器的信噪比损失最小。
在离散情况下,假设信号s =(s 1,s 2,!,s m )T ,滤波器加权因子U =(u 1,u 2,!,u n )T ,输出信号 =(!1,!2,!,!N )T =SU,N=m+n-1,其中S为N(n矩阵S=s100!0s2s10!0s3s2s1!0!!!!!s m s m-1s m-2!s m-n0s m s m-1!s m-n+1!!!000!s m04则式(7)~(9)可以写成minimiz&U&2=U)U(10)条件是|SU-d|∃e,即S U-d∃e(11)-S U+d∃e(12)其中 d=(e++e-)/2=(d1,d2,!,d N)T e=(e+-e-)/2=(e1,e2,!,e N)T 文献[1]详细推导了这个问题的解决过程。
本文只引用了最后的结论。
根据对偶理论,当U*=-S)∀*/2时,即可满足式(10)。
其中∀*的递推公式为∀k+1=∀k+a I(∀k)(13)式中I(∀k)为N(1矩阵;k为叠代次数, k=1,2,!,∗;I(∀k)=(l1,l2,!,l N),其中l(∀k)=-[(SS)∀)i/2+d1+e i]i+#+,Q+-[(SS)∀)i/2+d1-e i]i+#-,Q-i+Q07(14)i=1,2,!,N 其中:#+={i: ∀>0}#-={i: ∀<0}Q+={i+#0:-[(SS)∀)i/2+d1+ e i]>0}#0={i: ∀=0}Q-={i+#0:-[(SS)∀)i/2+d1-e i]<0}Q0={#0-Q+-Q-}则有U k=-S)∀k/2, k=1,2,!(15)且&U k-U+&−0 k−∗(16)这样,我们就可以得到EC滤波器的加权系数。
设计EC滤波器的出发点就是使噪声能量最小。
尽管文献中是以白噪声为例进行推算的[1],但可以很容易地推广到色噪声的情况。
另外,也可以不采用&U&准则而很方便地采用别的准则(如最小均方准则)来设计滤波器。
和LS法一样,我们选用了20kH z的多普勒频移补偿点,码元时宽T=0.7 s。
此时输入信号是有20kHz多普勒频移的13位巴克码的匹配滤波后的模序列{A l}。
设上包络e+=(0.01,0.01,!,0.01,10, 10,10,0.01,!,0.01)T;下包络e-=(-0.01,-0.01,!, -0.01,9.9,-0.01,!,-0.01)T。
运用式(13)进行叠代运算,再用式(15)进行计算,即可得到EC滤波器的权系数。
当k=200,a=1.5时,通过计算机编程运算,得到最大值归一化的权系数如表3所示。
表3 EC补偿式旁瓣抑制滤波器的权系数根据得到的权系数,运用计算机进行模拟仿真,可以得到以上三种补偿滤波器的性能,如图3和表4所示。
图3 四种滤波器的多普勒容限带宽表4 三种补偿滤波器的比较频偏f d /kH z 最大旁瓣/dB第一旁瓣/dB积累旁瓣/dB 信噪比损失/dB 0-31.80-29.16-32.49 -12.91 -16.39 -14.04-20.82-18.86-21.31 1.2540.5630.6685-32.11-31.13-32.39-11.51-14.19-12.41-22.48-20.41-22.870.6650.5370.58910-32.05-33.89-32.31-10.35-12.48-11.08-23.77-21.98-24.010.5830.5190.52415-32.01-36.81-32.23-9.357-11.06 -9.964-24.46-23.35-24.490.5610.5550.52420-31.96-35.47-32.16-8.492-9.852-8.993-24.34-24.00-24.280.5960.6580.58125-31.92-33.94-32.09-7.72-8.786-8.129-23.48-23.40-23.160.6910.8160.69530-31.88-32.26-32.02-7.018-7.823-7.345-22.09-21.73-21.690.8431.0340.86935-31.77-30.69-31.67-6.353-6.925-6.607-20.56-19.90-20.15 1.0611.3171.10640-30.87-29.63-29.96-5.544-5.946-5.739-19.52-20.07-19.351.3871.7111.453注:表中的数据分别为LS ,LP 和EC 补偿滤波器的数据。