CN104280621B - 一种区域电磁环境射线传播路径并行搜索方法 - Google Patents
一种区域电磁环境射线传播路径并行搜索方法 Download PDFInfo
- Publication number
- CN104280621B CN104280621B CN201410592670.7A CN201410592670A CN104280621B CN 104280621 B CN104280621 B CN 104280621B CN 201410592670 A CN201410592670 A CN 201410592670A CN 104280621 B CN104280621 B CN 104280621B
- Authority
- CN
- China
- Prior art keywords
- ray
- region
- scan
- propagation
- paths
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Measuring Magnetic Variables (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种区域电磁环境射线传播路径并行搜索方法,是将计算的电磁环境区域进行栅格化处理,同时并行搜索从辐射源出发的所有射线传播路径,统计到达栅格的所有路径,获得整个区域的场强分布。本发明的积极效果是:在完成区域中一个小栅格点场强计算的同时,完成对整个区域场强分布的计算,大大降低了计算运算量;对区域中所有射线的路径搜索可以同时并发执行,非常适合于通过cuda、OpenCL、openMP、Mpi、MapReduce等并行计算架构进行加速;将对整个区域场强分布的计算量简化为对区域中一个小栅格的场强计算量,扩大区域大小所带来的附加运算量将大大降低,不再成为限制区域电磁计算的瓶颈。
Description
技术领域
本发明属于电磁计算技术领域,具体涉及一种区域电磁环境射线传播路径并行搜索方法。
背景技术
在预测诸如城市小区、电子靶场等区域的电波传播特性方面,应用于区域射线传播路径搜索的射线跟踪法得到了广泛的应用,它将电磁波在路径上的传播简化为直射、反射和绕射,确立每条传播路径上的衰减损耗,再根据各路径对场的贡献得出总的场强。
为获得整个区域的电磁场强分布,传统的射线跟踪法对全区域所有射线传播路径进行一次搜索仅能获得区域中一个点的场强,即对区域中一个点的场强计算是通过遍历搜索完整个区域的射线路径而确定的,要获得整个区域的场强分布就要对整个区域中的射线传播路径进行多次串行遍历搜索,这将产生很大的运算量,不能满足对区域电磁环境中射线传播路径快速搜索的需求。
发明内容
为了克服现有技术的上述缺点,本发明提供了一种区域电磁环境射线传播路径并行搜索方法,在完成区域中一个点的场强计算的同时,并行完成对整个区域场强分布的计算,将对整个区域电磁环境计算的运算量缩减为对区域中一个点场强计算的运算量,实现对区域电磁环境中射线传播路径的快速搜索,降低计算整个区域场强分布的运算量。
本发明解决其技术问题所采用的技术方案是:一种区域电磁环境射线传播路径并行搜索方法,包括如下步骤:
步骤一、对电磁环境区域进行栅格化处理;
步骤二、并行搜索从辐射源出发的所有射线传播路径;
步骤三、统计到达栅格的所有路径,将到达每个栅格点的每一条路径相应的场强值进行叠加,得到该栅格点的场强,进而获得整个区域场强分布。
步骤二所述并行搜索从辐射源出发的所有射线传播路径的方法为:
(1)确定射线传播方向:从辐射源到区域边界的一个栅格点作为一条电波传播射线;
(2)同步并行地对所有射线传播路径进行第一轮跟踪搜索,直至满足搜索终止条件时,结束第一轮跟踪搜索,并相应地记录跟踪射线到达每个栅格点的路径衰减值;
如果第一轮跟踪搜索产生了多条射线分支则执行第二轮跟踪搜索,以至第i(i=3,4,5,…)轮跟踪搜索,每一轮跟踪搜索过程同第一轮跟踪搜索,直至不再产生多条射线分支,则完成对整个区域电磁波传播射线的路径搜索;其中,第二轮到第i轮都是针对上一轮新产生的射线分支传播路径进行跟踪搜索的。
所述搜索终止条件是指下述三个条件中的任意一个条件:
1)射线传播路径衰减为最小阈值;
2)射线抵达区域边界;
3)射线碰到障碍物产生多条射线分支。
与现有技术相比,本发明的积极效果是:
1)计算效率高。本发明在完成区域中一个小栅格点场强计算的同时,完成对整个区域场强分布的计算,大大降低了计算运算量,提高了计算效率。
2)算法并行度高。本发明对区域中所有射线的路径搜索可以同时并发执行,非常适合于通过cuda、OpenCL、openMP、Mpi、MapReduce等并行计算架构进行加速。
3)可扩展性好。本发明将对整个区域场强分布的计算量简化为对区域中一个小栅格的场强计算量,扩大区域大小所带来的附加运算量将大大降低,不再成为限制区域电磁计算的瓶颈。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是本发明方法的流程框图。
具体实施方式
本发明提出的一种区域电磁环境射线传播路径并行搜索方法,是将计算的电磁环境区域进行栅格化处理,同时并行搜索从辐射源出发的所有射线传播路径,统计到达栅格的所有路径,获得整个区域的场强分布。具体方法如图1所示,包括如下步骤:
步骤一、区域栅格化处理:
根据研究区域大小、区域环境(地形、建筑物等)、计算能力等综合因素确定基准栅格尺寸,并按照基准栅格尺寸对整个区域进行栅格划分。
步骤二、并行搜索从辐射源出发的所有射线传播路径:
(1)确定射线传播方向:
从辐射源出发,到研究区域的一个栅格接收点为一个射线传播方向,选取区域边界栅格点作为从辐射源出发的所有射线传播的终点,即从辐射源到区域边界的一个栅格点确立一条电波传播射线。
(2)同步并行地对所有射线传播路径进行第一轮跟踪搜索,直至:
1)射线传播路径衰减为最小阈值;
2)射线抵达区域边界;
3)射线碰到障碍物产生多条(至少两条)射线分支;
上述三个条件中任意一个条件满足时,结束第一轮跟踪搜索,并相应地记录跟踪射线到达每个栅格点的路径衰减值。
如果第一轮跟踪搜索产生了多条射线分支则执行第二轮跟踪搜索,以至第i轮(i=3,4,5,…)跟踪搜索,每一轮跟踪搜索过程同第一轮跟踪搜索,直至不再产生多条射线分支,则完成对整个区域电磁波传播射线的路径搜索。其中,第二轮到第i轮计算都是对上一轮新产生的射线分支进行计算。
假定整个研究区域被划分为M个小栅格,传统的区域电磁环境计算对区域中的所有栅格不是同步并行计算的,遍历区域中的传播射线一次仅能确立一个小栅格的场强值,要获得整个区域的场强分布,就要执行M次对区域中传播射线的遍历,若完成一个小栅格场强计算的运算量为N,则对整个区域电磁场强的运算量为M×N;而本发明的区域电磁环境路径并行搜索计算方法,通过对整个区域中的射线传播路径进行搜索遍历,在得到区域中一个小栅格的场强值的同时,也获得了整个区域的场强分布,将运算量缩小为原来的1/M,大大提升了区域电磁环境计算效率。
步骤三、统计到达栅格的所有路径,获得整个区域的场强分布:
统计到达每个栅格点的所有路径,将到达每个栅格点的每一条路径相应的场强值进行叠加,得到该栅格点的场强,进而获得整个区域场强分布。对每一个栅格点场强的计算公式如下:
其中,Ei是到达当前计算栅格点的第i条射线确立的电场强度。
Claims (1)
1.一种区域电磁环境射线传播路径并行搜索方法,其特征在于:包括如下步骤:
步骤一、对电磁环境区域进行栅格化处理;
步骤二、并行搜索从辐射源出发的所有射线传播路径:
(1)确定射线传播方向:从辐射源到区域边界的一个栅格点作为一条电波传播射线;
(2)同步并行地对所有射线传播路径进行第一轮跟踪搜索,直至满足如下任一搜索终止条件:
1)射线传播路径衰减为最小阈值;
2)射线抵达区域边界;
3)射线碰到障碍物产生多条射线分支;
(3)结束第一轮跟踪搜索,并相应地记录跟踪射线到达每个栅格点的路径衰减值;
如果第一轮跟踪搜索产生了多条射线分支则执行第二轮跟踪搜索,以至第i(i = 3,4,5,…)轮跟踪搜索,每一轮跟踪搜索过程同第一轮跟踪搜索,直至不再产生多条射线分支,则完成对整个区域电磁波传播射线的路径搜索;其中,第二轮到第i轮都是针对上一轮新产生的射线分支传播路径进行跟踪搜索的;
步骤三、统计到达栅格的所有路径,将到达每个栅格点的每一条路径相应的场强值进行叠加,得到该栅格点的场强,进而获得整个区域场强分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410592670.7A CN104280621B (zh) | 2014-10-29 | 2014-10-29 | 一种区域电磁环境射线传播路径并行搜索方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410592670.7A CN104280621B (zh) | 2014-10-29 | 2014-10-29 | 一种区域电磁环境射线传播路径并行搜索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104280621A CN104280621A (zh) | 2015-01-14 |
CN104280621B true CN104280621B (zh) | 2017-03-15 |
Family
ID=52255696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410592670.7A Active CN104280621B (zh) | 2014-10-29 | 2014-10-29 | 一种区域电磁环境射线传播路径并行搜索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104280621B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222410B (zh) * | 2019-05-30 | 2020-12-11 | 北京理工大学 | 一种基于Hadoop MapReduce的电磁环境仿真方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592690A (zh) * | 2009-05-05 | 2009-12-02 | 上海大学 | 基于射线跟踪法的电磁波传播预测方法 |
JP4676157B2 (ja) * | 2004-03-31 | 2011-04-27 | 古河電気工業株式会社 | 数値計算方法 |
CN102546039A (zh) * | 2010-12-20 | 2012-07-04 | ***通信集团北京有限公司 | 一种无线电磁波传播预测方法及装置 |
CN102981064A (zh) * | 2012-10-09 | 2013-03-20 | 中国人民解放军63892部队 | 飞机外部射频电磁环境预测方法及预测*** |
CN103149457A (zh) * | 2013-03-01 | 2013-06-12 | 西安电子科技大学 | 空间电磁强度分布分析方法 |
CN103987077A (zh) * | 2013-02-08 | 2014-08-13 | ***通信集团广东有限公司 | 无线网络仿真方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005227202A (ja) * | 2004-02-16 | 2005-08-25 | Mitsubishi Electric Corp | 電磁界特性解析装置、電磁界特性解析方法および電磁界特性解析プログラム |
WO2012159330A1 (zh) * | 2011-07-15 | 2012-11-29 | 华为技术有限公司 | 电磁波传播特性的分析方法和装置 |
-
2014
- 2014-10-29 CN CN201410592670.7A patent/CN104280621B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4676157B2 (ja) * | 2004-03-31 | 2011-04-27 | 古河電気工業株式会社 | 数値計算方法 |
CN101592690A (zh) * | 2009-05-05 | 2009-12-02 | 上海大学 | 基于射线跟踪法的电磁波传播预测方法 |
CN102546039A (zh) * | 2010-12-20 | 2012-07-04 | ***通信集团北京有限公司 | 一种无线电磁波传播预测方法及装置 |
CN102981064A (zh) * | 2012-10-09 | 2013-03-20 | 中国人民解放军63892部队 | 飞机外部射频电磁环境预测方法及预测*** |
CN103987077A (zh) * | 2013-02-08 | 2014-08-13 | ***通信集团广东有限公司 | 无线网络仿真方法及装置 |
CN103149457A (zh) * | 2013-03-01 | 2013-06-12 | 西安电子科技大学 | 空间电磁强度分布分析方法 |
Non-Patent Citations (3)
Title |
---|
一种基于三维射线弹跳法的并行电磁传播预测算法;杨锦辉 等;《装备指挥技术学院学报》;20101231;第21卷(第6期);第91-96页 * |
基于地形网格剖分的射线追踪法的时间优化;刘艳梅 等;《环境技术》;20140430(第S1期);第95-98页 * |
无线通信中电波传播和场强预测的研究;刘海涛;《万方学位论文》;20060918;第97-104,122页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104280621A (zh) | 2015-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102997928B (zh) | 一种城际路网索引和匹配方法 | |
CN107610487B (zh) | 基于动态随机车流相位差协调机制的区域交通控制***和方法 | |
CN105116913A (zh) | 植保无人机作业航线规划方法及装置 | |
CN104200042B (zh) | 一种定位数据轨迹跟踪方法 | |
CN103236169B (zh) | 城市干线公交站间双向分段绿波信号设置方法 | |
CN109190813A (zh) | 一种基于双层规划的共享自行车区域投放规划方法 | |
CN108648444B (zh) | 一种基于网格模型的信号交叉口运行评价方法 | |
CN103747262B (zh) | 一种基于gpu的运动估计方法 | |
CN103268014A (zh) | 一种具有四个方向隐身效果的柱状电磁波隐身器件 | |
CN104280621B (zh) | 一种区域电磁环境射线传播路径并行搜索方法 | |
CN107341284A (zh) | 高精度预测低频电波传播特性的双向抛物方程方法 | |
CN104007762B (zh) | 一种电力巡检机器人的导航方法 | |
CN104820859A (zh) | 基于预处理和自适应遗传模拟退火算法的相位解缠方法 | |
CN105974175A (zh) | 一种布置埋地管道局部接地点的方法及装置 | |
CN105528753A (zh) | 一种基于网络病毒传播模型的空气污染溯源方法 | |
CN110188430A (zh) | 交流330kV悬垂塔分段式复合绝缘子均压结构优化设计方法 | |
CN103513350A (zh) | 一种基于粒子群算法的阵列波导器件对准耦合方法与装置 | |
Tian et al. | Effectiveness of lead–lag phasing on progression bandwidth | |
Tarek et al. | Fog computing for optimized traffic control strategy | |
CN104809307A (zh) | 面向曼哈顿空间的发动机直角管路遗传算法规划方法 | |
CN105205299B (zh) | 电大目标电磁散射特性快速降维分析方法 | |
CN103198485B (zh) | 一种基于方向权重的红外弱小目标检测前跟踪方法 | |
CN106842183B (zh) | 一种基于全局最近邻法的波位间点迹合并的方法及*** | |
CN105303022A (zh) | 快速获取目标电磁散射特性的高斯波束方法 | |
Zhao et al. | Improved ant colony optimization algorithm and its application for path planning of mobile robot in 3-D space |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |