Longley-Rice无线电波传输模型

1、目录目 录 i第一章 绪 论 11.1 研究背景 11.2 模型分类及参数 11.3 Longley-Rice 模型传输损耗 21.4 Longley-Rice 模型的适用围 3第二章 传播模型分析及建模 . 52.1 传播模型的分析与选择 52.2 Longley-Rice 模型建模 82.2.1 衍射传播损耗 9222视距(LOS)传播损耗 102.2.3 向散射传播损耗 102.3 仿真分析 102.3.2 SRTM 高程数据的抽取过程 112.3.3 接收机高度描述 122.3.4 结果 13第三章 Longley-rice 等效模型的建立 153.1 longley-rice 的现状

2、和不足 153.2 Longley-Rice 模型的改进 15第四章 利用 longley-rice 模型设计的可视化电磁环境 174.1 电磁环境可视化 174.2 三维电磁环境体数据生成 184.2.1 不规则地形剖面提取 184.2.3 电磁环境体数据计算 194.3 在三维数字地球上体绘制电磁环境 204.3.1 硬件加速的直接体绘制. 204.3.2 电磁环境体数据包围盒. 204.3.3 体数据纹理坐标转换 . 21第五章用 Longley-Rice 大气折射修法的改进 255.1 对 Okumura-Hata 和 Longley-Rice 的比较 255.2 Longley-Ri

3、ce 模型中对大气折射误差的修正 255.3 大气折射修法的改进 26第六章 展望 29参考文献 31第一章 绪 论1.1 研究背景Lon gley-Rice无线电波传输模型是由Lon gley和Riee提出的无线电波传输模型， 该模型为统计模型， 它是以无线电波传播理论为依据， 并结合了许多实际环境下的测量 数据，所以该模型被称为半经验预测模型。 正是因为该模型是以无线电波传播理论为依 据，有加上了大量的实际测量数据， 因此该模型得到了广泛的应用。 该模型被称为不规 则地面模型，可以用该模型预测自由空间中由地形非规则性所造成的中值传输衰落。1.2 模型分类及参数Longley-rice 模型

4、有 :2 种模式。当地形路径易据很详细时，特定路径参数就很容 易被确定。这不预测方式为“点到点预测” 。如果地形数据不够训确，可以利用 Longley-Rice 模型估计特定参数的值这种预测方式为“区域预测” 。Longley-rice 模型有 4 种变化模式，分别为单信号模式、单天线模式、移动模式和广 播模式。在 longley-rice 模型的早期程序中，所有点对点预测的计算都是基于单天线 模式，这里天线的位置是不变的。后来，由于对计算精度需求的提高，人们才引入其他 模式。在各种变化模式中，变化的主要是时间、位置和情景 3 个参量，或者说是一个三 维变量。目前， Longley-Rice

5、无线电波模型已有相关的计算机仿真程序，可以用来对无线 电波传输的损耗进行计算。 当无线电波传输路径已知时， 计算机的仿真程序可以通过无 线电波传播路径的长度、极化方向、无线电波频率、地面有效半径、收发天线高度、地 面导电常数以及表面绕射率等参数确定无线电波传输损耗的大小。Longley-Rice 预测模型主要有以下参数 : 天线极化方式 : 可以采用水平极化或者垂直极化。 Longley-Rice 模型中假设发射 天线和接受天线具有相同的极化方式 ; 折射率 : 空气的折射率决定了无线电波的“弯曲”程度。在一般的模型中，空气 折射率用地面有效曲率来代替，通常取 1.333 。 介电常数 : 地

6、面的相对介电常数和电导率1.3 Longley-Rice 模型传输损耗根据无线电波的传播围的不同， Longley-Rice 模型的传输损耗可大致分为三种情 况，它们分别为 :(1) 视距传播损耗(2) 绕射传播损耗(3) 散射传播损耗。当无线电波传播于视距围时， 无线电波传播方式主要以反射传播为主。 通过对地貌 地形的路径及对流层的绕射特点进行分析， 预测出无线电波传输中值传输衰落， 将双线 地面反射模型用来模拟地平线以的传输场强。视距传播模型的适用围为 : dmin d dLS 。(2) 绕射传播损耗 当无线电波传播于超视距围时， 无线电波传播方式主要以绕射传播为主。 但是，当 无线电波的

7、传播环境为不规则的地形时， 会有两种理论用于计算绕射损耗。 其中一种适 用于预测非球形但光滑的地面上无线电波的传播， 而另外一种则适用于非常不规则的地 面上无线电波的传播。采用 Fresnel-Kirchoff 刃形模型来预测无线电波传播的绕射损 耗，其计算结果是上述两种理论结果的加权和。适用围：dLs d dx为绕射传播距离，单位为km(3) 散射传播损耗 当无线电波的传输距离为大大超出地平线的远距离无线电波的传输时， 无线电波传 播方式主要以散射传播为主，预测理论以绕射理论为主。适用围 ： dx d 为散射传播距离，单位为 km。综上所述，Longley-Rice模型传输损耗Lb为：LbL

8、refL free(3-15)其中 :Lfree32.45 20lg d 20lg f( 3-16)Lrefmax( 0, L be k1d k2lgd),dmin d dLsLbed md ,dLs ddxLbes msd,dx d(3-17)式中:d为传播距离，单位为km;f为无线电波频率，单位为MHz; dLs为光滑地面距离；dx表示此处的绕射损耗和散射损耗相等；Lbe、Lbed、Lbes分别表示自由空间下视距、绕射和散射时的传播损耗值；ki和k2为传播损耗系数；md和ms分别为绕射和散射损耗系数。dmin d dLs为视距传播距离，dLs d dx为绕射传播距离，dx d为散射传播距

9、离。1.4 Longley-Rice 模型的适用围Longley-Rice 模型中的实测数据大多数取自10-1 000 MHz的频率围，其中20-100MHZ的数据涉及5-50km的距离和1-9 m的收、发信天线高度；较高频段的数据涉及5-1 000 km的距离，10-1500 m的发射天线高度和3-9m的接收天线高度。数据来源于世 界各地， 但主要是美国， 多数为移动记录结果。 2 Longley-Rice 模型给出了参考衰减值 的计算公式及不同环境下相关修正因子的详细说明，公式中所使用的参数包括 : 不规则 地形参数、频率、收发信机天线高度和表面折射率等 3 。同时还引入了反映介质特性的

10、2 个参数 : 介电常数和导电率。以传播理论为依据，加上极其丰富的实测数据，使得 Longley-Rice 模型使用围得到了拓展，其适用围如下:1) 频率 f : 20-40000MHz；2) 收、发信机天线高度 :0.5-3000m；3) 覆盖半径 :1-2000km； 44）表面折射率：250-400Ns 。表1.1地面的相对介电常数和导电率相对介电常数导电率/（S/M）一般的地面156 005导电性差的地面40.001导电性好的地血250. 020纯水810. 010815. 000表1.2气候类型和相应地表折射率气険类里览表折疥扌率M （K-units）1近赤迪气幄ciih刚果）360

11、2亚热带大陆性气解cftii；苏丹）3203亚热笊沟洋性鰹（如】非洲西海岸）370r 4:汐潦气贋（如：就哈拉沙潼）2R05温带大陆性气候3016温带海洋性Hit地上】（Sa：英国）320在Longley-Rice模型中，温带大陆性气候为温地区大片陆地上的典型气候，其典 型特征为显著f气温变化和四季交替。在中纬度沿海地区，强大旷海风为大陆带来了 湿润的空气，因此这里主要是温带海洋性气候 。英国、美国西海岸和欧洲部分地区就 是这种气候的典型代表。对于小于100 km的播路径而言，温带大陆性气候和温带海洋性气候造成的差别微乎其微。但是对于更长的路径而言，带海洋性气候带来了更多的折射，使得在约10%

12、寸间其场强大于温带大陆性气候。第二章 传播模型分析及建模2.1 传播模型的分析与选择飞行器从起飞到飞临目标上空，一般情况下可能会途经平原、丘陵、高山、河流甚 至是海洋等不规则地形，对通信信道损耗的预测需要考虑不同的天然地形环境的影响。 同时还要考虑树木、 建筑物和其他遮挡物等人为因素的影响。 电波传播预测模型大体可 分为两类 : 一类是基于电磁波传播理论 7 ，根据具体的适用环境，确定电磁环境的边界 条件， 求解麦克斯韦电磁波方程式， 进而确定出电磁波的传播路径和传播场强值， 该类 模型通常适用于计算近区场电磁传播， 而对远区场而言边界条件难以确定， 需要考虑的 因素增多，计算相当复杂 ; 另

13、一类是利用数理统计方法，通过将大量数据筛选后进行统 计分析，并结合部分电磁理论来确定对电磁波传播损耗影响较大的因素， 再利用数据拟 合等方法得到电磁波的传播预测模型， 属半经验模型， 对远区场的电磁波预测大都使用 该类模型。通过长期的测试、研究，人们总结归纳了多种适用于远距离的电波传播预测模型，如 Okumur Hata 模型、 Egli 模型和 Longley-Rice 模型等。Okumura 模型以准平坦地形大城市区的中值场强或路径损耗作为参考， 在工程实际 中多用于市区、郊区和开阔地等地形起伏不大的地区 8 。对于起伏较大的不规则地形， 如丘陵地形、水陆混合地形和 孤立山峰， 其传播损耗

14、应在准平坦地形的中值传播损耗的基础上， 加上适当的修正因子 进行校正。 Okumura 模型以曲线图的形式给出，不便于快速的仿真，而 Okumura-Hata 模型是Hata在Okumura曲线图的基础上，通过曲线拟合所作的经验公式：L 69.55 26.16 lg f 13.82 lg ht A(hr) (44.9 6.55 lg ht )(lg d)Loss' (21)式中：f为电波频率，单位MHz;d为通信距离，单位km;ht、hr为收发天线高度，单位 m; Loss'为地形修正因子，A(hJ是移动天线有效高度修正因子，单位dB;为距离修正 因子9。Egli模型是通过在V

15、HF频段和UHF频段对不规则地形上得到的大量实测数据综合分析的基础上提出的一种经验模型，以地形起伏和障碍物高度不超过15m为准，对于地形起伏和障碍物超过15 m的，运用修正因子加以修正。该模型仅适用于视距围。Lon gley-Rice模型被称为不规则地形模型，以电波传播理论为依据，结合丰富的实测数据，用以预测在自由空间中由地形的非规则性造成的中值传播损耗。该模型具有两种预测模式。当能够获取详细的地形剖而数据时10，可以采用点对点模式，如果没有 地形数据，预测模型需要估算与路径相关的参数，需采用区域模式。表2.1不同传播预测模型的适用围适用频 率 / M Hz即离范S /km收发大线 高度/ t

16、n适用地形Edi40- 1 000视即.30/3起伏较低的 不规则地Jg(Ik':murn-Hnia150- 1 9201 10030-200/1-10郊冋、IF禍 地、丘陵等Longley Rice20 - Id 0001-2 0000.5-3 000不规则地形表2.1描述了以上三种模型的具体适用围。 从表中可以看出:Egli模型的适用频率 围较窄，距离围仅为视距，Egli模型不适用于地形高度起伏太大的山区，而 Okumura-Hata模型和Longley-Rice模型均可用于飞行器通信仿真系统。 对这两种预测 模型在开阔地和起伏较大的丘陵地区进行仿真，如图 1和图2所示。10103

17、040506070«090100通信RtA/kmo o o o O8 6 4 2 0i> i- It Al 1图2.1开阔地传输损耗对比图2.1中，Longley-Rice模型的预测值比Okumura-Hata模型的预测值明显偏低。Okumura-Hata模型在两种地形下的预测值变化较小，丘陵地区仅比开阔地偏高20 dB左右，而Longley-Rice模型的变化较大，丘陵地区比开阔地高 30 dB左右。从仿真可 以看:Longley-Rice 模型比Okumura-Hata模型对地形的变化更加敏感，特别是图2.2中Longley-Rice模型的点对点模式能够实时地反映地形对电磁

18、波传播的影响，比区域 模式更加适用于传播地形复杂的飞行器通信信道预测。文献和巨中均使用了 Longley-Rice模型作为海而电波传播模型，但由于应用环境 是海洋，地形不规则度较小，因此，使用的是区域模式;11文献中提出使用Longley-Rice 模型作为地而和导弹通信信道模型，在区域模式下仿真了频率、地形、气候类型和天线 位置对电波传播衰减的影响，但未给出点对点模式下地形影响的仿真结果。本文通过抽取传播路径地形高程值，从以飞行器高度作为接收天线高度的角度，使用点对点模式对 电波传播衰减进行了仿真。2.2 Lon gley-Rice模型建模Lon gley-Rice模型引入了电磁波频率f.、

19、收发天线有效高度ht、hr：及位置、极化方向、地形不规则度厶h、地球表而折射率Ns、地而电导率 和相对介电常数r等因 素，在考虑电波本身特性的基础上，同时兼顾了传播环境的电气特性。 不同路径长度的 传播损耗参考中值Lref的计算如下：max(0, Lbe kid k? Ig d),dmin d dlsLrefLbed md , d Ls d dx(2 2)Lbes msd,dx d式中:d为传播距离，单位为km;f为无线电波频率，单位为MHz; dLs为光滑地面距离；dx 表示此处的绕射损耗和散射损耗相等；Lbe、Lbed、Lbes分别表示自由空间下视距、绕 射和散射时的传播损耗值；&

20、和k2为传播损耗系数；md和叫分别为绕射和散射损耗系 数。dmin d dLS为视距传播距离，dLs d dx为绕射传播距离，dx d为散射传播距离。通过式(2 2)可以分别计算视距传播损耗、衍射传播损耗和散射传播损耗。同时, 再考虑到自由空间传播损耗，整个传播路径上的总体损耗为LbLrefL free(2 3)其中:Lfree 32.45 20lg d 20lg f(24)式中:d为传播距离，单位为km;f为无线电波频率，单位为 MHz衍射传播损耗dLs ddx为衍射传播距离，单位:km.不规则地形中的衍射传播损耗通过结合基于菲涅耳 适用于光滑地而的Vogler修正模型来计算。k f / f

21、0, e a(1 0.04665eNs/N1)dLs. 2hej Te, j 1,2dLs dLs1 dLs2, dL d L1dL2(2hej为收发天线有效高度，单位：m。Xac (k e2) 1/3为地球有效曲率-基尔霍夫理论的双刃峰模型和(2 5)(26)7)mdLdiff (max( dLs,dL 1.38Xae)2.76X ae) Ldiff (maX d Ls,d L 1 .38X ae)2.76Xae(2-8)LedLdiff (max( dLs,dL 1.38Xae)mid(maXdLs，dL 1.38Xae)(29)Ldiff (d)(1 w)Lk wLr Lfo(210)L

22、fomin(15, glg(1 矶血(ls)(211)式(2 10)根据参数二确定双刃峰衰落Lk和圆形地球衰落Lr在衰落计算中的比式(2 11)为杂波干扰衰落。222视距(LOS)传播损耗dmin d dLs为视距12传播距离，单位km.在LOS以反射传播机制为主，采 用双线地而反射模型计算。定义 d2 dLsL2 Led mdd2如果 Led 0，那么 d0 min (1/2dL,4 10 5he1he2 f)，否则do max( 1/2dL,min( Led /md ,dl 2)(211)3 1d1-do -dL(212)4 4Lo LLos(do)，L1 LlosG)(213)式中，Ll

23、°s (1 w)Ld wLt根据参数w确定衍射之外的损耗和双线理论损耗的比重。向散射传播损耗d x d为散射传播距离，单位km计算过程为d5 dL Ds( 214)d6 dL 2Ds( 215)ms (Lscat(d5)Lscat (d6)/Ds( 216)dx maXdLs,dL Xaelg( f), (L5 Led msd5)/(md ms)(217)式中 Lscat(s) 10lg(kH 4) F( s, Ns) H。(218)2.3仿真分析地形剖而数据的获取应用Longley-Rice模型的点对点模式进行计算时，需要获取收发信机之间详细的地形剖而数据。在仿真过程中采用质量较高

24、的航天飞机雷达地形测绘使命高程数据SRTM,分辨率为90m,SRTM数据有多种存储格式，此处使用ASCII格式存储的数据，通过读取ASCII文件头获取存储高程文件的基本信息 冋，如行列数目Nr和Ne、起始经纬 度X。和Yo以及数据元大小Sceii等，然后计算目标位置点相对起始点的偏移量厶 n就能得到该点的高程数据。如果该点不能与文件中的位置相对应， 则使用插值方法，根据若干 相邻点的高程值求出此点的高程值。2.3.2 SRTM高程数据的抽取过程将ASCII文件中除去基本信息的实际高程数据网格化，网格数目为Nr*Ne；网格的起始点坐标0(0,0)，终点坐标E(Nr, Ne).根据收发点的经纬度确

25、定收发点在网格中的坐标(Xt,Yt)(H,ScellScell(Xr,Yr)(严ScelYrScell),向下取整(1) 根据收发点的坐标计算采样点数N max(XT Xr,Yt Yr)和采样点间距d=90 ，单位为m其中YtYrXt Xr 采样点满足直线方程：y x Yr Xr，由此得到采样点的坐标为 (Xi,Y) (x,yj) (Xi, Xi ), i (0,1,2,3.,N)进而求得所有采样点的高程值，也就得到了收发点之间的地形剖而数据。在实验过程中，选取发射点坐标 T (115.9,42.7)，接收点坐标R(117,43.5)，采样点数N =1321,d=111.3 m，收发点间距d=

26、146.9 km.抽取后的地形剖而如图2.3所示。图2.3从SRTM由取的地形剖面图接收机高度描述在Longley-Rice模型的点对点模式14中，需要明确信号发射机的结构高度以及接 收机的结构高度，此处也就是飞行器距离地而的垂直距离。仿真时，可以假设飞行器在 飞行过程中，飞行的海拔高度在一段时间变化不大。但飞行器接收机的结构高度却时刻 在变化着。如图2.4所示，飞行器在位置1和位置2时的飞行高度相同，均为hm,.当在位置1 时，飞行器处于山体的最顶端上空，此时的接收机结构高度为 hg2 hm hj，且接收机 处于发射机的视距围。当在位置2时，飞行器飞临山体的低海拔处，此时接收机的结构 高度变

27、为hg2'，由于山体的阻挡，电波的传输以衍射为主。图2.4飞行器在不同位置时的接收机高度描述234结果在仿真计算中各工作参数设定为：电波工作频率f =1 500 MHz，全向天线，垂直极 化。发射机天线高度hg1 = 100 m飞行器飞行海拔高度hm =1 665 m地而电导率 =0.005,相对介电常数。r=15,地球表而折射率Ns = 320N.发射机坐标为(E115. 9,N42. 7)，飞行器接收机的坐标随飞行器以固定的飞行海拔高度向着坐标为(E117. 0,N43.5)的点沿直线飞行而不断变化，该地区的气候类型为亚热带大陆性气候。仿真结果如图2.5所示。左坐标轴表示地形高程值

28、，单位 m其中最左端标“ ” 处为发射机所在地；右坐标轴表示电波传播衰减，单位a;图上端带箭头的虚线表示飞行器的飞行轨迹。从图2.5中可以看到，在视距最远端点之前的传播衰减值缓慢增长， 且曲线较为光滑，这是因为飞行器飞临视距最远端点之前一直处于发射机的视距围，可以认为是自由空间传播。在视距最远端点处，飞行器和发射机之间恰好有山峰阻挡，故而此处的传播衰减值发生陡变，从 158dB剧增为188dB.此点后的电波传播就以绕射衍 射为主。值得注意的是当飞行器飞临地形中最高的山峰时，即在多重衍射区处，电波传 播衰减不但没有增加，反而有所减少。这种现象就是波在衍射区域遇到阻挡物被反射回 来时会进一步加强原

29、有波的强度。随着飞行器越飞越远进入散射区域后，传播衰减值呈 缓慢增长的趋势，地形的变化己经不是影响电波传播的主要因素，因此，不会引起传播衰减太大的波动。图2.5点对点模式下的传播衰减仿真从图 2.5 的仿真结果可以看出 : 在散射区域之前，电波的传播损耗随地形的高低起 伏发生明显的变化，说明应用 Longley-Rice 模型的点对点模式可以较精确地预测电波 传播衰减。第三章Longley-rice 等效模型的建立3.1 Ion gley-rice的现状和不足Lon gley-Rice 模型是由Lon gley和Rice提出的著名模型，它是一种统计模型，以传播理论为依据，同时结合了数千组实测数

30、据，因此称其为半经验预测模型。由于该模型以传播理论为依据，加上及其丰富的实测数据15，使其得到了广泛的应用。但是Lon gley-Rice模型没有考虑接收机附近的因素以及多径的影响，所以预测精度存在一 定的误差。这里旨在建立一种模型，让接收机附近的因素包含在 Longley-Rice模型中, 使模型更加完善，精度得到响应的提高。3.2 Lon gley-Rice模型的改进由于Longley-Rice模型不能反映接收机附近的路径损耗情况，为了使模型更加完善，提高预测的精度，作者对 Longley-Rice模型做了改进，用等效散射 模型描述了接收机附近的路径损耗情况。由于接收机天线的高度通常很低，

31、电波在传播过程中会遇到各种障碍物、树木以及起伏的地形，引起了电波的反射、折射和绕射，于是到达接收机的电 波是上述电磁波的叠加，如图3-1所示。笈射机图3-1电波的传播模型采用等效的散射体来描述接收机附近的电波传播情况，在该模型中，有n个等效散射体分布在接收机附近以r为半径的圆上，其中有一个散射体在发射机与接收机的视线传播路径上，如图 3.2所示。图3.2等效散射模型第n条路径的电波到达角度为：r 2i sin( i ),i0,1,2,3 .n 1。d n式中，d为发射机和接收机之间的距离；由i , d和r可以确定第n条路径的长度为：dn f ( i ,r,d),i0,1,2.n 1.于是将各个

32、路径的损耗叠加可以求得径的总损耗为：n 1Lp Li, i 0,1,2.n 1.i 1式中，Li为第i条路径的损耗。对Longley-Rice模型做了改进，考虑接收机附近的因素以及多径的影响，建立了 接收机附近的散射模型，使得Longley-Rice模型更加完善，减小了电波传播损耗计算的误差，提高了电波传播预测的精度。第四章 利用 longley-rice 模型设计的可视化电磁环境4.1 电磁环境可视化电磁环境是指存在于给定场所的所有电磁现象的总和， 包括自然环境因素和人为环 境因素，本文主要研究人为电磁设备辐射构成的三维电磁环境 . 由于各种电磁设备数量 众多、使用频段交错、辐射功率大，构成

33、的整个电磁环境异常复杂，相互间的干扰时有 发生，尤其在目前高技术条件下，敌我双方的电磁对抗更加重了电磁环境的复杂性 . 由 于电磁环境不可见， 指挥员不能直观地了解战场中的电磁态势， 所以要快速准确地做出 决策存在巨大的困难 . 计算机图形技术能形象直观地展示数字信息，通过数据可视化， 使用户能直观了解数据隐含的信息，带给用户 16 强烈的视觉感受 .因此，电磁环境可视 化在一定程度上能减少指挥人员对复杂电磁环境掌握的盲目性， 使指挥员对所处环境有 一个直观的认识，为快速准确地决策提供支持 .目前在虚拟战场环境中， 电磁环境的可视化多采用二维的表现方式， 不能体现三维 电磁环境情况 . 本文结

34、合虚拟现实技术，在三维数字地球上展示三维虚拟电磁环境，为 指挥员提供详细的三维电磁态势， 指挥员还能动态调整电磁设备参数， 并直观查看电磁 态势相应的变化， 通过交互更加直接了解各种条件下的复杂电磁环境， 为其快速准确地 决策提供直观的辅助支持 . 然而，把复杂的三维电磁环境高效且准确地展示到数字地球 上是非常困难的 . 为了可视化电磁环境，需采用合适的电波传播模型计算电磁设备传播 数据.虽然时域有限差分 (FDTD) 方法能精确描述电波传播， 但由于其需耗费大量的计算 时间，难以满足实时动态可视化三维电磁环境的要求 .Longley-rice 电波传播模型是一 种基于统计分析和电磁理论的电波

35、传播计算模型 17 ，它考虑了地形和大气的影响 . 通过 数字地形模型，可采用 Longley-rice 。点到点的预测模式计算电波传播，它能精细考 虑不规则地形对电波传播的影响，且该模式相比 FDTD能较快计算电波传播数据，适合 在数字地球上动态展示三维电磁环境 . 直接体绘制方法是一种可视化分析体数据的有效 工具，由于图形硬件能力的提高和功能的增强，采用图形硬件图形处理器(GPU)加速的直接体绘制方法可在主流图形硬件条件下实时绘制，并能得到高质量的绘制结果. 目前对电磁环境中雷达探测围的可视化己有一些研究， 主要采用了单层等值面的方法展示雷 达三维探测围，能直观展示雷达最大探测围的包络 ;

36、 文献采用等值面提取的方法对高压 电线周围的电场进行可视化， 18 可清晰展示特定电位大小的电场覆盖围， 但由于等值面 方法对表现体数据细节的不足，尚不能展示电磁环境部的细节信息 . 多层半透明等值面 绘制方法通过展示体数据的多层信息， 能在一定条件下弥补单层等值面的不足， 但没有 对体数据进行信息分析的切片、切割能力 . 直接体绘制技术能详细展示体数据信息，还 能从多个角度切割分析数据， 如文献采用直接体绘制方法绘制了电磁态势， 但没有采用 GPU加速的方法，不能实时展示电磁态势体数据.因此本文扩展了文献中的硬件加速直 接体绘制方法，通过坐标转换把三维电磁环境直接实时绘制到数字地球上。4.2

37、 三维电磁环境体数据生成为了对虚拟电磁环境进行展示， 首先计算虚拟环境中设置的各个电磁设备向三 维空间中辐射电磁波的情况 . 本文采用 Longleyice 电波传播模型，即 不规则地形模型来预测不同频段的电波传播，它能较准确地估计点波传播损失. 通过Longleyice 电波传播模型得到各个 19 电磁设备在三维空间中的辐射损失值， 进而计算 得到电磁波三维空间功率密度数据，以此作为三维电磁环境表现的基础 .4.2.1 不规则地形剖面提取不规则地形模型需要地形剖面数据来精确计算地形对电波传播损失的影响， 数字地 形高程数据用来生成地形剖面数据 . 如图 4.1 (a) 所示，电磁设备向 四周

38、发射的线条表示电波传播路径，背景线条是数字高程模型地形网格 . 为得到电波传 播路径上地形高程数据， 采用双线性插值方法， 按照数字地形高程数据的分辨率大小等 间隔在数字地形上采样生成路径剖面高程数据20，图4.1 (b)示出了采样得到的地形剖 面示意图。何采样地形出)地形剖面数据图4.1地形剖面选取423电磁环境体数据计算由式P Ft L Ls可计算得到三维电磁环境功率密度值，为减少存储空间和提高绘制效率，把电磁数据离散为一定间隔的网格数据，如图4.2 (a)所示.首先，按照高度方向把电磁空间分成n层，再把每层按照经纬度方向离散为m*k个网格点，用式F Ft L Ls计算每个离散网格点上的功

39、率密度值采用上述离散方法计算生成的电磁 环境数据大小为m*k*n,图4.2 (b)示出了离散电磁数据在三维环境中的强度情况，分 别用不同灰度表示强度大小.电磁数括"二细环览中的眾应图4.2电磁环境体数据4.3在三维数字地球上体绘制电磁环境431硬件加速的直接体绘制基于GPU勺光线投射体绘制把三维体数据沿光线方向投射到二维图像上， 采用体数 据包围盒作为绘制代理网格，并把包围盒的每个顶点颜色设置为三维纹理坐标， 采用多 遍绘制技术绘制体数据信息.传统基于硬件加速的直接体绘制方法采用规则包围盒绘制体数据，且体数据是规则网格.但是生成的电磁环境体数据并不是规则体数据，如图4.2 (b)所示

40、，电磁环境体数据采用球坐标系进行组织，因此传统基于GPU勺光线投射直接体绘制方法不能直接应 用到电磁环境体数据的绘制.下面将详细介绍扩展传统基于GPU的光线投射直接体绘制 方法，巧妙通过坐标转换实现在三维数字地球上展示电磁环境体数据.电磁环境体数据包围盒图4.3 (a)示出了在数字地球上表示的电磁环境体数据不规则包围盒的线框图， 包围盒底部的4个顶点与其上部对应的4个顶点经纬度值相同21，而海拔不同，但包围 盒的经纬度与三维电磁环境体数据计算区域的经纬度相同.假设三维电磁环境体数据计 算区域是从位置V°( 0，0, 0)至位置Vi( i, i, i)位置点均采用球坐标，其中n, n,

41、 n,(n0,1)分别为纬度、经度和海拔.图4.3(b)用填充色标示了三维电磁环境体数据的剖面，不规则包围盒垂直截面ABCD用粗实线表示.根据几何关系可明显得到包围 盒的底部B点和C点海拔为°，顶部A点和D点的海拔为H R1 R(4-1 )cos其中，R为地球半径；B为max 1°, 1。)/2.从图3 (b)可见，电磁环境体数据包围盒大于体数据的计算区域围，这样可得到体数据包围盒海拔高度：底部为°，顶部为H.包帀盒(坊包用盒裁面图4.3电磁体数据包围盒体数据纹理坐标转换在传统基于GPI光线投射体绘制方法中，包围盒顶点颜色值设置为三维体数据纹理 的纹理坐标，通过绘

42、制包围盒背面到纹理作为投射光线终止点.计算得到的三维电磁环境体数据是不规则的，难于计算其包围盒上的体数据纹理坐标， 且由于包围盒中其纹理 坐标系统不是均匀变化的，不能直接把三维体数据纹理坐标设置成顶点坐标， 通过线性 插值得到其他纹理坐标.本文提出一种方法，不用三维体数据纹理坐标作为顶点颜色值， 而是采用顶点的球坐标位置，通过坐标转换把球坐标转换为纹理坐标.首先，把球坐标转换为笛卡儿坐标x r cos cosy r cos sin( 4-2)z r si n其中，r, a，B为球坐标；x,y,z为笛卡儿坐标.用式(4-2)把包围盒的8个顶点坐标转 换为笛卡儿坐标，然后设置其颜色值.通过GPU孚

43、点纹理的支持，包围盒浮点颜色值被 绘制到了浮点纹理上，片段着色程序查询浮点纹理就能得到包围盒的顶点位置值，即得到投射光线的终止点.绘制包围盒正面作为投射光线的起始点22，因此投射光线方向就是终止点与起始点 的矢量差，沿着光线方向就能采样得到该光线穿过的体数据值.但该光线上的位置是在 笛卡儿坐标系下的坐标，不能直接从三维体数据纹理中得到体数据，需把笛卡儿坐标转 化为纹理坐标.由于体数据计算围为V°( °, °, °)至V,( 1, 1, J，首先笛卡儿坐标(x,y,z)转换为球坐标(r, a , B )，arcsi n?(4-3)ry arcta n 丄x

44、通过球坐标插值可得到纹理坐标为v 0(4-4)i or osi o其中体数据纹理坐标u,v,s的有效围为0,1，即在体数据计算围.因此如果计算得 到的纹理坐标不在该围，可采用空白区域剔除方法，不绘制该区域处的体数据，可加速 体数据绘制.采用OpenGL图形接口和Cg语言实现了硬件加速的电磁环境直接体绘制方法 实验 中在三维数字地球上设置了 20部电磁设备，假设为全向天线，地形数据采用 90 m分辨 率.电磁设备参数见表1.Longley-ric 。电波传播参数设置如下：折射率为320,介电常 数和地面电导率分别为15 F/m和0.005 S / m.地形剖面采样间隔为100 m,实验用个 人计

45、算机(PC)配置为 Intel Core 2 Duo 2.8 GHz ，显卡为 GeForee 8600 GT.图4.4示出了电磁环境在三维数字地球上的情况，电磁环境体数据分辨率为100 X 100 X 30,经度和纬度等分间隔都为0.O5度，高度等分间隔为200 m覆盖围约为500 km X 500 km,每个采样点数据采采用16位浮点数存储.如果体数据围较大，可把体数据划分为多个分块区域，每块采用同样的方法绘制，然后按照分块从后往前的顺序融合到 三维环境中，因此本文方法适合较大数据量的电磁环境体数据绘制釆用Longley-rice 电波传播模型计算实验用电磁环境体数据耗时8.836 s，窗

46、口大小为1024 X 768时绘制平均帧率为12帧/s.相比文献，本文方法可实时直接体绘制 电磁环境，且各个电磁设备参数可交互改变，系统将重新计算生成新的电磁环境体数据. 电磁环境体数据通过颜色映射工具(见图4.4(c)为不同的功率密度设置不同的颜色和 透明度，三维电磁环境体绘制效果如图4.4(a)所示.电磁环境在三维环境中的情况可直观地动表4.1电磁设备参数編号频率/MHz编号频率/MHz功率八W145tl2 5001121HM)220-J2125 IM)1PD】3仙创)0斗350I 5001430)2 500他2 non154X0SOOh2001 5(10162<X)I 000730

47、01 200】丁3S04S0S4001 2UU183502 500g3002(X)0194001 SOO101025002050)1 200态展示，而且可交互改变电磁设备参数，能极大辅助电磁设备设计和规划.图4.4(b)示出了关闭图4.4(a)中一部电磁设备后的情况，可清楚看到该电磁 设备影响情况.电磁环境体绘制(11动态关闭屯磁设备M颜色映射工貝-图4.4在三维地球上绘制的结果图4.5通过体数据切割技术示出了电磁环境在不同方向上的辐射情况，相比文献的等值面方法，本文方法能展示电磁体数据部细节.其中图4.5(a)是500 m海拔上电磁 环境切片，图4.5 (b)是5 km海拔上的电磁环境切片，

48、图4.5 (c)是沿经度、纬度方向 切割电磁环境的情况.( a) 5CX)m 髙度(b) 5 kni 髙度切块操作图4.5体数据切割这是一种传统体绘制方法的扩展方法，巧妙通过坐标转换把电磁环境体数据及时展 示到了三维数字地球上.采用Lon gley-rice 电波传播模型凶计算电磁环境数据，该模型 考虑了不规则地形和大气的影响，能较准确地预测电波传播.最后在数字地球上实现了三维电磁环境绘制.结果表明，本文方法能及时展示虚拟电磁环境，而且可动态调整电 磁设备参数(如频率、功率等)，及时展示不同电磁参数下电磁环境， 可用于描述大围三 维电磁态势.虽然采用了 Longley-rice 传播模型来生成

49、电磁环境体数据，但本文方法也 适合其他方法生成的体数据的及时可视化，因此下一步计划采用更为准确的电波传播模 型，并采用并行计算的方式提高电磁数据计算效率，增强电磁环境体数据准确性.第五章用 Longley-Rice 大气折射修法的改进5.1 对 Okumura-Hata 和 Longley-Rice 的比较 信道模型分为大尺度模型和小尺度模型两种， 大尺度型刻画信号传输过程中经历的 路径损耗中值，小尺度模型刻画小围 （几个波长 ）信号的幅度和相位上的快速衰落。 在大 尺度传播衰落模型研究方面先后建立了以下著名的信道模型 : Longley-Rice 模型、在 几公里到几百公里围 Okumura

50、-Hata 被广泛用来预测接收信号的场强，但海上移动信道 路径损耗预测利用 Longley-Rice 模型更合适， Longley-Rice 模型考虑了更多的与地形 有关的因素，包括面折射率、地面导电率、介电常数以及地面粗糙度等，还考虑了不同 的气候类型和天线的位置标准等。文献通过对 Okumura-Hata 和 Longley-Rice 模仿真结果与实测数据的比较，指出 预测海上传播损耗时后更确切。 Longley-Rice 模型适合于传输距离为 1-2000 公里， 频率围为2-2000MHZ,天线高度0-10000米，支持多种气候类型、多种地形、多种反 射面传播特性、三种天线位置标准、天

51、线激化方式等参数。该模型以传播理论为依据， 同时结合了数千组实测数据，因此称其为半经验预测模型。 Longley-Rice 模型从 1967 年提出以来在当时航海等大区制通信场强预测模型中得到了广泛的应用， 并不断被深人 研究和改进， Hufford 等人根据模型提供的理论和原始数据给出了模型的实用算法。由于大气折射的影响， 电波在大气中传播的实际路径是一条曲线， 与收发信机之间 的实际距离存在着误差，预测信道传播路径损耗时需要对实际距离加以修正。5.2 Longley-Rice 模型中对大气折射误差的修正Longley-Rice 模型中对大气折射误差的修正采用等效地球半径法， 等效地球半径

52、法计算简单， 但是精度不高， 该方法只利用了大气地面折射率的水平差异， 未利用垂直 剖面信息， 而实际大气折射率垂直差异， 尤其在近地面 lkm 的垂直差异比较大， 本文对Longley-Rice 模型大气折射误差修法加以改进。改进方法充分利用大气折射率垂直剖 面信息，利用电波实际传播路径和收发天线地面水平距离之间的映射关系， 将模型输人 的收发天线之间的实际距离映射为电波实际传播距离， 能更精确反映传播距离对电波传 播损耗的影响。改进算法预测结果比原有方法更接近实测结果，具有较大的实用价值。5.3 大气折射修法的改进Longley-Rice 模型参考衰减为距离的分段函数 max(0,Lbe

53、k1d k2 lg d ), dmin d dLsLrefLbed md,dLs d dx( 5-1)Lbes msd,dx d定义的三段分别称为视距区域， 绕射区域和散射区域。 模型控制置信度的参数由统 计量给出， L(qT,qL,qS )作为时间，位置，形势分量的函数，衰减不会超出这个值。也 就是说，在必的情况下，至少有私的地点衰减不会超过 L(qT,qL,qS )，并且这个时间至 少为 qT。对于光滑地球地平线dLsb的距离，是基于平坦地面的二径反射理论和一个绕射损耗 的外插值计算的；对恰好超过从dLs到dx,的地平线距离(为绕射损耗和散射损耗相等的 地点) ，参考衰减值是刃峰绕射和光滑

54、地球绕射计算的加权平均。加权因子为频率、地 形不规则因子和天线高度的函数。 对于很不规则地形， 从终端看地平线障碍物可看做是 锐利山脊，绕射损耗是基于 Epstein-Peterson 近似按双刃峰路径计算的。对围大于的 超视距路径，参考衰减按按绕射损耗计算和正向散射损耗计算中取较小者。由于大气层密度分布的不均匀性， 导致电波传播产生折射现象， 使得电波传播距离 与实际距离不符。在 Longley-Rice 模型中，传输距离参数给出的是发射机和接收机之 间的地面水平距离， 但影响传输损耗的距离是电波实际传播的距离 (视在距离) ，并且该 模型预测损耗值是传播距离的分段函数， 传播距离直接影响传

55、输衰落曲线的转折点， 因 此必须对给出的传播距离参数加以修正。 Longley-Rice 模型采用等效地球半径法修正 大气折射误差，等效地球半径定义为散射梯度N或地球表面折射值 N的函数。定义N=(n-1) x 10A6其中，。为地球表面大气折射指数。该方法对折射率剖面的计算没有利 用折射率负梯度的值在不同气候类型之间有较大差异的特性，实际中这种差异比较明显，而折射率剖面直接影响到电波传播的实际路径，从而影响预测精度。下面提出一种改进的折射修法，以改进Longley-Rice电波传播衰落预测精度。电波折射示意图如图1所示，T为发射天线相对平均海平面的高度，T'为T在海 平面上的投影，R为接收天线相对海平面的高度，R'为R在海平面上的投影，0为地心 原点，OT'.OR'为地球半径，P为射线上的某点，B为射线在 P点处的仰角，dl为射线 在P点处的一小段弧线，dr为地球径向延长线在P点处的线段小量，d为射线上P点 对应的地心角小量，R'T'为收发天线之间的地面水平距离，ds为R'T'在P点与地心连 线和地球表面的交点处水平距离小量， TT'为发射天线高度，单位 m RR'为接