资源描述
电磁场与电磁波实验报告题目：校园无线信号场强特性的研究 姓名班级学号序号 目 录一、实验目的2二、实验内容2三、实验原理5四、实验步骤51、实验对象选取52、数据采集5五、实验数据21、原始数据录入72、数据处理流程7六、实验结果与分析81、主楼周边电磁场信号强度分析82、主楼室内不同楼层楼道信号强度分析11七、问题分析与解决151、Matlab 仿真问题研究与解决232、场强分布的研究23 3、模型拟合.24八、分工安排及心得体会25附录I：原始数据26附录II：源代码30一 实验目的1. 掌握在室内环境下场强的正确测试方法，理解建筑物穿透损耗的概念；2. 通过实地测量，分析建筑物穿透损耗随频率的变化关系；3. 研究建筑物穿透损耗与建筑材料的关系。4. 掌握在移动环境下阴影衰落的概念以及正确测试方法。 二实验内容 利用DS1131场强仪和拉杆天线，实地测量信号场强。1. 研究具体现实环境下阴影衰落分布规律，以及具体的分布参数如何；2. 研究在校园内电波传播规律与现有模型的吻合程度，测试值与模型预测值的预测误差如何；3. 研究建筑物穿透损耗的变化规律三实验原理无线通信系统是由发射机、发射天线、无线信道、接收机、接收天线所组成。对于接收者，只有处在发射信号覆盖的区域内，才能保证接收机正常接收信号，此时，电波场强大于等于接收机的灵敏度。因此，基站的覆盖区的大小，是无线工程师所关心的。决定覆盖区大小的因素主要有：发射功率、馈线及接头损耗、天线增益、天线架设高度、路径损耗、衰落、接收机高度、人体效应、接收机灵敏度、建筑物的穿透损耗、同播、同频干扰。【阴影衰落】阴影衰落是电磁波在空间传播时受到地形起伏、高达建筑物群的阻挡，在这些障碍物后面会产生电磁场的阴影，造成场强中值的变化，从而引起信号衰减。阴影衰落的信号电平起伏是相对缓慢的，又称为慢衰落，其特点是衰落与无线电传播地形和地物的分布、高度有关。在无线信道里，造成慢衰落的最主要原因是建筑物或其他物体对电波的遮挡。在测量过程中，不同测量位置遇到的建筑物遮挡情况不同，因此接收功率不同，这样就会观察到衰落现象。在阴影衰落的情况下，移动台被建筑物遮挡，它所收到的信号是各种绕射、反射、散射波的合成。阴影衰落一般表示为电波传播距离r的m次幂与阴影损耗的正态对数分量的乘积。移动用户和基站之间的距离为r时，传播路径损耗和阴影衰落可以表示为：其中，为0均值的高斯分布随机变量，单位为，标准偏差为，单位也。对数正态分布描述了在传播路径上，具有相同的T-R距离时，不同的随机阴影效应。这样利用高斯分布可以方便的分析阴影的随机效应。它的概率密度函数是：应用于阴影衰落时，上式中的x表示某一次测量得到的接受功率，m表示以表示的接收功率的均值或中值，表示接收功率的标准差，单位为。阴影衰落的标准差同地形、建筑物类型、建筑物密度等有关，在市区的频段其典型值是。阴影衰落分布的标准差（dB）频率（MHz）准平坦地形不规则地形（米）城市郊区501503001503.55.5479111345067.51115189006.58141821 【大尺度路径衰落】大尺度衰落是由移动通信信道路径上的固定障碍物（建筑物、山丘、树林等）的阴影引起的，衰落特性一般服从律，平均信号衰落和关于平均衰落的变化具有对数正态分布的特征。利用不同测试环境下的移动通信信道的衰落中值计算公式，可以计算移动通信系统的业务覆盖区域。实际上，大尺度衰落不仅与时间有关，还与距离和载波频率有关。为了表达方便，上式中省略了距离因子d和载频。基于理论和测试的传播模型指出，无论室内还是室外信道，平均接收信号功距离的对数而衰减。 (2.2)或 (2.3)式中，n为路径损耗指数，表明路径损耗随距离增长的速率；是近地参考距离，由测试决定；d为发射机和接收机距离。在自由空间传播时，n为2，当有障碍物时，n变大。在移动通信系统中，路径损耗是影响通信质量的一个重要因素。大尺度平均路径损耗：用于测量发射机和接收机之间信号的平均衰落，定义为有效发射功率和平均接受功率之间的（dB）差值，根据理论和测试的传播模型，无论室内或室外信道，平均接受信号功率随距离对数衰减，这种模型已被广泛的使用。对任意的传播距离，大尺度平均路径损耗表示为：即平均接收功率为：其中，为路径损耗指数，表明路径损耗随距离增长的速度；为近地参考距离；为发射机与接收机之间的距离。测试表明，对于任意d，特定位置的路径损耗又服从随机正态分布，对数正态分布描述了在传播路径上，具有相同距离时，不同的随机阴影效应。这种现象叫对数正态阴影。 决定路径损耗大小的首要因素是距离，此外，它还与接收点的电波传播条件密切相关。为此，我们引进路径损耗中值的概念。中值是使实测数据中一半大于它而另一半小于它的一个数值（对于正态分布中值就是均值）。人们根据不同的地形地貌条件，归纳总结出各种电波传播模型。1. 自由空间模型2. 布灵顿模型3. EgLi模型4. Hata-Okumura模型 【建筑物的穿透损耗】1. 定义建筑物的穿透损耗大小对于研究室内无线信道具有重要意义。穿透损耗又称大楼效应，一般指建筑物一楼内的中值电场强度和室外附近街道上中值电场强度之差。发射机位于室外，接收机位于室内，电波从室外进入到室内，产生建筑物的穿透损耗，由于建筑物存在屏蔽和吸收作用，室内场强一定小于室外的场强，造成传输损耗。室外至室内建筑物的穿透损耗定义为：室外测量的信号平均场强减去在同一位置室内测量的信号平均场强，用公式表示为：是穿透损耗，单位，是在室内所测的每一点的功率，单位，共M个点，是在室外所测的每一点的功率，单位，共N个点。电磁波透射墙体示意图2.影响因素无线信号对透射能力是频率及高度的函数，天线对信号透射也有非常重要的影响。以室外发射机的天线高度远小于建筑物本身高度考虑，在测试中显示，随着频率的增加透射损耗（净穿透损耗）减小。例如下表数据：（卢东如2014.楼宇信号强度及穿透损耗研究）地点（状态）频率（MHz）透射损耗（dB）建筑物底层（第一组比较）44116.4896.511.614007.6建筑物底层（第二组比较）90014.2180013.4230012.8窗前比没有窗户小6dB透射损耗即净穿透损耗是指不考虑绕射等因素综合影响，墙体和玻璃等物体的穿透损耗，在绕射和反射等几乎没有或不明显的环境下测量得到的，近似穿透损耗。3.不同建筑材料各类墙体净穿透损耗：（卢东如 2014.楼宇信号强度及穿透损耗研究）墙体类型损耗（dB）普通窗户玻璃（35mm）1落地玻璃墙（2030mm）3.6建筑物外墙（钢筋混泥土）1926建筑物内墙（砖石）710钢筋混凝土类建筑贯穿损耗中值：数据自（通信人家园 论坛|中国第一通信社区）建筑物类型贯穿损耗中值（dB）标准偏差（dB）小城市市区钢筋混凝土类107.3小城市郊区钢筋混凝土类5.88.7大城市市区钢筋混凝土类187.7大城市郊区钢筋混凝土类13.19.5大城市市区金属框架类27-由于我国的城市环境与国外有很大的不同，贯穿中值一般比国外同类名称高8-10dB。四实验步骤1. 实验对象选取这次实验数据采集地点我们选择了主楼（图），主要测量三个部分的信号强度：1. 靠近主楼测量其周边一圈信号强度，由西北角开始顺时针绕主楼测量（图）。测量时分为东南西北四个方向，其中东西侧较长，南北侧较短，近似为一个矩形；2. 测量1,2,4,6,10,13楼室内信号强度，由西北角开始顺时针绕室内走廊测量，记录特殊结构造成的信号强度波动。由于主楼每层建筑结构大致相同，变量单一便于对比不同楼层高度对信号强度的影响；3. 测量电梯内部信号强度，电梯从上而下运行过程中记录测量数据。 北邮主楼卫星图北邮主楼3D指示图同时，经过查证，电视信号发射塔位于测量地点的西南方向约6.5公里处（图）。信号发射塔方位图在选频方面，频点选用的是低频190.758MHz（CH8，中央2套伴音频道频点），因为这个频点信号接收很好，且本班另一测量小组将在教二室内测量相同频点的信号强度，我们两组将进行钢筋混泥土建筑结构和砖石建筑结构的对比。2. 数据采集本组频点选择：190.758MHz（CH8，中央2套伴音频道频点），利用场强仪DS1131测量无线信号的强度（单位），估测频点相应半波长伟0.78米，每走1步约为0.8米（恰好为楼道地砖2格）读一次数并进行信号电平记录。共测量的近1000组数据进行分析比较，分别是主楼外北侧、东侧、南侧、西侧以及主楼一层、二层、四层、六层、十层、十三层，电梯内部。在测量建筑物外的信号强度时，采用自西北角起顺时针方向围绕该建筑物一周的测量方法。在测量室内信号时，采用自西北角起顺时针方向围绕该楼层楼道一周的测量方法。（图）。测量当天天气情况：多云转晴、微风。场强仪倾斜45使用，天线仅拔出第一节。主楼内部结构示意图五实验数据1.原始数据录入将测量得到的数据录入Excel表格，分别以东、南、西、北以及各个楼层方向为单独的sheet，便于数据处理时进行MATLAB的数导入和获取。记录时间、地点、电平值的同时记录测量时周围的可能影响电平值因素，便于后期数据分析(表)。 数据记录表格原始实验数据见附件2. 数据处理流程采集到的数据有近1000组，需要对数据进行细致的处理以便得到明确的结论。下图所示为数据处理的流程图：3.数据处理MATLAB关键语句本次实验数据的处理主要利用EXCEL和MATLAB两个软件完成，其中MATLAB中的关键处理语句主要分为以下几个部分：1. 原始数据读入MATLABdata=xlsread(data.xlsx,10);Data=data; %将数据转置为行向量 Data=-1.*Data;L=length(Data); %数据长度功能：录入数据，行列转置后各数据取相反数并计算数据个数。2. 作信号强度散点图x=(1:1:L);plot(x,Data,Linewidth,1.3, Color,0.3,0.3,1);meandata=mean(Data);hold on;plot(x,meandata,r., Color,1,0.4,0.6);功能：根据原始数据作出各信号强度的散点分布和平均值曲线。3. 作信号强度频率分布直方图y,x=hist(Data,10)%画样本数目直方图y=y/L/mean(diff(x);h1=bar(x,y,1);axis(-80,-45,0,0.20);grid on功能：根据原始数据划分十个信号强度分组，计算落于每一信号强度范围中的测量数据个数，作出频率分布直方图。4. 作统计概率曲线和拟合正太分布曲线y=normpdf(x,mean1,std1);%在相同均值标准差下，画正态分布图plot(x,y,r, Linewidth,1.8, Color,1,0.4,0.6)axis(-80,-45,0,0.20);n,m=hist(Data);%在同一坐标系下，画出统计概率图plot(m,n/L/mean(diff(m), Linewidth,1.8, Color,0.3,0.3,1 )功能：根据原始数据的频率分布直方图作出统计概率曲线，并利用测量数据的均值和方差作出拟合的正态分布曲线。5. 作信号强度表面图Data2=Data,zeros(1,L),(1:L);%转换矩阵Data2=reshape(Data2,L,3);surf(Data2);%画surf图，表示统计分布功能：根据原始数据用不同颜色表示不同点的信号强度，可更加直观的展示出信号强度的变化和分布规律。6. 作信号强度概率累积分布图cdfplot(Data)功能：根据原始数据作出信号强度概率累积分布图。六实验结果与分析 1. 主楼外围电磁场强度及其分析 （1）主楼外围北侧图 主楼外围北侧 最大值 最小值 平均值 标准差 -56.1dbmw -68.7dbmw -60.98dbmw 2.86由主楼外围北侧信号强度散点图（图）中标记的A点的电平值可知，信号强度突然降低，参见原始数据表格备注知该处在测量数据时有大量自行车及机动车停放于测试点周围，影响信号强度。（2）主楼外围东侧 最大值 最小值 平均值 标准差 -52.6dbmw -69dbmw -60.22dbmw 3.23（3）主楼外围南侧 最大值 最小值 平均值 标准差 -48.7dbmw -67.4dbmw -56.57dbmw 4.16 （4）主楼外围西侧 最大值 最小值 平均值 标准差 -59dbmw -67.8dbmw -58.22dbmw 4.02由信号强度散点图中标记的AB段的电平值可知，信号强度突然降低，参见原始数据表格备注知该处在测量数据时突然刮起一阵强风，可能影响信号强度。（5）主楼外围总体a.散点图、直方图及正太拟合 最大值 最小值 平均值 标准差 -48.7dbmw -69dbmw -58.96dbmw 3.98 由实测的主楼外围信号强度的平均值及标准差大致判断本次测量选取的CH8频道（190.758MHz）信号强度较好，适宜测量。b.概率累积分布图由概率累积分布曲线（图）可知，主楼外围信号强度概率分布较为接近S型曲线，说明此次测量时主楼外围总体信号环境较好，突发干扰较少。c场强空间分布图主楼外围方位 最大值(dbmw) 最小值(dbmw) 均值(dbmw) 标准差 北侧-56.1-68.7-60.982.86 东侧-52.6-69-60.223.23 南侧-59-67.8-58.224.02 西侧-56.1-68.7-60.982.86总体-48.7-69-58.963.98（6）主楼外围周边信号强度分析：l 不同方位信号强度对比：可以看出南侧的信号强度最高，因为主楼南侧的遮挡物较少、距离教二楼较远且建筑物较低，比较空旷，因此信号衰减小，信号强度高；主楼北侧的区域有教一楼遮挡，东侧有科学会堂遮挡，因此受到阴影衰落影响较严重，信号强度较差。同时根据场强空间分布图知西南角信号强度最好，经考证该方位为发射塔方向。因此此次测量的数据与实际情况相符。l 高低频信号强度对比： 结合同班另一测量小组（组员：王文雪，张玉洁），所测量的教二楼外部频道2（57.75MHz）的信号强度数据，进行高低频的对比： 信号电平 西侧 南侧 东侧 北侧 高频（535.25MHz） -70.1778 -65.8644 -76.0391 -69.4128低频（57.75MHz） -61.0018 -58.9465 -61.0298 -58.8750可以看出，无论是哪个方位，高频信号的衰减明显要大于低频，这与理论相符。同时低频也是南北侧信号强度高于东西侧。 2.主楼各楼层楼道电磁场强度及其分析（1）主楼一楼走廊 a. 散点图、直方图及正太拟合 最大值 最小值 平均值 标准差 -58.7dbmw -75dbmw -70.13dbmw 3.25b. 概率累积分布图c场强空间分布图 由场强空间分布图可知主楼一层西侧大厅中央处信号强度最好，根据实测数据表格备注可知此处为大厅开门处，且大门材料为玻璃，穿透损耗较低，符合不同材料对穿透损耗影响的相关规律。（2）主楼二楼走廊 a. 散点图、直方图及正太拟合 最大值 最小值 平均值 标准差 -66.5dbmw -75.4dbmw -72.11dbmw 1.81比较一楼和二楼的电平特征值，二楼电平均值更低，平均信号更差，但是标准差明显小于一楼，可排除一楼大厅开门处引起的信号强度分布不均，波动较大的影响，数据更具有代表性。故在分析楼内信号数据，并与外围数据进行比较分析穿透损耗时舍弃一层数据。b. 概率累积分布图（3）主楼四楼走廊（以四层为例分析楼层内各处影响信号强度因素）a.直方图及正太拟合 最大值 最小值 平均值 标准差 -53dbmw -74.9dbmw -68.35dbmw 4.33b. 概率累积分布图c信号强度散点图参考实验数据表格备注分析数据，现对信号强度散点图中标注的峰值进行说明分析。l AB段左侧有厕所，信号非常规波动，且每层楼均有此现象，目前原因未知有待探究；l C点左侧有铁门，引起信号强度陡然降低，反映出铁质建筑材料相比钢筋混凝土建筑材料对信号衰减影响更严重，符合不同建筑材料对穿透损耗影响的相关规律；l D点左侧有门可通外侧楼梯，信号强度较好；l E点右侧有一开门实验室，实验室内部四面墙壁附有特殊材料，为蓝色泡沫质感，且表面呈内凹蜂窝状，此处信号强度陡增，明显优于周围信号强度，推断为该材料使反射信号较强，对信号有明显的聚集效果。l F点左侧一外围实验室向着室内方向开着门，室外信号可穿透玻璃窗进入室内，使此处信号强度较好，明显高于均值。d场强空间分布图 由场强空间分布图配合散点图（对应处已标注相同字母），可以更直观比较出楼层内各处影响信号强度因素。如E处因有开门的特殊实验室信号强度明显优于周围各处。（4）主楼六楼走廊 a. 散点图、直方图及正太拟合 最大值 最小值 平均值 标准差 -47.5dbmw -74.9dbmw -63.41dbmw 5.65b. 概率累积分布图与低楼层信号强度均值相比，六层信号强度明显增强，由于主楼周围并没有高层建筑遮挡，此处信号衰减主要由建筑物穿透损耗引起，可忽略阴影衰落对信号强度的影响。（5）主楼十楼走廊a. 散点图、直方图及正太拟合 最大值 最小值 平均值 标准差 -60dbmw -73.4dbmw -63.35dbmw 5.87 b. 概率累积分布图c场强空间分布图由场强空间分布图可直观比较得出，西南侧信号强度明显好于东北侧，由于十层楼层较高无其他建筑物遮挡，仅考虑穿透损耗造成的影响。推断为电视信号发射塔位于主楼西南侧造成此信号分布不均现象。（6）主楼十三层楼道a. 散点图、直方图及正太拟合 最大值最小值平均值标准差-50.9dbmw-73.4dbmw-62.60dbmw4.84 b. 概率累积分布图3.主楼各楼层总体信号强度综合分析1 不同楼层信号强度对比：主楼楼层最大值(dBmw)最小值(dBmw)均值(dBmw)标准差一层-58.7-75-70.133.25二层-66.5-75.4-72.111.81四层-53-74.9-68.354.33六层-47.5-74.9-63.415.65十层-60-73.4-63.355.87十三层-50.9-73.4-62.604.84可以看出，除去受大厅开门处影响较大的一层数据后，随着楼层的升高，信号强度逐渐增强，在数据中显示到达6层及其以上楼层信号强度明显增加，分析其原因是高楼层对应的室外比较空旷，没有周围高层建筑物的阻挡，因此衰减很小。（图）同时发现室内楼道测量时数据的标准差很大，这是因为每层楼两侧部分空间狭小，信号较弱衰减很大，而中心楼梯口处一般比较空旷，信号强度较强，因此差别较大。从散点图中也可以看出，楼道两侧的信号衰减较大，而中间部分信号衰减较小。 在各测量部分中，电梯内部信号衰减最为明显，因为电梯位于主楼中央天井处，由金属材料封闭四周，穿透损耗值最高，引起严重的信号衰减。在电梯内部单频点测量时，只能听到严重的噪声，且随楼层高度变化，信号强度没有明显变化，测量过程中若电梯开门，信号强度明显增强。2.高低频信号强度对比： 结合同班另一测量小组（组员：王文雪，张玉洁），所测量的教二楼内部楼道频道2（57.75MHz）的信号强度数据，进行高低频的对比： 信号电平(dbmw) 一层 二层 三层 四层 五层高频（535.25MHz） -75.3429 -73.3209 -68.0047 -63.7930 -60.3256低频（57.75MHz） -72.4524 -62.5128 -59.0394 -57.1768 -57.9962 可以看出，同样是高频的衰减大于低频的衰减，与理论相符。且基本上也是随着楼层的上升衰减下降，信号强度增强。4.主楼东南西北各方位电磁场强度及其分析信号电平（dBmw）北侧东侧南侧西侧频率（190.75MHz）-67.7714-67.9756-62.5490-64.7054图 各楼层各方位信号强度均值分布由雷达图直观看出各个楼层，各个方位信号强度分布。观测由中心到四个方位走向的信号分布，比较四个方位总体西侧和南侧信号强度较好，其中南侧最佳。不同楼层曲线比较得出，除去受大厅开门影响的一楼，其余各层随楼层高度增加，信号强度增强，符合阴影衰落对信号强度影响的规律。5.主楼总体a. 直方图及正太拟合最大值 （dBmw）最小值 （dBmw）平均值 （dBmw）标准差-47.5-75.4-66.615.84b. 概率累积分布图地点最大值（dBmw）最小值（dBmw）平均值（dBmw）标准差主楼室内总体-47.5-75.4d-66.615.84主楼外围总体-48.7-69-58.963.98二楼室内-66.5-75.4-72.111.81（2）穿透损耗从控制变量的思想出发，考虑主楼外围受周围建筑物阴影衰落的影响较大，而主楼高层可忽略阴影衰落影响，故不应以主楼室内总体信号强度与外围比较估算穿透损耗。与此同时一楼大厅开门处引起的信号强度分布不均，波动较大的影响，数据不具代表性。故取二层室内数据与主楼外围总体相比，估算钢筋混凝土材料的穿透损耗，反映出建筑物对信号的屏蔽和吸收作用。公式如下：P= 1 Ni=1NPi(outside)-1Mj=1MPj(inside)由主楼外围总体测量信号均值减去二层室内总体信号均值，可得穿透损耗：P=(-58.96)-( -72.11)=13.15dBmw分析说明：（1）阴影衰落分布规律随机对数正态分布 从实验结果中各区域实测数据的概率分布柱状图与拟合正态分布，及累积概率曲线和对应正态分布函数曲线(均值、方差相同)对比可得，无线信号强度分布符合随机对数正态分布。实际测量主楼外围信号强度，根据所测数据分析得出南侧的信号强度最高，因为主楼南侧的遮挡物较少、距离教二楼较远且建筑物较低，比较空旷，信号衰减小，信号强度高；主楼北侧的区域有教一楼遮挡，东侧有科学会堂遮挡，信号衰减影响较严重，信号强度较差。以上均符合阴影衰落原理（阴影衰落是电磁波在空间传播时受到地形起伏、高达建筑物群的阻挡，在这些障碍物后面会产生电磁场的阴影，造成场强中值的变化，从而引起信号衰减）。由雷达图直观看出各个楼层，各个方位信号强度分布。观测由中心到四个方位走向的信号分布，比较四个方位总体西侧和南侧信号强度较好，其中南侧最佳。与此同时我们也在本次实验中发现，除了周围建筑物，周围地方的物品停放，人流，天气等等都对当时的信号有影响。例如本次实验主楼外围北侧的大量自行车停放，西侧的大风等等。 （2）空间结构对信号强度的影响所选的主楼的空间结构比较复杂,但是可以从信号强度谱图中看出,楼道中靠近楼梯和窗口的地方信号强度都比较好,厕所附近信号有不规律的随机波动，而被建筑物包围的内侧信号相对于外侧信号而言有了明显的衰减,尤其是在拐角的阴影处,信号更差,而越是测量地点越是开阔的地方,信号强度则越强.与此同时也可比较得出不同建筑物对穿透损耗的影响，例如有铁门，引起信号强度陡然降低，反映出铁质建筑材料相比钢筋混凝土建筑材料对信号衰减影响更严重；相反的一楼西侧大厅玻璃门处测量信号强度明显优于同层其他地区，优于二层西侧地区信号强度；电梯内部信号衰减最为明显，因为电梯位于主楼中央天井处，由金属材料封闭四周，穿透损耗值最高，引起严重的信号衰减。以上均符合不同建筑材料对穿透损耗影响的相关规律。（3）楼层与信号强度关系由本次实验每层楼信号强度均值做比较以及雷达图上不同楼层曲线得出，除去受大厅开门影响的一楼，其余各层随楼层高度增加，信号强度增强，符合阴影衰落对信号强度影响的规律，在数据中显示到达6层及其以上楼层信号强度明显增加，分析其原因是高楼层对应的室外比较空旷，没有周围高层建筑物的阻挡，因此衰减很小。六问题分析与解决1. Matlab仿真问题研究与解决获取的数据需要进行matlab数据分析，绘出散点图、概率累积分布图和正太分布拟合图，并进行分析。主要解决的问题有：1.从excel表格里获取每个地方的数据列向量：直接从excel里导入的是矩阵，而需要分析的是每个地点的数据列向量，于是用从矩阵中获取列向量的方法（例如：data1:54,1）2.绘制信号强度散点图：用matlab的plot函数绘制散点图，需要生成一个和数据点数相同的行向量作为横坐标，进行绘制。3.绘制概率累积分布图：用matlab的cdfplot函数进行概率累积分布图的绘制4.绘制正太分布拟合图：统计不同区间内数据的分布直方图，再在根据平均值和标准差进行正太分布的拟合。（具体代码请见附录）2. 场强分布的研究根据matlab的数据分析和绘图，我们进行了教二楼地点的信号强度分布特性。通过对比东、西、南、北侧的信号强度特性，发现教二楼南北侧的信号强度要由于东西侧。南侧由于测量时距离教二楼较远，因此信号衰减较小，而北侧距离教二楼较近，因此受到阴影衰落比较严重，因此虽然较空旷但信号衰减略大。通过对比不同楼层的信号强度，发现随着楼层的上升，信号强度逐渐增强；且每个楼道的两侧信号强度较小，中心信号强度较大。通过对比高低频的信号的信号强度，发信低频信号的衰减明显要小于高频信号。这也与自由空间传播的损耗：空间损耗=20lg（F）+20lg（D）+32.4；F为频率，单位：MHz；D为距离，单位：Km；频率越高虽好越大相符。 3.模型拟合对教二楼周边进行模型拟合：取频率f=535.25MHZ，设发射天线在主楼，取发射天线到移动台水平距离为=150m，取天线高度（即主楼高度）=40m，取移动台高度为=1.5m,得到各模型下的路径损耗计算值为：自由空间模型：55.6930 布灵顿模型：51.4806市区Hata-Okumura模型：71.1620 开阔地Hata-Okumura模型：56.7704比较可知教二楼周边较符合市区Hata-Okumura模型。 1、2、七分工安排测量过程和数据整理：由王雅萱读数，杨淑铃记录matlab仿真：由王雅萱完成讨论分析结果：共同完成论文整理：由杨淑铃完成八心得体会杨淑铃： 经过这次的实验，我了解到了主楼外部和内部的场强分布，观察场强分布的变化可以学习到影响场强大小的各种原因以及其原理，像是阴影衰落和穿透损耗等相关的知识。于是我在移动通信课程中学到的理论知识通过本次实地测量更加深入理解。此次实验中得知在测量衰减的时候，无论是一排自行车或是一道门都会对场强衰减度产生影响，在高处空旷处得体场强相对就会稳定，高出的地方信号也会变好，衰减也会小一些。不禁想起那些生活中爬到楼顶去打电话的实际例子，原来都是有科学的依据和相应的数学模型。同时测量的过程也是给我留下了深刻的印象，从下午测到晚上天黑，在主楼人员稀疏的环境中两个女生也真的是为科学事业勇往直前呀，相信很多年后也是一笔谈资，一个值得怀念的实验过程。张诗悦： 这次实验我收获了很多： 首先我了解了分析电磁场强度的方法，通过实地测量电磁场强度，我们了解了学校教二楼周边的电磁场强度分布，以及我们教二楼每一层楼道的电磁场强度分布。 其次，我学会了使用matlab处理数据。散点图的绘制，概率累积分布图的绘制，以及直方图和正太分布的拟合。 最后，我们通过与经典模型的对比得出了教二楼周边符合的基本模型。 通过这次实验，我们第一次做了一个实地测量分析的实验，有趣有有意义，让我们亲身经历了实践当中电磁场分布的分析，收获了很多也学到了很多。 谢谢老师的指导和同学的帮助！附录I：原始数据1.教二楼周边数据教二楼西侧电平教二楼南侧电平教二楼东侧电平教二楼北侧电平-71.7-67.7-73.1-69.2-72-68.3-70.5-67.6-69.2-63.6-76.4-68.2-73.2-69.2-75.5-66.7-70.5-66.2-80.5-70.2-72-67.2-71.4-68.1-68.4-66.2-71.6-70.9-71.5-65.2-75.5-69.7-64.4-64.3-72.5-68.8-65.8-70.4-75.4-62.5-71.2-72.3-76.8-66.4-70.3-65.2-72.8-75.8-73-70.8-71.3-68.7-68.2-68.8-73.5-68.8-65-65.7-71.7-68-64.7-67.2-73.5-74.9-68-71.4-79.6-70.2-69.4-64.5-82.2-66.9-70.3-62.4-77.9-67.3-70.8-65-80.2-64.3-70.7-68.3-79-74.2-71.5-63.7-80.6-63.4-69.2-68.9-78.2-73.3-69.9-71.2-81.3-72.2-74.3-66-82.1-73.3-73.1-72-73.8-69.6-68.9-63.5-81.9-65.4-70.3-66-76-69.4-71.8-64.8-77.4-73.4-71.4-72.3-78.8-65.2-67.9-68.9-78.2-72-67.4-63.8-77.2-67.9-73.7-63.2-77.2-68.1-79.2-62.7-76.8-64.2-73.7-63-80.3-69.7-70.9-67.9-72.9-68.8-72.5-63.8-73.2-71.9-66.8-64.1-78.3-69-71.9-66.8-77.9-67.4-74.5-67.8-76.2-69.8-71.8-66.3-81-66.8-73-64.8-77.2-69.4-66.9-67.5-69-66.7-76.2-63.7-69.9-67.6-66.2-63.7-72.1-69-73.9-66.2-69.4-68.7-72.3-69-65.4-67.2-62.5-69.7-69.2-64-68.8-65.4-63.4-67.4-68.7-71.6-73-68.8-67.7-65.9-65.9-65.9-73.1-64.9-65.5-70.9-60.9-61.2-61.1-65.1-65-67.2-63.3-64.9-66.3-69.1-65-67.2-69.4-71.4-60.1-70.9-64.4-66.9-66.5-65.9-65.2-68.4-70.9-67.5-62.8-63.3-70.9-68.8-67.2-67.3-63.2-72.2-65.2-70.3-60.4-75.9-65-73.4-61.2-76.3-66.2-74.4-67.3-71.5-60.6-70-67.7-71.2-67.1-78.9-69.2-71.2-60.2-70.6-61.2-70.9-60.8-74.9-62.2-73.5-66.5-71.9-65.8-73.5-64.8-69-68.5-64.62.教二楼每个楼层楼道数据教二楼一层电平教二楼二层电平教二楼三层电平教二楼四层电平教二楼五层电平-82.4-77.8-70.5-66-65.8-77.9-73.3-73.1-72.5-66.6-80.4-78.1-78.5-66.5-67-81.3-77.4-77.9-68.7-63.7-78.3-81.2-74-73.9-59.8-74.7-74.7-72.1-71.7-54.8-80.4-79.9-73.8-66.5-48.6-77.2-75.8-72.2-72.9-65.9-75.9-71.4-65.1-65.8-58.2-78.5-74.2-72.9-63.4-53.2-75.7-75.5-63-70.4-47.6-77.1-72.6-66.6-63.4-54.8-74.7-75.4-68.4-60.1-59.9-71.4-77.6-64.7-61.2-56.8-76.5-72.2-69.9-64.7-54.9-74.7-70.7-63.9-58.1-63.7-70.8-68.6-66.2-65.5-68-72.4-68.8-59.8-57.8-66.9-67.9-70.4-60.1-64.1-63.4-67.7-72.6-58.4-62.4-59.4-66.6-70.6-57.7-64-63.3-71.4-67.3-59.9-63-55.4-66.4-61.5-57.3-54.8-62.7-70.3-63.7-60.5-57.2-58.7-71.4-66.3-64.5-64.5-57.9-71.2-66.7-63.2-60-54.6-74.9-74.3-67.3-60.5-60.4-70.3-69.6-65.4-62.6-60.1-76-71.8-70.1-60.7-58.9-69.3-68.6-67.5-63.8-54.6-66.8-70.1-66.4-54.6-52.1-71-70.8-68.7-55.8-56.7-76.4-69.3-69.3-55.9-48.4-73.9-70.3-69-58.1-57.6-74.8-76.6-70.3-68.1-63.4-74.2-78.6-71.2-67.3-64.5-80.2-78.9-73.4-59.2-59.8-81.6-74.1-70-60.2-60.4-72.4-80.5-69.3-69.8-65.4-75.9-79.1-68.9-70.9-68.4-82.1-84.7-74.2-71-66.7-76.9-74-73.4-60.5-67-76.8-77.2-75.6-65-68-79.6-80.4-81.5-84.6-79.3-79.7附录II：源代码%录入原始数据data=xlsread(west.xlsx,1);Data=data; %将数据转置为行向量 Data=-1.*Data;L=length(Data); %数据长度%作出原始数据散点图 figure(11)subplot(1,2,1);x=(1:1:L);plot(x,Data,Linewidth,1.3, Color,0.3,0.3,1);meandata=mean(Data);hold on;plot(x,meandata,r., Color,1,0.4,0.6);legend(信号电平,平均值(-62.5959);title(十三楼信号强度分布（CH8频道，f=190.758MHz）);xlabel(信号个数);ylabel(信号电平值（dBmw）);%作出样本频率分布直方图、概率分布曲线和拟合正太分布曲线subplot(1,2,2);y,x=hist(Data,10)%画样本数目直方图y=y/L/mean(diff(x);h1=bar(x,y,1);axis(-80,-45,0,0.20);grid onh = get(gca,Children); set(h,FaceColor,.85 .85 .85) %设置颜色min1=min(Data)%最小值max1=max(Data)%最大值mean1=mean(Data)%均值std1=std(Data)%标准差hold on x=-80:0.5:-45.4;y=normpdf(x,mean1,std1);%在相同均值标准差下，画正态分布图plot(x,y,r, Linewidth,1.8, Color,1,0.4,0.6)axis(-80,-45,0,0.20);hold onn,m=hist(Data);%在同一坐标系下，画出统计概率图plot(m,n/L/mean(diff(m), Linewidth,1.8, Color,0.3,0.3,1 )legend(实际样本分布 ,，相同的正态分布,样本概率分布)%显示图例grid ontext(-60,0.12,最大值= ,num2str(max1);%在图中标明最值text(-60,0.11,最小值= ,num2str(min1);text(-60,0.10,平均值= ,num2str(mean1,%.2fn);text(-60,0.09,标准差= ,num2str(std1,%.2fn);str=主楼总体;信号电平概率分布直方图与正态分布拟合曲线; title(str);xlabel(电平值(dBmw);ylabel(概率p(x);%作出原始数据累积概率分布曲线figure(2)cdfplot(Data)title(主楼外围西侧信号概率累积分布曲线（CH8频道，f=190.758MHz）);ylabel(信号概率);xlabel(信号电平值（dBmw）);%作出原始数据信号强度二维表面图Data2=Data,zeros(1,L),(1:L);%转换矩阵Data2=reshape(Data2,L,3);figure(12)surf(Data2);%画surf图，表示统计分布title(主楼十三楼走廊自西北角顺时针信号电平分布);xlabel(主楼十三楼走廊北至东侧-主楼十三楼走廊南至西侧);axis(1,L,1,2);caxis(-80,-50);colorbar(horiz);
展开阅读全文