TR-R2多站雷达系统的近程应用分析与仿真

2011-04-18 23:48:34来源: 互联网

本文将宽带FMCW发射波形引入近程TR-R2多站雷达系统,利用FMCW雷达信号的高分辨率特点,分析了近程TR-R2系统的目标定位,提出了利用回波信号频谱估计细柱状目标长度的方法,并导出了计算目标长度的解析公式.速度矢量在目标探测与分析中,尤其是在近程目标的识别中具有重要的意义.本文首次分析了TR-R2系统中目标速度的测量方法并给出了算法公式.文中对系统性能进行了较全面的分析与仿真.
  关键词:多基地雷达;调频连续波;系统分析

Analysis and Simulation for the Short Range Application of TR-R2 Multistatic Radar System

WEI Chong-yu,XU Shan-jia,WANG Dong-jin
(University of Science & Techndogy of China,Hefei 230027,China)

  Abstract:The FMCW is introduced into short range TR-R2 multistatic radar system and an analysis is given,for locating the geometric center of a thin cylindrical target in short range using the high resolution of the wideband FMCW signal.A method and its analytical expressions for calculating the length of the target using echo spectrum are presented.Target velocity measurement is of important significance for target detection and analysis,especially for recognizing a target in short range.In this paper,an algorithm for measuring target velocity with TR-R2 system is analyzed for the first time,and formulae for the calculation of velocity are shown.System performances are discussed and simulated with the present algorithm.
  Key words:multistatic radar;FMCW;system analysis;simulation

一、引  言
  多站雷达具有许多适应现代战争环境的独特优点,所以近年来重新获得人们的重视.然而,以往这方面的研究[1~3]大都是针对远程应用.对多站雷达近程应用的研究很少.对于那些目标密度低,目标形状简单(一般为细柱状),但要求系统反应速度快的近程应用场合(作用距离在3km以内),怎样高速地对目标几何中心定位,分析目标特性并估计目标尺寸;如何分析系统对目标的探测精度;多站系统该如何布局等等,都是尚未很好研究的实际问题.本文将分析TR-R2多基地雷达系统的近程应用.
  高线性度FMCW发射波形具有高分辨率特点,雷达系统采用宽带高线性度FMCW发射波形时能够实现对目标的高精度定位.但单纯的FMCW信号会发生严重的时延和频移耦合,从而使FMCW信号不能发挥应有的效能.同时这种信号的匹配滤波也会因多卜勒失真而变差.当目标速度v<0.1c/TB(c为光速,T为信号时宽,B为信号带宽)时,这种影响可以忽略[4],然而,对于大时宽带宽信号(如时宽1ms,带宽500MHz),且目标速度较高时,上述条件一般不再满足.另外,单频脉冲是一种简单的发射波形,其信号处理容易,频移测量精度高且与时延不发生耦合,很适合于用来测量目标速度.为此本文将FMCW信号与单频脉冲信号两者的优点结合起来,采用两种波形交替发射,对目标回波综合处理的办法完成时延与频移的解耦,从而构成了一种能充分发挥FMCW信号效能,同时实现高精度定位与测速的新方法.系统发射信号的频率关系如图1所示.

t56.gif (1198 bytes)

图1 发射频率变化规律

  在0-T1时段内系统发射单频脉冲,完成多卜勒频移测量.在T1-T时段内,系统发射高线性度FMCW信号,接收机根据0-T1时段的频移测量,在回波信号中扣除频移影响,完成时延的高精度测量,实现目标几何中心定位,而后利用回波谱宽估算细柱状目标长度,并求出目标速度矢量.
  文中假定测量周期T在ms量级,在这段时间内近似认为目标速度是不变的.同时考虑系统的实时性要求,本文将不讨论目标细节的识别,而是根据近程环境特点仅仅估计细柱状目标长度.文中所采用的算法也都是能一次完成计算的解析算法.

二、几何中心定位算法
  多站系统的几何布局如图2所示.在T1-T时段内,t时刻的发射频率为

νT=ν0+kt (1)

t时刻经目标几何中心到达Si的信号频率为

g57-1.gif (555 bytes) (2)

式中RT,Ri分别为目标几何中心与ST及Si的距离,k为FMCW线性调频的速率.t时刻直达Si的信号频率为

g57-2.gif (417 bytes) (3)

RTSi为发射机直达Si的距离.假设已根据0-T1时段测出的目标多卜勒频移,并在式(2)、(3)中扣除了该频移的影响.这时,Si对νTSi与νTOSi进行相关处理,求得频差νdi

g57-3.gif (790 bytes) (4)

t57-1.gif (1447 bytes)

图2 多站雷达系统布局

t57-2.gif (1103 bytes)

图3 多站系统的收发频率关系

假设目标坐标为(x,y,z),发射机布在S3处.并设di为LTi与LT3之差,即

di=LTi-LT3=Ri-R3,i=1,2,3 (7)

  根据距离关系,将站坐标Si(xi,yi,zi)和目标坐标代入式(7)展开整理并求解可得

g57-4.gif (611 bytes)
rh=[x y]T=W-1(β-d)
ra=[x y z]T=[f1(η) f2(η) f3(η)]T (8)

式中 a=1+dTWTW-1d;
b=2rTh3W-1d-2dTW-TW-1β-2z3;
c=βTW-TW-1β-2rTh3W-1β-β3,
g57-5.gif (1345 bytes)
x3j=x3-xj,y3j=y3-yj,z3j=z3-zj,j=1,2
β3j=0.5(d2j+2djR3-ρ2j+ρ23),j=1,2
这里ρi与Si与原点的距离.
η=[d1 d2 R3]T=[η1 η2 η3]T.只要三个雷达站没有布局在同一直线上,W就是非奇异的,式(8)有确定的解.式中正负号的选取可参考文献[3],但通过合理选择坐标系可以将目标置于上半球或下半球,从而使符号选取简化.

三、定位性能分析
  1.定位坐标偏移量
  当TR-R系统对目标定位时,直接测得量的d1,d2,R3,设测量误差矢量为

δη=[δd1 δd2 δR3]T=[δη1 δη2 δη3]T (9)

目标定位误差矢量为

δra=ra-r0=[δx δy δz]T (10)

其中r0为目标的真实位置矢量,由式(7)得

δηi=δLTi-δLT3=δLTi-2δR3 (11)

上式表明,δηi之间不独立.
  设各雷达站对频率的测量相互独立,频率测量误差符合零均值正态分布.则各雷达站对距离的测量相互独立,且测量误差δLTi及δR3都符合零均值正态分布.而δηi与δLTi及δR3是线性关系,因而δηi也符号零均值正态分布,这时有

E[δηi]=E[δLTi]-2E[δR3]=0 (12)

但是,由于目标几何中心定位坐标与ηi的函数关系是非线性的,因此δηi将使定位坐标产生偏移.即使E[δra]≠0,偏移量的期望值可通过式(8)估算.
  为求δra的数学期望,将ra=f(η)在r0附近展为η的Taylor级数,保留到二阶偏微分项,并考虑到δηi的相关性,可得定位坐标偏移量期望值的表达式

g57-6.gif (2500 bytes) (13)

式中σ2为单站雷达测距方差,b1,b2,b3分别表示坐标偏移量bx,by,bz.
  2.目标定位误差的GDOP因子
  一阶近似情况下,δra的协方差矩阵P为
p=E{(δra-E[δra][δra-E[δra]]T)}=FSFT-bbbTb (14)
其中  bb=[bx by bz]T,

g57-7.gif (3209 bytes)

P的对角线元素为定位坐标方差σ2x,σ2y,σ2z.目标定位误差的GDOP因子为

g57-8.gif (626 bytes) (15)

四、目标速度测量
  在0-Ti时段系统发射单频脉冲,已经测得运动目标回波的多卜勒频移,由此可求得目标速度在相应方向的投影,将这些投影看做空间矢量,就可估计目标速度.
  在TR-R2多站雷达系统中,S3构成单站雷达的同时还分别与S1,S2构成双基地系统.根据单站与双基地系统多卜勒关系得到方程

g57-9.gif (1796 bytes) (16)

vi为目标速度在Si——目标视线上的投影,νdi为Si中测得的目标多卜勒频移,ts58.gif (92 bytes)i=vi+v3.由式(16)解得vi,并以矢量

v=[q1 q2 q3]T (17)

表示目标速度,利用Si——目标视线的方向余弦(cosαi,cosβi,cosγi),将vi表示为

vi=q1cosαi+q2cosβi+q3cosγi,i=1,2,3 (18)

整理后写成矩阵形式解得

V=Φ-1μ (19)

式中

g58-1.gif (1833 bytes)

Φ为已知,因此可求得速度v,速度的数值由下式计算

g58-2.gif (461 bytes) (20)

五、速度测量性能分析
  当各雷达站的频率测量相互独立,测量误差符合具有相同方差的零均值正态分布时,由于ts58.gif (92 bytes)1与νdi是线性关系,则ts58.gif (92 bytes)i也是正态分布的

g58-3.gif (717 bytes) (21)

式中σ2v为ts58.gif (92 bytes)1的方差,ts58.gif (92 bytes)0i为ts58.gif (92 bytes)i的均值,式(16)表明vi与ts58.gif (92 bytes)i是线性关系,故目标在各雷达站方向的速度投影vi也符合正态分布.同样式(19)表明qi与vi也是线性关系,因此qi也符合正态分布,其概率密度公式为

g58-4.gif (927 bytes) (22)

其中,g58-5.gif (772 bytes),为qi的均值,σ2qi=[k21i+k22i+1/4(k1i+k2i-k3i)2]σ2v为qi的方差,而kij=Δ′ij/Δ′,Δ′=det(Φ),Δ′ij为Δ的代数余子式.
  根据式(22)写出q1,q2,q3的联合概率密度函数

g58-6.gif (1218 bytes) (23)

式中B为qi的协方差矩阵,其元素为

B(i,j)=E[(qi-q0i)(qj-q0j)],μt=[q01 q02 q03]T

为求v的概率密度函数,对式(23)作如下的变量代换,令

q1=ξcosθcosφ,q2=ξcosθsinφ,q3=ξsinθ.

变换后的变量取值范围相应变为

g58-7.gif (806 bytes)

此变换的Jacobian行列式为J=ξ2cosθ,这样ξ即v的概率密度函数可以写成

g58-8.gif (1621 bytes) (24)

目标速度期望值及其方差分别为

g58-9.gif (833 bytes) (25)
g58-10.gif (1053 bytes) (26)

六、定位误差对速度测量的影响
  目标速度估算是在定位之后完成的.由于速度估算公式中要用到Si——目标视线的方向余弦,这些方向余弦的精度取决于目标几何中心定位误差.因此,定位精度将直接影响速度测量精度.为便于分析,下面在考虑定位误差影响时,不考虑其它因素.
  令δv为因定位误差所产生的速度误差,将速度v展为定位坐标的Taylor级数,在一阶近似条件下求得其偏移量计算式为

g58-11.gif (1003 bytes) (27)

速度估计方差为

g58-12.gif (3004 bytes) (28)

其中ε2i=2σ2

g58-13.gif (5218 bytes)

七、细柱状目标长度估计
  雷达发射FMCW波形时,体目标回波频谱将占据一定的带宽,对回波信号作频谱分析则可以提取到此频宽信息.
  设细柱状目标的长度为m,接收机Si的接收信号谱宽为Δνi,目标在Si方向的投影为mi,则有如下关系成立

g58-14.gif (1995 bytes) (29)

这里,目标投影及频谱宽度均应有正负号,符号的选取由目标相对各雷达站的位置关系及目标的空间取

[1] [2]

关键字:分析  仿真

编辑:神话 引用地址:http://www.eeworld.com.cn/mndz/2011/0418/article_7233.html
本网站转载的所有的文章、图片、音频视频文件等资料的版权归版权所有人所有,本站采用的非本站原创文章及图片等内容无法一一联系确认版权者。如果本网所选内容的文章作者及编辑认为其作品不宜公开自由传播,或不应无偿使用,请及时通过电子邮件或电话通知我们,以迅速采取适当措施,避免给双方造成不必要的经济损失。
论坛活动 E手掌握
微信扫一扫加关注
论坛活动 E手掌握
芯片资讯 锐利解读
微信扫一扫加关注
芯片资讯 锐利解读
推荐阅读
全部
分析
仿真

小广播

独家专题更多

富士通铁电随机存储器FRAM主题展馆
富士通铁电随机存储器FRAM主题展馆
馆内包含了 纵览FRAM、独立FRAM存储器专区、FRAM内置LSI专区三大部分内容。 
走,跟Molex一起去看《中国电子消费品趋势》!
走,跟Molex一起去看《中国电子消费品趋势》!
 
带你走进LED王国——Microchip LED应用专题
带你走进LED王国——Microchip LED应用专题
 
电子工程世界版权所有 京ICP证060456号 京ICP备10001474号 电信业务审批[2006]字第258号函 京公海网安备110108001534 Copyright © 2005-2016 EEWORLD.com.cn, Inc. All rights reserved