北邮无线场强实验报告.docx
《北邮无线场强实验报告.docx》由会员分享，可在线阅读，更多相关《北邮无线场强实验报告.docx（35页珍藏版）》请在冰豆网上搜索。
北邮无线场强实验报告
电磁场与电磁波实验报告
校园无线信号场强特性的研究
班级：
2011211116
学号：
2011210466，2011210467
班内序号：
21，22
姓名：
李亚东，钱晋
一、实验目的
1.掌握在移动环境下阴影衰落的概念以及正确的测试方法；
2.研究校园内各种不同环境下阴影衰落的分布规律；
3.掌握在室内环境下场强的正确测量方法，理解建筑物穿透损耗的概念；
4.通过实地测量，分析建筑物穿透损耗随频率的变化关系；
5.研究建筑物穿透损耗与建筑材料的关系。
二、实验原理
1.电磁波的传播方式
无线通信系统是由发射机、发射天线、无线信道、接收机、接收天线所组成。
对于接受者，只有处在发射信号的覆盖区内，才能保证接收机正常接受信号，此时，电波场强大于等于接收机的灵敏度。
因此基站的覆盖区的大小，是无线工程师所关心的。
决定覆盖区的大小的主要因素有：
发射功率，馈线及接头损耗，天线增益，天线架设高度，路径损耗，衰落，接收机高度，人体效应，接收机灵敏度，建筑物的穿透损耗，同播，同频干扰等。
电磁场在空间中的传输方式主要有反射﹑绕射﹑散射三种模式。
当电磁波传播遇到比波长大很多的物体时，发生反射。
当接收机和发射机之间无线路径被尖锐物体阻挡时发生绕射。
当电波传播空间中存在物理尺寸小于电波波长的物体﹑且这些物体的分布较密集时，产生散射。
散射波产生于粗糙表面，如小物体或其它不规则物体﹑树叶﹑街道﹑标志﹑灯柱。
2.尺度路径损耗
在移动通信系统中，路径损耗是影响通信质量的一个重要因素。
大尺度平均路径损耗：
用于测量发射机与接收机之间信号的平均衰落，即定义为有效发射功率和平均接受功率之间的（dB）差值，根据理论和测试的传播模型，无论室内或室外信道，平均接受信号功率随距离对数衰减，这种模型已被广泛的使用。
对任意的传播距离，大尺度平均路径损耗表示为：
即平均接收功率为：
其中，定义n为路径损耗指数，表明路径损耗随距离增长的速度，d0为近地参考距离，d为发射机与接收机之间的距离。
公式中的横杠表示给定值d的所有可能路径损耗的综合平均。
坐标为对数－对数时，平均路径损耗或平均接收功率可以表示为斜率10ndB/10倍程的直线。
n依赖于特定的传播环境，例如在自由空间，n为2；当有阻挡物时，n比2大。
决定路径损耗大小的首要因素是距离，此外，它与接受点的电波传播条件密切相关。
为此，我们引进路径损耗中值的概念，中值是使实验数据中一半大于它而另一半小于它的一个数值（对于正态分布中值就是均值）。
人们根据不同放入地形地貌条件，归纳总结出各种电波传播模型。
下边介绍几种常用的描述大尺度衰落的模型。
常用的电波传播模型：
1）自由空间模型
自由空间模型假定发射天线和接收台都处在自由空间。
我们所说的自由空间一是指真空，二是指发射天线与接收台之间不存在任何可能影响电波传播的物体，电波是以直射线的方式到达移动台的。
自由空间模型计算路径损耗的公式是：
其中Lp是以dB为单位的路径损耗，d是以公里为单位的移动台与基站之间的距离，f是以MHz为单位的移动工作频点或工作频段的频率。
空气的特性可近似为真空，因此当发射天线与移动台距离地面都较高时，可以近似使用自由空间模型来估计路径损耗。
2）布灵顿模型
布灵顿模型假设发射天线和移动台之间的地面是理想平面大地，并且两者之间的距离d远大于发射天线的高度ht，或移动台的高度hr，此时的路径损耗计算公式为：
其中距离d的单位是公里，天线高度ht及hr的单位是米，路径损耗Lp的单位是dB。
3）EgLi模型
前述的自由空间模型及布灵顿模型都是基于理论分析得出的计算公式。
EgLi公式则是从大量实测结果中归纳出来的中值预测公式，属于经验模型，其计算式为：
其中路径损耗Lp的单位是dB，距离d的单位是公里，天线高度ht及hr的单位是米，工作频率f的单位是MHz，地形修正因子G的单位是dB。
G反应了地形因素对路径损耗的影响。
EgLi模型认为路径损耗同接收点的地形起伏程度Δh有关，地形起伏越大，则路径损耗也越大。
当Δh用米来测量时，可按下式近似的估计地形的影响：
若将移动台的经典高度值hr=1.5m代入EgLi模型则有：
4）Hata-Okumura模型
Hata-Okumura模型也是根据实测数据建立的模型。
当移动台的高度为典型值hr=1.5m时，按Hata-Okumura模型计算路径损耗的公式为：
市区内的Hata模型为：
简化后为：
5）CCIR模型
CCIR给出了反映自由空间路径损耗和地形引入的路径损耗联合效果的经验公式。
具体公式请参见模型计算部分。
3.阴影衰落
在无线信道里，造成慢衰落的最主要原因是建筑物或其它物体对电波的遮挡。
在测量过程中，不同位置遇到的建筑物遮挡情况不同，因此接收功率也不同，这样就会观察到衰落现象。
由于这种原因造成的衰落也叫“阴影效应”或“阴影衰落”。
在阴影衰落的情况下，移动台被建筑物所遮挡，它收到的信号是各种绕射反射，散射波的合成。
所以，在距基站距离相同的地方，由于阴影效应的不同，它们收到的信号功率有可能相差很大，理论和测试表明，对任意的d值，特定位置的接受功率为随机对数正态分布即：
其中，Xσ为0均值的高斯分布随机变量，单位dB；标准偏差σ，单位dB。
对数正态分布描述了在传播路径上，具有相同T-R距离时，不同的随机阴影效应。
这样利用高斯分布可以方便地分析阴影的随机效应。
正态分布，也叫高斯分布，概率密度函数为：
应用于阴影衰落时，上式中的
表示某一次测量得到的接收功率，
表示以dB表示的接收功率的均值或中值，表示接收功率的标准差，单位是dB。
阴影衰落的标准差同地形，建筑物类型，建筑物密度等有关，在市区的150MHz频段其典型值是5dB。
除了阴影效应外，大气变化也会导致阴影衰落。
比如一天中的白天，夜晚，一年中的春夏秋冬，天晴时，下雨时，即使在同一个地点上，也会观察到路径损耗的变化。
但在测量的无线信道中，大气变化造成的影响要比阴影效应小的多。
下面是阴影衰落分布的标准差，其中
（dB）是阴影效应的标准差。
（dB）
频率（MHz）
准平坦地形
不规则地形（米）
城市
郊区
50
150
300
150
3.5~5.5
4~7
9
11
13
450
6
7.5
11
15
18
900
6.5
8
14
18
21
表2.1阴影衰落分布的标准差
（dB）
4.建筑物的穿透损耗的定义
建筑物穿透损耗的大小对于研究室内无线信道具有重要意义。
穿透损耗又称大楼效应，一般指建筑物一楼内的中值电场强度和室外附近街道上中值电场强度dB之差。
发射机位于室外，接收机位于室内，电波从室外进入到室内，产生建筑物的穿透损耗，由于建筑物存在屏蔽和吸收作用，室内场强一定小于室外的场强，造成传输损耗。
室外至室内建筑物的穿透损耗定义为：
室外测量的信号平均场强减去同一位置室内测量的信号平均场强。
用公式表示为：
是穿透损耗，单位是dB；
是在室内所测的每一点的功率，单位是
，共
个点；
是在室外所测的每一点的功率，单位是
，共
个点。
三、实验内容
利用DS1131场强仪，实地测量信号场强。
1）研究具体现实环境下阴影衰落分布规律，以及具体的分布参数如何。
2）研究在校园内电波传播规律与现有模型的吻合程度，测试值与模型预测值的预测误差如何。
3）研究建筑物穿透损耗的变化规律。
四、实验步骤
1、实验对象的选择
由于对衰落的现象比较感兴趣，我们选择了室外作为实验场所，我们围绕教三一周进行测量。
并为了研究穿透损耗现象，我们测量了教三一层内部的场强。
最后，为了研究开阔地段的信号衰落，我们还从教三北侧主楼西侧从毛主席像测量到西门。
在选频方面，我们决定采用频道15，频率487.25MHz，此时的波长约为0.62m，我们大约1步读取一个数据。
2、实验路线如红线所示
图1
3、数据采集及录入
我们将测量得到的数据录入Excel表格，得到一个由6张分表构成的总表：
教三西侧、教三北侧、教三东侧、教三南侧、教三北侧主楼西侧从毛主席像到西门、教三一层内部。
每个地点的数据占据一个sheet表格的一列，共6个sheet，450个数据。
如下图所示。
图2，数据表格截图
4、数据处理流程
采集到的数据很多，需要对数据进行细致的处理以便得到明确的结论。
下图所示为数据处理的流程图。
图3
5、matlab程序代码
clearall;
closeall;
north=-1.*[62.1,49.2,57.1,55.5,60.0,69.1,55.2,65.1,65.6,58.5,65.2,67.7,55.2,66.1,64.3,74.1,58.3,70.0,63.1,66.1,58.1,68.1,61.2,70.0,64.0,60.1,58.0,61.2,66.3,57.3,62.3,61.1,59.0,72.2,67.2,65.3,63.1,63.1,65.2,72.7,63.8,72.0,70.4,59.2,65.2,58.3,71.2,57.3,60.0,69.3,60.2,69.4,75.3,62.3,65.0,67.0,59.2,72.1,64.3,60.3,60.0,65.4,69.2,62.2,63.2,62.2,61.0,64.0,60.0,65.0,66.3,60.2,64.0,74.2,63.2,59.5,59.6,66.4,60.0,67.3,65.4,65.4,70.0,60.0,62.3,64.3,70.4,59.5,69.3,64.0,70.6,69.5,78.6,61.3,70.5,60.0,67.4,70.4,71.0,60.4,63.3,72.3,71.0,72.2,66.4,75.4,68.3,64.3,70.4,58.0,58.4,57.3,70.4,60.4,61.2,65.3,68.4];
west=-1.*[67.1,73.2,56.2,62.3,64.2,62.4,71.2,55.2,57.3,62.3,55.2,70.4,59.2,64.2,64.0,60.2,58.2,68.1,64.3,64.2,60.0,58.2,65.3,61.2,60.3,69.1,61.0,61.2,63.2,70.1];
south=-1.*[64.2,70.2,65.1,53.2,67.2,64.2,57.2,60.0,59.2,65.3,62.1,59.1,66.0,81.2,62.2,58.2,74.2,64.1,59.2,60.0,65.3,65.7,78.3,76.8,59.2,68.2,59.3,72.1,59.2,59.0,63.2,66.8,63.2,74.8,66.0,65.3,70.0,61.1,58.3,61.2,57.3,54.2,55.5,54.3,66.2,67.3,58.0,55.1,65.2,68.7,71.2,59.2,63.3,73.2,67.1,58.7,66.1,77.2,79.1,66.2,75.0,68.0,60.1,60.1,64.2,71.2,70.2,68.2,67.2,74.5,57.2,69.2,65.1,69.2,73.3,64.5,61.7,62.0,62.0,64.2,70.0,69.4,57.3,72.3,60.2,69.5,63.2,52.3,61.6,64.3,60.9,59.5,66.6,67.3];
east=-1.*[65.2,63.4,59.1,55.0,63.0,66.1,70.4,67.2,72.3,58.0,61.1,60.7,56.3,65.0,73.2,60.2,61.2,65.3,59.2,65.2,62.0,59.4,54.0,58.1,57.1,56.2,58.2,56.4,54.0,57.0,61.2,56.3,69.2,61.1,53.2,63.2,54.2,54.1,50.8,56.2];
zhu=-1.*[62.9,77.5,66.7,71.9,76.0,60.7,62.7,64.0,64.1,61.9,73.1,56.6,69.6,60.2,59.5,60.5,65.1,64.5,68.4,60.7,65.0,62.0,58.4,61.2,58.0,56.7,71.8,69.0,63.8,71.2,56.4,72.6,70.2,63.2,62.7,62.7,60.7,63.4,64.5,62.8,76.4,65.6,59.8,66.1,53.2,56.0,56.5,65.9,53.5,60.2,72.5,56.1,58.4,53.4,57.6,50.5,56.7,60.0,53.2,62.6,59.3,75.9,60.0,54.7,50.3,57.4,59.5,54.3,67.0,55.8,57.4,58.4,68.0,60.0,58.9,61.1,52.9,54.2,59.2,55.4,62.9,59.8,61.0,64.2,66.1,60.8,59.0,66.5,68.4,60.2,63.1,64.9,72.9,66.4];
in=-1.*[69.0,76.5,79.2,75,72,70,68,68,71,74,75,59,69,74,65,77,69,71,65,78,70,66,78,84,74,63,58,61,77,80,74,72,62,58,70,69,71,72,86,70,59,70,68,60,62,72,66,74,67,68,73,80,73,73,72,64,78,70,63,75,70,62,76,78,68,68,73,59,68,58,63,58,64,72,77];
%-------?
三北侧------%
figure（11）
subplot（1,2,1）;%所画图形的行数，列数，序号
hist（north）%?
?
?
?
?
?
?
?
axis（[-80,-50,0,30]）;%设定坐标轴的范围
gridon%使绘制的图有网格，便于观察横纵坐标
str={'?
三北侧';'?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;%标记横坐标
ylabel（'?
?
?
?
（?
）'）;
legend（'?
?
?
?
?
?
'）;%添加图解注释
h=get（gca,'Children'）;%?
?
?
?
set（h,'FaceColor',[.8.81]）
min1=min（north）%?
?
?
max1=max（north）%?
?
?
mean1=mean（north）%?
?
std1=std（north）%?
?
?
subplot（1,2,2）
x=-80:
0.5:
-50;
y=normpdf（x,mean1,std1）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（x,y,'r'）%r代表曲线为红色
axis（[-80,-50,0,0.7]）;
holdon%此后绘制的图形都将添加在这个图形的基础上
[n,m]=hist（north）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（m,n/46）
legend（'µ,s?
?
?
?
?
?
?
','?
?
?
?
?
?
'）%?
?
?
?
gridon
text（-60,0.15,['?
?
?
=',num2str（max1）]）;%?
?
?
?
?
?
?
text（-60,0.2,['?
?
?
=',num2str（min1）]）;
text（-60,0.25,['?
?
?
=',num2str（mean1,'%.2f\n'）]）;
text（-60,0.3,['?
?
?
=',num2str（std1,'%.2f\n'）]）;
str={'教三北侧';'?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;
ylabel（'?
?
p（x）'）;
tb1etow2=[north,zeros（1,117）,（1:
117）];%?
?
?
?
，生成1*46全零矩阵
tb1etow2=reshape（tb1etow2,117,3）;
figure（12）
surf（tb1etow2'）;%?
surf?
?
?
?
?
?
?
title（'?
三北侧'）;
xlabel（'?
三北侧'）;
axis（[1,117,1,2]）;
caxis（[-80,-50]）;
colorbar（'horiz'）;%格式水平
%-------?
三西侧------%
figure（11）
subplot（1,2,1）;%所画图形的行数，列数，序号
hist（west）%?
?
?
?
?
?
?
?
axis（[-75,-50,0,10]）;%设定坐标轴的范围
gridon%使绘制的图有网格，便于观察横纵坐标
str={'?
三西侧';'?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;%标记横坐标
ylabel（'?
?
?
?
（?
）'）;
legend（'?
?
?
?
?
?
'）;%添加图解注释
h=get（gca,'Children'）;%?
?
?
?
set（h,'FaceColor',[.8.81]）
min1=min（west）%?
?
?
max1=max（west）%?
?
?
mean1=mean（west）%?
?
std1=std（west）%?
?
?
subplot（1,2,2）
x=-75:
0.5:
-50;
y=normpdf（x,mean1,std1）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（x,y,'r'）%r代表曲线为红色
axis（[-75,-50,0,0.15]）;
holdon%此后绘制的图形都将添加在这个图形的基础上
[n,m]=hist（west）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（m,n/46）
legend（'µ,s?
?
?
?
?
?
?
','?
?
?
?
?
?
'）%?
?
?
?
gridon
text（-60,0.15,['?
?
?
=',num2str（max1）]）;%?
?
?
?
?
?
?
text（-60,0.2,['?
?
?
=',num2str（min1）]）;
text（-60,0.25,['?
?
?
=',num2str（mean1,'%.2f\n'）]）;
text（-60,0.3,['?
?
?
=',num2str（std1,'%.2f\n'）]）;
str={'教三西侧';'?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;
ylabel（'?
?
p（x）'）;
tb1etow2=[west,zeros（1,30）,（1:
30）];%?
?
?
?
，生成1*46全零矩阵
tb1etow2=reshape（tb1etow2,31,3）;
figure（12）
surf（tb1etow2'）;%?
surf?
?
?
?
?
?
?
title（'?
三西侧'）;
xlabel（'?
三西侧'）;
axis（[1,30,1,2]）;
caxis（[-75,-50]）;
colorbar（'horiz'）;%格式水平
%-------?
三南侧------%
figure（11）
subplot（1,2,1）;%所画图形的行数，列数，序号
hist（south）%?
?
?
?
?
?
?
?
axis（[-85,-50,0,25]）;%设定坐标轴的范围
gridon%使绘制的图有网格，便于观察横纵坐标
str={'?
三南侧';'?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;%标记横坐标
ylabel（'?
?
?
?
（?
）'）;
legend（'?
?
?
?
?
?
'）;%添加图解注释
h=get（gca,'Children'）;%?
?
?
?
set（h,'FaceColor',[.8.81]）
min1=min（south）%?
?
?
max1=max（south）%?
?
?
mean1=mean（south）%?
?
std1=std（south）%?
?
?
subplot（1,2,2）
x=-85:
0.5:
-50;
y=normpdf（x,mean1,std1）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（x,y,'r'）%r代表曲线为红色
axis（[-85,-50,0,0.5]）;
holdon%此后绘制的图形都将添加在这个图形的基础上
[n,m]=hist（south）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（m,n/46）
legend（'µ,s?
?
?
?
?
?
?
','?
?
?
?
?
?
'）%?
?
?
?
gridon
text（-60,0.15,['?
?
?
=',num2str（max1）]）;%?
?
?
?
?
?
?
text（-60,0.2,['?
?
?
=',num2str（min1）]）;
text（-60,0.25,['?
?
?
=',num2str（mean1,'%.2f\n'）]）;
text（-60,0.3,['?
?
?
=',num2str（std1,'%.2f\n'）]）;
str={'教三南侧';'?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;
ylabel（'?
?
p（x）'）;
tb1etow2=[south,zeros（1,94）,（1:
94）];%?
?
?
?
，生成1*46全零矩阵
tb1etow2=reshape（tb1etow2,94,3）;
figure（12）
surf（tb1etow2'）;%?
surf?
?
?
?
?
?
?
title（'?
三南侧'）;
xlabel（'?
三南侧'）;
axis（[1,94,1,2]）;
caxis（[-85,-50]）;
colorbar（'horiz'）;%格式水平
%-------?
三东侧------%
figure（11）
subplot（1,2,1）;%所画图形的行数，列数，序号
hist（east）%?
?
?
?
?
?
?
?
axis（[-75,-50,0,10]）;%设定坐标轴的范围
gridon%使绘制的图有网格，便于观察横纵坐标
str={'?
三东侧';'?
?
?
?
?
?
?
?
?
?
?
'};
title（str）;
xlabel（'?
?
?
（dBmw）'）;%标记横坐标
ylabel（'?
?
?
?
（?
）'）;
legend（'?
?
?
?
?
?
'）;%添加图解注释
h=get（gca,'Children'）;%?
?
?
?
set（h,'FaceColor',[.8.81]）
min1=min（east）%?
?
?
max1=max（east）%?
?
?
mean1=mean（east）%?
?
std1=std（east）%?
?
?
subplot（1,2,2）
x=-75:
0.5:
-50;
y=normpdf（x,mean1,std1）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（x,y,'r'）%r代表曲线为红色
axis（[-75,-50,0,0.2]）;
holdon%此后绘制的图形都将添加在这个图形的基础上
[n,m]=hist（east）;%?
?
?
?
?
?
?
?
?
?
?
?
?
?
plot（m,n/46
升级会员 