[发明专利]一种天地波混合雷达回波谱仿真方法

摘要

本发明公开了一种天地波混合雷达回波谱仿真方法。本方法它包括以下步骤天地波雷达目标定位计算；传播衰减计算；噪声背景下雷达回波谱仿真；海杂波背景下雷达回波谱仿真。本发明从目标定位模型出发，完整地表示出了天地波混合雷达完整过程自天波雷达发射高频电磁波直至地波接收站接收海面回波以及在这个过程中所遇到的上述噪声、衰减、杂波，最后的回波谱则表明了该新型体制雷达真的能够检测出某一海面上给定的目标；分析了该类型雷达的探测性能，为以后研究新的处理算法，指导阵地建设提供一定的参考。

主权项

一种天地波混合雷达回波谱仿真方法，其特征在于，它包括以下步骤：步骤一：天地波雷达目标定位计算；步骤二：传播衰减计算；步骤三：噪声背景下雷达回波谱仿真；步骤四：海杂波背景下雷达回波谱仿真；在步骤一中，雷达目标定位计算采用的公式是：P3＝P′‑2P1其中，P3为目标到雷达接收机的距离；P′为雷达发射的电磁波经目标反射后被雷达接收机接收时所传播的距离，其为可测量；P1为雷达发射机到电离层等效反射点的距离，其值不可直接测得，但可由下面方程联立得出：P12＝h2+(D/2)2D2＝P32+D02‑2P3D0cosθ其中各个参数的含义：h为电离层等效的反射点，可以实际测得；D为雷达发射机与目标连线距离，此量未知，但其满足式D2＝P32+D02‑2P3D0cosθ；D0为雷达发射机和接收机之间连线距离，为已知量；θ为目标和雷达接收机连线与直线D0之间的夹角，此量为已知量，这样就可以根据式上述三个公式完成目标距离P3的测量，实现目标的定位；在步骤二中，传播衰减计算包括电离层吸收衰减和海浪附加衰减；对于电离层吸收衰减的计算采用的公式是：其中各个参数的含义：为在100Km高处的入射角；f为工作频率(MHz)；fH为在高为100km处的磁旋频率(MHz)；I为吸收指数，其值可由下式求出：I＝(1+0.0037R)(cos0.881χ)1.3，其中R为太阳黑子数在12个月的流动平均值；χ为100km吸收处所对应太阳天顶角；考虑到电离层存在的法拉第旋转效应，根据电离层吸收衰减的计算公式计算的结果还有加上3dB；对于海浪附加衰减的计算方法是：由于海浪附加衰减公式比较复杂，往往采用利用已有数据读值和线性插值计算的方式近似得出，若想得到工作频率f0，海况为s0的附加衰减Lw，首先根据读图确定工作频率f1海况为s0的附加衰减Lw1以及工作频率f2海况为s0的附加衰减Lw2，使得f0∈(f1,f2)，则所求频率的海浪附加衰减值为在步骤三中噪声背景下雷达回波谱仿真包括目标信噪比计算以及回波谱生成；目标信噪比计算公式为：SNR＝PtGtGrλ2σTa/[(4π)3P2P32kT0FaLaLwLs]其中Pt为天波发射站发射功率；P为天波路径P＝2P1；Gt为发射机天线增益；La为电离层吸收损耗；σ的目标后向散射面积；P3为地波传播路径；Lw为地波传播损耗因子；Gr为地波接收站天线增益；λ为雷达工作波长；Ls为系统损耗因子；k为波尔兹曼常数；Fa为环境噪声系数；T0为噪声温度温度；Ta为相干积累时间；目标回波谱生成要得到目标的时域信号，第n个天线阵元接收到的目标相对于第m个发射脉冲的回波信号为表达式为：其中，τ为回波信号相对于发射信号的时延；fT为运动目标的Doppler频率；λ为电波波长；d为天线阵元间距；R为目标距离雷达接收机距离；θ为目标初始位置相对于接收天线阵列的法线方向的夹角；vR为目标径向速度；noisen,m(t)为均值为0方差为1的高斯白噪声；a为目标信号幅度a＝SNR/(B×M×N)；B、M和N分别为做方位变换、距离变换和速度变换的点数，对于Sn,m(t)表达式的信号在方位维做傅里叶变换，之后再进行距离‑速度的二维快速傅里叶变换，就可以得到某个方向某个距离单元上的回波数据；在步骤四中主要计算海杂波回波谱，包括一阶海杂波和二阶海杂波；对于一阶海杂波，其所在位置的计算公式：其中±号表示朝向及背离雷达波束的谐振海浪所产生的正负一阶海杂波，f0为雷达工作频率，单位为MHz，fB的单位为Hz，β为雷达波束与水平面的夹角；一阶海杂波中能量较大的一阶回波可以等效为散射面积为几万平米的目标，而两个一阶海杂波峰能量比r是雷达波束方向与海面风向夹角θ的函数，其理论公式如下：式中s为经验值，取值范围为2～6；对于二阶海杂波的仿真较为复杂，首先根据理论公式得到理想的一二阶海杂波谱，为了节省计算量，设置能量阈值来去除能量较小的二阶海杂波；为了获得连续的二阶谱区域，舍弃间隔不大于频率分辨率的谱线，最后利用极值点的频率和幅度得到仿真出的回波谱。