北邮电磁场实验,场强仪测校园磁场分布

1、无线信号场强特性的研究2009211111班孙瑞&祁剑飞 20132014学年第二学期 电磁场与电磁波实验报告 专 业 通信工程 班 级 2011211109 姓 名 唐天宇、廖浩帆 学 号 2011210258、2011210259 班内序号 15、16 开课系室 信通院实验中心 报告日期 2014年5月 1、 实验目的1. 掌握在移动环境下阴影衰落的概念以及正确测试方法。 2. 研究校园内各种不同环境下阴影衰落的分布规律。 3. 掌握在室内环境下的场强的正确测试方法，理解建筑物穿透损耗概念。 4. 通过实地测量，分析建筑物穿透损耗随频率的变化关系。 2、 实验原理无线通信系统是由发

2、射机，发射天线，无线信道，接收机，接收天线所组成。对于接收者，只有处在发射信号的覆盖区内，才能保证接收机正常接收信号， 此时，电波场强大于等于接收机的灵敏度。因此，基站的覆盖区的大小，是无线工程师所关心的。决定覆盖区的大小的主要因素有：发射功率，馈线及接头损耗， 天线增益，天线架设高度，路径损耗，衰落，接收机高度，人体效应，接收机灵敏度，建筑物的穿透损耗，同播，同频干扰。（1） 大尺度路径损耗在移动通信系统中，路径损耗是影响通信质量的一个重要因素。 大尺度路径损耗：用于测量发射机与接收机之间信号的平均衰落，即定义为有效发射功率和平均接收功率之间的(dB)差值，根据理论和测试的传播模型，无论室内

3、或室外信道，平均接收信号功率随距离对数衰减，这种模型已被广泛地采用。对任意的传播距离，大尺度平均路径损耗表示式为：PLddB=PLd0+10nlog(d/d0)即平均接收功功率为：PrddBm=PtdBm-PLd0-10nlogdd0=Prd0dBm-10nlogdd0其中，n为路径损耗指数，表明路径损耗随距离增长的速度；d0为近地参考距离； d为发射机与接收机(T-R)之间的距离。横杠表示给定值d的所有可能路径损耗的综合平均。坐标为对数-对数时，平均路径损耗或平均接收功率可表示为斜率10ndB/10倍程的直线。n值取决于特定的传播环境。例如自由空间，n为2，当有阻挡物时，n比2大。决定路径损

4、耗大小的首要因素是距离，此外，它还与接收点的电波传播条件密切相关。为此，我们引进路径损耗中值的概念。中值是使实测数据中一半大于它而另一半小于它的一个数值(对于正态分布中值就是均值)。人们根据不同的地形地貌条件，归纳总结出各种电波传播模型：1) 自由空间模型； 2) 双径模型；3) Hata模型；4) Hat-cost231模型；5) Okumura 模型。（2） 阴影衰落在无线信道里，造成慢衰落的最主要原因是建筑物或其它物体对电波的遮挡。在测量过程中，不同位置遇到的建筑物遮挡情况不同，因此接收功率也不同，这样就会观察到衰落现象。由于这种原因造成的衰落也叫“阴影效应”或“阴影衰落”。在阴影衰落的

5、情况下，移动台被建筑物所遮挡，它收到的信号是各种绕射、反射、散射波的合成。所以，在距基站距离相同的地方，由于阴影效应的不同，它们收到的信号功率有可能相差很大，理论和测试表明，对任意的d值，特定位置的接受功率为随机对数正态分布即：PrddBm=PrdBm+X=Prd0dBm-10nlogdd0+X其中，X为0均值的高斯分布随机变量，单位为dB，标准偏差为，单位也是dB。 对数正态分布描述了在传播路径上，具有相同T-R 距离时，不同的随机阴影效应。这样利用高斯分布可以方便地分析阴影的随机效应。正态分布，也叫高斯分布，它的概率密度函数是：fx=122e-(x-m)222应用于阴影衰落时，上式中的x

6、表示某一次测量得到的接收功率，m 表示以dB 表示的接收功率的均值或中值，s 表示接收功率的标准差，单位是dB。阴影衰落的标准差同地形，建筑物类型，建筑物密度等有关，在市区的150MHz 频段其典型值是5dB。除了阴影效应外，大气变化也会导致阴影衰落。比如一天中的白天，夜晚，一年中的春夏秋冬，天晴时，下雨时，即使在同一个地点上，也会观察到路径损耗的变化。但在测量的无线信道中，大气变化造成的影响要比阴影效应小的多。下面是阴影衰落分布的标准差，其中s(dB)是阴影效应的标准差。s s (dB)频率（MHz）准平坦地形不规则地形h(m)城市郊区501503001503.5-5.54-79111345

7、067.51115189006.58141821（3） 建筑物的穿透损耗的定义建筑物穿透损耗的大小对于研究室内无线信道具有重要意义。穿透损耗又称大楼效应，一般指建筑物一楼内的中值电场强度和室外附近街道上中值电场强度dB 之差。 发射机位于室外，接收机位于室内，电波从室外进入到室内，产生建筑物的穿透损耗，由于建筑物存在屏蔽和吸收作用，室内场强一定小于室外的场强，造成传输损耗。室外至室内建筑物的穿透损耗定义为：室外测量的信号平均场强减去同一位置室内测量的信号平均场强。用公式表示为：P= 1 Ni=1NPi(outside)-1Mj=1MPj(inside)P是穿透损耗，单位是dB，Pj是在室内所测

8、的每一点的功率，单位是dBv，共M 个点，Pi是在室外所测的每一点的功率，单位是dBv，共N个点。3、 实验内容利用DS1131 场强仪，实地测量学六内外信号场强 （1）研究学六内外不同环境下阴影衰落分布规律，以及具体的分布参数如何； （2）研究建筑物穿透损耗的变化规律。4、 实验步骤（1）数据采集频点486494MHz (中央3)地点学六外四周及楼内4层楼道测量步长1个数据/步长说明：=cf=0.8m，每步长平均约60cm，相当于一个波长内测量1.3个数据。（2）数据录入（原始数据见附录1） 将原始数据录入Excel表格：在不同工作表下依次录入不同区域的实测数据，单列存储，便于读取。（3）数

9、据处理 1）用Matlab对数据进行编程处理绘制相关曲线或图形：实测数据曲线、信号强度图谱、概率分布柱状图、累积分布曲线等；计算数据特征量：均值、标准差、中值等。2）得出测试区域阴影衰落的分布规律，研究路径损耗与楼层关系，并计算学六内外的穿透损耗。5、 程序代码 程序代码根据所要实现的内容进行模块化编程，编写了5个函数文件和3个命令文件（代码详见附录2）。函数文件：函数名函数功能datain读入数据datacurve绘制实测数据曲线，并叠加均值线colorscheme绘制信号强度图谱histogram绘制概率分布柱状图，并叠加拟合的正态分布概率密度曲线distribution绘制累积概率分布图

10、及相对应的正态分布概率累积曲线命令文件：文件名函数功能No1绘制不同区域实测数据线、概率分布柱状图、累积概率分布图No2绘制学五四周信号强度图谱No3同坐标显示4楼楼道实测数据曲线、绘制楼层信号强度图谱、绘制各楼层数据均值与楼层关系曲线6、 实验结果（1）学五楼外不同区域阴影衰落 1）东侧（NS）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量） 2）北内侧（EW）l 实测数据曲线 l 样本分布柱状图l 累积概率分布图（含数据特征量）3）南内侧（WE）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量）4）北外侧（WE）l 实测数据曲线l 样本分布柱状图l 累

11、积概率分布图（含数据特征量）5）西内侧（NS）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量）6）5楼北侧楼道（EW）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量）7）5楼西侧楼道（NS）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量）8）5楼南侧楼道（WE）l 实测数据曲线l 样本分布柱状图l 累积概率分布图（含数据特征量）5）学六四周信号强度图谱（2）学六内5楼楼道信号强度图谱7、 结果分析各区域数据特征：最小值(dBmw)最大值(dBmw)均值(dBmw)标准差外_北侧-80.2-59.9-69.158623.958087外_

12、东侧-85.4-67.2-74.513734.008712内_北侧-80-63-71.576274.224923内_南侧-86-66.5-75.866674.110818内_西侧-86.7-69.2-75.7253.7990665楼_北侧-87.2-72.3-82.391782.7458495楼_南侧-88.7-79.4-84.266152.6291645楼_西侧-89.1-81.2-85.076922.090013实测数据说明：l 学六院外侧测试数据全部在区间-60,-86dBmw内，学六院内侧测试数据区间为-63,-87；且外侧平均值均大于院内侧平均值。l 一层北侧的信号（不管是内侧还是外

13、侧）都比其他方向较为好；l 学六五层的测试数据区间为-72，-89，楼道内的信号和宿舍楼外的信号相比明显差了很多，测试数据有所不足的是楼道内的数据测试的是五层的而不是一层，因为学六的一层结构更为复杂，并且信号太差，所以测试的是宿舍楼的最高层。l 学六五层和一层的数据相比，稳定了很多，标准差明显较小，这与周围建筑遮挡，学六本身墙面遮挡有关。 分析说明： （1）阴影衰落分布规律随机对数正态分布 从实验结果中各区域实测数据的概率分布柱状图与拟合正态分布，及累积概率曲线和对应正态分布函数曲线(均值、方差相同)对比可得，无线信号强度为随机对数正态分布。 （2）穿透损耗由学六内外测试数据均值对比可知，学六

14、外的信号强度比学六内的信号强度大，这反映出建筑物对信号的屏蔽和吸收作用，即建筑物的穿透损耗。 应用公式： P= 1 Ni=1NPi(outside)-1Mj=1MPj(inside) 由学六外测量信号均值减去五楼楼道信号均值，可得穿透损耗： P=(-73.368058)-( -83.911617)=10.543559dBmw （3）空间结构对信号强度的影响所选的学六宿舍楼的空间结构比较复杂,但是可以从信号强度谱图中看出,楼道中靠近楼梯和窗口的地方信号强度都比较好,而被建筑物包围的内侧信号相对于外侧信号而言有了明显的衰减,尤其是在拐角的阴影处,信号更差,而越是测量地点越是开阔的地方,信号强度则越

15、强.当然,在测试楼道数据时,因为测量时有些宿舍的门窗是打开的,一定程度上会影响所测量的数据.（4）楼层与信号强度关系 因为时间和学六楼层结构的缘故,我们并没有测量一层楼道内的数据,所测量的数据并不具有楼层与信号强度的对比性,不过可想而知的是,在一定程度范围内(楼层高度相对有限)的情况下,越是底层的信号的穿透损耗越高,越是高层的楼层的空间结构对信号强度的影响越小实验总结 本次实验中，我们小组合理分工，团结协作，通力完成实验要求内容。在选择所要测量的地点时，我们选择了和我们平时生活息息相关并且空间结构更为复杂的学六宿舍楼。从之前的图片上也可以了解到，学六宿舍楼不仅结构上和学校内的所有其他建筑都不一

16、样，并且也处于校园的边缘。在原先的设想中，我们还想测量学六一层楼道、三个不同位置上的楼梯、内外侧拐角宿舍内外、内外侧正常普通宿舍内外以及相同宿舍内部，不同床位上的信号强度。因为在我们平时的生活中，往往是靠窗边的同学才能收得到手机信号（移动为例），而靠门的同学则根本收不到手机的信号。因为时间精力有限，并且我们的场强仪的测量信号频率的上限够不到中国移动联通的下行频段，所以只选择了中央3频道在学六宿舍楼内外侧的几个方向的测量数据。同时，当我们测试已经进行了一半之后我们才发现学六宿舍附近因为施工的原因，宿舍楼的南外侧和西外侧无法测量数据，给本次实验带来不小的麻烦。 总的来说，通过本次实验，我们学会了如

17、何使用场强仪，并且对于信号的阴影衰减以及穿透损耗等概念有了个清晰的认识，还对于通信系统的统计性和随机性有了一定的认知，对于理论课上所学到的知识更加深刻。附录附录1 原始数据登记表(见纸板报告)附录2 Matlab文件 函数文件： datain.m function Data,L=datain(s)data=xlsread('data.xls',s); %从Excel表格中读取工作表s的数据Data=data' %将数据转置为行向量 L=length(Data); %数据长度end datacurve.m function datacurve(Data,L)n=1:1:L

18、;m=ones(1,L)*mean(Data); %数据均值向量plot(n,Data,n,m) %同坐标绘制实测数据曲线及其均值线legend('电平值曲线','电平值均值')xlabel('测量点')ylabel('电平值x(dBmw)')gridend colorscheme.m function colorscheme(Data,L,az,el)x=1:1:L; %x为测量数据点向量，维数为1*Ly=linspace(0,1,2); %y为绘图辅助向量，维数为1*2(向量取值范围和列数也可设置为其他值)x,y=meshgr

19、id(x,y); %生成xy平面网格坐标矩阵z=Data;Data; %z的维为2*Lsurf(x,y,z); %绘制三维曲面图xlabel('测量点')axis(1,L,0,1)view(az,el) %观察视角end histogram.m function histogram(Data)histfit(Data,30); legend('测量样本分布','拟合正态分布');xlabel('电平值(dBmw)');ylabel('样本数量(个)');grid;end distribution.m functio

20、n distribution(Data) h,stats=cdfplot(Data);%绘制累积概率分布,stats为包含数据特征的结构体hold on %图形保持 %添加与原始数据相对应(均值&标准差)的正态分布的分布函数曲线x=(stats.min-5) : 0.1 : (stats.max+5); mu=stats.mean;sigma=stats.std;f=normcdf(x,mu,sigma); %根据原始数据特征(均值&标准差)生成对应正态分布plot(x,f,'r') %绘制正态分布的分布函数曲线xlabel('电平值x(dBmw)

21、9;)legend('电平值分布函数','对应正态分布函数') %添加数据特征说明text(stats.median,0.45,'最大值=',num2str(stats.max),'dBmw','FontSize',14)text(stats.median,0.35,'最小值=',num2str(stats.min),'dBmw','FontSize',14)text(stats.median,0.25,均 值=',num2str(stats.mean),'dBmw','FontSiz