基于海杂波的低空电磁波传播特性估计方法技术

本发明专利技术涉及一种基于海杂波功率谱的低空电磁波传播特性估计方法，属于海洋环境特征及电磁波传播特性实时估计技术领域。主要用于解决海上大气折射率剖面参数估计与电磁波传播特性评估问题。区别于基于海杂波功率估计雷达探测效能的方法，该方法主要采用基于对实测海杂波功率谱估计和通过预设大气折射率剖面模型参数的电磁波传播扩散功率谱变化求解比对，反演实际大气折射率剖面、电磁波传播衰减量、求解大气折射率剖面参数后验估计值，得到大气折射率剖面图和电磁波传播衰减图。本发明专利技术工程实现简单，具有军事隐蔽性好、预报误差小、成本低、装备可嵌入加载等优点，主要适用于军事上和海洋监测中海上大气环境和电磁波传播衰减变化的测量与估计。

【技术实现步骤摘要】

本专利技术涉及的基于海杂波功率谱估计的低空大气折射率特征参数及电磁波传播特性评估实现方法，属于海洋环境特征及电磁波传播特性实时估计

技术介绍
对海洋大气特征及电磁波传播特性的检测与评估是实现对海监测与海上目标探测效能判断的重要条件。目前为了弥补正演方式通过测量海面一定高度处的实时气象参数，基于各大气折射率模型递推得出大气折射率剖面，再结合电磁波传播方程求解得出电磁波在空间的传播因子的方法装备实现困难、成本高、误差大等不足，提出了通过海杂波的功率反演电磁波传播特性的方法，该方法基于对实测海杂波的强度与距离数据变化分析通过遗传算法全局优化搜索大气折射率剖面特征参数，进而求解出电磁波传播衰减量。但是受海情、雷达工作方式等参量的影响海杂波随机性较强，同时影响海杂波的多参量之间存在复杂的相互依赖性，因此通过建立海杂波模型来获取具有一定可信度和精度的海杂波强度，需要经过较长的试验验证过程。因此，寻求其它实用、简便、并能够实时准确估计低空电磁波传播特性的技术成为日益重要的发展方向。
技术实现思路
本专利技术的目的在于针对目前基于气象数据、海杂波功率估计海洋环境参数及电磁波传播衰减特性工程实现困难、成本高、误差大等问题，提供了一种基于海杂波功率谱估计的采用正演加反演方式估计海洋环境特性参数及电磁波传播衰减特性的技术及其实现方法。该方法主要通过对海雷达回波数据经分离优解得到海杂波数据及其特征参数，采用相关图法进行功率谱估计；然后假设一种大气折射率模型，将该模型中的随机环境参数代入电磁波传播模型，结合电磁波距离扩散方程及海杂波模型计算出海杂波功率谱变化，模型计算出的海杂波功率谱与实际试验中测得的海杂波功率谱的谱峰值、谱峰位置及谱宽进行比对后，利用粒子群算法反演出实际的大气折射率剖面；最后对大气折射率剖面参数进行后验估计，验证反演准确率，并绘制出大气折射率剖面图和电磁波传播衰减曲线图。其特征在于I)基于雷达回波的海杂波特征参数分离优解技术特征在于对采集的雷达回波信号进行平滑滤波，以抑制雷达系统噪声；采用了通过傅里叶变换计算多普勒频移，并结合对雷达回波的空间分布、幅度统计特性、时间相关性、频谱特性以及时频二维特性的分析， 将海杂波数据从雷达回波数据中进行分离。对分离出的海杂波数据，本专利技术采用数字正交滤波、脉冲压缩和相位检波，形成海杂波视频数据，通过相关图法估计出海杂波功率随频率的变化。2)基于海杂波功率谱估计的海洋大气折射率剖面参数匹配寻优估计方法其特征在于根据实测的海杂波功率谱估计，结合电磁波传播模型、海洋环境物理模型及海杂波模型，利用粒子群算法对大气折射率剖面参数进行优化搜索，并实现后验估计，解决了基于气象数据、海杂波功率进行大气折射率剖面和电磁波传播衰减估计的多种不足。建立大气折射率模型时,为反映海洋环境的动态特性,对大气折射率模型在空间范围内进行细化，距离步进为200km,垂直高度为100m,高度步进为Im;然后,根据电磁波传播模型，利用随机环境参数计算出电磁波传播的衰减因子，再根据海洋环境物理模型及海杂波模型，估计出海杂波功率谱变化，并与实测的随距离离散化的海杂波功率谱进行比对， 在环境参数的取值范围内采用神经网络方法进行循环搜索，直至搜索出环境参数的最优解。其中本专利技术为了解决海杂波的随机性对反演准确度的影响，每30分钟利用实际采集的雷达回波数据，对海杂波模型进行一次修正；最后，通过奇异值分解对大气折射率剖面模型进行后验估计，检查基于最优解的大气折射率剖面与实测的大气折射率剖面的符合程度， 若两者误差较大，修正大气折射率剖面模型，重新反演。本专利技术技术实现软件平台为Visual Studio 2005,硬件平台为双核处理器、主频 2.8GHz、3GB内存，操作系统为Windows XP。在软件设计中，为了提高故障排查隔离率，算法可替换性以及工程的稳定性，本专利技术采用了模块化设计，对各模型算法进行了封装，以.DLL 动态数据链接库的方式进行调用。本专利技术与现有技术相比，其显著优点为能够实时、准确地估计海洋探测环境特性参数及电磁波传播特性；能够实现在对海雷达上进行嵌入式加载，除数据采集器外，无其它额外设备。该方法工程实现简单，在双核处理器上可以实现实时处理，军事隐蔽性好，成本低，且工程软件实现时高度集成化，故障排除能力强，稳定性较高。下面结合附图对本专利技术作进一步详细描述。附图说明图I是海杂波特征参数分离优解处理流程图2是大气折射率剖面参数反演流程图。具体实施方式下面结合附图对本评估技术的实现方法及具体实施步骤进行详细说明。I)实测一定方位范围内的海杂波功率谱变化Pcx(W)；本专利技术中在估计海杂波功率谱时，首先需要确定待反演的方位范围，并利用数据采集器采集雷达回波数据；其次需要将海杂波特征参数从采集的雷达回波数据中进行分离优解。本专利技术设计了一种有多个系数的低线性相位有限脉冲相应(FIR)滤波器，对原始雷达回波数据进行平滑滤波。特别注意的是在设计该滤波器时，需确保滤波器的带宽大得足以不衰减海杂波数据。然后本专利技术通过傅里叶变换计算多普勒频移，结合对杂波的空间分布、幅度统计特性、时间相关性、频域特性以及时频二维特性分析的方法，将该方位范围内的海杂波数据从各种信号数据中分离出来。其中本专利技术采用双线性变换维格纳-威利分布方法进行时频分析；最后对海杂波数据进行数字正交滤波、脉冲压缩和相位检波，形成海杂波视频数据，通过相关图法估计出各个方位上海杂波功率谱。2)建立大气折射率模型，将环境参数m映射成不同高度处的修正折射率M ；目前对于环境模型的描述主要有两种方式三段式折射率剖面模型(适用于表面波导)和四段式折射率剖面模型(适用于贴地波导和悬空波导)。由于海上大气波导一般为蒸发波导，所以本专利技术中环境模型采用四段式折射率剖面模型。本专利技术在建立该模型时， 通过“莫宁-奥布霍夫”相似定理，一共采用11个参数对大气折射率模型进行描述，其中5 个参数用来描述大气折射率剖面垂直结构，6个参数用来描述大气折射率剖面水平结构。包括蒸发波导高度、反演层的厚度、陷获层基高、混合层斜率和反演层的M增量，对应环境参数m= ，将大气折射率剖面在垂直结构上分为4层。3)利用修正折射率M，结合电磁波传播模型和海洋环境物理模型，计算雷达回波功率随距离的变化P ( )，代入海杂波后向散射系数，通过相关图法计算海杂波的功率谱；目前典型的研究电磁波传播衰减的方法有三种，分别是抛物线方程法、几何光学理论以及波导模理论。本专利技术中将抛物线方程法与其它两种方法有机结合，弥补抛物线方程法忽略后向散射系数，对大于与水平线成20°的传播夹角，解的相位不够精确等不足。求解时需限定垂直空间的最大高度和仰角的最大度数0_。利用傅里叶分步积分对电磁波传播时空间场分布进行求解，求解时在距离上以Ax为步长步进，Ax按雷达距离分辨率的倍数进行选取，一般取30 100m，取值越小，计算电磁波的衰减因子时间越长，同时计算结果越精确。本专利技术中海杂波后向散射系数的计算采用对GIT经验模型进行修正的方法，针对海杂波的随机性，每30分钟采集一次回波数据，对模型进行修正。其中有效浪高采用信噪比结合浮标测得的有效浪高拟合得出。本专利技术中计算海杂波功率谱时，对海杂波功率随距离的变化进行自相关序列估计，再进行DFT从本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：王玲玲，王向敏，徐超，
申请(专利权)人：中国船舶重工集团公司第七二四研究所，
类型：发明
国别省市：