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衍射传播损耗92.2.2视距(LOS)传播损耗102.2.3向散射传播损耗102.3仿真分析102.3.2 SRTM高程数据的抽取过程112.3.3接收机高度描述122.3.4 结果13第三章Longley-rice等效模型的建立153.1 longley-rice的现状和不足153.2 Longley-Rice模型的改进15第四章利用lon

2、gley-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-Rice模型中对大气折射误差的修正255.3大气折射修正方法的改进26第六章展望30参考文献31第一章绪论1.1研究背景Longley-Ric

3、e无线电波传输模型1是由Longley和Riee提出的无线电波传输模型，该模型为统计模型，它是以无线电波传播理论为依据，并结合了许多实际环境下的测量数据，所以该模型被称为半经验预测模型。正是因为该模型是以无线电波传播理论为依据,有加上了大量的实际测量数据，因此该模型得到了广泛的应用。该模型被称为不规则地面模型，可以用该模型预测自由空间中由地形非规则性所造成的中值传输衰落。1.2模型分类及参数Longley-rice模型有:2种模式。当地形路径易据很详细时，特定路径参数就很容易被确定。这不预测方式为“点到点预测”。如果地形数据不够训确，可以利用Longley-Rice模型估计特定参数的值这种预测

4、方式为“区域预测”。Longley-rice模型有4种变化模式，分别为单信号模式、单天线模式、移动模式和广播模式。在longley-rice模型的早期程序中，所有点对点预测的计算都是基于单天线模式，这里天线的位置是不变的。后来，由于对计算精度需求的提高，人们才引入其他模式。在各种变化模式中，变化的主要是时间、位置和情景3个参量，或者说是一个三维变量。目前，Longley-Rice无线电波模型已有相关的计算机仿真程序，可以用来对无线电波传输的损耗进行计算。当无线电波传输路径已知时，计算机的仿真程序可以通过无线电波传播路径的长度、极化方向、无线电波频率、地面有效半径、收发天线高度、地面导电常数以及

5、表面绕射率等参数确定无线电波传输损耗的大小。Longley-Rice预测模型主要有以下参数： 天线极化方式:可以采用水平极化或者垂直极化。Longley-Rice模型中假设发射天线和接受天线具有相同的极化方式； 折射率:空气的折射率决定了无线电波的“弯曲”程度。在一般的模型中，空气折射率用地面有效曲率来代替，通常取1.333。 介电常数:地面的相对介电常数和电导率1.3Longley-Rice模型传输损耗根据无线电波的传播范围的不同，Longley-Rice模型的传输损耗可大致分为三种情况，它们分别为：(1) 视距传播损耗(2) 绕射传播损耗(3) 散射传播损耗。当无线电波传播于视距范围内时，

6、无线电波传播方式主要以反射传播为主。通过对地貌地形的路径及对流层的绕射特点进行分析，预测出无线电波传输中值传输衰落，将双线地面反射模型用来模拟地平线以内的传输场强。视距传播模型的适用范围为:d<d<d。minLS(2) 绕射传播损耗当无线电波传播于超视距范围内时，无线电波传播方式主要以绕射传播为主。但是,当无线电波的传播环境为不规则的地形时，会有两种理论用于计算绕射损耗。其中一种适用于预测非球形但光滑的地面上无线电波的传播,而另外一种则适用于非常不规则的地面上无线电波的传播。采用Fresnel-Kirchoff刃形模型来预测无线电波传播的绕射损耗，其计算结果是上述两种理论结果的加权

7、和。适用范围:d<d<d为绕射传播距离，单位为km。Lsx(3) 散射传播损耗当无线电波的传输距离为大大超出地平线的远距离无线电波的传输时，无线电波传播方式主要以散射传播为主，预测理论以绕射理论为主。适用范围:d<d为散射传播距离，单位为km。x综上所述，Longley-Rice模型传输损耗L为：bL二L+L(3-15)breffree其中：L二32.45+20lgd+20lgf(3-16)freemax(0,L+kd+klgd),d<d<dbe12minLsLL+m,d<d<d(3-17)refbeddLsxL+md,d<dbessx式中:d为

8、传播距离，单位为km;f为无线电波频率,单位为MHz;d为光滑地面距离;dLsx表示此处的绕射损耗和散射损耗相等;L、L、L分别表示自由空间下视距、绕bebedbes射和散射时的传播损耗值;k和k为传播损耗系数;m和m分别为绕射和散射损耗系12ds数。d<d<d为视距传播距离，d<d<d为绕射传播距离，d<d为散射传播距minLSLsxX离。1.4Longley-Rice模型的适用范围Longley-Rice模型中的实测数据大多数取自10-1000MHz的频率范围，其中20-100MHz的数据涉及5-50km的距离和1-9m的收、发信天线高度;较高频段的数据涉及5

9、-1000km的距离，10-1500m的发射天线高度和3-9m的接收天线高度。数据来源于世界各地，但主要是美国，多数为移动记录结果。©Longley-Rice模型给出了参考衰减值的计算公式及不同环境下相关修正因子的详细说明，公式中所使用的参数包括：不规则地形参数、频率、收发信机天线高度和表面折射率等3。同时还引入了反映介质特性的2个参数:介电常数和导电率。以传播理论为依据，加上极其丰富的实测数据，使得Longley-Rice模型使用范围得到了拓展，其适用范围如下：1) 频率f:20-40000MHz;2) 收、发信机天线高度:0.5-3000m;3) 覆盖半径:1-2000km;4）

10、表面折射率:250-400Ns。表1.1地面的相对介电常数和导电率相对介电常数导电率/（S/m）一般的地面150.005导电性差的地面40.001导电性好的地面250.020纯水810.010海水815.000表1.2气候类型和相应地表折射率气候类型地表折射率（N-units）1近赤道气候（如：刚果）3602亚热带大陆性气候（如：苏丹）3203亚热带海洋性气候（如：非洲西海岸）3704沙漠气候（如：撒哈拉沙漠）2805温带大陆性气候3016温带海洋性气候（陆地上）（如：英国）320在Longley-Rice模型中，温带大陆性气候为温地区大片陆地上的典型气候，其典型特征为显著f气温变化和四季交替

11、。在中纬度沿海地区，强大旷海风为大陆带来了湿润的空气，因此这里主要是温带海洋性气候。英国、美国西海岸和欧洲部分地区就是这种气候的典型代表。对于小于100km的播路径而言，温带大陆性气候和温带海洋性气候造成的差别微乎其微。但是对于更长的路径而言，带海洋性气候带来了更多的折射，使得在约10%时间内其场强大于温带大陆性气候。第二章传播模型分析及建模2.1传播模型的分析与选择飞行器从起飞到飞临目标上空，一般情况下可能会途经平原、丘陵、高山、河流甚至是海洋等不规则地形，对通信信道损耗的预测需要考虑不同的天然地形环境的影响。同时还要考虑树木、建筑物和其他遮挡物等人为因素的影响。电波传播预测模型大体可分为两

12、类:一类是基于电磁波传播理论7，根据具体的适用环境，确定电磁环境的边界条件，求解麦克斯韦电磁波方程式，进而确定出电磁波的传播路径和传播场强值，该类模型通常适用于计算近区场电磁传播，而对远区场而言边界条件难以确定，需要考虑的因素增多，计算相当复杂;另一类是利用数理统计方法，通过将大量数据筛选后进行统计分析，并结合部分电磁理论来确定对电磁波传播损耗影响较大的因素，再利用数据拟合等方法得到电磁波的传播预测模型，属半经验模型，对远区场的电磁波预测大都使用该类模型。通过长期的测试、研究，人们总结归纳了多种适用于远距离的电波传播预测模型，如OkumurHata模型、Egli模型和Longley-Rice模

13、型等。Okumura模型以准平坦地形大城市区的中值场强或路径损耗作为参考，在工程实际中多用于市区、郊区和开阔地等地形起伏不大的地区。对于起伏较大的不规则地形，如丘陵地形、水陆混合地形和孤立山峰，其传播损耗应在准平坦地形的中值传播损耗的基础上，加上适当的修正因子进行校正。Okumura模型以曲线图的形式给出，不便于快速的仿真，而Okumura-Hata模型是Hata在Okumura曲线图的基础上，通过曲线拟合所作的经验公式：L二69.55+26.16lgf-13.82lgh-A(h)+(44.9-6.55lgh)(lgd)y+Loss'(21)trt式中:f为电波频率，单位MHz;d为通

14、信距离，单位km;h、h为收发天线高度，单位trm;Loss'为地形修正因子，A(h)是移动天线有效高度修正因子，单位dB;Y为距离修正r因子9。Egli模型是通过在VHF频段和UHF频段对不规则地形上得到的大量实测数据综合分析的基础上提出的一种经验模型，以地形起伏和障碍物高度不超过15m为准，对于地形起伏和障碍物超过15m的，运用修正因子加以修正。该模型仅适用于视距范围内。Longley-Rice模型被称为不规则地形模型，以电波传播理论为依据，结合丰富的实测数据，用以预测在自由空间中由地形的非规则性造成的中值传播损耗。该模型具有两种预测模式。当能够获取详细的地形剖而数据时10，可以采

15、用点对点模式，如果没有地形数据，预测模型需要估算与路径相关的参数，需采用区域模式。表2.1不同传播预测模型的适用范围模型名称适用频率/MHz距离范韦ktn收发天线高度加适用地形Egli40-1000视距30/3起伏较低的不规则地形()kumura-Hata15019201-10030-200,1-10郊区、开阔地、丘陵等Loiighy-Ricr20-40000170000.5-3000不规则地形表2.1描述了以上三种模型的具体适用范围。从表中可以看出:Egli模型的适用频率范围较窄，距离范围仅为视距，Egli模型不适用于地形高度起伏太大的山区，而Okumura-Hata模型和Longley-R

16、ice模型均可用于飞行器通信仿真系统。对这两种预测模型在开阔地和起伏较大的丘陵地区进行仿真，如图1和图2所示。Okumura-HadLongIey-RicM4K式FreeSpace100/B°i)IQ3Q<3040SO6070W90100通信Ki.Sl/km图2.1开阔地传输损耗对比图2.1中，Longley-Rice模型的预测值比Okumura-Hata模型的预测值明显偏低。Okumura-Hata模型在两种地形下的预测值变化较小，丘陵地区仅比开阔地偏高20dB左右，而Longley-Rice模型的变化较大，丘陵地区比开阔地高30dB左右。从仿真可以看:Longley-Ric

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

18、点模式对电波传播衰减进行了仿真。6021020旳40506G70fiO901Q0通信距离/km图2.2丘陵地区传播损耗对比2.2 Longley-Rice模型建模Longley-Rice模型引入了电磁波频率f.、收发天线有效高度h、h:及位置、极化方tr向、地形不规则度Ah、地球表而折射率N、地而电导率b和相对介电常数e等因素,sr在考虑电波本身特性的基础上，同时兼顾了传播环境的电气特性。不同路径长度的传播损耗参考中值L的计算如下:refmax(0,L+kd+klgd),d<d<d(22)be12minLsLL+m,d<d<drefbeddLsxL+md,d<db

19、essx式中:d为传播距离，单位为km;f为无线电波频率,单位为MHz;d为光滑地面距离;dLsx表示此处的绕射损耗和散射损耗相等;L、L、L分别表示自由空间下视距、绕bebedbes射和散射时的传播损耗值；k和k为传播损耗系数;m和m分别为绕射和散射损耗系12ds数。d<d<d为视距传播距离，d<d<d为绕射传播距离，d<d为散射传播距minLSLsxX离。通过式(22)可以分别计算视距传播损耗、衍射传播损耗和散射传播损耗。同时，再考虑到自由空间传播损耗，整个传播路径上的总体损耗为L二L+L(23)breffree其中：L二32.45+20lgd+20lgf(2

20、4)free式中:d为传播距离，单位为km;f为无线电波频率，单位为MHz。2.2.1衍射传播损耗d<d<d为衍射传播距离，单位:km.Lsx不规则地形中的衍射传播损耗通过结合基于菲涅耳-基尔霍夫理论的双刃峰模型和适用于光滑地而的Vogler修正模型来计算。k=f/f丫=丫(10.04665eNs/nJ0,ea(25)d=、；2h/y,j=1,2Lseje(26)d=d+d,d=d+dLsLs1Ls2LL1L2(27)h.为收发天线有效咼度，单位:m。eJX=(ky2)-1/3为地球有效曲率acem=L(max(d,d+1.38X)+2.76X)一LddiffLsLaeaediff

21、(max(d,d+1.38X)1(2-8)LsLae2.76XaeL二L(max(d,d+1.38X)m(max(d,d+1.38X)(2一9)eddiffLsLaedLsLaeL(d)二(1w)L+wL+L(210)diffkrfoL=min(15,glg(1+akhha(d)(2一11)fog1g2hLs式(210)根据参数二确定双刃峰衰落L和圆形地球衰落L在衰落计算中的比kr重。式(211)为杂波干扰衰落。2.2.2视距(LOS)传播损耗d<d<d为视距12传播距离，单位km.在LOS内，以反射传播机制为主,minLS采用双线地而反射模型计算。定乂d=dL=L+md2Ls2e

22、dd2如果L>0，那么d=min(1/2d,4x10-5hhf)，否则ed0Lele2d=max(1/2d,min(-L/m,d2)(2一11)0LeddL31d=d+d(212)14o4lL=L(d)，L=L(d)(213)0Los01Los1式中，L=(1w)L+wL根据参数w确定衍射之外的损耗和双线理论损耗的比重。Losdt2.2.3向散射传播损耗d<d为散射传播距离，单位km。计算过程为xd=d+D5Ls(214)d=d+2D(215)6Lsm=(L(d)L(d)/D(216)sscat5scat6sd=maxd,d+Xlg(f),(LLmd)/(mm)xLsLae5eds

23、5ds(217)式中L(s)=10lg(kHO4)+F(Os,N)+H(218)scats02.3仿真分析2.3.1地形剖而数据的获取应用Longley-Rice模型的点对点模式进行计算时，需要获取收发信机之间详细的地形剖而数据。在仿真过程中采用质量较高的航天飞机雷达地形测绘使命高程数据SRTM，分辨率为90m,SRTM数据有多种存储格式，此处使用ASCII格式存储的数据,通过读取ASCII文件头获取存储高程文件的基本信息13,如行列数目N和N、起始经re纬度X和Y以及数据元大小S等，然后计算目标位置点相对起始点的偏移量厶n就能00cell得到该点的高程数据。如果该点不能与文件中的位置相对应，

24、则使用内插值方法，根据若干相邻点的高程值求出此点的高程值。2.3.2 SRTM高程数据的抽取过程将ASCII文件中除去基本信息的实际高程数据网格化，网格数目为N*N网格的re起始点坐标O(0,0)，终点坐标E(N,N).根据收发点的经纬度确定收发点在网格中的坐re标xT7yTSLcellS-cellxy44SSLcell-Lcell),),向下取整。(Xt，勺=(Xr，Yr)二(根据收发点的坐标计算采样点数N=max(|X-X|,Y-Y|)和采样点间距TRTRd=90九，单位为m，其中X=广r.X-XTR(2)米样点满足直线方程：y=Xx+Y-XX，由此得到米样点的坐标为RR(X,Y)=(x,

25、y)=(xXx),ie(0,1,2,3.,N)进而求得所有采样点的高程值，也就得到了iiiii,i收发点之间的地形剖而数据。在实验过程中，选取发射点坐标T(115.9,42.7),接收点坐标R(117,43.5),采样点数N=1321,d=111.3m,收发点间距d=146.9km.抽取后的地形剖而如图2.3所示。中抽取地形剖面数据-<左鎖笔<全猴图2.3从SRTM抽取的地形剖面图2.3.3接收机高度描述在Longley-Rice模型的点对点模式】4中，需要明确信号发射机的结构高度以及接收机的结构高度，此处也就是飞行器距离地而的垂直距离。仿真时，可以假设飞行器在飞行过程中，飞行的海

26、拔高度在一段时间内变化不大。但飞行器接收机的结构高度却时刻在变化着。如图2.4所示，飞行器在位置1和位置2时的飞行高度相同，均为h，当在位置1m时，飞行器处于山体的最顶端上空，此时的接收机结构高度为h二h-h，且接收机g2m1处于发射机的视距范围内。当在位置2时，飞行器飞临山体的低海拔处，此时接收机的结构高度变为h'，由于山体的阻挡，电波的传输以衍射为主。g2接收机位戏2图2.4飞行器在不同位置时的接收机高度描述2.3.4结果在仿真计算中各工作参数设定为:电波工作频率f=1500MHz，全向天线，垂直极化。发射机天线高度h=100m，飞行器飞行海拔高度h=1665m，地而电导率q=0.

27、g1m005,相对介电常数。e=15,地球表而折射率N=320N.发射机坐标为(E115.9,N42.7)，rs飞行器接收机的坐标随飞行器以固定的飞行海拔高度向着坐标为(E117.0,N43.5)的点沿直线飞行而不断变化，该地区的气候类型为亚热带大陆性气候。仿真结果如图2.5所示。左坐标轴表示地形高程值，单位m，其中最左端标“.处为发射机所在地;右坐标轴表示电波传播衰减，单位a;图上端带箭头的虚线表示飞行器的飞行轨迹。从图2.5中可以看到，在视距最远端点之前的传播衰减值缓慢增长，且曲线较为光滑，这是因为飞行器飞临视距最远端点之前一直处于发射机的视距范围内，可以认为是自由空间传播。在视距最远端点

28、处，飞行器和发射机之间恰好有山峰阻挡，故而此处的传播衰减值发生陡变，从158dB剧增为188dB.此点后的电波传播就以绕射衍射为主。值得注意的是当飞行器飞临地形中最高的山峰时，即在多重衍射区处，电波传播衰减不但没有增加，反而有所减少。这种现象就是波在衍射区域遇到阻挡物被反射回来时会进一步加强原有波的强度。随着飞行器越飞越远进入散射区域后，传播衰减值呈缓慢增长的趋势，地形的变化己经不是影响电波传播的主要因素，因此，不会引起传播衰减太大的波动。70050LOO图2.5点对点模式下的传播衰减仿真从图2.5的仿真结果可以看出:在散射区域之前，电波的传播损耗随地形的高低起伏发生明显的变化，说明应用Lon

29、gley-Rice模型的点对点模式可以较精确地预测电波传播衰减。第三章Longley-rice等效模型的建立3.1longley-rice的现状和不足Longley-Rice模型是由Longley和Rice提出的著名模型，它是一种统计模型，以传播理论为依据，同时结合了数千组实测数据，因此称其为半经验预测模型。由于该模型以传播理论为依据，加上及其丰富的实测数据15，使其得到了广泛的应用。但是Longley-Rice模型没有考虑接收机附近的因素以及多径的影响，所以预测精度存在一定的误差。这里旨在建立一种模型，让接收机附近的因素包含在Longley-Rice模型中，使模型更加完善，精度得到响应的提高

30、。3.2 Longley-Rice模型的改进由于Longley-Rice模型不能反映接收机附近的路径损耗情况，为了使模型更加完善，提高预测的精度，作者对Longley-Rice模型做了改进，用等效散射模型描述了接收机附近的路径损耗情况。由于接收机天线的高度通常很低，电波在传播过程中会遇到各种障碍物、树木以及起伏的地形，引起了电波的反射、折射和绕射，于是到达接收机的电波是上述电磁波的叠加，如图3-1所示。发射机图3-1电波的传播模型采用等效的散射体来描述接收机附近的电波传播情况，在该模型中，有n个等效散射体分布在接收机附近以r为半径的圆上，其中有一个散射体在发射机与接收机的视线传播路径上，如图3

31、.2所示。图3.2等效散射模型第n条路径的电波到达角度为：r2兀&.=sin(i),i=0,l,2,3.n-1。1dn式中，d为发射机和接收机之间的距离;由9，d和r可以确定第n条路径的长度为:id=f(9,r,d),i二0,1,2.n1.ni于是将各个路径的损耗叠加可以求得径的总损耗为：L=£L,i=0,1,2n1.pii=1式中，L为第i条路径的损耗。i对Longley-Rice模型做了改进，考虑接收机附近的因素以及多径的影响，建立了接收机附近的散射模型，使得Longley-Rice模型更加完善，减小了电波传播损耗计算的误差，提高了电波传播预测的精度。第四章利用longl

32、ey-rice模型设计的可视化电磁环境4.1电磁环境可视化电磁环境是指存在于给定场所的所有电磁现象的总和，包括自然环境因素和人为环境因素，本文主要研究人为电磁设备辐射构成的三维电磁环境.由于各种电磁设备数量众多、使用频段交错、辐射功率大，构成的整个电磁环境异常复杂，相互间的干扰时有发生，尤其在目前高技术条件下，敌我双方的电磁对抗更加重了电磁环境的复杂性.由于电磁环境不可见，指挥员不能直观地了解战场中的电磁态势，所以要快速准确地做出决策存在巨大的困难.计算机图形技术能形象直观地展示数字信息，通过数据可视化，使用户能直观了解数据隐含的信息，带给用户161强烈的视觉感受.因此，电磁环境可视化在一定程

33、度上能减少指挥人员对复杂电磁环境掌握的盲目性，使指挥员对所处环境有一个直观的认识，为快速准确地决策提供支持.目前在虚拟战场环境中，电磁环境的可视化多采用二维的表现方式，不能体现三维电磁环境情况.本文结合虚拟现实技术，在三维数字地球上展示三维虚拟电磁环境，为指挥员提供详细的三维电磁态势，指挥员还能动态调整电磁设备参数，并直观查看电磁态势相应的变化，通过交互更加直接了解各种条件下的复杂电磁环境，为其快速准确地决策提供直观的辅助支持然而，把复杂的三维电磁环境高效且准确地展示到数字地球上是非常困难的为了可视化电磁环境，需采用合适的电波传播模型计算电磁设备传播数据虽然时域有限差分(FDTD)方法能精确描

34、述电波传播，但由于其需耗费大量的计算时间，难以满足实时动态可视化三维电磁环境的要求.Longley-rice电波传播模型是一种基于统计分析和电磁理论的电波传播计算模型17，它考虑了地形和大气的影响.通过数字地形模型，可采用Longley-rice。点到点的预测模式计算电波传播，它能精细考虑不规则地形对电波传播的影响，且该模式相比FDTD能较快计算电波传播数据，适合在数字地球上动态展示三维电磁环境.直接体绘制方法是一种可视化分析体数据的有效工具，由于图形硬件能力的提高和功能的增强，采用图形硬件图形处理器(GPU)加速的直接体绘制方法可在主流图形硬件条件下实时绘制，并能得到高质量的绘制结果目前对电

35、磁环境中雷达探测范围的可视化己有一些研究，主要采用了单层等值面的方法展示雷达三维探测范围，能直观展示雷达最大探测范围的包络;文献采用等值面提取的方法对高压电线周围的电场进行可视化，181可清晰展示特定电位大小的电场覆盖范围，但由于等值面方法对表现体数据细节的不足，尚不能展示电磁环境内部的细节信息多层半透明等值面绘制方法通过展示体数据的多层信息，能在一定条件下弥补单层等值面的不足，但没有对体数据进行信息分析的切片、切割能力直接体绘制技术能详细展示体数据信息，还能从多个角度切割分析数据，如文献采用直接体绘制方法绘制了电磁态势，但没有采用GPU加速的方法，不能实时展示电磁态势体数据因此本文扩展了文献

36、中的硬件加速直接体绘制方法，通过坐标转换把三维电磁环境直接实时绘制到数字地球上。4.2三维电磁环境体数据生成为了对虚拟电磁环境进行展示，首先计算虚拟环境中设置的各个电磁设备向三维空间中辐射电磁波的情况.本文采用Longleyice电波传播模型，即不规则地形模型来预测不同频段的电波传播，它能较准确地估计点波传播损失通过Longleyice电波传播模型得到各个19电磁设备在三维空间中的辐射损失值，进而计算得到电磁波三维空间功率密度数据，以此作为三维电磁环境表现的基础.4.2.1不规则地形剖面提取不规则地形模型需要地形剖面数据来精确计算地形对电波传播损失的影响，数字地形高程数据用来生成地形剖面数据如

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

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

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

40、111a,p,y,(n二0,1)分别为纬度、经度和海拔.图4.3(b)用填充色标示了三维电磁环境体数nnn据的剖面，不规则包围盒垂直截面ABCD用粗实线表示根据几何关系可明显得到包围盒的底部B点和C点海拔为y，顶部A点和D点的海拔为0(4-1)H二耳-Rcos0其中，R为地球半径；陽max(a-a,p-p)/2.从图3(b)可见，电磁环境体数据1010包围盒大于体数据的计算区域范围，这样可得到体数据包围盒海拔高度:底部为y，顶0部为H.ChyjifLI金诫祈图4.3电磁体数据包围盒4.3.3体数据纹理坐标转换在传统基于GPU光线投射体绘制方法中，包围盒顶点颜色值设置为三维体数据纹理的纹理坐标，

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

42、理的支持，包围盒浮点颜色值被绘制到了浮点纹理上，片段着色程序查询浮点纹理就能得到包围盒的顶点位置值，即得到投射光线的终止点.绘制包围盒正面作为投射光线的起始点22,因此投射光线方向就是终止点与起始点的矢量差，沿着光线方向就能采样得到该光线穿过的体数据值.但该光线上的位置是在笛卡儿坐标系下的坐标，不能直接从三维体数据纹理中得到体数据，需把笛卡儿坐标转化为纹理坐标.由于体数据计算范围为V(a,p,y)至V(a,p,y)，首先笛卡儿坐标00001111(x,y,z)转换为球坐标(r,a,P)，r=Qx2+y2+z2(4-3)za=arcsmrP=arctanx通过球坐标插值可得到纹理坐标为a-a0-

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

44、e2Duo2.8GHz，显卡为GeForee8600GT.图4.4示出了电磁环境在三维数字地球上的情况，电磁环境体数据分辨率为100x100x30,经度和纬度等分间隔都为0.O5度，高度等分间隔为200m覆盖范围约为500kmx500km，每个采样点数据采采用16位浮点数存储如果体数据范围较大，可把体数据划分为多个分块区域，每块采用同样的方法绘制，然后按照分块从后往前的顺序融合到三维环境中，因此本文方法适合较大数据量的电磁环境体数据绘制.釆用Longley-rice电波传播模型计算实验用电磁环境体数据耗时8.836s，窗口大小为1024x768时绘制平均帧率为12帧/s.相比文献，本文方法可实

45、时直接体绘制电磁环境，且各个电磁设备参数可交互改变，系统将重新计算生成新的电磁环境体数据.电磁环境体数据通过颜色映射工具(见图4.4(c)为不同的功率密度设置不同的颜色和透明度，三维电磁环境体绘制效果如图4.4(a)所示.电磁环境在三维环境中的情况可直观地动表4.1电磁设备参数编号频率/MHz功率编匕频率/MHz功H)HH)220025001235050035003000350600斗.15()15001斗30025005.32()2(HH)1548050062001500162(M)10(H)7500120017350450S40012001S35025009500

46、2)194(M)1500101025002()5(M)1200态展示，而且可交互改变电磁设备参数，能极大辅助电磁设备设计和规划.图4.4(b)示出了关闭图4.4(a)中一部电磁设备后的情况，可清楚看到该电磁设备影响情况.同电磁环境休绘制佝动态关闭电磁设备(C)颜色映射L鳧图4.4在三维地球上绘制的结果图4.5通过体数据切割技术示出了电磁环境在不同方向上的辐射情况，相比文献的等值面方法，本文方法能展示电磁体数据内部细节.其中图4.5(a)是500m海拔上电磁环境切片，图4.5(b)是5km海拔上的电磁环境切片，图4.5(c)是沿经度、纬度方向切割电磁环境的情况.(a)500ni髙度旳,5kni高

47、度（讨切块操作图4.5体数据切割这是一种传统体绘制方法的扩展方法，巧妙通过坐标转换把电磁环境体数据及时展示到了三维数字地球上.采用Longley-rice电波传播模型23计算电磁环境数据，该模型考虑了不规则地形和大气的影响，能较准确地预测电波传播.最后在数字地球上实现了三维电磁环境绘制.结果表明，本文方法能及时展示虚拟电磁环境，而且可动态调整电磁设备参数（如频率、功率等），及时展示不同电磁参数下电磁环境，可用于描述大范围三维电磁态势.虽然采用了Longley-rice传播模型来生成电磁环境体数据，但本文方法也适合其他方法生成的体数据的及时可视化，因此下一步计划采用更为准确的电波传播模型，并采用

48、并行计算的方式提高电磁数据计算效率，增强电磁环境体数据准确性.第五章用Longley-Rice大气折射修正方法的改进5.1 对Okumura-Hata和Longley-Rice的比较信道模型分为大尺度模型和小尺度模型两种，大尺度型刻画信号传输过程中经历的路径损耗中值，小尺度模型刻画小范围内（几个波长）信号的幅度和相位上的快速衰落。在大尺度传播衰落模型研究方面先后建立了以下著名的信道模型：Longley-Rice模型、在几公里到几百公里围内Okumura-Hata被广泛用来预测接收信号的场强，但海上移动信道路径损耗预测利用Longley-Rice模型更合适，Longley-Rice模型考虑了更多

49、的与地形有关的因素，包括面折射率、地面导电率、介电常数以及地面粗糙度等，还考虑了不同的气候类型和天线的位置标准等。文献通过对Okumura-Hata和Longley-Rice模仿真结果与实测数据的比较,指出预测海上传播损耗时后更确切。Longley-Rice模型适合于传输距离为1-2000公里，频率范围为2-2000MHZ，天线高度0-10000米，支持多种气候类型、多种地形、多种反射面传播特性、三种天线位置标准、天线激化方式等参数。该模型以传播理论为依据，同时结合了数千组实测数据，因此称其为半经验预测模型。Longley-Rice模型从1967年提出以来在当时航海等大区制通信场强预测模型中得

50、到了广泛的应用，并不断被深人研究和改进，Hufford等人根据模型提供的理论和原始数据给出了模型的实用算法。由于大气折射的影响，电波在大气中传播的实际路径是一条曲线，与收发信机之间的实际距离存在着误差，预测信道传播路径损耗时需要对实际距离加以修正。5.2 Longley-Rice模型中对大气折射误差的修正Longley-Rice模型中对大气折射误差的修正采用等效地球半径法，等效地球半径法计算简单，但是精度不高，该方法只利用了大气地面折射率的水平差异，未利用垂直剖面信息，而实际大气折射率垂直差异，尤其在近地面lkm内的垂直差异比较大，本文对Longley-Rice模型大气折射误差修正方法加以改进

51、。改进方法充分利用大气折射率垂直剖面信息，利用电波实际传播路径和收发天线地面水平距离之间的映射关系，将模型输人的收发天线之间的实际距离映射为电波实际传播距离，能更精确反映传播距离对电波传播损耗的影响。改进算法预测结果比原有方法更接近实测结果，具有较大的实用价值。5.3大气折射修正方法的改进Longley-Rice模型参考衰减为距离的分段函数max(0,L+kd+klgd),d<d<dbe12minLsL=<L+m,d<d<d(5-1)refbeddLsxL+md,d<dbessx定义的三段分别称为视距区域，绕射区域和散射区域。模型控制置信度的参数由统计量给出

52、，L(qT,qL,qS)作为时间，位置，形势分量的函数，衰减不会超出这个值。也就是说，在必的情况下，至少有私的地点衰减不会超过L(qT,qL,qS),并且这个时间至少为qT。对于光滑地球地平线d的距离，是基于平坦地面的二径反射理论和一个绕射损耗Lsb的外插值计算的;对恰好超过从d到d，的地平线距离(为绕射损耗和散射损耗相等的地Lsx点)，参考衰减值是刃峰绕射和光滑地球绕射计算的加权平均。加权因子为频率、地形不规则因子和天线高度的函数。对于很不规则地形，从终端看地平线障碍物可看做是锐利山脊，绕射损耗是基于Epstein-Peterson近似按双刃峰路径计算的。对范围大于的超视距路径，参考衰减按按

53、绕射损耗计算和正向散射损耗计算中取较小者。由于大气层密度分布的不均匀性，导致电波传播产生折射现象，使得电波传播距离与实际距离不符。在Longley-Rice模型中，传输距离参数给出的是发射机和接收机之间的地面水平距离，但影响传输损耗的距离是电波实际传播的距离(视在距离)，并且该模型预测损耗值是传播距离的分段函数，传播距离直接影响传输衰落曲线的转折点，因此必须对给出的传播距离参数加以修正。Longley-Rice模型采用等效地球半径法修正大气折射误差，等效地球半径定义为散射梯度N或地球表面折射值N的函数。定义N=(n-1)x10A6其中，。为地球表面大气折射指数。该方法对折射率剖面的计算没有利用

54、折射率负梯度的值在不同气候类型之间有较大差异的特性，实际中这种差异比较明显，而折射率剖面直接影响到电波传播的实际路径，从而影响预测精度。下面提出一种改进的折射修正方法，以改进Longley-Rice电波传播衰落预测精度。电波折射示意图如图1所示，T为发射天线相对平均海平面的高度，T'为T在海平面上的投影，R为接收天线相对海平面的高度，R'为R在海平面上的投影，0为地心原点，OTOR为地球半径，P为射线上的某点，0为射线在P点处的仰角，dl为射线在P点处的一小段弧线，dr为地球径向延长线在P点处的线段小量，d申为射线上P点对应的地心张角小量，R'T'为收发天线之间

55、的地面水平距离，ds为R'T'在P点与地心连线和地球表面的交点处水平距离小量，TT'为发射天线高度，单位m。RR'为接收天线高度，单位m,益为电波传输视在距离，单位km,n为P点处的电波折射率，h=OP,9°射线在T点处的初始仰角。图1电波传播折射几何图利用Longley-rice模型计算电波传输路径损耗时，收发天线间的距离参数给出的是地面水平距离，即R'T'，而实际电波传播距离为TR,下面根据Snell定理，利用射线描迹法推导R'T'和TR之间的关系式。首先推导地面水平距离RT计算公式:由图2知:drds=OT'dq,dq=cot0-h从而,drds=OT'cot0h进而得,(5-2)T'R=S=TOT'cot01drhTT'由SNELL定理可导出，cot0=nr0-0cos0nr0阻cos0nrcos00zQ:n2r2n2r2cos201000(5-3)nr从而，mRR'OTrcos0TR=J&=o0drn2r2n2r2cos20TT'000下面推导实际电波传播距离TR的计算公式:(5-4)设光速为c,电波从发射天线T传到R的时间为t,则1ORn,ct,t=drc一sin0TR