CN115238737B 一种基于合作式航海雷达信号的蒸发波导剖面反演方法 （西北工业大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN115238737B(45)授权公告日2025.07.08

(21)申请号202210779365.3

(22)申请日2022.07.01

(65)同一申请的已公布的文献号申请公布号CN115238737A

(43)申请公布日2022.10.25

(73)专利权人西北工业大学

地址710072陕西省西安市友谊西路

(72)发明人杨坤德杨帆史阳张皓王淑文段顺利

(74)专利代理机构西安凯多思知识产权代理事务所(普通合伙)61290

专利代理师王鲜凯

(51)Int.CI.

GO6F18/24(2023.01)

GO6F17/12(2006.01)

G01S13/95(2006.01)

GO6F123/00(2023.01)

(56)对比文件

US2013279561A1,2013.10.24CN106707256A,2017.05.24

审查员吴海旋

权利要求书2页说明书6页附图4页

(54)发明名称

一种基于合作式航海雷达信号的蒸发波导剖面反演方法

(57)摘要

CN115238737B本发明涉及一种基于合作式航海雷达信号的蒸发波导剖面反演方法，采用时间累积、方向锁定、波束定时的方法，在已知目标舰船航海雷达参数的基础上，获取监测链路上的路径损失数据，通过建立监测链路上水平非均匀的蒸发波导剖面模型，并结合电磁波高级传播模型、目标函数、全局寻优算法、插值同化算法等，以水平均匀环境反演结果为初值，设置目标函数搜索阈值，反演得到监测链路上非均匀蒸发波导分布。合理利用舰船航海雷达的已知信息，实时被动接收射频信号，反演蒸发波导剖面，降低反演难度，提高反演效率，实现蒸发波导剖面的被动信号反演，

CN115238737B

发射端参数

发射端参数

目标船

导航市达信号

发射端参数

高级传播模型监测点

拷贝场

匹配场

目标函数

全局寻优算法循环次数=1

循环次数1

否

满足终止条件

输出非均匀蒸发波导剖面

初始M蒸发波导模型一

CN115238737B权利要求书1/2页

2

1.一种基于合作式航海雷达信号的蒸发波导剖面反演方法，其特征在于步骤如下：

步骤1、建立拷贝场数据库：以设定间隔，计算蒸发波导高度h=zED的大气修正折射率分布：

其中：M?、M(z)分别为海面和高度z处的大气修正折射率，M?取值330;h=zED为蒸发波导高度，粗糙度长度z?=1.5×10??,△M=M(h)-M?为蒸发波导强度，M(h?)为蒸发波导高度处的修正折射率；

以上得到的数据作为建立水平均匀的蒸发波导环境数据；

步骤2、路径损失仿真：依据目标舰船信息设定发射源仿真参数，代入到高级传播模型中，计算对应蒸发波导剖面的路径损失：

PL=[PL-20logFl

PL是自由空间损耗，单位是dB:PLs=32.45+20logf+20logr;

其中：f为频率，单位是MHz,r为发射天线和接收天线之间的距离，单位是km;

所述F是传播因子：F=lu(x,z)|√x

其中：u(x,z)|表示利用抛物方程模型计算得到的蒸发波导环境下的电磁场场强分布；

x表示传播距离，单位为m;

步骤3、计算实测路径损失：根据航海雷达自转特性，提取一个雷达自转周期内的最大监测信号电平及其出现的时间，结合目标船提供的位置信息、发射功率、天线高度、天线增益，计算监测链路的路径损失：

PL=[P+G′+G+G?-P-IL]

其中：t表示时间累积对应的时间参数，t=1,2,…,N;PL为实际测得的路径损失，单位为dB;P为发射信号功率，单位为dBm;G为发射天线增益，单位为dBi;G为接收天线增益，单位为dBi;G为接收信号放大器增益，单位为dB;P为接收信号的电平值，单位为dBm;IL为整个系统的等效系统插损，单位为dB;

步骤4、建立反演目标函数φ(p):

其中，p表示蒸发波导剖面模型参数，N在拷贝场数据库的搜索次数；

利用全局优化算法，搜索求解使φ(p)最小的一组最优解，反演得到单链路上的均匀蒸发波导修正折射率剖面；

步骤5、计算拷贝场向量：将反演结果即单链路上的均匀蒸发波导修正折射率剖面作为初值，在监测链路上建立水平不均匀的蒸发波导