Loading [MathJax]/jax/output/SVG/jax.js

基于凸优化方法的认知雷达低峰均比波形设计

郝天铎 崔琛 龚阳 孙从易

李航舰, 王宇, 邓云凯, 王伟, 张衡. 分布式星载SAR系统时间同步和波束同步误差分析[J]. 雷达学报, 2018, 7(2): 244-253. doi: 10.12000/JR17028
引用本文: 郝天铎, 崔琛, 龚阳, 孙从易. 基于凸优化方法的认知雷达低峰均比波形设计[J]. 雷达学报, 2018, 7(4): 498-506. doi: 10.12000/JR18002
Li Hangjian, Wang Robert, Deng Yunkai, Wang Wei, Zhang Heng. Analysis of Time and Beam Synchronization Errors for Distributed Spaceborne SAR System[J]. Journal of Radars, 2018, 7(2): 244-253. doi: 10.12000/JR17028
Citation: Hao Tianduo, Cui Chen, Gong Yang, Sun Congyi. Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization[J]. Journal of Radars, 2018, 7(4): 498-506. doi: 10.12000/JR18002

基于凸优化方法的认知雷达低峰均比波形设计

DOI: 10.12000/JR18002
基金项目: 国家部委基金
详细信息
    作者简介:

    郝天铎(1989–),男,内蒙古呼和浩特人,国防科技大学电子对抗学院博士生,主要研究方向为认知雷达信号处理和最优化理论。E-mail: haotianduo0423@126.com

    崔 琛(1962–),男,河北易县人,国防科技大学电子对抗学院教授,博士生导师,主要研究方向为雷达信号处理及雷达对抗技术。E-mail: kycuichen@163.com

    龚 阳(1992–),男,湖北随州人,国防科技大学电子对抗学院博士生,主要研究方向为雷达信号处理。E-mail: 13156527915@163.com

    孙从易(1992–),男,山东青岛人,国防科技大学电子对抗学院硕士生,现就职于96630部队,主要研究方向为MIMO雷达信号处理。E-mail: 1337325128@163.com

    通讯作者:

    郝天铎   haotianduo0423@126.com

Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization

Funds: The National Ministries Foundation
  • 摘要: 为了提高雷达发射波形的检测性能,同时使发射机发挥其最大效能,以发射波形的低峰均比(PAR)为约束条件,该文提出了一种信号相关杂波背景下的认知雷达发射波形和接收机滤波器联合优化方法。首先,面向距离扩展目标检测问题,构建关于雷达输出信干噪比(SINR)的优化模型;然后将该模型转化为Rayleigh商形式,给出了接收机权值的解析表达式;在此基础上,通过半正定松弛,将关于发射波形半正定矩阵的非凸问题转化为凸问题,求得发射波形的最优矩阵解;最后,将秩1近似法和最近邻方法相结合,从最优矩阵解中提取出发射波形的最优向量解。该方法在给定PAR取值范围内可使波形的输出SINR达到最大,PAR=2时波形的SINR值与能量约束下优化波形的SINR值相同,并且比PAR=1时所得波形高出约0.5 dB。仿真结果验证了所提方法的有效性。

     

  • 合成孔径雷达(Synthetic Aperture Radar, SAR)可以实现全天时、全天候的高分辨率成像[1],星载SAR系统为全球地理遥感、地形测绘、形变监测等应用提供了丰富的高分辨数据。分布式星载SAR系统是一种特殊模式的合成孔径雷达系统,雷达发射机和接收机位于不同平台,这种收发分置的特点带来了很多优势,比如灵活的基线配置,多样的接收模式,可以提供更精细的地物分类信息等,具有单基情况无可比拟的优势。但是收发分置的特点也带来了很多新的挑战,尤其是时间同步、相位同步和波束同步等三大同步问题[2,3]。由于收发雷达不共频率源,导致了时间同步和相位同步问题[4],时间同步问题是指,发射雷达和接收雷达使用不同的晶振产生脉冲信号,导致了触发脉冲不一致,进而导致信号发射时刻与接收窗起始时刻不一致,使距离向的录取回波减少,时间同步误差会在回波中引入相位误差,并且晶振偏差会随时间累积,对回波数据产生较大影响,严重降低雷达图像质量[5,6]。为达到时间同步的要求,Tandem-X SAR系统采用频率一致性指标为10–8量级的晶振,并采用跳跃PRI的方法,确保600 s内,时间同步误差小于5 μs[7]。相位同步问题是指,由于晶振频率的差异,导致了发射信号调制频率和接收信号解调频率不一致,造成解调回波中存在沿方位时间变化的相位误差,该误差在成像后传递到图像域,导致方位向位置偏移、主瓣展宽、积分旁瓣比恶化等问题,并且相位误差最终会传递到干涉相位,造成干涉测高精度损失。TanDEM-X SAR系统采用脉冲交替对传的同步方案解决相位同步的问题。波束同步问题是指,由于收发天线的控制精度不同,使得收发天线无法在同一时刻完全覆盖地面同一区域,导致双基合成天线方向图的增益降低,同时接收天线无法获得全部的目标散射信号,造成有效观测幅宽损失。为保证系统系能,TanDEM-X SAR系统的波束指向精度控制在0.01°,最差情况下导致的天线增益损失小于0.1 dB。

    为了消除相位同步带来的影响,可以建立专门的同步链路,文献[8]提出了卫星间连续波对传的同步方法,文献[9]和文献[10]进一步提出了脉冲交替对传的方法,文献[11]提出一种基于信号交换的调频连续波合成孔径雷达相位同步方法。无论采用哪种同步链路方式,目的都是获得相位误差信息,最后在进行数据处理时提取并补偿相位误差,这也是解决分布式卫星SAR相位同步的优选方法[12]。文献[13]建立了一种时间同步误差模型,分析了时间同步误差对成像和干涉测高性能的影响,但是没有分析时间同步和相位同步的联系。时间同步和相位同步来自同一误差源,都是由于收发平台采用了独立的晶振系统来校准系统时钟[14],所以两者之间存在一定的关系。波束同步的研究主要集中在波束同步控制方法,以及波束同步误差对雷达成像性能的影响,文献[15]分析了几种可行的飞行控制策略,同时分析了波束同步对收发平台和姿态控制精度的要求,文献[16]分析了系统波束控制精度需要满足的要求。但是,在不同的卫星姿态控制精度和波束指向精度情况下,波束同步误差对分布式星载SAR系统的综合影响还需要进一步分析。

    本文结构安排如下。第2节介绍了时间同步误差模型和波束同步误差模型。第3节基于误差模型分析了时间同步和相位同步导致的相位误差之间的关系,并分析晶振稳定度需要满足的条件。第4节分析了波束同步误差对合成天线方向图、双基图像信噪比、多普勒去相干以及波束覆盖偏差的具体影响。第5节,进行仿真验证,为分布式SAR系统设计提供参考。

    收发雷达之间的时间同步误差模型可以表示为常数项、线性项和随机项3部分[13],随机项由频率源的短期稳定性决定,由于随机的特点,该项影响结果不具备积累性,一般只有在频率源稳定度较低或者所需相干积累时间较长时有较大影响[1720],因此该模型可以表达为:

    ε(t)=Δt+βt (1)

    其中, $\Delta t$ 为常数项,表示初始状态主辅星间雷达累积的定时信号之差; $\beta t$ 为线性项,由主辅星晶振频率偏差导致, $\beta $ 可以表示为:

    β=ΔPRI1ΔPRI2PRI=Δf1f1Δf2f2Δf1Δf2f0 (2)

    其中,f0表示标称频率, $\Delta {f_1}$ $\Delta {f_2}$ 分别表示发射雷达和接收雷达的频率偏移。当存在同步误差时,在以发射信号为参考的时间坐标系中,不同脉冲的回波在接收雷达采样窗中的位置会发生偏移[21]

    卫星姿态误差会影响卫星和地面目标的波束传播距离,同时影响回波多普勒频率,如图1所示,X轴从星下点指向卫星方向,Y轴指向方位向方向,Z轴指向距离向方向,卫星姿态误差可以分解为偏航角、俯仰角和横滚角,其中Rs表示斜距,R0表示最近斜距,P0表示姿态导引后期望的波束中心地面落点,P表示波束中心实际落点, ${\theta _{\rm{y}}}$ 表示偏航角, ${\theta _{\rm{p}}}$ 表示俯仰角, ${\theta _{\rm{r}}}$ 表示横滚角。可见,偏航角和俯仰角分别导致波束中心线绕X轴和Z轴旋转,进而偏离零多普勒面,均会导致斜视角。横滚角导致波束中心线绕Y轴旋转,所以横滚角引起波束在距离向移动,在距离向产生波束指向误差。

    图  1  卫星姿态对系统的影响
    Figure  1.  Impacts of yaw angle, pitch angle and roll angle

    根据几何关系可以得到,偏航角和俯仰角引起的斜视角分别为:

    θsy=arctan(tanθysinθl) (3)
    θsp=arctan(tanθpcosθl) (4)

    其中, ${\theta _{\rm{l}}}$ 为下视角。横滚角引起的距离向波束指向误差为:

    θsr=θr (5)

    主星发射的线性调频脉冲信号可以表示为:

    s(τ,t)=rect(τTp)exp(jπkτ2)exp(jπf0τ) (6)

    其中, ${\mathop{\rm rect}\nolimits} \left[ \cdot \right]$ 为矩形窗函数, $\tau $ 为距离向快时间,t为方位向慢时间,Tp为脉冲宽度,f0为发射信号载频,k为调频率。将时间同步误差模型代入发射信号,经辅星雷达接收解调后含有时间同步误差的回波信号为:

    sD(τ,t)=rect(τRt/cε(t)Tp)rect(ttc)exp(jπk(τRt/cε(t))2)exp(j2πRt/λ)exp(j2πf0ε(t)) (7)

    其中,c为光速, $\lambda $ 为波长,tc为波束中心偏离时间,Rt为雷达和卫星之间的波束传播距离,

    R(t)=R12+V21(tt01)2+R22+V22(tt02)2 (8)

    其中,R1R2分别表示收发最近斜距,t01t02分别表示波束中心穿越时刻,V1V2分别表示收发卫星的速度。为简化分析,下文中,假设两星同时穿越波束中心时刻且以该时刻为方位时间零点,即t01=t02=0。

    以主星发射信号,辅星接收信号为例,采用距离多普勒聚焦算法,推导时间同步误差对SAR成像的影响,假设主星雷达晶振为标准频率源,辅星雷达晶振含有误差项,即时间同步误差全部表现在辅星获得的图像上。

    辅星雷达只存在固定时间误差的情况下,即

    ε(t)=Δt (9)

    成像处理之后的辅星图像将存在时间同步误差项,具体可以表示为:

    σ(τ,t)=sinc(B(τ(R1+R2)/cΔt))sinc(Bdt)exp(j2πR1+R2λ)exp(j2πfact)exp(j2πf0Δt) (10)

    其中,Bd是多普勒带宽,fac是多普勒中心频率。由式(10)可以得到,固定时间同步误差在辅星图像中引入距离向主瓣偏移和固定相位误差,固定相位误差在经过干涉处理之后会产生常量干涉相位误差。主瓣偏移量和干涉相位误差分别为:

    Δr=Δtc2 (11)
    φe=2πf0Δt (12)

    常量相位误差会影响绝对测高精度。

    辅星雷达只存在线性时间误差的情况下,即

    ε(t)=βt (13)

    对式(7)中带有线性时间误差的相位项进一步推导可以得到

    φlinear=2πβf0t=2πΔf1Δf2f0f0t=2π(Δf1Δf2)t (14)

    由于晶振频率存在偏移,t时刻主星的晶振频率为 ${f_1} = {f_0} + \Delta {f_1}$ ,辅星的晶振频率为 ${f_2} = {f_0} + \Delta {f_2}$ ,根据文献[9]中对脉冲对传相位同步方案的分析,辅星需要补偿的相位偏差为:

    φc(t)=2π(Δf2Δf1)t+π(Δf2Δf1)(τ12+τsys)πfdτsys+12(nφ2(t+τsys)+nφ2(t+τ)nφ1(t)nφ1(t+τ+τsys))+12(φSNR1(t)φSNR2(t))+Δφsys(t)+Δφant(t) (15)

    其中,fd是多普勒频率, ${\tau _{{\rm{sys}}}}$ 是同步脉冲的发射延迟, ${n_{\varphi 1}}$ ${n_{\varphi 2}}$ 分别代表晶振的相位噪声, ${\varphi _{{\rm{SNR1}}}}$ ${\varphi _{{\rm{SNR2}}}}$ 分别代表主星和辅星的接收机噪声, $\Delta {\varphi _{{\rm{sys}}}}\left( t \right)$ 代表通道链路的相位噪声, $\Delta {\varphi _{{\rm{ant}}}}\left( t \right)$ 代表天线方向图引入的相位噪声。

    式(15)的第1项 $2{{π}} \left( {\Delta {f_2} - \Delta {f_1}} \right)$ ,等于式(7)中的线性偏差 ${\varphi _{{\rm{linear}}}}$ ,即随方位时间线性变化的相位误差部分可以通过相位同步方案补偿,补偿后的回波信号经过距离徙动校正后在方位频域的表达形式为:

    SRcom(τ,fa)=sinc(B(τR1+R2c+βfaKa))Wa(fafac)exp(j2π(R1+R2)/λ)exp(jπf2aKa) (16)

    其中,Ka为方位向调频率。式(16)中 $\beta \left( {{{{f_{\rm{a}}}} \mathord{\left/ {\vphantom {{{f_{\rm{a}}}} {{K_{\rm{a}}}}}} \right. } {{K_{\rm{a}}}}}} \right)$ 项会造成距离徙动校正误差,使方位信号弥散在若干距离单元上,导致主瓣展宽。如果要求距离徙动校正误差造成的主瓣展宽小于2%,那么要求未被校正的徙动量小于半个分辨率单元[22],即

    max(βKafa)<0.50.886B (17)

    根据时频关系可以得到

    βKaKaTsar<0.50.886Bβ<0.443BTsar (18)

    其中Tsar为合成孔径时间。

    满足式(18)的情况下,经方位压缩后的SAR信号为:

    σ(τ,t)=sinc(B(τR1+R2c))sinc(Bdt)exp(j2πR1+R2λ) exp(j2πfact) (19)

    此时主要导致主瓣展宽,严重情况可能导致散焦。

    理想情况下,SAR系统处于匀速运动状态,但是在实际工作情况下,卫星的运动受到控制系统精度的影响,会产生姿态误差,所以必须保证分布式系统的姿态控制精度[23]。波束同步误差可以分为两个方面,一方面是卫星姿态引起的误差,另一方面是波束指向误差。波束同步误差会导致天线照射范围不一致,导致信噪比损失,多普勒带宽偏移,方位向天线增益降低等问题[24]。同时,由于照射区域的偏差,还会导致距离向重叠幅宽减小。

    本节首先引入文献[7]中介绍的基于椭圆轨道的双星一发双收工作模式下的姿态导引方式,在此基础上,分析波束同步误差对多普勒去相干、天线增益以及双基图像信噪比损失的影响。

    双星有两种姿态导引方式[7],一种是只针对主星进行零多普勒姿态导引,此时,场景内的回波多普勒中心频率不为0。另一种是主辅双星均进行零多普勒姿态导引,此时辅星与主星照射的共同区域等效为正侧视,即场景内多普勒中心频率为0。图2图3给出了在主辅星顺轨基线从0到3 km情况下,不同姿态导引方式下方位向合成天线方向图的情况,其它仿真参数如表1所示。

    图  2  仅辅星进行多普勒导引
    Figure  2.  The first attitude steering strategy
    图  3  主星辅星分别进行多普勒导引
    Figure  3.  The second attitude steering strategy

    图2中,由于仅主星进行零多普勒姿态导引,导致合成后的方位天线方向图指向等效于斜视,即场景回波的多普勒中心频率不为0,产生方位多普勒去相干,导致相干系数发生变化,如图4所示。

    图  4  顺轨基线对相干系数的影响
    Figure  4.  Contributions to correlation coefficient versus along-track baseline

    考虑到顺轨基线和波束同步误差的综合影响,相干系数的计算方法可以表示为:

    γdc=1|ΔfDop1+ΔfDop2Bd| (20)

    $\Delta {f_{{\rm{Dop1}}}}$ 表示由顺轨基线导致的多普勒频率偏移,计算公式为[7]

    ΔfDop1vdalongλr (21)

    其中,dalong为基线长度,r为平均斜距。

    $\Delta {f_{{\rm{Dop2}}}}$ 表示由波束指向误差和卫星姿态误差导致的多普勒频率偏移,

    ΔfDop2=2vsin(θsq+Δθa)λ2vsin(θsq)λ (22)
    Δθa=θerr+θsy+θsp (23)

    其中, ${\theta _{{\rm{sq}}}}$ 表示斜视角, $\Delta {\theta _{\rm{a}}}$ 表示方位向波束同步误差,包括方位向波束指向误差以及偏航角和俯仰角对斜视角的影响部分。

    图3中,由于对主辅双星均进行零多普勒姿态导引,合成后的方位天线方向图等效于正侧视,此时对多普勒去相干影响很小。考虑到主辅双星之间顺轨基线的影响,合成后的天线方向图幅度会有一定的下降,从仿真中可以看出,在顺轨基线小于3 km的情况下,导致的幅度衰减小于0.05 dB,可以忽略。在不同顺轨基线情况下,主辅星同时进行零多普勒导引的合成天线方向图对频谱加权的计算公式为[7]

    Wbistat(fDop;ΔfDop)=sinc(l2v(fDopΔfDop))sinc(l2v(fDop+ΔfDop)) (24)

    其中,l为天线长度,v为卫星速度, $\Delta {f_{{\rm{Dop}}}}$ 是由顺轨基线导致的多普勒频率偏移。

    在此基础上,加入波束指向误差和卫星姿态误差的影响,顺轨基线和波束同步误差的综合加权公式可以表示为:

    Wbistat(fDop;ΔfDop)=sinc(l2v(fDopΔfDop1ΔfDop2))sinc(l2v(fDop+ΔfDop1+ΔfDop2)) (25)

    由于波束同步误差,导致成像后辅星图像信噪比降低,信噪比损失量为[25]

    SNRloss=θbw1/2θbw1/2(sinc2(lsinφ/λ)sinc2(lsin(φΔθa)/λ))dφθbw1/2θbw1/2(sinc2(lsinφ/λ)sinc2(lsinφ/λ))dφ (26)

    其中, $\varphi $ 为方位向离轴角, ${\theta _{{\rm{bw1}}}}$ 为方位向波束宽度。

    针对距离向,波束控制误差主要影响主、辅星的重叠幅宽,当主、辅星距离向同时存在相反方向的指向偏差时,距离向的重叠区域(图5图6中黄色区域)会大幅减小,影响后续的干涉处理。根据几何模型,可以推导出覆盖区域偏差的计算公式:

    图  5  无误差的波束覆盖情况
    Figure  5.  No beam pointing error case
    图  6  有误差的波束覆盖情况
    Figure  6.  Beam pointing error exists case
    ΔW=(H+Re)cos(θl+θbw2/2Δθr)(H+Re)2cos2(θl+θbw2/2Δθr)2HRe(H+Re)cos(θl+θbw2/2+Δθr)+(H+Re)2cos2(θl+θbw2/2+Δθr)2HRe (27)

    其中,H为卫星高度,Re为地球半径, ${\theta _{\rm{l}}}$ 为下视角, ${\theta _{{\rm{bw2}}}}$ 为俯仰向波束宽度, $\Delta {\theta _{\rm{r}}}$ 为距离向波束同步误差。

    仿真参数如表1所示。

    表  1  仿真参数表
    Table  1.  Simulation parameters
    参数 数值
    载频(GHz) 1.30
    卫星高度(km) 700
    方位向天线长度(m) 6
    脉冲宽度(μs) 80
    距离向带宽(MHz) 100
    PRF(Hz) 1600
    卫星有效速度(m/s) 7500
    距离向波束指向精度(°) 0.06
    方位向波束指向精度(°) 0.05
    姿态控制精度(°) 0.02
    顺轨基线(km) 3
    下载: 导出CSV 
    | 显示表格

    仿真中加入参数为 $\Delta t = 100 \ {\rm{ ns}}$ 的固定时间误差,点目标分析结果如图7所示,通过计算可以得到理论上的距离向位置偏移为 $\Delta r = \Delta t \cdot c/2{\rm{ = 15 \ m}}$ ,仿真得到的对应偏差为15 m,如图7(b)所示。

    图  7  固定时间误差影响仿真
    Figure  7.  Constant time error simulations

    时间同步误差会引起主瓣展宽,如式(18)所示,仿真参数下,当 $\beta < 1.2 \times {10^{ - 9}}$ 时,引起的展宽可以忽略。当 $\beta = 5 \times {10^{ - 8}}$ 时,距离向和方位向剖面图如图8所示,由仿真结果可以看到明显的主瓣展宽。

    图  8  回波域相位补偿后的结果
    Figure  8.  Time error simulations after compensation in echo domain

    本节的仿真给出了当顺轨基线为3 km时,波束指向误差和卫星姿态误差对系统性能的影响。图9图10分别比较了方位向波束指向精度在0.01°到0.05°情况下,对合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.05°波束指向误差造成方位天线增益下降小于0.05 dB,造成相干系数减小0.02。

    图  9  波束指向引起的合成天线方向图增益损失
    Figure  9.  Gain loss versus beam pointing accuracy
    图  10  波束指向引起的相干系数变化
    Figure  10.  Coherent coefficient versus beam pointing accuracy

    图11图12分别给出了卫星姿态控制精度在0.005°到0.020°的情况下,对方位向合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.020°的姿态控制误差造成方位天线增益下降小于0.02 dB,造成相干系数减小0.01。

    图  11  卫星姿态引起的合成天线方向图增益损失
    Figure  11.  Gain loss versus satellite attitude accuracy
    图  12  卫星姿态引起的合成天线方向图增益损失
    Figure  12.  Coherent coefficient versus satellite attitude accuracy

    不同波束指向和卫星姿态误差情况下,对合成天线方向图及相干系数的综合影响结果如图13图14所示。波束指向精度为0.05°,卫星姿态控制精度0.02°的情况下,和无误差情况对比,造成方位天线增益下降小于0.1 dB,造成相干系数减小0.03。

    图  13  波束指向和卫星姿态对天线增益的综合影响
    Figure  13.  Overall influence on antenna gain versus beam pointing accuracy and attitude accuracy
    图  14  波束指向和卫星姿态对相干系数的综合影响
    Figure  14.  Overall influence on coherent coefficient versus beam pointing accuracy and attitude accuracy

    不同波束指向和卫星姿态误差情况下,辅星图像信噪比损失情况如表2所示,0.05°的波束指向误差和0.02°的姿态误差导致的辅星图像信噪比损失小于0.1 dB。

    表  2  波束同步误差造成的辅星图像信噪比损失(dB)
    Table  2.  Influence of beam synchronization error on SNR (dB)
    波束指向精度(°) 姿态控制精度(°)
    0.01 0.02
    0.01 0.021 0.025
    0.02 0.028 0.032
    0.03 0.036 0.039
    0.04 0.046 0.048
    0.05 0.055 0.057
    下载: 导出CSV 
    | 显示表格

    0°~0.06°的波束指向误差对距离向重叠幅宽的影响如图15所示。横滚角会在距离向产生波束控制误差,对距离向覆盖偏差产生影响,0.02°的姿态控制误差对重叠幅宽的影响如图16所示。横滚角从0°步进至0.02°,距离向波束指向精度从0°步进至0.06°情况下,两者对距离向重叠幅宽的综合影响如图17所示。上述结果表明,0.02°的姿态控制误差对距离向重叠幅宽的影响为368 m, 0.06°的波束指向误差对距离向重叠幅宽的影响为1100 m,由于波束指向和姿态控制导致的误差均为随机误差,所以综合影响采用几何平均数的形式,两者综合之后对距离向重叠幅宽的影响为1164 m。考虑到轨道控制精度的影响,在重轨工作模式下,两次测量中卫星的位置会在一定范围内偏移,在卫星姿态误差和波束指向误差基础上,系统设计中应综合考虑卫星位置偏移的影响。

    图  15  波束指向对重叠幅宽的影响
    Figure  15.  Overlapping swath error versus beam pointing accuracy
    图  16  横滚角对重叠幅宽的影响
    Figure  16.  Overlapping swath error versus roll angle
    图  17  横滚角和波束指向对重叠幅宽的影响
    Figure  17.  Overlapping swath error versus roll angle and beam pointing error

    由上述分析可知,在仿真精度情况下,波束同步精度对多普勒去相干、天线增益以及双基图像信噪比损失的影响较小,但对距离向重叠幅宽影响较大,因此需要在波位设计中保留一定的幅宽余量,或者根据幅宽需求,对波束同步精度提出要求。不同卫星姿态控制精度和波束指向精度对重叠幅宽的综合影响总结在表3可供系统设计参考。针对不同的卫星系统,可以根据系统参数结合式(27)求得对波束同步误差的限制。

    表  3  波束同步误差对重叠幅宽的影响(m)
    Table  3.  Influence of beam synchronization error to overlapping area (m)
    波束指向精度(°) 姿态控制精度(°)
    0.01 0.02
    0.01 260 411
    0.02 411 520
    0.03 582 663
    0.04 759 823
    0.05 938 991
    0.06 1119 1164
    下载: 导出CSV 
    | 显示表格

    本文分析了分布式SAR系统中时间同步误差和波束同步误差的影响,根据固定项和线性项的时间同步误差模型,推导时间同步和相位同步的关系,并分析线性项需要满足的条件。波束同步部分,本文分析了波束指向误差和卫星姿态误差对方位向合成天线方向图、多普勒去相干以及双基图像信噪比损失情况的影响,并根据几何关系推导距离向重叠幅宽的损失量。根据仿真参数,给出实验结果,以便在系统设计中充分考虑同步问题的影响,提出指标要求,为系统设计提供指导。

  • 图  1  相关杂波下的信号模型

    Figure  1.  Signal model

    图  2  确定TIR和随机CIR

    Figure  2.  Determinate TIR and random CIR

    图  3  随迭代次数的算法性能分析

    Figure  3.  Algorithm performance analysis with iterations

    图  4  不同算法信干噪比随着迭代次数变化的曲线

    Figure  4.  Comparison of the SINR for different algorithms with iterations

    图  5  不同算法信干噪比随着CNR变化的曲线

    Figure  5.  Comparison of the SINR for different algorithms with different CNR

    图  6  本文算法不同PAR下的信干噪比性能对比

    Figure  6.  Comparison of the SINR vs different PAR

    图  7  不同PAR波形的实部和虚部

    Figure  7.  Real and imaginary parts of different PAR

  • [1] Haykin S. Cognitive radar: A way of the future[J]. IEEE Signal Processing Magazine, 2006, 23(1): 30–40. DOI: 10.1109/MSP.2006.1593335
    [2] Guerci J R. Cognitive Radar: The Knowledge-Aided Fully Adaptive Approach[M]. Norwood, MA: Artech House, Inc., 2010: 13–30
    [3] Wang Z J, Babu P, and Palomar D P. Design of PAR-constrained sequences for MIMO channel estimation via majorization-minimization[J]. IEEE Transactions on Signal Processing, 2016, 64(23): 6132–6144. DOI: 10.1109/TSP.2016.2607181
    [4] Stoica P, He H, and Li J. Optimization of the receive filter and transmit sequence for active sensing[J]. IEEE Transactions on Signal Processing, 2012, 60(4): 1730–1740. DOI: 10.1109/TSP.2011.2179652
    [5] 唐波, 张玉, 李科, 等. 杂波中MIMO雷达恒模波形及接收机联合优化算法研究[J]. 电子学报, 2014, 42(9): 1705–1711. DOI: 10.3969/j.issn.0372-2112.2014.09.007

    Tang Bo, Zhang Yu, Li Ke, et al. Joint constant-envelope waveform and receiver design for MIMO radar in the presence of clutter[J]. Acta Electronica Sinica, 2014, 42(9): 1705–1711. DOI: 10.3969/j.issn.0372-2112.2014.09.007
    [6] 郝天铎, 周青松, 孙从易, 等. 非准确先验知识下认知雷达低峰均比稳健波形设计[J]. 电子与信息学报, 2018, 40(3): 532–540. DOI: 10.11999/JEIT170560

    Hao Tian-duo, Zhou Qing-song, Sun Cong-yi, et al. Low-PAR robust waveform design for cognitive radar with imprecise prior knowledge[J]. Journal of Electronics&Information Technology, 2018, 40(3): 532–540. DOI: 10.11999/JEIT170560
    [7] Sen S. Characterizations of PAPR-constrained radar waveforms for optimal target detection[J]. IEEE Sensors Journal, 2014, 14(5): 1647–1653. DOI: 10.1109/JSEN.2014.2299283
    [8] Bell M R. Information theory and radar waveform design[J]. IEEE Transactions on Information Theory, 1993, 39(5): 1578–1597. DOI: 10.1109/18.259642
    [9] Tang B and Tang J. Robust waveform design of wideband cognitive radar for extended target detection[C]. IEEE International Conference on Acoustics, Speech and Signal Processing, Shanghai, China, 2016: 3096–3100. DOI: 10.1109/ICASSP.2016.7472247
    [10] Chen C Y and Vaidyanathan P P. MIMO radar waveform optimization with prior information of the extended target and clutter[J]. IEEE Transactions on Signal Processing, 2009, 57(9): 3533–3544. DOI: 10.1109/TSP.2009.2021632
    [11] Chen X X, Deng X B, and Hao Z M. Waveform design for extended target detection under a peak to average power ratio constraint[C]. Proceedings of 2016 CIE International Conference on Radar, Guangzhou, China, 2016: 1–4. DOI: 10.1109/RADAR.2016.8059566
    [12] Tang Y H, Zhang Y D, Amin M G, et al. Wideband multiple-input multiple-output radar waveform design with low peak-to-average ratio constraint[J]. IET Radar,Sonar&Navigation, 2016, 10(2): 325–332. DOI: 10.1049/iet-rsn.2015.0189
    [13] Gorji A A and Adve R S. Waveform optimization for random-phase radar signals with PAPR constraints[C]. IEEE International Radar Conference, Lille, France, 2014: 1–5. DOI: 10.1109/RADAR.2014.7060340
    [14] Yue W Z, Zhang Y, Liu Y M, et al. Radar constant-modulus waveform design with prior information of the extended target and clutter[J]. Sensors, 2016, 16(6): 889. DOI: 10.3390/s16060889
    [15] 唐波. 宽带认知雷达低峰均比波形快速设计算法[J]. 航空学报, 2016, 37(2): 688–694. DOI: 10.7527/S1000-6893.2015.0125

    Tang B. Efficient design algorithm of low PAR waveform for wideband cognitive radar[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(2): 688–694. DOI: 10.7527/S1000-6893.2015.0125
    [16] Romero R A, Bae J, and Goodman N A. Theory and application of SNR and mutual information matched illumination waveforms[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(2): 912–927. DOI: 10.1109/TAES.2011.5751234
    [17] Chen P, Qi C H, Wu L N, et al. Estimation of extended targets based on compressed sensing in cognitive radar system[J]. IEEE Transactions on Vehicular Technology, 2017, 66(2): 941–951. DOI: 10.1109/TVT.2016.2565518
    [18] 邹鲲, 廖桂生, 李军. 复合高斯杂波中知识辅助检测器的先验信息感知方法[J]. 中国科学: 信息科学, 2014, 44(8): 993–1003. DOI: 10.1360/N112013-00116

    Zou Kun, Liao Gui-sheng, and Li Jun. Prior information cognitive method for knowledge aided detector in compound Gaussian clutter[J]. Scientia Sinica Informationis, 2014, 44(8): 993–1003. DOI: 10.1360/N112013-00116
    [19] 张贤达. 矩阵分析与应用[M]. 第2版, 北京: 清华大学出版社, 2013: 447–450

    Zhang Xian-da. Matrix Analysis and Application[M]. Second Edition, Beijing: Tsinghua University Press, 2013: 447–450
    [20] Luo Z Q, Ma W K, So A M C, et al. Semidefinite relaxation of quadratic optimization problems[J]. IEEE Signal Processing Magazine, 2010, 27(3): 20–34. DOI: 10.1109/MSP.2010.936019
    [21] Yue W Z, Zhang Y, and Xie J W. Radar constant-modulus waveform design for multiple extended targets[J]. IEICE Transactions on Fundamentals of Electronics,Communications and Computer Sciences, 2017, E100-A(3): 888–892. DOI: 10.1587/transfun.E100.A.888
    [22] Grant M and Boyd S. CVX: Matlab software for disciplined convex programming (Web page and software) 2008[EB/OL]. http://cvxr.com/cvx/
    [23] Leshem A, Naparstek O, and Nehorai A. Information theoretic adaptive radar waveform design for multiple extended targets[J]. IEEE Journal of Selected Topics in Signal Processing, 2007, 1(1): 42–55. DOI: 10.1109/JSTSP.2007.897047
    [24] Yang Y and Blum R S. MIMO radar waveform design based on mutual information and minimum mean-square error estimation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2007, 43(1): 330–343. DOI: 10.1109/TAES.2007.357137
  • 加载中
图(7)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-02
  • 修回日期:  2018-04-23
  • 网络出版日期:  2018-08-28

目录

/

返回文章
返回