CN102967855A - 异地连续定标rcs测量中改进背景相减技术的信号处理方法 - Google Patents
异地连续定标rcs测量中改进背景相减技术的信号处理方法 Download PDFInfo
- Publication number
- CN102967855A CN102967855A CN2012105358165A CN201210535816A CN102967855A CN 102967855 A CN102967855 A CN 102967855A CN 2012105358165 A CN2012105358165 A CN 2012105358165A CN 201210535816 A CN201210535816 A CN 201210535816A CN 102967855 A CN102967855 A CN 102967855A
- Authority
- CN
- China
- Prior art keywords
- calibration
- background
- target
- rcs
- measurement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,该方法的步骤如下:步骤1:t1时刻,测量包含定标支架的测试场背景回波SBC(f,t1);步骤2:t2时刻,安装定标体,并测量定标体回波SC(f,t2);步骤3:t3时刻,测量包含目标支架的测试场背景回波SBT(f,t3);步骤4:t4时刻,安装目标,并测量目标回波ST(f,t4);步骤5:提取时变传递函数参数、幅度相位校正、背景相量相减和目标RCS定标处理。有别于传统RCS测量与定标处理模型,本发明提出的新模型为RCS测量与定标处理提供了新的理论依据,具有以下优点:引入目标宽带复散射函数的概念,用于指导确定测试与处理的流程、形成背景杂波提取和精确背景相减处理方法和算法,也用于分析影响外场测量不确定度的主要因素,具有普适性。
Description
技术领域
本发明涉及雷达散射截面(RCS)测量以及雷达信号处理的技术领域,特别涉及一种异地连续定标RCS测量中改进背景相减技术的信号处理方法。
背景技术
在雷达散射截面(RCS)外场测量中,通常采用异地连续定标测量技术,此时定标体和目标位于不同的雷达距离上,且一般采用低散射金属支架作为支撑定标体和被测目标的支架,其几何关系示意图如图1所示。
由图1,定标体和待测目标的接收回波功率均满足雷达方程(参见许小剑、黄培康,《雷达系统及其信息处理》,电子工业出版社,2010):
式中,Pr和Pt分别为雷达接收和发射功率;G为雷达天线增益;σ为目标的RCS;R为雷达到目标的距离;λ为雷达波长;L表示各种损耗的综合影响。
若目标RCS为σt,定标体RCS为σc,则两者之间的关系为:
式中Rc和Rt分别为定标体和待测目标的测量距离,为系统传输损耗项,故常数可通过测试中系统和几何关系参数直接计算得到;σc为定标体的精确RCS值,为已知量;Prc和Prt分别为定标体和目标回波的功率,有:
式中St和Sc分别为一次测量采样中,测量目标和测量定标体时的雷达复回波信号。
对于宽带扫频RCS测量,有:
St(f)=T(f)+Bt(f)+N (4)
和
Sc(f)=C(f)+Bc(f)+N (5)
式中,St(f)和Sc(f)分别表示测目标和测定标体时雷达接收到的回波信号;T(f)表示目标真实回波;C(f)表示定标体真实回波;Bt(f)和Bc(f)分别表示测目标和测定标体时的真实背景回波;N表示外部和接收机内部噪声影响;上述回波信号采样均为复相量。
由式(2)~(5),有:
由式(4)~(6)可见,仅当背景回波Bt(f)和Bc(f)、噪声N均为零时,按照(6)式定标处理后得到的目标RCS测量值才代表其真值:
因此,为了提高RCS测量精度,必须消除背景和噪声的影响。在实际测量中,一般采用大信噪比测量,也即噪声对回波的影响可以忽略,有:
St(f)≈T(f)+Bt(f) (8)
和
Sc(f)≈C(f)+Bc(f) (9)
因此,可采用背景向量相减技术提高RCS测量精度,此时的定标方程为:
比较式(10)和(7),并由式(8)和(9)可知,当可精确测得背景向量和时,即 时,由式(10)所得到的目标RCS测量值趋于目标RCS真值。因此,能否精确测得目标的RCS,在很大程度上取决于能否精确地完成背景向量相减。
在RCS外场测试中,通常采用上述图1所示的异地连续定标测试技术,即将定标体和目标置于不同的径向距离上,测量雷达设有两个硬件距离门,一个对定标体选通,一个对目标选通,并在任何测量周期内均同时对目标和定标体回波进行采样,从而实现RCS测量连续定标。此时,RCS测量和定标的步骤为:
(2)安装定标体和目标,并测量定标体回波Sc(f)和目标回波St(f);
上述现有技术的技术方案的缺点在于:研究表明,在测试场环境不受外界影响且测量背景、定标体和目标之间的时间间隔较短时,采用上述测量步骤可以完成目标的精确RCS测量。但是,在外场测试条件下,由于环境受天气等条件影响大,且对大型目标测量中目标吊装等需要耗费很长时间,即使采用异地定标体连续定标可以系统漂移造成的定标体回波随时间的变化,由于背景测量并不是连续测量的,而且它同目标测量之间的时间间隔长,由于测量雷达系统发生幅度与相位漂移和测试场环境因素发生变化,将导致按照上述测量与处理步骤不能很好地完成背景向量相减处理,进而导致目标RCS测试误差大大增加。
发明内容
本发明所要解决的技术问题:本发明提出一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,通过采用该信号处理方法,可以补偿雷达系统漂移和测试场环境因素带来的影响,大大提高背景相减的精度,进而提高RCS测量精度。
本发明采用的技术方案为:一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,该方法的步骤如下:
步骤-1:t1时刻,测量包含定标支架的测试场背景回波SBC(f,t1);
步骤-2:t2时刻,安装定标体,并测量定标体回波SC(f,t2);
步骤-3:t3时刻,测量包含目标支架的测试场背景回波SBT(f,t3);
步骤-4:t4时刻,安装目标,并测量目标回波ST(f,t4);
步骤-5:利用定标体测量数据提取时变传递函数的3个幅度相位补偿参数、进行测量数据的幅度与相位校正、背景相量相减和目标RCS定标处理,具体处理方式如下:
式中,Ei(f)和Es(f)分别表示目标处雷达入射场和雷达天线处目标散射场;它同RCS之间的关系为
假定RCS测量满足高信噪比测量条件,即噪声的影响可忽略,若引入测量系统-测试场的时变传递函数对雷达接收回波的影响,则对于宽带扫频幅相测量,测目标和定标体时各自的雷达回波信号可分别表示为:
ST(f,t)=HT(f,t)·[T(f)+BT(f,t)] (12)
和
SC(f,t)=HC(f,t)·[C(f)+BC(f,t)] (13)
式中,HT(f,t)和HC(f,t)分别表示测目标和测定标体时测量系统-测试场的传递函数,是随时间慢变化的,反映了测量雷达系统漂移和测试场类型及其电参数特性随时间的变化特性;ST(f,t)和SC(f,t)分别表示雷达接收到的目标回波和定标体回波,两者均受到测量雷达-测试场时变传递函数的影响;T(f)和C(f)分别表示目标和定标体在给定姿态下的真实散射场,在给定姿态和雷达参数条件下不随时间而变化;BT(f,t)和BC(f,t)分别表示测目标和测定标体时的固有背景散射,两者均随时间t变化,表示由于目标/定标体同支架之间的耦合散射或背景状态的改变,造成放置目标/支架和空支架测量时,固有背景散射的变化,注意到上述散射信号均表示为复数相量;
这样,在外场异地定标RCS测量中,背景相减和RCS定标处理的过程数学上可表示为:
式中t1、t2、t3、t4分别表示测量定标背景、测定标体、测目标背景、测目标的时间段,SBT(f,t)和SBC(f,t)分别表示安放目标支架时测得的背景回波和安放定标体支架时测得的背景回波,两者可表示为:
SBT(f,t)=HT(f,t)·BT(f,t) (15)
SBC(f,t)=HC(f,t)·BC(f,t) (16)
这样,由式(12)~(16),有:
式中,ΔC(f,Δt12)和ΔT(f,Δt34)为经背景抵消处理后的剩余背景误差,
式(17)~(19)构成了异地定标RCS测量中测试与定标处理的基本数学模型。
其中,为了作进一步分析,假定:
精确的RCS测量和定标要求:
其中A0(f)为幅频特性,φ0(f)为相频特性;且
ΔT(f,Δt34)=0,ΔC(f,Δt12)=0 (23)
即剩余背景抵消误差为0。
其中,若测试场和定标与目标支架均经过仔细设计和优化,一般测试过程中定标区和目标区的背景是不随时间变化的,即有BC(f,t2)=BC(f,t1)=BC(f),BT(f,t4)=BT(f,t3)=BT(f),则此时式(18)和(19)变为:
因此,若要求剩余背景误差为0,则要求:
HT(f,t4)=HT(f,t3) (26)
HC(f,t2)=HC(f,t1) (27)
综合式(20)~(22)、(26)和(27),有:
AT(f,t4)=AT(f,t3),φT(f,t4)=φT(f,t3) (28)
AC(f,t2)=AC(f,t1),φC(f,t2)=φC(f,t1) (29)
式(28)~(30)构成了对精确的RCS测量、背景抵消和定标处理的基本要求,但是,由于真实测量雷达的参数总是随时间漂移的,也即上述三式的要求是不可能得到满足的,除非在测量过程中采用必要的辅助测量和信号处理方法对时变的传递函数幅度和相位作补偿处理,使得经过补偿处理后式(28)~(30)同时得到满足。
本发明与现有技术相比的优点在于:
与现有传统的RCS测量与定标处理模型不同,本发明所提出的新模型为RCS测量与定标处理提供了新的理论依据,且具有以下优点:
(1)新模型引入目标宽带复散射函数的概念,并将测量雷达系统和测试场随时间变化的传递函数引入电磁散射测量的背景相量相减和定标处理,进而建立用于地平场条件下目标宽带标电磁散射测量背景相量相减、定标处理的完整数学模型,解决了外场测试中诸多时变因素导致的测量定标误差大的问题;
(2)新模型可用于指导确定测试与处理的流程、形成背景杂波提取和精确背景相减处理方法和算法,同时也用于分析影响外场测量不确定度的主要因素;
(3)新模型具有普适性,不但用于地面平面场,也适用于采用自由空间场的室内、外测试场等,具体应用中只需根据测试场的不同设计对数学模型作适当简化即可。
附图说明
图1为现有技术中异地定标RCS测量几何关系示意图;
图2为传递函数参数估计处理流程图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
为了便于讨论宽带RCS幅度和相位测量,定义宽带复散射函数为:
式中,Ei(f)和Es(f)分别表示雷达入射场(目标处)和目标散射场(雷达天线处);它同RCS之间的关系为
假定RCS测量满足高信噪比测量条件,即噪声的影响可忽略。若引入测量系统-测试场的时变传递函数对雷达接收回波的影响,则对于宽带扫频幅相测量,测目标和定标体时各自的雷达回波信号可分别表示为:
ST(f,t)=HT(f,t)·[T(f)+BT(f,t)] (12)
和
SC(f,t)=HC(f,t)·[C(f)+BC(f,t)] (13)
式中,HT(f,t)和HC(f,t)分别表示测目标和测定标体时测量系统-测试场的传递函数,是随时间慢变化的,反映了测量雷达系统漂移和测试场类型及其电参数特性随时间的变化特性;ST(f,t)和SC(f,t)分别表示雷达接收到的目标回波和定标体回波,两者均受到测量雷达-测试场时变传递函数的影响;T(f)和C(f)分别表示目标和定标体在给定姿态下的真实散射场,在给定姿态和雷达参数条件下不随时间而变化;BT(f,t)和BC(f,t)分别表示测目标和测定标体时的固有背景散射,两者均随时间t变化,表示由于目标/定标体同支架之间的耦合散射或背景状态的改变,造成放置目标/支架和空支架测量时,固有背景散射的变化。注意到上述散射信号均表示为复数相量。
这样,在外场异地定标RCS测量中,背景相减和RCS定标处理的过程数学上可表示为:
式中t1、t2、t3、t4分别表示测量定标背景、测定标体、测目标背景、测目标的时间段SBT(f,t)和SBC(f,t)分别表示安放目标支架时测得的背景回波和安放定标体支架时测得的背景回波,两者可表示为:
SBT(f,t)=HT(f,t)·BT(f,t) (15)
SBC(f,t)=HC(f,t)·BC(f,t) (16)
这样,由式(12)~(16),有:
式中,ΔC(f,Δt12)和ΔT(f,Δt34)为经背景抵消处理后的剩余背景误差,
式(17)~(19)构成了异地定标RCS测量中测试与定标处理的基本数学模型。
为了作进一步分析,假定:
精确的RCS测量和定标要求:
其中A0(f)为幅频特性,φ0(f)为相频特性;且
ΔT(f,Δt34)=0,ΔC(f,Δt12)=0 (23)
即剩余背景抵消误差为0。
现根据图1中的异地定标测量的几何关系作进一步分析。若测试场和定标与目标支架均经过仔细设计和优化,一般测试过程中定标区和目标区的背景是不随时间变化的,即有BC(f,t2)=BC(f,t1)=BC(f),BT(f,t4)=BT(f,t3)=BT(f),则此时式(18)和(19)变为:
因此,若要求剩余背景误差为0,则要求:
HT(f,t4)=HT(f,t3) (26)
HC(f,t2)=HC(f,t1) (27)
综合式(20)~(22)、(26)和(27),有:
AT(f,t4)=AT(f,t3),φT(f,t4)=φT(f,t3) (28)
AC(f,t2)=AC(f,t1),φC(f,t2)=φC(f,t1) (29)
式(28)~(30)构成了对精确的RCS测量、背景抵消和定标处理的基本要求。但是,由于真实测量雷达的参数总是随时间漂移的,也即上述三式的要求是不可能得到满足的,除非在测量过程中采用必要的辅助测量和信号处理方法对时变的传递函数幅度和相位作补偿处理,使得经过补偿处理后式(28)~(30)同时得到满足。这正是本项发明的基本出发点所在。
根据图1所示几何关系,在单站测量且收发共用天线的情况下,根据雷达方程(1),有:
其中,αT(t)和αC(t)为随频率线性变化的时变相位因子,βT(t)和βC(t)分别为不随频率变化的相位因子;c为传播速度,f为雷达频率。假设天线增益项同时包含了雷达发射机和发射天线增益随时间的漂移、接收机增益和接收天线增益随时间漂移的影响。
注意到如图1所示,对于技术设计良好的测试场,
AC(f,t2)=AC(f,t1),φC(f,t2)=φC(f,t1) (36)
(1)G(f,t)包含了发射机功率和相位漂移造成的系统增益影响;
(2)RT(t)和RC(t)为距离因子,对系统参数随时间的变化不敏感,故可认为是不随时间变化的常数,可记为RT(t)=RT,RC(t)=RC;
(3)LT(f,t)、LC(f,t)为衰减因子,对时间变化一般不敏感,但随雷达工作频率是变化的,可简记为LT(f)和LC(f);
(4)ΔR(t)的微小变化造成随频率线性变化的相位发生改变,相当于目标参考中心的微小平移,不会对目标RCS测量结果产生实质性影响,故也可假定为不随时间变化的常数,不会影响对后续问题的分析和讨论,可记为ΔR(t)=ΔR。
这样,式(20)和(21)可写为:
有
式中H0(f)为传递函数中可通过测量几何关系和雷达参数确定的幅度和线性相位因子,不随时间推移而变化,
以下讨论两种常见的RCS测量情况。
1)同地定标情况:目标和定标体放置在同一支架上,在时间上按先后顺序测量。若在RCS测试中采用同地定标,有RC=RT,LC(f)=LT(f),此时有A(f)=1,ΔR=0,故有:
H0(f)=1 (41)
2)异地连续定标情况:定标体与目标放置在不同距离上,利用两个距离门实现定标体、目标回波信号同时采集和连续定标测量。此时,一般测试过程中定标体与目标之间的距离较近,且雷达系统的相对带宽较小(否则难以满足地面平面场条件),故可认为近似等于常数,因此有:
由于在RCS定标处理中,可根据测量几何关系和雷达参数,并依据雷达方程精确计算出H0(f),故该项传递函数一般不会引入明显的定标处理误差。
这样,实际RCS测量和处理中可对定标公式(14)作如下修正:
注意到上式中仅第一个分式是随时间变化的。记Δtmn=tn-tm,表示tm和tn两次不同测量之间的时间间隔(m,n=1,,2,3,4)。这样,式(17)可重写为:
式(44)中,Hm(f,Δt24)表示在测量定标体(t2时刻)与测量目标(t4时刻)期间,测试场-测量系统传递函数随测试时间变化引起的测量不确定性,ΔT(f,ΔtT)和ΔC(f,ΔtC)分别表示测目标和测定标体过程中,测试场与雷达系统传递函数随时间变化产生的剩余背景误差所带来的测量不确定性。
因此,对于任何RCS测试场,均可通过式(44)来分析影响RCS测试不确定度的主要因素,并采取具体措施尽可能消除这些不确定因素,提高RCS测量精度。例如:
例-1:测量雷达系统漂移和测试场环境电特性变化的影响及其消除方法:
在式(44)中,
式中
为测量雷达系统(含天线)随时间和频率漂移引起的不确定性,通常天线增益变化影响很小,不确定性主要是发射机和接收机漂移引起的。
为测试场在测目标和测定标体时两者之间传播因子不一致引起的不确定性,主要是由于测试场电参数随时间和气象等条件变化引起的。
如果测试过程中通过测试流程设计或辅助测量手段可以导出传递函数Hm(f,Δt24),则在定标处理中用1/Hm(f,Δt24)进行补偿处理,测量雷达系统漂移和测试场环境电特性变化的影响得以消除。
例-2:背景相减剩余误差的影响及其消除方法
式(44)中剩余背景误差所带来的测量不确定性可表示如下:
式中,
为测目标路径传播因子引起的不确定性;
为测定标体路径传播因子变化引起的不确定性。
可见,背景抵消剩余误差受到雷达系统漂移HG(f,Δt34)(测目标支架背景与测目标期间的雷达系统漂移)、HG(f,Δt12)(测定标支架背景与测定标体期间的雷达系统漂移)、HFT(f,Δt34)(测目标支架背景与测目标期间的路径传播因子变化)、以及HFC(f,Δt12)(测定标支架背景与测定标体期间的路径传播因子变化)四者的影响。
同例-1中一样,若在测试中通过测试流程设计或辅助测量手段可导出HG(f,Δt12)·HFC(f,Δt12)和HG(f,Δt34)·HFT(f,Δt34),则在背景相减处理前可对背景数据进行补偿处理,从而消除剩余背景误差。
由此可见,本发明所提出的RCS测试定标处理数学模型完整表达了RCS测试和处理过程中影响测试不确定度的随时间变化和不随时间变化的各种因素。
假设定标区背景和目标区背景测量分别为和从背景测量到某次定标体和目标测量的过程中,随时间变化的测量系统传递函数变化为H(f)。测量系统漂移和测试场环境变化对测量回波带来的影响可由三个参数来表征:幅度漂移造成的增益变化A、随测试频率线性变化的相移和不随频率变化的固定相移φ0=ejβ。因此,总的传递函数为:
式中c为传播速度,f为测量雷达频率。
因此,如果不作幅相补偿,直接用已测量的背景数据进行背景相减和定标处理,由式(10)有: (55)
可见,由于传递函数H(f)的影响,将产生背景相减误差,进而影响RCS测量精度。传递函数H(f)越偏离1,背景相减误差越大。
为此,需要通过对测量数据处理求出H(f)的3个参数,即可完成幅度增益和相移特性补偿,达到提高背景相减效果,进而提高RCS测量精度的目的。其中:
1、具体的测量和幅相补偿处理流程如下:
(1)背景杂波测量:在未安装定标体和目标的条件下,对仅包含定标支架和金属目标支架的测试场背景进行宽带扫频测量,得到定标区背景回波Bmc(f)和目标区背景回波Bmt(f);
(2)定标体参考回波测量:安装标准定标体并进行宽带扫频测量,得到包含定标体真实回波和定标区背景回波的定标体测量参考回波:
Sc0(f)=C(f)+Bc0(f) (56)
由于定标体属于简单小目标,测中可以快速安装,故从背景杂波测量到定标体参考回波测量之间的系统漂移和环境变化可以忽略不计,因而此时近似有Bmc(f)≈Bc0(f);
(3)定标体和目标测量:安装目标并进行宽带扫频测量,得到包含定标体真实回波和定标区背景回波的定标体测量回波:
Sc(f)=C(f)+Bc(f) (57)
以及包含目标真实回波与目标区背景的目标测量回波:
St(f)=T(f)+Bt(f) (58)
Bc(f)=Bmc(f)·H(f) (59)
Bt(f)=Bmt(f)·H(f) (60)
注意到目标与定标体是同时测量的,故此时两者的传递函数相同。
(5)背景相减与定标:根据式(10)背景相减和定标方程有:
故可完成精确的背景相减和定标处理;由此可获得精确的目标RCS扫频测量σt(f)。
可见,上述处理的关键是要从测量数据中准确估计出H(f)的3个参数
2、传递函数参数估计方法:
(1)参考定标体测量回波的AR模型:对参考定标体的测量回波进行背景相减处理,并建立其AR模型M0(fn),
其中,n=0,…,N-1表示信号的样本点数;P、a0k和p0k分别代表模型的阶数、复幅度及极点;阶数P也代表散射中心的个数;复幅度a0k包含了信号的幅度及固定相位信息;极点p0k中包含着信号的线性相位信息,即散射中心的相对距离和随频率衰减(或增大)的情况。
注意:(a)由于定标体一般采用金属圆柱、金属球等简单体,其主要散射中心通常只有1~3个,故模型阶数可以低至2~3;(b)实际测试场设计中,一般要求测量中所选择的定标体其RCS高出定标区背景电平20dB以上,故此处一般也可不作背景相减处理;(c)对于定标体测量数据,即使此处不做背景相减处理,由于AR模型的阶数很低,反映了数据中主要散射中心的回波特性,而这些主要散射中心正是定标体本身的散射,故AR模型M0(fn)可以精确地反映定标体的回波特性。
(2)定标体测量回波的AR模型:对某次特定的目标与定标体连续测量,建立定标体测量回波的AR模型Mc(fn),
其中,n=0,…,N-1表示每组合成信号的样本点数;P、ack和pck分别代表模型的阶数、复幅度及极点。各参数的物理意义与参考定标体AR模型的类似。
模型阶数P选择与参考定标测量模型阶数完全相同。同样,AR模型Mc(fn)精确地反映了定标体在当前时刻的回波特性。
(3)传递函数参数提取
理想情况下,如果不存在系统漂移和环境影响,应该有M0(fn)=Mc(fn)。现由于非理想条件,M0(fn)与Mc(fn)不等,两者之间存在传递函数H(f)的影响,可表示为Mc(fn)=M0(fn)·H(f)。因此,可利用不同时间点测得的定标体数据建立AR模型M0(fn)与Mc(fn),并通过建立以下代价函数求得H(f)的3个参数:
(a)模型幅度参数估计算法
一旦确定AR模型的阶数,估计AR模型幅度可采用全局最小均方误差(TLS)方法,因为其中需要用到奇异值分解(SVD),所以也称之为SVD-TLS算法。研究表明,TLS算法是目前AR谱估计中抗噪性能最好的一种算法,已广泛应用于各个领域的信号分辨与估计中(参见:张贤达,现代信号处理,第二版,北京:清华大学出版社,2005)。
(b)模型极点估计算法
模型极点的估计一般可采用Barabell提出的root-MUSIC方法(参见:A.J.Barabell,Improving the resolution performance of eigenstructure based direction-finding algorithms,Proc.ICASS’83,336-339)或改进的root-MUSIC算法(参见:张贤达,现代信号处理,第二版,北京:清华大学出版社,2005)。
Claims (4)
1.一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,其特征在于:该方法的步骤如下:
步骤-1:t1时刻,测量包含定标支架的测试场背景回波SBC(f,t1);
步骤-2:t2时刻,安装定标体,并测量定标体回波SC(f,t2);
步骤-3:t3时刻,测量包含目标支架的测试场背景回波SBT(f,t3);
步骤-4:t4时刻,安装目标,并测量目标回波ST(f,t4);
步骤-5:利用定标体测量数据提取时变传递函数的3个幅度相位补偿参数、进行测量数据的幅度与相位校正、背景相量相减和目标RCS定标处理,具体处理方式如下:
式中,Ei(f)和Es(f)分别表示目标处雷达入射场和雷达天线处目标散射场;它同RCS之间的关系为
假定RCS测量满足高信噪比测量条件,即噪声的影响可忽略,若引入测量系统-测试场的时变传递函数对雷达接收回波的影响,则对于宽带扫频幅相测量,测目标和定标体时各自的雷达回波信号可分别表示为:
ST(f,t)=HT(f,t)·[T(f)+BT(f,t)] (12)
和
SC(f,t)=HC(f,t)·[C(f)+BC(f,t)] (13)
式中,HT(f,t)和HC(f,t)分别表示测目标和测定标体时测量系统-测试场的传递函数,是随时间慢变化的,反映了测量雷达系统漂移和测试场类型及其电参数特性随时间的变化特性;ST(f,t)和SC(f,t)分别表示雷达接收到的目标回波和定标体回波,两者均受到测量雷达-测试场时变传递函数的影响;T(f)和C(f)分别表示目标和定标体在给定姿态下的真实散射场,在给定姿态和雷达参数条件下不随时间而变化;BT(f,t)和BC(f,t)分别表示测目标和测定标体时的固有背景散射,两者均随时间t变化,表示由于目标/定标体同支架之间的耦合散射或背景状态的改变,造成放置目标/支架和空支架测量时,固有背景散射的变化,注意到上述散射信号均表示为复数相量;
这样,在外场异地定标RCS测量中,背景相减和RCS定标处理的过程数学上可表示为:
式中t1、t2、t3、t4分别表示测量定标背景、测定标体、测目标背景、测目标的时间段,SBT(f,t)和SBC(f,t)分别表示安放目标支架时测得的背景回波和安放定标体支架时测得的背景回波,两者可表示为:
SBT(f,t)=HT(f,t)·BT(f,t) (15)
SBC(f,t)=HC(f,t)·BC(f,t) (16)
这样,由式(12)~(16),有:
式中,ΔC(f,Δt12)和ΔT(f,Δt34)为经背景抵消处理后的剩余背景误差,
式(17)~(19)构成了异地定标RCS测量中测试与定标处理的基本数学模型。
2.根据权利要求1所述的一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,其特征在于:为了作进一步分析,假定:
精确的RCS测量和定标要求:
其中A0(f)为幅频特性,φ0(f)为相频特性;且
ΔT(f,Δt34)=0,ΔC(f,Δt12)=0 (23)
即剩余背景抵消误差为0。
3.根据权利要求2所述的一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,其特征在于:若测试场和定标与目标支架均经过仔细设计和优化,一般测试过程中定标区和目标区的背景是不随时间变化的,即有BC(f,t2)=BC(f,t1)=BC(f),BT(f,t4)=BT(f,t3)=BT(f),则此时式(18)和(19)变为:
因此,若要求剩余背景误差为0,则要求:
HT(f,t4)=HT(f,t3) (26)
HC(f,t2)=HC(f,t1) (27)
综合式(20)~(22)、(26)和(27),有:
AT(f,t4)=AT(f,t3),φT(f,t4)=φT(f,t3) (28)
AC(f,t2)=AC(f,t1),φC(f,t2)=φC(f,t1) (29)
式(28)~(30)构成了对精确的RCS测量、背景抵消和定标处理的基本要求,但是,由于真实测量雷达的参数总是随时间漂移的,也即上述三式的要求是不可能得到满足的,除非在测量过程中采用必要的辅助测量和信号处理方法对时变的传递函数幅度和相位作补偿处理,使得经过补偿处理后式(28)~(30)同时得到满足。
4.根据权利要求3所述的一种异地连续定标RCS测量中改进背景相减技术的信号处理方法,其特征在于:利用定标体在不同时间点的测量数据作为辅助测量,通过建立不同时间点定标体测量数据的自回归(AR)模型,并采用优化算法求得时变传递函数幅度与相位补偿的3个参数,进而完成幅度相位校正,使得校正处理后的测量数据满足式(28)~(30)所构成的对精确RCS测量、背景抵消和定标处理的基本要求,从而达到精确背景抵消和精确RCS定标的目的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210535816.5A CN102967855B (zh) | 2012-12-12 | 2012-12-12 | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210535816.5A CN102967855B (zh) | 2012-12-12 | 2012-12-12 | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102967855A true CN102967855A (zh) | 2013-03-13 |
CN102967855B CN102967855B (zh) | 2014-04-09 |
Family
ID=47798129
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210535816.5A Active CN102967855B (zh) | 2012-12-12 | 2012-12-12 | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102967855B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376446A (zh) * | 2013-07-09 | 2013-10-30 | 电子科技大学 | 一种rcs测试背景对消系统及方法 |
CN105182329A (zh) * | 2015-10-09 | 2015-12-23 | 上海无线电设备研究所 | 一种小双站角复合反射特性时域测量方法 |
CN105891795A (zh) * | 2016-05-23 | 2016-08-24 | 北京环境特性研究所 | 一种双站rcs测量定标方法 |
CN106443611A (zh) * | 2015-11-11 | 2017-02-22 | 北京航空航天大学 | 一种弱散射目标的rcs测量方法 |
CN106646383A (zh) * | 2016-09-28 | 2017-05-10 | 中国科学院电子学研究所 | 一种对月壤结构探测仪数据去除背景的方法 |
CN107102304A (zh) * | 2017-04-17 | 2017-08-29 | 北京环境特性研究所 | 基于地平场的非均匀照射不确定度分析方法 |
CN107153179A (zh) * | 2017-05-26 | 2017-09-12 | 中国电子科技集团公司第四十研究所 | 一种雷达目标rcs与散射中心同步测试方法 |
CN108717181A (zh) * | 2018-06-27 | 2018-10-30 | 成都飞机工业(集团)有限责任公司 | 一种新型室外场静态rcs前置联动标校测试方法 |
CN110441746A (zh) * | 2019-08-20 | 2019-11-12 | 北京环境特性研究所 | 一种时域门变换方法和装置 |
CN110441745B (zh) * | 2019-08-16 | 2021-04-30 | 北京环境特性研究所 | 一种基于宽带雷达俯视测量目标rcs的方法和系统 |
CN112964280A (zh) * | 2016-06-17 | 2021-06-15 | 亚德诺半导体集团 | 提供关于传递函数的不确定性的测量结果的设备和方法 |
CN113050061A (zh) * | 2021-03-30 | 2021-06-29 | 北京环境特性研究所 | 一种rcs测量背景电平对消处理方法及装置 |
CN113050055A (zh) * | 2021-03-25 | 2021-06-29 | 北京环境特性研究所 | 一种有源定标设备等效rcs标定方法及系统 |
CN115079112A (zh) * | 2022-07-21 | 2022-09-20 | 中国航发四川燃气涡轮研究院 | 一种航空发动机近地动态rcs试验测试系统及测试方法 |
CN116577785A (zh) * | 2023-07-12 | 2023-08-11 | 南京纳特通信电子有限公司 | 外场对地目标rcs成像与测量设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253376A (zh) * | 2011-04-14 | 2011-11-23 | 西北工业大学 | 一种基于二维微波成像的低散射共形天线rcs测试方法 |
CN102401893A (zh) * | 2010-09-10 | 2012-04-04 | 上海无线电设备研究所 | 一种近场目标对称极化rcs测试方法 |
CN102062857B (zh) * | 2010-12-10 | 2012-08-22 | 中国舰船研究设计中心 | 矩形平板近场雷达散射截面“平顶”效应预估方法 |
-
2012
- 2012-12-12 CN CN201210535816.5A patent/CN102967855B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102401893A (zh) * | 2010-09-10 | 2012-04-04 | 上海无线电设备研究所 | 一种近场目标对称极化rcs测试方法 |
CN102062857B (zh) * | 2010-12-10 | 2012-08-22 | 中国舰船研究设计中心 | 矩形平板近场雷达散射截面“平顶”效应预估方法 |
CN102253376A (zh) * | 2011-04-14 | 2011-11-23 | 西北工业大学 | 一种基于二维微波成像的低散射共形天线rcs测试方法 |
Non-Patent Citations (1)
Title |
---|
吴鹏飞 等: "地面平面场RCS测量异地定标误差分析", 《雷达学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376446A (zh) * | 2013-07-09 | 2013-10-30 | 电子科技大学 | 一种rcs测试背景对消系统及方法 |
CN105182329A (zh) * | 2015-10-09 | 2015-12-23 | 上海无线电设备研究所 | 一种小双站角复合反射特性时域测量方法 |
CN106443611A (zh) * | 2015-11-11 | 2017-02-22 | 北京航空航天大学 | 一种弱散射目标的rcs测量方法 |
CN106443611B (zh) * | 2015-11-11 | 2020-05-05 | 北京航空航天大学 | 一种弱散射目标的rcs测量方法 |
CN105891795A (zh) * | 2016-05-23 | 2016-08-24 | 北京环境特性研究所 | 一种双站rcs测量定标方法 |
CN105891795B (zh) * | 2016-05-23 | 2018-01-16 | 北京环境特性研究所 | 一种双站rcs测量定标方法 |
CN112964280B (zh) * | 2016-06-17 | 2023-11-24 | 亚德诺半导体国际无限责任公司 | 提供关于传递函数的不确定性的测量结果的设备和方法 |
CN112964280A (zh) * | 2016-06-17 | 2021-06-15 | 亚德诺半导体集团 | 提供关于传递函数的不确定性的测量结果的设备和方法 |
CN106646383B (zh) * | 2016-09-28 | 2019-09-06 | 中国科学院电子学研究所 | 一种对月壤结构探测仪数据去除背景的方法 |
CN106646383A (zh) * | 2016-09-28 | 2017-05-10 | 中国科学院电子学研究所 | 一种对月壤结构探测仪数据去除背景的方法 |
CN107102304B (zh) * | 2017-04-17 | 2019-06-25 | 北京环境特性研究所 | 基于地平场的非均匀照射不确定度分析方法 |
CN107102304A (zh) * | 2017-04-17 | 2017-08-29 | 北京环境特性研究所 | 基于地平场的非均匀照射不确定度分析方法 |
CN107153179B (zh) * | 2017-05-26 | 2019-11-22 | 中国电子科技集团公司第四十一研究所 | 一种雷达目标rcs与散射中心同步测试方法 |
CN107153179A (zh) * | 2017-05-26 | 2017-09-12 | 中国电子科技集团公司第四十研究所 | 一种雷达目标rcs与散射中心同步测试方法 |
CN108717181B (zh) * | 2018-06-27 | 2022-05-06 | 成都飞机工业(集团)有限责任公司 | 一种新型室外场静态rcs前置联动标校测试方法 |
CN108717181A (zh) * | 2018-06-27 | 2018-10-30 | 成都飞机工业(集团)有限责任公司 | 一种新型室外场静态rcs前置联动标校测试方法 |
CN110441745B (zh) * | 2019-08-16 | 2021-04-30 | 北京环境特性研究所 | 一种基于宽带雷达俯视测量目标rcs的方法和系统 |
CN110441746A (zh) * | 2019-08-20 | 2019-11-12 | 北京环境特性研究所 | 一种时域门变换方法和装置 |
CN113050055A (zh) * | 2021-03-25 | 2021-06-29 | 北京环境特性研究所 | 一种有源定标设备等效rcs标定方法及系统 |
CN113050055B (zh) * | 2021-03-25 | 2023-02-24 | 北京环境特性研究所 | 一种有源定标设备等效rcs标定方法及系统 |
CN113050061A (zh) * | 2021-03-30 | 2021-06-29 | 北京环境特性研究所 | 一种rcs测量背景电平对消处理方法及装置 |
CN115079112A (zh) * | 2022-07-21 | 2022-09-20 | 中国航发四川燃气涡轮研究院 | 一种航空发动机近地动态rcs试验测试系统及测试方法 |
CN116577785A (zh) * | 2023-07-12 | 2023-08-11 | 南京纳特通信电子有限公司 | 外场对地目标rcs成像与测量设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102967855B (zh) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102967855B (zh) | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 | |
CN102998665B (zh) | 目标雷达散射截面测量与定标处理的方法 | |
US9037414B1 (en) | Methods and apparatus for electromagnetic signal polarimetry sensing | |
CN104635222A (zh) | 一种低可探测目标rcs测量中的背景提取与抵消处理方法 | |
CN101738601B (zh) | 基于雷达近场回波功率谱特征的机车速度测量系统及方法 | |
CN105372691A (zh) | 一种模糊度固定的长基线卫星编队gnss相对定位方法 | |
CN103792531A (zh) | 基于天线阵的多目标微位移测量方法 | |
Mohanakumar et al. | Technical details of a novel wind profiler radar at 205 MHz | |
CN107300561A (zh) | 基于多遥感器联合探测的海洋盐度卫星 | |
CN101900692A (zh) | 大面积土壤湿度测量方法 | |
CN114397425B (zh) | 基于广义延拓逼近的gnss-ir土壤湿度反演方法 | |
Park et al. | Retracking considerations in spaceborne GNSS-R altimetry | |
CN101561499A (zh) | 一种单站多普勒测距定位方法 | |
CN105301570A (zh) | 一种机载顺轨干涉sar系统的外场定标方法 | |
CN105572635A (zh) | 基于最小二乘法的单站无源快速定位方法 | |
CN103760546A (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN104950307A (zh) | 一种机载三通道sar-gmti的精确定位方法 | |
Rao et al. | Advanced meteor radar installed at Tirupati: System details and comparison with different radars | |
CN111337548A (zh) | 高精度信噪比拟合模型及基于该模型的土壤湿度反演方法 | |
Kumar et al. | Initial results from SKiYMET meteor radar at Thumba (8.5 N, 77 E): 1. Comparison of wind measurements with MF spaced antenna radar system | |
Zhang et al. | An evaluation of the Chinese HY-2B satellite’s microwave scatterometer instrument | |
Yun et al. | CFOSAT Rotating Fan‐Beam Scatterometer Backscatter Measurement Processing | |
Fuhrmann et al. | Investigating GNSS multipath effects induced by co-located Radar Corner Reflectors | |
Lozinsky et al. | ICEBEAR-3D: A low elevation imaging radar using a non-uniform coplanar receiver array for E region observations | |
CN103792519A (zh) | 一种基于有源定标器的星载雷达时钟漂移在轨校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |