CN107728117A - 机载双基地sar一发两收杂波抑制方法 - Google Patents
机载双基地sar一发两收杂波抑制方法 Download PDFInfo
- Publication number
- CN107728117A CN107728117A CN201710858178.3A CN201710858178A CN107728117A CN 107728117 A CN107728117 A CN 107728117A CN 201710858178 A CN201710858178 A CN 201710858178A CN 107728117 A CN107728117 A CN 107728117A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- eta
- tau
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种机载双基地SAR一发两收杂波抑制方法,本发明的方法利用—发两收的天线配置,通过在回波域去除与方位时间有关的距离历史的一次和二次项,消除发射站的影响,在图像域补偿固定相位,消除由于两接收天线水平排布位置不同导致的固定距离历史误差,从而使得两通道对静止目标的回波完全相同,所以在相减之后杂波被抑制,而运动目标由于其速度的影响在相减后仍有剩余,从而检测出运动目标,解决了传统DPCA方法不能够适用于双基地的问题。本发明的优点是该方法与现有的其他双基雷达杂波抑制方法STAP相比,处理更简单,运算量更小,且不需要考虑双基雷达的杂波分布模型。
Description
技术领域
本发明属于雷达技术领域,具体涉及机载双基地雷达对于地面运动目标的检测技术。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种全天候、高分辨率微波遥感成像雷达,它利用雷达天线和目标区域间的相对运动来获得空间的高分辨。目前在地形测绘,灾区探测等领域得到广泛应用。
偏移相位中心天线(Displaced Phase Center Antenna,DPCA)技术最初是为补偿雷达平台的运动效应而提出来的,其实质是补偿由于雷达平台运动造成的杂波谱展宽。传统的DPCA方法采用两个移位相位中心,由前置天线发射和接收第一个脉冲,在后置天线到达同一空间位置时,发射和接收下一个脉冲,即前后两个脉冲是在空间同一位置发射和接收的,从而使天线相位中心在相继的发射脉冲里相对于地面静止不动,这样采用两脉冲对消就可以抑制杂波。
而对于双基地系统,尤其是收发平台同时运动的双基地系统,在一发两收的模式下,由于发射通道与接收通道位于不同的平台上,固定时间间隔后发射通道位置将发生偏移,从而无法满足固定时间间隔内收-发通道等效相位中心位于相同的位置的条件,造成传统的DPCA方法无法适用。
在文献:“合成孔径雷达运动目标检测和成像研究[D]”郑明洁,中国科学院电子学研究所,2003中研究了针对单基地的基于DPCA技术的多天线SAR动目标检测方法,通过时移和相位补偿,使相邻通道的相位中心重合;文献:“基于二天线的双基地DPCA技术[J]”陈娟,王盛利,电子与信息学报,2007,29[7]:1687-1690中研究了针对双基地的二天线DPCA杂波抑制方法,但是地面接收系统是固定不动的,在固定的时间间隔后仍然能够满足两通道回波的距离历史相同的条件;文献:“GEO星机双基地SAR动目标检测与成像方法研究[D]”,廉濛,哈尔滨工业大学,2013中研究了GEO星机双基地SAR的动目标检测方法,由于发射平台运动,无法保持等效相位中心位置不变,通过计算出针对某一目标的距离历史在一个脉冲发射周期内的变化,对其中一路信号补偿这一差异来达到等效相位中心不变的目的,但却需要在不同的方位时间上分别进行补偿。以上文献都没有提出有效的适用于双基地一发两收构型的DPCA运动目标检测方法。
发明内容
本发明的目的是针对现有技术存在的缺陷,提出一种机载双基地SAR一发两收杂波抑制方法,以克服传统DPAC方法在机载双基地构型下等效相位中心无法重合、杂波无法抑制的问题。
本发明的技术方案为:一种机载双基地SAR一发两收杂波抑制方法,具体包括如下步骤:
步骤一:建立机载双基地SAR一发两收的空间几何结构,并完成参数初始化;发射站系统和接收站系统分置于两个不同的机载平台上,发射站有一个发射天线,接收站包含两个接收天线,两个接收天线沿平台飞行方向水平排布,工作时,由发射天线发射信号,两个接收天线同时接收回波信号;
在直角坐标系中,O为坐标原点,设P(x,y,0)为地面目标,其中,x为该地面目标的X轴坐标,y为该地面目标的Y轴坐标;发射天线的位置坐标为(xT,yT,HT),其中,xT为发射天线的X轴坐标,yT为发射天线的Y轴坐标,HT为发射天线的Z轴坐标,第一接收天线的坐标为(xR1,yR1,HR1),其中,xR1为第一接收天线的X轴坐标,yR1为第一接收天线的Y轴坐标,HR1为第一接收天线的Z轴坐标;第二接收天线的坐标为(xR2,yR2,HR2),其中,xR2为第二接收天线的X轴坐标,yR2为第二接收天线的Y轴坐标,HR2为第二接收天线的Z轴坐标;发射平台的飞行速度为vT,接收平台的飞行速度为vR,两平台飞行方向都平行于Y轴;两接收天线的基线长度为d,RT,c为发射天线在波束中心时刻到目标的斜距,RR1,c为第二接收天线在波束中心时刻到目标的斜距,RR2,c为第二接收天线在波束中心时刻到目标的斜距,其中发射天线的斜视角θT为RT,c与Z轴方向的夹角,第一接收天线的斜视角θR1为RR1,c与Z轴方向的夹角,第二接收天线的斜视角θR2为RR2,c与Z轴方向的夹角。
步骤二:获取机载双基地SAR两通道到地面目标P(x,y,0)的回波,距离压缩后的表达式为:
其中,A为常数,Kr为距离向脉冲的调频率,Tp为距离向信号时宽,τ为距离向快时间变量,η为方位向慢时间变量,f0表示载频,c为光速,R1(η)为目标到第一接收天线的双基距离历史,R2(η)为目标到第二接收天线的双基距离历史;
步骤三:将双基距离历史在波束中心时刻η=0泰勒展开,则两通道的回波信号可以表示为:
步骤四:分别对以上两通道信号乘以两个补偿函数,函数分别为:
补偿后得到两通道的回波为:
步骤五:对两通道信号分别进行距离走动校正,并将通道2的信号方位时间移动得到:
步骤六:通过方位向傅里叶变换分别对两通道的信号成像,得到:
其中,Ts为合成孔径时间,fη为方位向多普勒频率。
步骤七:对通道2的图像进行相位补偿,补偿函数如下:
得到:
步骤八:对两通道的成像结果进行相减,得到:
I12(τ,fa)=I1(τ,fa)-I′2(τ,fa)=0
从上式可以看出,经过处理后两通道对静止目标的成像结果相同,所以在对消后结果为零,即实现了静止杂波目标的抑制。
本发明的有益效果:本发明的方法利用—发两收的天线配置,通过在回波域去除与方位时间有关的距离历史的一次和二次项,消除发射站的影响,在图像域补偿固定相位,消除由于两接收天线水平排布位置不同导致的固定距离历史误差,从而使得两通道对静止目标的回波完全相同,所以在相减之后杂波被抑制,而运动目标由于其速度的影响在相减后仍有剩余,从而检测出运动目标,解决了传统DPCA方法不能够适用于双基地的问题。本发明的优点是该方法与现有的其他双基雷达杂波抑制方法STAP相比,处理更简单,运算量更小,且不需要考虑双基雷达的杂波分布模型。本发明可以应用于机载双基地雷达动目标检测等领域。
附图说明
图1是本发明的机载双基地SAR一发两收杂波抑制方法的流程框图。
图2是本发明提供的机载双基SAR两通道空间几何结构。
图3是本发明具体实施方式采用的几何结构参数表。
图4是本发明具体实施方式中对消处理前运动目标与静止目标的成像结果。
图5是本发明具体实施方式中对消处理后运动目标与静止目标的成像结果。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2010上验证正确。
本发明的方法解决了传统DPCA方法由于收发分置、发射通道与接收通道位于不同的平台上,无法满足固定时间间隔内收-发通道等效相位中心对齐,杂波难以抑制的问题。它的特点是利用双基SAR的一发两收模式,建立机载双基SAR两通道的回波模型,推导两个通道距离历史差距的表达式,通过分别在回波域和图像域的相位补偿处理消除这一差距,从而使两个通道静止杂波的距离历史相同,有效解决了传统DPCA方法无法适用的问题。具体流程图如图1所示,包括如下步骤:
步骤一:建立机载双基SAR两通道的空间几何结构,并完成参数初始化;
如图2所示,在直角坐标系中,设发射站的位置坐标为(-10km,-13km,12km);接收站2的坐标为(0km,-10km,10km),则接收站1的坐标为(0km,-10km-d,10km);发射站的飞行速度为200m/s,飞行方向与Y轴平行;两接收站的飞行速度都为200m/s,飞行方向与Y轴平行;接收站通道间距为0.2m。具体参数如图3所示。
步骤二:设置两个仿真点,一个为静止目标,一个为运动目标。其中,静止目标位置为(300m,100m,0),运动目标位置为(500,80m,0),运动目标距离向速度为Vx=2m/s,方位向速度为Vy=1m/s。
步骤三:将双基距离历史在波束中心时刻η=0泰勒展开,则两通道的回波信号可以表示为:
步骤四:分别对两通道信号乘以两个补偿函数,函数分别为:
补偿后得到两通道的回波为:
步骤五:对两通道信号分别进行距离走动校正,并将通道2的信号方位时间移动得到:
步骤六:通过方位向傅里叶变换分别对两通道的信号成像,得到:
步骤七:对通道2的图像进行相位补偿,补偿函数如下:
得到:
步骤八:对两通道的成像结果进行相减,得到:
I12(τ,fa)=I1(τ,fa)-I′2(τ,fa)=0
经过上述步骤处理,即可完成对两通道回波的DPCA处理,对消处理前运动目标与静止目标的成像结果如图4所示,对消处理后运动目标与静止目标的成像结果如图5所示,从图中可以看出,在对消处理前,静止目标与运动目标的成像结果无法分辨,经过对消处理后,静止目标被消除,剩余了运动目标,且杂波抑制效果达到了22dB,从而验证了本方法的有效性。
Claims (1)
1.一种机载双基地SAR一发两收杂波抑制方法,具体包括如下步骤:
步骤一:建立机载双基地SAR一发两收的空间几何结构,并完成参数初始化;发射站系统和接收站系统分置于两个不同的机载平台上,发射站有一个发射天线,接收站包含两个接收天线,两个接收天线沿平台飞行方向水平排布,工作时,由发射天线发射信号,两个接收天线同时接收回波信号;
在直角坐标系中,O为坐标原点,设P(x,y,0)为地面目标,其中,x为该地面目标的X轴坐标,y为该地面目标的Y轴坐标;发射天线的位置坐标为(xT,yT,HT),其中,xT为发射天线的X轴坐标,yT为发射天线的Y轴坐标,HT为发射天线的Z轴坐标,第一接收天线的坐标为(xR1,yR1,HR1),其中,xR1为第一接收天线的X轴坐标,yR1为第一接收天线的Y轴坐标,HR1为第一接收天线的Z轴坐标;第二接收天线的坐标为(xR2,yR2,HR2),其中,xR2为第二接收天线的X轴坐标,yR2为第二接收天线的Y轴坐标,HR2为第二接收天线的Z轴坐标;发射平台的飞行速度为vT,接收平台的飞行速度为vR,两平台飞行方向都平行于Y轴;两接收天线的基线长度为d,RT,c为发射天线在波束中心时刻到目标的斜距,RR1,c为第二接收天线在波束中心时刻到目标的斜距,RR2,c为第二接收天线在波束中心时刻到目标的斜距,其中发射天线的斜视角θT为RT,c与Z轴方向的夹角,第一接收天线的斜视角θR1为RR1,c与Z轴方向的夹角,第二接收天线的斜视角θR2为RR2,c与Z轴方向的夹角。
步骤二:获取机载双基地SAR两通道到地面目标P(x,y,0)的回波,距离压缩后的表达式为:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>&lsqb;</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>&lsqb;</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
</mrow>
其中,A为常数,Kr为距离向脉冲的调频率,Tp为距离向信号时宽,τ为距离向快时间变量,η为方位向慢时间变量,f0表示载频,c为光速,R1(η)为目标到第一接收天线的双基距离历史,R2(η)为目标到第二接收天线的双基距离历史;
步骤三:将双基距离历史在波束中心时刻η=0泰勒展开,则两通道的回波信号可以表示为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>T</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>R</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&eta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>T</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>R</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&eta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<mi>&pi;</mi>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
<mo>}</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<mi>&pi;</mi>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
<mo>}</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤四:分别对以上两通道信号乘以两个补偿函数,函数分别为:
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mfrac>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>T</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>R</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&eta;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mfrac>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>T</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>R</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&eta;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>R</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
补偿后得到两通道的回波为:
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
步骤五:对两通道信号分别进行距离走动校正,并将通道2的信号方位时间移动得到:
<mrow>
<msubsup>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>&tau;</mi>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
<mrow>
<msubsup>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>2</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>&tau;</mi>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
步骤六:通过方位向傅里叶变换分别对两通道的信号成像,得到:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;f</mi>
<mi>&eta;</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mo>}</mo>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>&tau;</mi>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>s</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>(</mo>
<mrow>
<mi>&tau;</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>&eta;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;f</mi>
<mi>&eta;</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mo>}</mo>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>&tau;</mi>
<mo>}</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Ts为合成孔径时间,fη为方位向多普勒频率。
步骤七:对通道2的图像进行相位补偿,补偿函数如下:
<mrow>
<mi>G</mi>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
得到:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>I</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>G</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>A</mi>
<mo>|</mo>
<msub>
<mi>K</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;f</mi>
<mi>&eta;</mi>
</msub>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mo>}</mo>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mo>{</mo>
<msub>
<mi>&pi;K</mi>
<mi>r</mi>
</msub>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mi>&tau;</mi>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤八:对两通道的成像结果进行相减,得到:
I12(τ,fa)=I1(τ,fa)-I2′(τ,fa)=0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710858178.3A CN107728117B (zh) | 2017-09-21 | 2017-09-21 | 机载双基地sar一发两收杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710858178.3A CN107728117B (zh) | 2017-09-21 | 2017-09-21 | 机载双基地sar一发两收杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107728117A true CN107728117A (zh) | 2018-02-23 |
CN107728117B CN107728117B (zh) | 2020-08-04 |
Family
ID=61207774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710858178.3A Active CN107728117B (zh) | 2017-09-21 | 2017-09-21 | 机载双基地sar一发两收杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107728117B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710117A (zh) * | 2018-08-24 | 2018-10-26 | 电子科技大学 | 一种合成孔径雷达杂波对消方法 |
CN110109113A (zh) * | 2019-05-07 | 2019-08-09 | 电子科技大学 | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 |
CN110133654A (zh) * | 2019-05-31 | 2019-08-16 | 电子科技大学 | 一种高轨卫星sar动目标检测方法 |
CN110133646A (zh) * | 2019-04-19 | 2019-08-16 | 电子科技大学 | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 |
CN110568434A (zh) * | 2019-07-15 | 2019-12-13 | 电子科技大学 | 一种多通道匀加速轨迹sar动目标二维速度估计方法 |
CN111707996A (zh) * | 2020-05-26 | 2020-09-25 | 北京理工大学重庆创新中心 | 基于改进grft-stap的geo星机sar动目标检测方法 |
CN112946650A (zh) * | 2021-04-15 | 2021-06-11 | 中山大学 | 一站固定双站低频超宽带sar运动目标检测与成像方法 |
CN113435299A (zh) * | 2021-06-23 | 2021-09-24 | 电子科技大学 | 基于空-时匹配的双基前视sar杂波抑制方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5113193A (en) * | 1985-11-12 | 1992-05-12 | The United States Of America As Represented By The Secretary Of The Air Force | Autonomous synchronization of a bistatic synthetic aperture radar (SAR) system |
CN101458334A (zh) * | 2007-12-14 | 2009-06-17 | 电子科技大学 | 一种双基地合成孔径雷达成像的运动补偿方法 |
CN102707283A (zh) * | 2012-06-15 | 2012-10-03 | 电子科技大学 | 一种固定站双基地合成孔径雷达成像方法 |
CN103207388A (zh) * | 2013-03-26 | 2013-07-17 | 中国科学院电子学研究所 | 一种斜视条件下的机载干涉sar定标方法 |
CN103308914A (zh) * | 2012-12-06 | 2013-09-18 | 北京理工大学 | 一站固定双基地干涉sar处理方法 |
CN103728617A (zh) * | 2014-01-13 | 2014-04-16 | 电子科技大学 | 双基地合成孔径雷达时域快速成像方法 |
CN104678361A (zh) * | 2015-02-03 | 2015-06-03 | 电子科技大学 | 机载双基地mimo雷达杂波模型建模方法 |
CN105182340A (zh) * | 2015-09-28 | 2015-12-23 | 电子科技大学 | 一种双基地前视sar运动补偿方法 |
CN106990396A (zh) * | 2017-05-25 | 2017-07-28 | 电子科技大学 | 一种双基地合成孔径雷达定位误差计算方法 |
-
2017
- 2017-09-21 CN CN201710858178.3A patent/CN107728117B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5113193A (en) * | 1985-11-12 | 1992-05-12 | The United States Of America As Represented By The Secretary Of The Air Force | Autonomous synchronization of a bistatic synthetic aperture radar (SAR) system |
CN101458334A (zh) * | 2007-12-14 | 2009-06-17 | 电子科技大学 | 一种双基地合成孔径雷达成像的运动补偿方法 |
CN102707283A (zh) * | 2012-06-15 | 2012-10-03 | 电子科技大学 | 一种固定站双基地合成孔径雷达成像方法 |
CN103308914A (zh) * | 2012-12-06 | 2013-09-18 | 北京理工大学 | 一站固定双基地干涉sar处理方法 |
CN103207388A (zh) * | 2013-03-26 | 2013-07-17 | 中国科学院电子学研究所 | 一种斜视条件下的机载干涉sar定标方法 |
CN103728617A (zh) * | 2014-01-13 | 2014-04-16 | 电子科技大学 | 双基地合成孔径雷达时域快速成像方法 |
CN104678361A (zh) * | 2015-02-03 | 2015-06-03 | 电子科技大学 | 机载双基地mimo雷达杂波模型建模方法 |
CN105182340A (zh) * | 2015-09-28 | 2015-12-23 | 电子科技大学 | 一种双基地前视sar运动补偿方法 |
CN106990396A (zh) * | 2017-05-25 | 2017-07-28 | 电子科技大学 | 一种双基地合成孔径雷达定位误差计算方法 |
Non-Patent Citations (2)
Title |
---|
ZHONGYU LI ET AL.: "Ground-Moving Target Imaging and Velocity Estimation Based on Mismatched Compression for Bistatic Forward-Looking SAR", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
黄钰林 等: "机载双站SAR分辨率特性分析", 《电波科学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710117A (zh) * | 2018-08-24 | 2018-10-26 | 电子科技大学 | 一种合成孔径雷达杂波对消方法 |
CN108710117B (zh) * | 2018-08-24 | 2022-07-15 | 电子科技大学 | 一种合成孔径雷达杂波对消方法 |
CN110133646A (zh) * | 2019-04-19 | 2019-08-16 | 电子科技大学 | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 |
CN110109113A (zh) * | 2019-05-07 | 2019-08-09 | 电子科技大学 | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 |
CN110133654A (zh) * | 2019-05-31 | 2019-08-16 | 电子科技大学 | 一种高轨卫星sar动目标检测方法 |
CN110568434B (zh) * | 2019-07-15 | 2022-06-03 | 电子科技大学 | 一种多通道匀加速sar动目标二维速度估计方法 |
CN110568434A (zh) * | 2019-07-15 | 2019-12-13 | 电子科技大学 | 一种多通道匀加速轨迹sar动目标二维速度估计方法 |
CN111707996A (zh) * | 2020-05-26 | 2020-09-25 | 北京理工大学重庆创新中心 | 基于改进grft-stap的geo星机sar动目标检测方法 |
CN111707996B (zh) * | 2020-05-26 | 2023-04-07 | 北京理工大学重庆创新中心 | 基于改进grft-stap的geo星机sar动目标检测方法 |
CN112946650A (zh) * | 2021-04-15 | 2021-06-11 | 中山大学 | 一站固定双站低频超宽带sar运动目标检测与成像方法 |
CN112946650B (zh) * | 2021-04-15 | 2023-08-15 | 中山大学 | 一站固定双站低频超宽带sar运动目标检测与成像方法 |
CN113435299A (zh) * | 2021-06-23 | 2021-09-24 | 电子科技大学 | 基于空-时匹配的双基前视sar杂波抑制方法 |
CN113435299B (zh) * | 2021-06-23 | 2022-05-13 | 电子科技大学 | 基于空-时匹配的双基前视sar杂波抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107728117B (zh) | 2020-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107728117A (zh) | 机载双基地sar一发两收杂波抑制方法 | |
CN108761404B (zh) | 一种基于二次相位函数参数估计及补偿的改进算法 | |
CN107966688A (zh) | 基于相位干涉技术的宽带雷达目标速度解模糊方法 | |
US8134493B2 (en) | System and method for precision geolocation utilizing multiple sensing modalities | |
CN105158748B (zh) | 一种高速目标多通道补偿聚焦与tbd混合积累检测方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN105891828B (zh) | 一种机载cssar雷达动目标的检测方法 | |
CN106443671A (zh) | 基于调频连续波的sar雷达动目标检测与成像方法 | |
CN105403887B (zh) | 基于ins的机载sar实时运动补偿方法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
CN104316903A (zh) | 一种三站时差定位性能试验评估方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
JP4899937B2 (ja) | Fmcw方式合成開口レーダ、偏流角検出方法、プログラムおよび記憶媒体 | |
CN106597437A (zh) | 高频机载调频连续波sar运动补偿及成像处理方法 | |
CN105487074B (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN107942314A (zh) | 基于lasso特征提取的多普勒穿墙雷达定位方法 | |
CN104749570A (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN107271997B (zh) | 机载多通道cssar地面运动目标运动参数估计方法 | |
CN106526583A (zh) | 一种基于天线方向图信息的地面运动目标定位方法 | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 | |
CN104950307A (zh) | 一种机载三通道sar-gmti的精确定位方法 | |
CN106483495A (zh) | 一种室内运动标签定位和测速方法 | |
CN101943753B (zh) | 一种微波雷达波束中心向速度的测量方法 | |
CN106526553A (zh) | 一种通用的sar卫星方位模糊度性能精确分析方法 | |
CN109521418B (zh) | 基于干涉场的地基雷达测角方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |