一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法
【专利摘要】本发明公开了一种主副瓣地海杂波主分布计算方法，包含以下步骤：根据雷达下视尾后探测波束各方向增益强弱，确定主副瓣地海杂波主分布对应区域；在对应区域内选取相应的杂波点，并计算所有杂波点的位置坐标；根据雷达的位置坐标、速度矢量参数及所有杂波点的位置坐标，计算所有杂波点与雷达的相对距离和相对多普勒频率；结合雷达发射信号的最大不模糊距离、脉冲重复频率参数，计算所有杂波点的视在距离和视在多普勒频率；量化所有杂波点的视在距离和视在多普勒频率，并将所有杂波点的量化距离和量化多普勒频率映射至距离?多普勒域，得到主副瓣地海杂波主分布的准确位置。本发明计算实时性强，有效提升雷达下视尾后探测威力和跟踪性能。
【专利说明】
一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法
技术领域
[0001] 本发明涉及雷达下视尾后探测技术领域，具体涉及一种雷达下视尾后探测下主副 瓣地海杂波主分布计算方法。
【背景技术】
[0002] 地海杂波是雷达下视尾后探测和跟踪目标性能的最大制约因素。下视尾后探测空 中目标，目标回波信号落入副瓣杂波区，雷达发射中重复频率波形用于探测目标。相对高重 复频率波形，采用中重复频率波形能够减小副瓣杂波距离混叠程度，并进一步通过距离维 切割降低副瓣杂波功率，使雷达具备下视尾后探测能力。但由于副瓣杂波在雷达处理的距 离-多普勒域上并非均匀分布，而是存在强弱分布，当目标回波信号落入强副瓣杂波区域 时，一方面大大降低了检测信杂比，影响雷达探测距离;另一方面极大增加了雷达错误锁定 或者跟踪副瓣杂波的概率，导致丢失目标，这种情况亟需改进。
【发明内容】
[0003] 本发明的目的在于提供一种雷达下视尾后探测下主副瓣地海杂波主分布计算方 法，能够在线计算获取地海杂波在距离-多普勒域上的主分布准确位置，解决强副瓣杂波严 重影响雷达下视尾后探测和跟踪目标性能的问题，为雷达下视尾后探测目标时使目标回波 信号避开强副瓣杂波，避免落入强副瓣杂波区域提供必要参考条件，有效提升雷达下视尾 后探测威力和跟踪性能。
[0004] 为了达到上述目的，本发明通过以下技术方案实现：一种雷达下视尾后探测下主 副瓣地海杂波主分布计算方法，其特点是，包含以下步骤：
[0005] S1、根据雷达下视尾后探测波束各方向增益强弱，确定主副瓣地海杂波主分布对 应区域；
[0006] S2、在主副瓣地海杂波主分布对应区域内选取相应的杂波点，并计算所有杂波点 的位置坐标；
[0007] S3、根据雷达的位置坐标、雷达的速度矢量参数及所有杂波点的位置坐标，计算所 有杂波点与雷达的相对距离和相对多普勒频率；
[0008] S4、结合雷达发射信号的最大不模糊距离和雷达发射信号的脉冲重复频率参数， 计算所有杂波点的视在距离和视在多普勒频率；
[0009] S5、量化所有杂波点的视在距离和视在多普勒频率，并将所有杂波点的量化距离 和量化多普勒频率映射至距离-多普勒域，得到主副瓣地海杂波主分布的准确位置。
[0010] 所述的步骤S1包含：
[0011] S1.1、分别建立基准坐标系和平台坐标系；
[0012] S1.2、根据实际雷达天线方向图强弱分布特征，以平台坐标系下的指向角，设置主 副瓣地海杂波远距边界点和近距边界点与雷达天线相对指向，并根据雷达在基准坐标系下 的位置和姿态分别计算远距边界点和近距边界点在基准坐标系下的坐标；
[0013] SI .3、在基准坐标系下，由远距边界点和近距边界点的连线构成主副瓣地海杂波 主分布对应区域的计算区间。
[0014] 所述的步骤S1.1中基准坐标系的原点设定在地面/海面上，X轴为北向，Y轴为天 向，Z轴为东向，雷达所处位置定义为PR，在基准坐标系下的坐标为(x R，yR，zR)，按照右手定 则确定雷达姿态，雷达姿态定义为（Ψ，θ，γ )，Ψ表示偏航角，Θ表示俯仰角，γ表示滚转角。 [0015]所述的步骤S1.1中的平台坐标系是将基准坐标系原点平移至雷达所处位置，再将 其依次按照Ψ、θ、γ姿态角进行三次转动，得到平台坐标系P RX〃Y〃Z〃，指向角(α，β)表示雷 达天线波束、地面/海面位置点与雷达在平台坐标系中的相对位置关系，α表示方位角，β表 示俯仰角。
[0016] 所述的步骤S1.2中远距边界点在基准坐标系下的坐标计算方法包含：
[0017] 在平台坐标系下，由远距边界点指向角计算得至I」指向单位矢量；
[0018] 通过欧拉坐标转换，将平台坐标系下的指向单位矢量转换到基准坐标系中，计算 得到基准坐标系下指向单位矢量；
[00?9]通过基准坐标系下指向单位矢量、平台坐标计算得到远距边界点在基准坐标系下 的坐标。
[0020] 所述的步骤S1.2中近距边界点在基准坐标系下的坐标计算方法包含：
[0021] 在平台坐标系下，由近距边界点指向角计算得到指向单位矢量；
[0022] 通过欧拉坐标转换，将平台坐标系下的指向单位矢量转换到基准坐标系中，计算 得到基准坐标系下指向单位矢量；
[0023]通过基准坐标系下指向单位矢量、平台坐标计算得到近距边界点在基准坐标系下 的坐标。
[0024]所述的步骤S2包含：
[0025] 在基准坐标系下，根据远距边界点的坐标和近距边界点的坐标计算远距边界点和 近距边界点的相对距离；
[0026] 在基准坐标系下，通过远距边界点的坐标和近距边界点的坐标计算远距边界点到 近距边界点连线与Ζ轴的夹角；
[0027] 设定远距边界点和近距边界点连线上取点个数，依次计算各选取点相对于近距边 界点的距离及各选取点在基准坐标系中的(χ，ζ)坐标；
[0028] 依次根据各选取点在基准坐标系中的(Χ，Ζ)坐标，从高程数据中查询各选取点的Υ 坐标，以得到所有杂波点在基准坐标系下的位置坐标。
[0029] 所述的步骤S3中所有杂波点与雷达的相对距离的计算方法为：
[0030] 在基准坐标系下，利用杂波点的位置坐标和雷达的位置坐标，根据两点间距离计 算公式可得所有杂波点与雷达的相对距离。
[0031] 所述的步骤S3中所有杂波点与雷达的相对多普勒频率的计算方法为：
[0032] 在基准坐标系下，利用杂波点的位置坐标和雷达的位置坐标，计算得到指向矢量；
[0033] 利用指向矢量和雷达速度矢量计算得到各杂波点和雷达连线与速度方向的夹角；
[0034] 利用平台合成速度值计算得打杂波点与雷达的相对多普勒频率。
[0035] 所述的远距边界点和近距边界点连线上取点个数Ν，需满足以下条件：
[0037] 其中，Δ R表示相邻取点直接的距离间隔，RFN表示远距边界点和近距边界点直接的 相对距离，N表示取点个数，ARr表示雷达距离分辨率。
[0038] 本发明一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法与现有技术相 比具有以下优点：通过确定主副瓣杂波主分布与地面/海面区域的对应关系简化杂波区域 选择，缩小计算区域，计算区域内所有杂波单元与雷达的相对距离和相对多普勒频率，并映 射至距离-多普勒域，得到杂波主分布的准确位置，能够在雷达下视尾后探测目标时，快速 计算主副瓣杂波回波信号在距离-多普勒域主分布位置，为雷达利用目标信号与强副瓣杂 波信号相对距离和频率关系判断或者避免目标信号落入或者穿越强副瓣杂波区提供核心 参数和重要依据;本发明计算实时性强，能够满足机载、弹载雷达对处理实时性的要求。
【附图说明】
[0039] 图1为本发明一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法的流程 图；
[0040] 图2为雷达下视尾后探测的几何关系框图；
[0041] 图3为本发明实施例中主副瓣地海杂波主分布计算结果示意图。
【具体实施方式】
[0042] 以下结合附图，通过详细说明一个较佳的具体实施例，对本发明做进一步阐述。
[0043] 如图1所示，一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法，包含以下 步骤：
[0044] S1、根据雷达下视尾后探测波束各方向增益强弱，确定主副瓣地海杂波主分布对 应区域。
[0045] 步骤S1具体包含：
[0046] SI. 1、分别建立基准坐标系和平台坐标系，如图2所示；
[0047]其中，基准坐标系的原点设定在地面/海面上，X轴为北向，Y轴为天向，Z轴为东向， 雷达所处位置定义为Pr，在基准坐标系下的坐标为（仙，71?，2[0，按照右手定则确定雷达姿 态，雷达姿态定义为（Ψ，θ，γ )，ψ表示偏航角，Θ表示俯仰角，γ表示滚转角;平台坐标系是 将基准坐标系原点平移至雷达所处位置，再将其依次按照ψ、θ、γ姿态角进行三次转动，得 到平台坐标系P RX〃Y〃Z〃，指向角(α，β)表示雷达天线波束、地面/海面位置点与雷达在平台 坐标系中的相对位置关系，α表示方位角，β表示俯仰角；
[0048] S1.2、根据实际雷达天线方向图强弱分布特征，以平台坐标系下的指向角，设置主 副瓣地海杂波远距边界点Pi和近距边界点Ρν与雷达天线相对指向，远距边界点PjP近距边 界点Pn的指向角分别为(aFar，fcar)和(aNear，fo ear)，并根据雷达在基准坐标系下的位置和姿 态分别计算远距边界APi和近距边界点Pn在基准坐标系下的坐标；
[0049]远距边界点丹在基准坐标系下的坐标计算方法包含：
[0050]在平台坐标系下，由远距边界点指向角（aFar，0Far)计算得到指向单位矢量（χ Ρ1， ypi^zpi)；
[0051] XPl = COS(0Far) x COS(aFar)
[0052] ypi = COs(0Far) X sin(apar)
[0053] ZPl = sin(-0Far)
[0054] 通过欧拉坐标转换，将平台坐标系下的指向单位矢量(xpl，ypl，zpl)转换到基准坐 标系中，计算得到基准坐标系下指向单位矢量(xui，yui，zui);
[0055]
[005?]通过基准坐标系下指向单位矢量、平台坐标计算得到远距边界点Pi在基准坐标系 下的坐标(xi，yi，zi);
[0058] 同理，近距边界点Pn在基准坐标系下的坐标计算方法包含：
[0059] 在平台坐标系下，由近距边界点指向角计算得到指向单位矢量；
[0060] 通过欧拉坐标转换，将平台坐标系下的指向单位矢量转换到基准坐标系中，计算 得到基准坐标系下指向单位矢量；
[0061]通过基准坐标系下指向单位矢量、平台坐标计算得到近距边界点Pn在基准坐标系 下的坐标(XN，yN，ZN);
[0063] S1.3、在基准坐标系下，由远距边界APi和近距边界点Pn的连线构成主副瓣地海杂 波主分布对应区域的计算区间。
[0064] S2、在主副瓣地海杂波主分布对应区域内选取相应的杂波点，并计算所有杂波点 的位置坐标。
[0065] 步骤S2具体包含：
[0066] 在基准坐标系下，根据远距边界点的坐标(X1，yi， Z1)和近距边界点的坐标(xN，yN， zN)计算远距边界点朽和近距边界点Pn的相对距离RFN;
[0068] 在基准坐标系下，通过远距边界点的坐标(XI，y 1，z 1)和近距边界点的坐标(xn，yN， ZN)计算远距边界点丹到近距边界点Pn连线与Z轴的夹角ε ;
[0070] 设定远距边界APi和近距边界点Pn连线上取点个数Ν，依次计算各选取点相对于近 距边界点Pn的距离RnN及各选取点在基准坐标系中的(X，Z)坐标(Xn，zn);
[0071 ]其中，取点个数N需满足条件
[0073] RnN= (N-n+1) X Δ R
[0074] xn=RnNX sin( ε )+xn
[0075] Zn = RnNX cos( ε )+ZN
[0076] 依次根据各选取点在基准坐标系中的(X，Z)坐标(Xn，zn)，从高程数据中查询各选 取点的Y坐标y n，其中yn=HCn，HCn为选取点所在地面/海面高度，以得到所有杂波点在基准坐 标系下的位置坐标(x n，yn，ζη)，η由1至ijN取值。
[0077] S3、根据雷达的位置坐标、雷达的速度矢量参数及所有杂波点的位置坐标，计算所 有杂波点与雷达的相对距离和相对多普勒频率。
[0078] 步骤S3具体包含：
[0079]在基准坐标系下，利用杂波点的位置坐标(Xn，yn，Zn)和雷达的位置坐标( XR，yR， zr)，根据两点间距离计算公式可得所有杂波点与雷达的相对距离Rn;
[0081 ]在基准坐标系下，利用杂波点的位置坐标(Xn，yn，Zn)和雷达的位置坐标( XR，yR， zr)，计算得到指向矢量/yi;
[0082] 户= (? -xn，），R ―為,? -&)
[0083] 利用指向矢和雷达速度矢量€ (vRx:，vRy,VRz)计算得到各杂波点和雷达连线 与速度方向的夹角Φ η;
[0085]利用雷达合成速度值计算得到杂波点与雷达的相对多普勒频率Fn;
[0087]其中，λ为雷达发射信号波长，雷达速度矢量由惯性装置获得，机载雷达或者弹载 雷达即对应的平台的速度矢量。
[0088] S4、结合雷达发射信号的最大不模糊距离和雷达发射信号的脉冲重复频率参数， 计算所有杂波点的视在距离和视在多普勒频率。
[0089]步骤S4中视在距离Δ Rn和视在多普勒频率Δ Fn的计算公式分别为：
[0090] Δ Rn=mod(Rn，Rmax)
[0091] 其中，Rmax为发射信号最大不模糊距离；
[0092] AFn=mod(Fn,PRF)
[0093] 其中，PRF为发射信号的脉冲重复频率。
[0094] S5、量化所有杂波点的视在距离和视在多普勒频率，并将所有杂波点的量化距离 和量化多普勒频率映射至距离-多普勒域（A QRn，Δ QFn)，（n = 1，2···Ν)，得到主副瓣地海杂 波主分布的准确位置。
[0095]步骤S5中，量化所有杂波点的视在距离△ Rn和视在多普勒频率△ Fn，得到量化距离 Δ QRn和量化多普勒频率Δ QFn，其中
[0097]其中，ArR为雷达距离分辨率；
[0099] 其中，Δ Fr雷达频率分辨率。
[0100] 如图3所示，给出了雷达坐标(61310，14958,52736)，雷达姿态(-38，-22,0)，速度 矢量（709，-371，564)，远距边界点指向（-0 · 5，-10)，近距边界点指向（0，-80)，雷达发射信 号的脉冲重复频率41kHz条件下，一个主副瓣地海杂波主分布计算结果图。
[0101] 尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的 描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的 多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。
【主权项】
1. 一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征在于，包含W下 步骤： 51、 根据雷达下视尾后探测波束各方向增益强弱，确定主副瓣地海杂波主分布对应区 域； 52、 在主副瓣地海杂波主分布对应区域内选取相应的杂波点，并计算所有杂波点的位 置坐标； 53、 根据雷达的位置坐标、雷达的速度矢量参数及所有杂波点的位置坐标，计算所有杂 波点与雷达的相对距离和相对多普勒频率； 54、 结合雷达发射信号的最大不模糊距离和雷达发射信号的脉冲重复频率参数，计算 所有杂波点的视在距离和视在多普勒频率； 55、 量化所有杂波点的视在距离和视在多普勒频率，并将所有杂波点的量化距离和量 化多普勒频率映射至距离-多普勒域，得到主副瓣地海杂波主分布的准确位置。2. 如权利要求1所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S1包含： 51.1、 分别建立基准坐标系和平台坐标系； 51.2、 根据实际雷达天线方向图强弱分布特征，W平台坐标系下的指向角，设置主副瓣 地海杂波远距边界点和近距边界点与雷达天线相对指向，并根据雷达在基准坐标系下的位 置和姿态分别计算远距边界点和近距边界点在基准坐标系下的坐标； 51.3、 在基准坐标系下，由远距边界点和近距边界点的连线构成主副瓣地海杂波主分 布对应区域的计算区间。3. 如权利要求2所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S1.1中基准坐标系的原点设定在地面/海面上，X轴为北向，Y轴为天向，Z 轴为东向，雷达所处位置定义为Pr，在基准坐标系下的坐标为(XR，yR，ZR)，按照右手定则确 定雷达姿态，雷达姿态定义为（Ψ，θ，丫），Ψ表示偏航角，Θ表示俯仰角，丫表示滚转角。4. 如权利要求3所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S1.1中的平台坐标系是将基准坐标系原点平移至雷达所处位置，再将其 依次按照Ψ、9、丫姿态角进行Ξ次转动，得到平台坐标系扣Χ"Υ"^，指向角(α，β)表示雷达 天线波束、地面/海面位置点与雷达在平台坐标系中的相对位置关系，α表示方位角，β表示 俯仰角。5. 如权利要求4所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S1.2中远距边界点在基准坐标系下的坐标计算方法包含： 在平台坐标系下，由远距边界点指向角计算得到指向单位矢量. 通过欧拉坐标转换，将平台坐标系下的指向单位矢量转换到基准坐标系中，计算得到 基准坐标系下指向单位矢量； 通过基准坐标系下指向单位矢量、平台坐标计算得到远距边界点在基准坐标系下的坐 标。6. 如权利要求5所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S1.2中近距边界点在基准坐标系下的坐标计算方法包含： 在平台坐标系下，由近距边界点指向角计算得到指向单位矢量. 通过欧拉坐标转换，将平台坐标系下的指向单位矢量转换到基准坐标系中，计算得到 基准坐标系下指向单位矢量； 通过基准坐标系下指向单位矢量、平台坐标计算得到近距边界点在基准坐标系下的坐 标。7. 如权利要求6所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S2包含： 在基准坐标系下，根据远距边界点的坐标和近距边界点的坐标计算远距边界点和近距 边界点的相对距离； 在基准坐标系下，通过远距边界点的坐标和近距边界点的坐标计算远距边界点到近距 边界点连线与Z轴的夹角； 设定远距边界点和近距边界点连线上取点个数，依次计算各选取点相对于近距边界点 的距离及各选取点在基准坐标系中的(X，Z)坐标； 依次根据各选取点在基准坐标系中的（X，Z)坐标，从高程数据中查询各选取点的Y坐 标，W得到所有杂波点在基准坐标系下的位置坐标。8. 如权利要求7所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S3中所有杂波点与雷达的相对距离的计算方法为： 在基准坐标系下，利用杂波点的位置坐标和雷达的位置坐标，根据两点间距离计算公 式可得所有杂波点与雷达的相对距离。9. 如权利要求7所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的步骤S3中所有杂波点与雷达的相对多普勒频率的计算方法为： 在基准坐标系下，利用杂波点的位置坐标和雷达的位置坐标，计算得到指向矢量. 利用指向矢量和雷达速度矢量计算得到各杂波点和雷达连线与速度方向的夹角； 利用平台合成速度值计算得打杂波点与雷达的相对多普勒频率。10. 如权利要求7所述的雷达下视尾后探测下主副瓣地海杂波主分布计算方法，其特征 在于，所述的远距边界点和近距边界点连线上取点个数N，需满足W下条件：其中，Δ R表示相邻取点直接的距离间隔，化N表示远距边界点和近距边界点直接的相对 距离，N表示取点个数，Δ化表示雷达距离分辨率。
【文档编号】G01S7/41GK106093883SQ201610377418
【公开日】2016年11月9日
【申请日】2016年5月31日
【发明人】杜科, 易堃, 谭智, 马超, 唐文明, 沈军, 郑博
【申请人】上海无线电设备研究所