基于三维抛物方程的不规则地形电波传播因子预测方法技术

本发明专利技术提出了一种基于三维抛物方程的不规则地形电波传播因子预测方法，用于解决二维抛物方程无法考虑横向地形对电波传播特性的影响问题，实现步骤为：采用数字高程地图建立三维不规则地形的几何模型；采用麦克斯韦方程推导出以位函数表示的宽角抛物方程；由激励源的电流分布函数求取宽角抛物方程的初始场；采用改进的分步傅里叶变换方法求解宽角抛物方程，得到宽角抛物方程的总场；构建整个空间内的衰减函数，并对宽角抛物方程总场的幅度进行修正；采用考虑吸收边界的宽角抛物方程总场，计算不规则地形电波传播因子。本发明专利技术提高了三维抛物方程的预测精度，且拓展了数字高程地图的应用范围，可用于复杂环境下的无线通信和GPS定位等。

【技术实现步骤摘要】
基于三维抛物方程的不规则地形电波传播因子预测方法
本专利技术属于无线电
，涉及一种电波传播因子预测方法，具体涉及一种基于三维抛物方程的不规则地形电波传播因子预测方法，可用于复杂环境下的无线通信和GPS定位等。
技术介绍
随着科学技术的飞速发展，无线通信已经越来越广泛地应用于军事、工业、教育等各个领域，给人们带来的影响是无可争议的。众所周知，无线通信系统都是由一段或多段无线电波传播电路实现的，因此，在无线通信领域中，电波传播特性的预测一直以来都是研究热点，获得人们愈来愈多的关注。在对流层电波传播中，复杂环境是影响电波传播特性的主要因素之一。复杂环境中的电波传播通常表现出反射、折射、绕射和散射等不同的传播机理，而地形地貌的影响主要表现为直达波与地面反射波的多径干涉效应，起伏地表对电波的绕射效应等。考虑到地形地貌的几何形态和物理特征的多样性，准确预测不规则地形上的电波传播特性显得尤为重要。目前，对对流层中的电波传播因子的预测均是基于传播模型实现的，传播模型可分为两类：一类是经验模型，例如Hata模型和COST231模型，通过对大量实测数据进行拟合，得到拟合公式，快速地给出特定场景下电波传播特性的统计结果。但是，由于这些模型不考虑传播环境的具体信息，故预测结果的准确度通常比较低，而且不具有通用性。另一类是建立在电磁波理论的基础上的理论模型，这种理论模型不仅基础牢固，还可以充分考虑传播环境的几何及形态特征，预测结果更准确，且更具普适性。在理论模型中，抛物方程模型由于能够快速求解大区域电波传播问题而被广泛关注，也逐渐成为准确预测复杂环境下电波传播特性的重要手段。抛物方程模型从波动方程出发推导出抛物方程，并采用分步傅里叶变换方法对其进行求解，将初始场和边界条件考虑在内，最终得到传播因子的预测结果，例如授权公告号为CN102722618B，名称为“一种基于抛物方程的准三维电磁环境模型构建与并行方法”的专利技术专利，公开了一种基于抛物方程的准三维电磁环境模型构建与并行方法，该方法主要采用切片法，即以天线发射机为中心，按照一定角度将三维空间剖分成单独的二维切片，把求解复杂的三维电磁环境仿真问题简化成了二维仿真问题，接下来从天线方向图出发构建准三维方法的初始场，并对复杂地理环境和大气结构的复杂环境进行准三维剖分，然后将剖分后的二维切片采用粗粒度的并行计算方法对求解过程进行了加速。该方法不需要进行实际的测量就可以得到空间中的电磁场分布，其计算复杂度大大减少，计算范围有了较大程度的扩展。但是在该方法中，二维切片之间是独立的，且二维抛物方程模型由于无法考虑横向的地形对电波传播的影响而导致预测精度受限，因此，开展三维不规则地形条件下的三维抛物方程方法的研究十分必要，对复杂环境下的无线通信具有重要价值。
技术实现思路
本专利技术的目的在于克服上述现有技术存在的不足，提出了一种基于三维抛物方程的不规则地形电波传播因子预测方法，从矢量波动方程出发推导出三维抛物方程，并通过改进的分步傅里叶变换方法对其进行求解，充分考虑横向地形对电波传播特性的影响，实现了对三维不规则地形上电波传播因子的精准预测。为实现上述目的，本专利技术采取的技术方案包括如下步骤：(1)在预测区域内建立直角坐标系，采用数字高程地图建立三维不规则地形的几何模型，获取数字高程地图建模文件，实现步骤为：(1a)采用二维非均匀网格对三维不规则地形进行覆盖，得到包含各网格节点处高度值的三维不规则地形几何模型；(1b)按照非规则格网数字高程地图的文件存储格式，输出三维不规则地形几何模型中各网格节点处的高度值，形成数字高程地图建模文件；(2)采用麦克斯韦方程推导出以位函数表示的宽角抛物方程，实现步骤为：(2a)对麦克斯韦方程进行推导，得到三维波动方程；(2b)根据三维波动方程中的电场和磁场的场分量，设定位函数电矢位和磁矢位并将其代入三维波动方程，得到由位函数表示的矢量波动方程；(2c)用ψe表示磁矢位的任意直角分量，用ψm表示电矢位的任意直角分量，并设定ψ＝ψe+ψm为任意标量场分量；(2d)忽略后向传播，采用Feit-Fleck近似方法计算微分算子，推导出宽角抛物方程：其中，为偏微分算子，x、y和z分别为直角坐标系的三个分量，i为虚数单位，k0为自由空间传播常数，为横向拉普拉斯算子，n为传播介质的折射率；(3)将不规则地形表面看作PEC平面，采用镜像原理，由激励源的电流分布函数求取宽角抛物方程的初始场ψ(0,y,z)；(4)对分步傅里叶变换方法进行改进，并采用改进的分步傅里叶变换方法求解宽角抛物方程，得到宽角抛物方程的总场ψ(x,y,z)，实现步骤为：(4a)将不规则地形对电波传播的影响等效为一系列单刃峰的绕射，并将不规则地形表面以下区域内宽角抛物方程的总场ψ(x,y,z)设为零，得到改进的分步傅里叶变换方法；(4b)从数字高程地图建模文件中读取三维不规则地形几何模型；(4c)忽略地表上的爬行波，根据初始场ψ(0,y,z)，采用改进的分步傅里叶变换方法求解宽角抛物方程，得到宽角抛物方程总场ψ(x,y,z)的二维傅里叶变换；(4d)采用奇偶分解法，通过奇场和偶场的二维傅里叶变换表示宽角抛物方程的总场ψ(x,y,z)的二维傅里叶变换，并利用FFT求出其逆变换，得到宽角抛物方程总场ψ(x,y,z)；(5)构建整个空间内的衰减函数，并对宽角抛物方程总场ψ(x,y,z)的幅度进行修正，实现步骤为：(5a)在直角坐标系的y方向和z方向各设置一个Tukey窗函数，并将这两个Tukey窗函数相乘，得到宽角抛物方程在(y,z)平面内的衰减函数；(5b)采用衰减函数对整个空间内的宽角抛物方程总场ψ(x,y,z)的幅度进行修正，得到考虑吸收边界的宽角抛物方程总场ψs(x,y,z)；(6)采用考虑吸收边界的宽角抛物方程总场ψs(x,y,z)，计算不规则地形电波传播因子F。本专利技术与现有技术相比，具有如下优点：(1)本专利技术在预测三维不规则地形上电波传播因子时，通过改进的分步傅里叶变换方法对三维抛物方程进行求解，避免了现有技术中采用二维抛物方程模型无法考虑横向地形对电波传播的影响，有效地提高了电波传播因子的预测精度。(2)本专利技术采用数字高程地图对三维不规则地形进行几何建模，解决了传统数字地图无法描述不规则起伏平面的问题，拓展了抛物方程模型的应用范围。(3)本专利技术由于可以考虑任意形状的三维不规则地形，而且不影响宽角抛物方程的求解效率，具有更高的通用性和更广的工程应用范围。附图说明图1是本专利技术的实现流程图；图2是本专利技术在平坦地面上预测垂直面和水平面内电波传播因子的仿真图；图3是本专利技术采用数字高程地图建立的三维不规则地形的几何模型图；图4是本专利技术在三维高斯起伏地形上预测垂直面和水平面内电波传播因子的仿真图；图5是本专利技术与现有技术电波传播因子预测结果的仿真对比图。具体实施方式以下结合附图和具体实施例，对本专利技术作进一步详细说明。参照图1，基于三维抛物方程的不规则地形电波传播因子预测方法，包括如下步骤：步骤1)在预测区域内建立直角坐标系，由于传统数字地图无法描述不规则起伏平面，故采用数字高程地图建立三维不规则地形的几何模型，获取数字高程地图建模文件，实现步骤为：步骤1a)在三维抛物方程模型中，x方向和y方向的剖分步长相差甚远本文档来自技高网...

【技术保护点】
一种基于三维抛物方程的不规则地形电波传播因子预测方法，包括如下步骤：(1)在预测区域内建立直角坐标系，采用数字高程地图建立三维不规则地形的几何模型，获取数字高程地图建模文件，实现步骤为：(1a)采用二维非均匀网格对三维不规则地形进行覆盖，得到包含各网格节点处高度值的三维不规则地形几何模型；(1b)按照非规则格网数字高程地图的文件存储格式，输出三维不规则地形几何模型中各网格节点处的高度值，形成数字高程地图建模文件；(2)采用麦克斯韦方程推导出以位函数表示的宽角抛物方程，实现步骤为：(2a)对麦克斯韦方程进行推导，得到三维波动方程；(2b)根据三维波动方程中的电场和磁场的场分量，设定位函数电矢位

【技术特征摘要】
1.一种基于三维抛物方程的不规则地形电波传播因子预测方法，包括如下步骤：(1)在预测区域内建立直角坐标系，采用数字高程地图建立三维不规则地形的几何模型，获取数字高程地图建模文件，实现步骤为：(1a)采用二维非均匀网格对三维不规则地形进行覆盖，得到包含各网格节点处高度值的三维不规则地形几何模型；(1b)按照非规则格网数字高程地图的文件存储格式，输出三维不规则地形几何模型中各网格节点处的高度值，形成数字高程地图建模文件；(2)采用麦克斯韦方程推导出以位函数表示的宽角抛物方程，实现步骤为：(2a)对麦克斯韦方程进行推导，得到三维波动方程；(2b)根据三维波动方程中的电场和磁场的场分量，设定位函数电矢位和磁矢位并将其代入三维波动方程，得到由位函数表示的矢量波动方程；(2c)用ψe表示磁矢位的任意直角分量，用ψm表示电矢位的任意直角分量，并设定ψ＝ψe+ψm为任意标量场分量；(2d)忽略后向传播，采用Feit-Fleck近似方法计算微分算子，推导出宽角抛物方程：其中，为偏微分算子，x、y和z分别为直角坐标系的三个分量，i为虚数单位，k0为自由空间传播常数，为横向拉普拉斯算子，n为传播介质的折射率；(3)将不规则地形表面看作PEC平面，采用镜像原理，由激励源的电流分布函数求取宽角抛物方程的初始场ψ(0,y,z)；(4)对分步傅里叶变换方法进行改进，并采用改进的分步傅里叶变换方法求解宽角抛物方程，得到宽角抛物方程的总场ψ(x,y,z)，实现步骤为：(4a)将不规则地形对电波传播的影响等效为一系列单刃峰的绕射，并将不规则地形表面以下区域内宽角抛物方程的总场ψ(x,y,z)设为零，得到改进的分步傅里叶变换方法；(4b)从数字高程地图建模文件中读取三维不规则地形几何模型；(4c)忽略地表上的爬行波，根据初始场ψ(0,y,z)，采用改进的分步傅里叶变换方法求解宽角抛物方程，得到宽...

【专利技术属性】
技术研发人员：郭立新，关晓伟，王亚姣，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61