一种天地波混合雷达回波谱仿真方法
【专利摘要】本发明公开了一种天地波混合雷达回波谱仿真方法。本方法它包括以下步骤：天地波雷达目标定位计算；传播衰减计算；噪声背景下雷达回波谱仿真；海杂波背景下雷达回波谱仿真。本发明从目标定位模型出发，完整地表示出了天地波混合雷达完整过程：自天波雷达发射高频电磁波直至地波接收站接收海面回波以及在这个过程中所遇到的上述噪声、衰减、杂波，最后的回波谱则表明了该新型体制雷达真的能够检测出某一海面上给定的目标；分析了该类型雷达的探测性能，为以后研究新的处理算法，指导阵地建设提供一定的参考。
【专利说明】一种天地波混合雷达回波谱仿真方法
【技术领域】：
[0001] 本发明涉及的是一种天地波混合雷达回波谱仿真方法。
【背景技术】：
[0002] 高频超视距雷达（HF-0THR)按其电波传播模式可划分为天波超视距雷达和地波 超视距雷达。天波超视距雷达（简称天波雷达）虽然可实现1000?4000公里的探测，但 受到了发射天线仰角以及电离层电子浓度的限制，天波雷达会存在700?1000公里内的近 距盲区以及电离层扰动影响。地波超视距雷达（简称地波雷达）对空中/海面目标的作用 距离可达100?300公里，且传播路径比较稳定，但受到其工作机理的限制，发射机通常部 署在沿海地区，其安全性大大降低，抗摧毁能力也较差。为此，雷达研究者们提出了一种天 地波混合雷达，这种雷达将天波电离层反射传播模式和地波绕射传播模式相结合，兼有两 种体制雷达的优点。但其目标定位相对复杂，且传播模式多样，增加了其信号处理算法和处 理步骤的复杂程度。因此，需要通过计算机仿真雷达回波谱，以解决上述问题。
【发明内容】
：
[0003] 本发明的目的在于克服现有技术的不足，提供一种天地波混合雷达回波谱仿真方 法。根据仿真的回波谱，研究如何提高新体制雷达的性能和适应该雷达的信号处理算法。
[0004] 为了解决【背景技术】所存在的问题，本发明采用以下技术方案：
[0005] -种天地波混合雷达回波谱仿真方法，包括以下步骤：
[0006] 步骤一：天地波雷达目标定位计算；
[0007] 步骤二：传播衰减计算；
[0008] 步骤三：噪声背景下雷达回波谱仿真；
[0009] 步骤四：海杂波背景下雷达回波谱仿真。
[0010] 进一步的，在步骤一中，雷达目标定位计算采用的公式是：
[0011] P3 = P/ -2P!
[0012] 其中，匕为目标到雷达接收机的距离；P'为雷达发射的电磁波经目标反射后被雷 达接收机接收时所传播的距离，其为可测量A为雷达发射机到电离层等效反射点的距离， 其值不可直接测得，但可由下面方程联立得出：
[0013] P/ = h2+(D/2)2
[0014] D2 = P32+D〇2-2P3D〇cos θ
[0015] 其中各个参数的含义：h为电离层等效的反射点，可以实际测得；D为雷达发射机 与目标连线距离，此量未知，但其满足式D 2 = P^+D/lPp^os Θ为雷达发射机和接收机 之间连线距离，为已知量；Θ为目标和雷达接收机连线与直线％之间的夹角，此量为已知 量，这样就可以根据式上述三个公式完成目标距离P 3的测量，实现目标的定位。
[0016] 进一步的，在步骤二中，传播衰减计算包括电离层吸收衰减和海浪附加衰减；
[0017] 对于电离层吸收衰减的计算采用的公式是：
[0018]
【权利要求】
1. 一种天地波混合雷达回波谱仿真方法，其特征在于，它包括以下步骤： 步骤一：天地波雷达目标定位计算； 步骤二：传播衰减计算； 步骤三：噪声背景下雷达回波谱仿真； 步骤四：海杂波背景下雷达回波谱仿真。
2. 根据权利要求1所述的一种天地波混合雷达回波谱仿真方法，其特征在于，在步骤 一中，雷达目标定位计算采用的公式是： P3 = P/ -2P! 其中，P3为目标到雷达接收机的距离^为雷达发射的电磁波经目标反射后被雷达接 收机接收时所传播的距离，其为可测量A为雷达发射机到电离层等效反射点的距离，其值 不可直接测得，但可由下面方程联立得出： P/ = h2+(D/2)2 D2 = P32+D02-2P3D 0cos θ 其中各个参数的含义：h为电离层等效的反射点，可以实际测得；D为雷达发射机与目 标连线距离，此量未知，但其满足式D2 = P^+D/lPp^os Θ为雷达发射机和接收机之间 连线距离，为已知量；Θ为目标和雷达接收机连线与直线％之间的夹角，此量为已知量，这 样就可以根据式上述三个公式完成目标距离P 3的测量，实现目标的定位。
3. 根据权利要求1所述的一种天地波混合雷达回波谱仿真方法，其特征在于，在步骤 二中，传播衰减计算包括电离层吸收衰减和海浪附加衰减； 对于电离层吸收衰减的计算采用的公式是：
其中各个参数的含义：
为在l〇〇Km高处的入射角；f为工作频率（MHz) ;fH为在 高为100km处的磁旋频率（MHz) ;1为吸收指数，其值可由下式求出：I = (1+0. 0037R) (cosO. 881 X ) i3,其中R为太阳黑子数在12个月的流动平均值；X为100km吸收处所对应 太阳天顶角；考虑到电离层存在的法拉第旋转效应，根据电离层吸收衰减的计算公式计算 的结果还有加上3dB ; 对于海浪附加衰减的计算方法是：由于海浪附加衰减公式比较复杂，往往采用利 用已有数据读值和线性插值计算的方式近似得出，若想得到工作频率4，海况为^的 附加衰减1^，首先根据读图确定工作频率4海况为 S(l的附加衰减Lwl以及工作频 率4海况为S(l的附加衰减L w2，使得& e (4￡2)，则所求频率的海浪附加衰减值为
4. 根据权利要求1所述的一种天地波混合雷达回波谱仿真方法，其特征在于，在步骤 三中噪声背景下雷达回波谱仿真包括目标信噪比计算以及回波谱生成； 目标信噪比计算公式为： SNR = PtGtGr λ 2 σ Ta/ [ (4 π ) 3P2P32kT〇FaL aLwLs] 其中Pt为天波发射站发射功率；P为天波路径P = 2Pi = ;Gt为发射机天线增益；La为 电离层吸收损耗；σ的目标后向散射面积；P3为地波传播路径；Lw为地波传播损耗因子；4 为地波接收站天线增益；λ为雷达工作波长；LS为系统损耗因子；k为波尔兹曼常数；Fa为 环境噪声系数；?；为噪声温度温度；T a为相干积累时间； 目标回波谱生成要得到目标的时域信号，第η个天线阵元接收到的目标相对于第m个 发射脉冲的回波信号为表达式为：
其中，
τ为回波信号相对于发射信号的时延；fT为运动目标的 Doppler频率；λ为电波波长；d为天线阵元间距；R为目标距离雷达接收机距离；Θ为目标 初始位置相对于接收天线阵列的法线方向的夹角为目标径向速度；noise^a)为均值 为〇方差为1的高斯白噪声；a为目标信号幅度a = SNRABXMXN) ;B、M和N分别为做方 位变换、距离变换和速度变换的点数，对于式（8)的信号在方位维做傅里叶变换，之后再进 行距离-速度的二维快速傅里叶变换，就可以得到某个方向某个距离单元上的回波数据。
5.根据权利要求1所述的一种天地波混合雷达回波谱仿真方法，其特征在于，在步骤 四中主要计算海杂波回波谱，包括一阶海杂波和二阶海杂波； 对于一阶海杂波，其所在位置的计算公式：
，其中土号表 示朝向及背离雷达波束的谐振海浪所产生的正负一阶海杂波，为雷达工作频率，单位为 MHz，fB的单位为Ηζ，β为雷达波束与水平面的夹角； 一阶海杂波中能量较大的一阶回波可以等效为散射面积为几万平米的目标，而两个一 阶海杂波峰能量比r是雷达波束方向与海面风向夹角Θ的函数，其理论公式如下：
，式中s为经验值，取值范围为2?6,-般取s = 4 ; 对于二阶海杂波的仿真较为复杂，首先根据理论公式得到理想的一二阶海杂波谱，为 了节省计算量，设置能量阈值来去除能量较小的二阶海杂波；为了获得连续的二阶谱区域， 舍弃间隔不大于频率分辨率的谱线，最后利用极值点的频率和幅度得到仿真出的回波谱。
【文档编号】G01S7/40GK104215946SQ201410390643
【公开日】2014年12月17日 申请日期:2014年8月8日 优先权日:2014年8月8日
【发明者】侯成宇, 陈迪夏, 添李达 申请人:哈尔滨工业大学