《任意姿态机载 PD 雷达通用杂波算法研究》由会员分享，可在线阅读，更多相关《任意姿态机载 PD 雷达通用杂波算法研究（5页珍藏版）》请在金锄头文库上搜索。
1、精品论文任意姿态机载 PD 雷达通用杂波算法研究李艳玲1, 杨建宇2, 崔国龙21四川邮电职业技术学院 成都 6100672电子科技大学电子工程学院 成都 610054E-mail: 摘要：在载机和天线处于任意姿态的情况下，本文采用距离环地面散射单元划分方法，推导了计算机载 PD 雷达距离-多普勒单元地物杂波的通用封闭表达式，并着重分析了载机和天线在任意姿态下方位角变 化对机载 PD 雷达杂波功率谱的影响，做出了不同姿态情况下的距离多普勒二维杂波图。仿真结果与理论 分析吻合较好，证明了该算法的有效性。关键词：机载 PD 雷达；杂波；距离多普勒1.引言地面杂波对采用下视工作方式的机载 PD 雷达
2、的目标检测和识别具有重要影响，由于 “机载”和“下视”环境下的特殊性，机载 PD 雷达杂波建模的主要工作是杂波功率谱的计 算方法。机载 PD 雷达的杂波谱是多普勒频率和距离的函数，由于机载下视雷达与地面之间存在 相对运动，再加上雷达天线方向图的影响，使机载雷达地面杂波的频谱发生了显著变化，需 认真研究，以便采用有效的系统设计和信号处理方法，提高雷达整体性能，近年来国内外许 多专家学者对机载 PD 雷达杂波建模进行了大量研究13。文献4提出一种适用于机载 PD 雷达的杂波谱的快速计算方法，但任意姿态下机载雷达杂波模型及算法的研究却少有提 及。文献5研究了任意姿态下机载雷达杂波功率模型，但主要是对
3、地形的研究，没有给出 杂波功率谱的具体计算方法。文献6 研究了任意姿态下机载雷达相干杂波模型，但该模型 实现较为困难，缺乏通用性。本文在此基础上推导了任意姿态机载 PD 雷达通用杂波算法， 并对此算法进行仿真，验证了算法的有效性。2.机载 PD 雷达杂波的主要特点地面散射单元的划分通常采用距离-多普勒地面单元划分法，即将地面按等距离环和等 多普勒线划分为若干个杂波单元(见图 1)，将每一个杂波单元视做一个面目标，假设杂波是 均匀的，即来自不同散射单元的反向散射信号是统计独立的，在空间上没有相干性，对所有 距离环的杂波单元做功率元的积分后再求和，则可计算出杂波的功率谱。c- 1 -距离分辨率和多
4、普勒分辨率分别为：R =f = PRF / N FFT2不模糊距离和不模糊频率表示为： Ru=c2PRFfu = PRFRm = R0 + mRu；m = 0 , 1 , 2 fn = f0 + nfu； n = 0 , 1 , 2 R0 , f0 为回波对应距离门和多普勒频移，在同一距离门和多普勒带宽内满足距离模糊和 频率模糊 的目标和 地面单元 应满足如 下距离和 多普勒频 率条 件 ：Rm R Rm + R fn f fn + f图 1 机载雷达与地面杂波源的几何关图 2 惯性参考坐标系、载体坐标系、天线坐标系及地面散射单元的空间关系对单基地雷达来说，对给定的杂波类型、信号频率、极化，后
5、向散射系数 0 与入射角有关。对于机载环境，存在天线旁瓣波束近似垂直入射形成的高度线杂波，杂波后向散射系数计算应该采用修正的地面散射模型1od 0 = sin s + os e (90D s )2 / 2 od 、 os 和 适于散射有关的参量， s 是擦地角。右边第一项是漫反射分量，第二项是镜面反射分量，即构成机载雷达的高度线回波。3.任意姿态机载 PD 雷达杂波功率的计算3.1 单个散射单元参数确定由于载机处于任意姿态，因此涉及到多个坐标系7，包括：惯性参考坐标系、载体坐 标系、天线坐标系等，如图 2 所示。惯性参考坐标系原点为载机质心，OG ZG 轴指向地心，OG X G轴指向正东，OG
6、 X G轴指向正北；惯性参考坐标系依次绕 OG ZG 、OG X G 、OG X G 轴分别作逆时针旋转偏航角P 、倾斜角 P 、滚动角 P ，得到载体坐标系，且 OP X P 轴与载机速度 矢量方向一致；载体坐标系依次绕 OP Z P 、 OPYP 轴作逆时针旋转偏航角A ，倾斜角 A ，得到天线坐标系， OA X A 轴与天线波束指向一致。若杂波散射单元在惯性参考坐标系中的方位角为 、俯仰角为 ，则根据天线坐标系 和惯性坐标系的转换关系可以确定该散射单元在天线坐标系中的方位角T 和俯仰角T 。cos T cosT cos cos T T Y A Z A X P Y P Y P 且 cos
7、sin = R ( ) R ( ) R ( ) R ( ) R ( ) cos sin sin Tcossin 0sin cos 0 - sin RZ ( ) = sin cos 0 RY ( ) = 0 1 0 0 0 1sin0 cos精品论文1 0 0 RX ( ) = 0 cos0 -sinsin cos 式中：P 、 P 、 P 分别为载机的偏航角、俯仰角和横滚角； A 、 A 为天线的偏航角和俯仰角。假上式右边列向量可表示：a a a T ，则散射单元在天线坐标系下的方位角 和俯仰角T 表示如下1 2 3 T同理2 2T = arctan(a2 / a1 )T = arctan(a
8、3 /a2 + a1 ) a = arccos(cos P cos P cos cos + cos P sin P cos sin sin P sin )3.2 距离-多普勒杂波功率谱闭合计算算法距离-多普勒单元与坐标系定义如图 1 所示，定义水平方向矢量与飞机高度线所在平 面为 XZ 平面，飞机速度方向与 X 轴夹角定义为俯冲角 ，飞行高度为 h ，飞行速度为均速V 。散射单元在坐标系中的方位角和俯仰角分别为 和 ，其余飞机速度矢量方向的夹角为来自距离 R0 距离波门杂波回波信号在频率 f0 的功率谱密度 c (R0 , f0 ) 等于满足距离模糊和多普勒模糊多个距离-多普勒单元的功率谱密度
9、之和。P 2 c (R0 , f0 ) =G (Tn ,Tm ) + G (Tn ,) D(R, f )r 22(4 )3 Lfm n1 0Tm m n2 RD(Rm , fn ) = 4 dARm R Rm + R;fn f f0 + fPr 为平均发射功率； 为发射波长； G(T ,T ) 为天线增益；T 、T 为杂波散射单元在天线 坐标下的方位角和俯仰角； L 为系统损耗； 0 为杂波后向散射系数。 0为计算 D(Rm , fn ) ，定义 C(r, ) 为：C(r, ) = AdAh R r; f R4其物理含义为内层的灯多普勒曲线与外层的距离圆所相交的面积，如图 所示。 D(R ,
10、f ) = 1 C (R+ R, f ) C (R, f ) C (R+ R, f+ f ) + C (R , f+ f )m n 2m n m n m n m n文献5利用平面直角坐标系 ( X ,Y , Z ) 与球坐标系 (R, , ) 之间的转换关系,对公式进行推导给出了关于 的一维积分形式闭合解：rL( ) = 0 cos sin sin 1 cos2 cos1 cos (sin sin( ) (sin( + ) sin ) 22 cos cos+ sin sin 2 cos1 sin sin cos cos sin 从而 C(r, ) 为 L( A ) L(r ) 的结果。即可得杂
11、波功率谱的闭合解。4.程序实现及仿真结果分析用 Matlab 程序实现任意姿态机载 PD 雷达通用杂波闭合算法的仿真，仿真主要参数设定 雷达参数发射功率、入射波长、收发综合损耗、发射脉宽、脉冲重复周期、FFT 滤波器组字长、距离波门个数载机参数载机高度、载机速度大小、载机偏航角、载机俯仰角、载机横滚角天线模型天线最大电压增益、3dB 功率点波束宽度、天线偏航角、天线俯仰角地面散射模型参数采用修正的地面散射模型图 3俯仰角对杂波功率谱的影响频率不模糊时，雷达照射地面的最远距离为 R = h / sin(T 3dB / 2) 而T 与天线的俯仰角 A 、载机的俯仰角 P 、以及散射单元的俯仰角 有
12、关，当角度较大时距离小，不会产生距离模糊。如图 3（1）所示可清晰分出主瓣杂波、旁瓣杂波和高度线杂波。随着角度减小，距离越大出现部分模糊，如图 3（2）所示。当角度接近 3dB 波束宽度时，主波束照射范围 内完全距离模糊，如图 3（3）所示，每个距离波门内杂波功率都基本一致，都存在主波束 杂波功率。图 4偏航角和方位角对杂波功率谱的影响天线偏航角 A 、载机偏航角 P 、和散射单元方位角 的变化会引起主瓣杂波位置的变DD化，主瓣杂波中心频率 f0 = (2v / ) (cos P cosT sin T + sin P sin T ) 图 4 给出 A = 5 ； A = 5 、D D P = 5 ； = 5 时的仿真图，杂波主瓣杂波中心频率由原来的 0KHz 变为-23KHz；-46KHz；23KHz。从仿真图可以看出符合理论结果因此，无论是频率维上还是距离维上，仿真结果都与理论分析一致。参考文献
