北邮电磁场与电磁波实验报告.docx
《北邮电磁场与电磁波实验报告.docx》由会员分享，可在线阅读，更多相关《北邮电磁场与电磁波实验报告.docx（25页珍藏版）》请在冰豆网上搜索。
北邮电磁场与电磁波实验报告
信息与通信工程学院
电磁场与电磁波实验报告
题目：
校园信号场强特性的研究
姓名
班级
学号
序号
薛钦予
2011210496
2011210496
21
一、实验目的
1.掌握在移动环境下阴影衰落的概念以及正确的测试方法；
2.研究校园内各种不同环境下阴影衰落的分布规律；
3.掌握在室内环境下场强的正确测量方法，理解建筑物穿透损耗的概念；
4.通过实地测量，分析建筑物穿透损耗随频率的变化关系；
5.研究建筑物穿透损耗与建筑材料的关系。
二、实验原理
1、电磁波的传播方式
无线通信系统是由发射机、发射天线、无线信道、接收机、接收天线所组成。
对于接受者，只有处在发射信号的覆盖区内，才能保证接收机正常接受信号，此时，电波场强大于等于接收机的灵敏度。
因此基站的覆盖区的大小，是无线工程师所关心的。
决定覆盖区的大小的主要因素有：
发射功率，馈线及接头损耗，天线增益，天线架设高度，路径损耗，衰落，接收机高度，人体效应，接收机灵敏度，建筑物的穿透损耗，同播，同频干扰等。
电磁场在空间中的传输方式主要有反射﹑绕射﹑散射三种模式。
当电磁波传播遇到比波长大很多的物体时，发生反射。
当接收机和发射机之间无线路径被尖锐物体阻挡时发生绕射。
当电波传播空间中存在物理尺寸小于电波波长的物体﹑且这些物体的分布较密集时，产生散射。
散射波产生于粗糙表面，如小物体或其它不规则物体﹑树叶﹑街道﹑标志﹑灯柱。
2、尺度路径损耗
在移动通信系统中，路径损耗是影响通信质量的一个重要因素。
大尺度平均路径损耗：
用于测量发射机与接收机之间信号的平均衰落，即定义为有效发射功率和平均接受功率之间的（dB）差值，根据理论和测试的传播模型，无论室内或室外信道，平均接受信号功率随距离对数衰减，这种模型已被广泛的使用。
对任意的传播距离，大尺度平均路径损耗表示为：
（式1）
即平均接收功率为：
（式2）
其中，定义n为路径损耗指数，表明路径损耗随距离增长的速度，d0为近地参考距离，d为发射机与接收机之间的距离。
公式中的横杠表示给定值d的所有可能路径损耗的综合平均。
坐标为对数－对数时，平均路径损耗或平均接收功率可以表示为斜率10ndB/10倍程的直线。
n依赖于特定的传播环境，例如在自由空间，n为2；当有阻挡物时，n比2大。
决定路径损耗大小的首要因素是距离，此外，它与接受点的电波传播条件密切相关。
为此，我们引进路径损耗中值的概念，中值是使实验数据中一半大于它而另一半小于它的一个数值（对于正态分布中值就是均值）。
人们根据不同放入地形地貌条件，归纳总结出各种电波传播模型。
下边介绍几种常用的描述大尺度衰落的模型。
常用的电波传播模型：
1）自由空间模型
2）布灵顿模型
3）EgLi模型
4）Hata-Okumura模型
3、阴影衰落
在无线信道里，造成慢衰落的最主要原因是建筑物或其它物体对电波的遮挡。
在测量过程中，不同位置遇到的建筑物遮挡情况不同，因此接收功率也不同，这样就会观察到衰落现象。
由于这种原因造成的衰落也叫“阴影效应”或“阴影衰落”。
在阴影衰落的情况下，移动台被建筑物所遮挡，它收到的信号是各种绕射反射，散射波的合成。
所以，在距基站距离相同的地方，由于阴影效应的不同，它们收到的信号功率有可能相差很大，理论和测试表明，对任意的d值，特定位置的接受功率为随机对数正态分布即：
（式3）
其中，Xσ为0均值的高斯分布随机变量，单位dB；标准偏差σ，单位dB。
对数正态分布描述了在传播路径上，具有相同T-R距离时，不同的随机阴影效应。
这样利用高斯分布可以方便地分析阴影的随机效应。
正态分布，也叫高斯分布，概率密度函数为：
（式4）
应用于阴影衰落时，上式中的
表示某一次测量得到的接收功率，
表示以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
表1.阴影衰落分布的标准差
（dB）
4、建筑物的穿透损耗的定义
建筑物穿透损耗的大小对于研究室内无线信道具有重要意义。
穿透损耗又称大楼效应，一般指建筑物一楼内的中值电场强度和室外附近街道上中值电场强度dB之差。
发射机位于室外，接收机位于室内，电波从室外进入到室内，产生建筑物的穿透损耗，由于建筑物存在屏蔽和吸收作用，室内场强一定小于室外的场强，造成传输损耗。
室外至室内建筑物的穿透损耗定义为：
室外测量的信号平均场强减去同一位置室内测量的信号平均场强。
用公式表示为：
（式5）
是穿透损耗，单位是dB；
是在室内所测的每一点的功率，单位是
，共
个点；
是在室外所测的每一点的功率，单位是
，共
个点。
三、实验内容
利用DS1131场强仪，实地测量信号场强。
1）研究具体现实环境下阴影衰落分布规律，以及具体的分布参数如何。
2）研究在校园内电波传播规律与现有模型的吻合程度，测试值与模型预测值的预测误差如何。
3）研究建筑物穿透损耗的变化规律。
四、实验步骤
1、实验对象的选择
我们选择了北邮足球场，因为我们想看看，空旷的地方和人流、树木、建筑物、车辆等等密集的地方对信号强弱的影响，还能得到学校中心位置的信号强度分布是什么样的。
（测量路线会在报告的数据采集部分标注出）
对于选频，我们选了频道五，主要是因为这个频率接近100Mhz，经过查看，是90.0MHZ左右，经过计算，其波长入=C/F，接近3m，所以我测试了下，我的步子比较大，两步正好是1.5m，那么如果注意控制脚步的话，还是很好的。
因此，在刚开始的时候，我读取数据，再用表格记录，这样比较准确；读取数据和记录数据同步，因此不会有很大的误差，尽量把误差控制在随机误差的范围。
2、数据采集
仪器：
场强仪DS1131
采用的单位：
老师讲的dBmw
采集：
半个波长读取一次数据，记录一次；把每步的长度尽量的控制，不让出现很大的偏差，读数据的时候，尽量等数据基本稳定以后读取，基本都是取中间的数据，因此公差不会太大。
测
3、数据录入
将测量得到的数据录入Excel表格。
我独自完成，得到一个数据将近500的表格。
63.6
64.6
71.2
68.5
66.3
66.9
79.7
60.1
55.5
61.2
59.1
60
57.6
60.2
55.8
64.2
59.1
61.2
59.2
64.9
63.9
60.1
63.6
62.7
61.3
60.9
62.3
57.1
58.1
63.7
57.6
57.5
58.1
63.7
64.6
62.7
66.6
61
57.9
60.3
63.4
65.2
65.4
62.4
61.3
62.9
66.8
59.3
67
61.1
62.6
63.4
63.2
66.7
62.7
67.2
64.3
58、1
59、1
57.2
65.3
72.3
63.7
60.5
63.7
60.3
65.1
62.5
65.8
68.9
64.3
61.5
61.2
60.4
61.3
69.9
67.96
63.67
59.9
66.5
65.9
71.2
6.6
65.4
72.1
65.1
64.3
62
56.2
59
65.2
59
59.2
58.8
58.1
61.1
60.1
65.7
59.4
57.5
61.5
56.2
64.9
60。
0
58.2
59
53.1
49.1
52.1
50.1
55.1
58.8
57.1
62.7
57.3
53.3
53.9
54.1
50.1
57.4
64.5
58.9
58.。
9
58.9
57.9
63.9
59.3
57.。
3
60.2
69.2
65.9
59.9
55.2
58.8
61.4
69.2
63.1
58.7
60.2
65.4
65.4
59.1
60.2
71.6
右侧北师
61.2
64.7
59.2
59.6
65
69.2
64.2
69.2
57.2
6.2
62.1
58.2
72.3
58.2
59.2
60.4
58.2
65.2
65.3
64.4
64.3
68.2
69.3
69.4
70.2
62
62.8
69.3
65.4
65.3
65.3
66.2
66.3
67.3
70.2
65.3
72.3
63.7
60.5
63.7
60.3
65.1
62.5
65.8
68.9
64.3
61.5
61.2
60.4
61.3
69.9
67.1
63.67
59.9
66.5
65.9
71.2
6.6
65.4
72.1
65.1
64.3
62
56.2
59
65.2
59
59.2
58.8
58.1
61.1
60.1
65.7
59.4
57.5
61.5
56.2
64.2
60.1
58.2
67.1
53.1
49.1
52.1
50.1
55.1
58.8
57.1
62.7
57.3
53.3
53.9
54.1
50.1
57.4
64.5
58.9
57.1
58.9
57.9
63.9
59.3
52.3
60.2
69.2
56.2
57.1
64.9
60
58.2
59
53.3
49.4
52.1
50.1
55.1
58.8
57.1
62.7
57.3
53.5
53.9
54.1
50.1
53.1
57.4
64.4
58.9
58.9
58.2
58
57.1
64.9
59.2
57.3\
60.2
69.2
65.2
59.9
55.2
58.8
61.5
69.2
53.3
59.7
60.265.4
59
60.2
71.3
71.6
700
65.3
69.2
59.2
64.7
66.2
62.2
62.4
66.3
59.4
56.4
54.9
59.9
59.9
62.4
65.5
62.5
6.3
59.9
564
54.9
59.9
59.9
62.4
65.9
65.9
62.5
63.3
59.3
63.1
63.2
59.3
64.2
63.2
59.5
62.9
60.1
53.3
54.2
58.2
53.2
61.2
54.3
55.4
55.4
61.2
60
28.3
28.2
29.2
59.2
60.2
53.5
54.2
58.2
53.2
61.2
54
55.4
55.4
61.2
60.6
58.2
58.8
59.2
60.2
65.3
62.8
60.9
68.3
53.2
61.2
54.3
55.4
61.2
60
58.2
58.8
59.2
60.2
60.2
65.3
62.8
60.9
67.6
68.1
68.2
68
64.9
64.2
6.2
58.7
67.2
65.3
60.2
70.2
63.3
62.1
61
67.3
67
67.6
68.1
69
64
69.3
69.2
66.3
61.3
59.3
63.1
63.5
65.2
57.8
55.8、
55.8
55.3
57.3
53.8
51.8
51.1
52.9
59.2
61.8
65.8
55.2
55.8
58.2
61.3
62.3
60.3
56
56.3
58.2
58.4
57.2
61.2
56
56.3
58.3
58.4
57.3
61.2
54.2
59.2
4、数据处理流程
数据。
数据处理的流程图。
5、实验结果与分析
h1=
229.0038
s1=
min:
39.3000
max:
61.7000
mean:
47.1760
median:
45.4500
std:
5.4407
h1=
671.0024
s1=
min:
37.5000
max:
58
mean:
45.8845
median:
45.5000
std:
4.5477
h1=
1.1130e+003
s1=
min:
30.4000
max:
58
mean:
43.6837
median:
44
std:
6.2097
h1=
1.5550e+003
s1=
min:
35.6000
max:
62.1000
mean:
49.8205
median:
49.5000
std:
6.4593
h1=
1.9970e+003
s1=
min:
43.8000
max:
66.4000
mean:
53.2384
median:
52.8000
std:
3.9029
westV=
47.1765.4407
southV=
45.88454.5477
eastV=
43.68376.2097
northV=
49.82056.4593
代码：
clearall;
closeall;
%-------------读取文件---------------%
w2e=xlsread（'data1.xlsx','sheet2'）;
s2n=xlsread（'data1.xlsx','sheet1'）;
e2w=xlsread（'data1.xlsx','sheet4'）;
n2s=xlsread（'data1.xlsx','sheet5'）;
e2wout=xlsread（'data1.xlsx','sheet3'）;
%-------------转换成矩阵------------%
w2e2=reshape（w2e,1,50）;
s2n2=reshape（s2n,1,58）;
e2w2=reshape（e2w,1,147）;
n2s2=reshape（n2s,1,73）;
e2wout2=reshape（e2wout,1,159）;
%---------为画平面场强图作准备----------%
w2e3=[w2e2,zeros（1,50）,[1:
50]];
w2e3=reshape（w2e3,50,3）;
s2n3=[s2n2,zeros（1,58）,[1:
58]];
s2n3=reshape（s2n3,58,3）;
e2w3=[e2w2,zeros（1,147）,[1:
147]];
e2w3=reshape（e2w3,147,3）;
n2s3=[n2s2,zeros（1,73）,[1:
73]];
n2s3=reshape（n2s3,73,3）;
e2wout3=[e2wout2,zeros（1,159）,[1:
159]];
e2wout3=reshape（e2wout3,159,3）;
%-----------科学会堂正面由西向东-----------%
figure（11）
subplot（1,2,1）;
histfit（w2e2）;%画柱状图
axis（[30,70,0,30]）;
gridon;
str={'科学会堂正面由西向东';'信号电平概率分布'};
title（str）;
xlabel（'电平值（-dBmw）'）;
ylabel（'样本数量（个）'）;
legend（'实际样本分布','理想概率分布线'）;
subplot（1,2,2）;
[h1,s1]=cdfplot（w2e2）%画累积概率分布图
axis（[30,70,0,1]）;
holdon;
w2emean=num2str（s1.mean）;
w2estd=num2str（s1.std）;
text（56,0.23,['最小值=',num2str（s1.min）]）;
text（56,0.18,['最大值=',num2str（s1.max）]）;
text（56,0.13,['均值=',num2str（s1.mean）]）;
text（56,0.08,['中值=',num2str（s1.median）]）;
text（56,0.03,['标准差=',num2str（s1.std）]）;
title（'对应累积概率分布'）;
figure（12）
surf（w2e3'）;%画衰落强度图
title（'科学会堂正面由西向东信号电平分布图'）;
xlabel（'科学会堂正面由西向东'）;
axis（[1,50,1,2]）;
caxis（[3070]）;
colorbar（'horiz'）;
%------------科学会堂正面由东向西-----------%
figure（21）
subplot（1,2,1）;
histfit（s2n2）;
axis（[30,70,0,30]）;
gridon;
title（'科学会堂正面由东向西信号电平概率分布'）;
xlabel（'电平值（-dBmw）'）;
ylabel（'样本数量（个）'）;
legend（'实际样本分布','理想概率分布线'）;
subplot（1,2,2）;
[h1,s1]=cdfplot（s2n2）
axis（[30,70,0,1]）;
holdon;
s2nmean=num2str（s1.mean）;
s2nstd=num2str（s1.std）;
text（56,0.23,['最小值=',num2str（s1.min）]）;
text（56,0.18,['最大值=',num2str（s1.max）]）;
text（56,0.13,['均值=',num2str（s1.mean）]）;
text（56,0.08,['中值=',num2str（s1.median）]）;
text（56,0.03,['标准差=',num2str（s1.std）]）;
title（'对应累积概率分布'）;
figure（22）
surf（s2n3'）;
title（'科学会堂正面由东向西信号电平分布图'）;
xlabel（'科学会堂正面由东向西'）;
axis（[1,58,1,2]）;
caxis（[3070]）;
colorbar（'horiz'）;
%-------------科学会堂西侧面由北向南-----------%
figure（31）
subplot（1,2,1）;
histfit（e2w2）;
axis（[30,70,0,100]）;
gridon;
str={'科学会堂西侧面由北向南';'信号电平概率分布'};
title（str）;
xlabel（'电平值（-dBmw）'）;
ylabel（'样本数量（个）'）;
legend（'实际样本分布','理想概率分布线'）;
subplot（1,2,2）;
[h1,s1]=cdfplot（e2w2）
axis（[30,70,0,1]）;
holdon;
e2wmean=num2str（s1.mean）;
e2wstd=num2str（s1.std）;
text（56,0.23,['最小值=',num2str（s1.min）]）;
text（56,0.18,['最大值=',num2str（s1.max）]）;
text（56,0.13,['均值=',num2str（s1.mean）]）;
text（56,0.08,['中值=',num2str（s1.median）]）;
text（56,0.03,['标准差=',num2str（s1.std）]）;
title（'对应累积概率分布'）;
figure（32）
surf（e2w3'）;
title（'科学会堂西侧面由北向南信号电平分布图'）;
xlabel（'科学会堂西侧面由北向南'）;
axis（[1,147,1,2]）;
caxis（[3070]）;
colorbar（'horiz'）;
%-----------足球场后面由西向东-----------%
figure（41）
subplot（1,2,1）;
histfit（n2s2）;
axis（[30,70,0,40]）;
gridon;
title（'足球场后面由西向东信号电平概率分布'）;
xlabel（'电平值（-dBmw）'）;
ylabel（'样本数量（个）'）;
legend（'实际样本分布','理想概率分布线'）;
subplot（1,2,2）;
[h1,s1]=cdfplot（n2s2）
axis（[30,70,0,1]）;
holdon;
n2smean=num2str（s1.mean）;
n2sstd=num2str（s1.std）;
text（56,0.23,['最小值=',num2str（s1.min）]）;
text（56,0.18,['最大值=',num2str（s1.max）]）;
text（56,0.13,['均值=',num2str（s1.mean）]）;
text（56,0.08,['中值=',num2str（s1.median）]）;
text（56,0.03,['标准差=',num2str（s1.std）]）;
title（'对应累积概率分布'）;
figure（42）
surf（n2s3'）;
title（'足球场后面由西向东信号电平分布图'）;
xlabel（'足球场后面由西向东'）;
axis（[1,70,1,2]）;
caxis（[3070]）;
colorbar（'horiz'）;
%-------足球场东侧面由南向北--------%
figure（61）
subplot（1,2,1）;
histfit（e2wout2）;
axis（[30,70
升级会员 