实验三_Okumura-Hata方法计算机仿真
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
姓名：123学号：321杭州电子科技大学
实验一Okumura-Hata方法计算机仿真
【实验目的】
⏹加深对奥村模型的理解；
⏹能够使用Matlab利用Okumura-Hata方法计算基本传输损耗；
⏹比较奥村模型和Okumura-Hata方法获得的基本传输损耗的差异，分析
Okumura-Hata方法的误差。
【实验内容】
⏹使用Matlab利用Okumura-Hata方法计算基本传输损耗；
⏹分析Okumura-Hata方法的误差；
【实验设备】
⏹一台PC机
【实验步骤】
1.采用Okumura-Hata方法分别计算大城市市区地区准平滑地形、郊区和开阔区，基站天线高度
是200米，手机天线高度是3米情况下，不同传播距离和不同载波频率条件下的传播损耗中值。画出相应的曲线。
2.将计算结果和通过奥村模型实测测得的结果进行比较，验证计算结果的正确性。
3.分析Okumura-Hata方法在何距离以及何频率范围内存在较大的误差。
【实验内容】
1.大城市
clear;
hb=200;
hm=3;
for d=[125103060100]
f1=150:0.1:300;
Lb11=69.55+26.16*log10(f1)-13.82*log10(hb)-(8.29*log10(1.54*hm).^2-1.1) +(44.9-6.55*log10(hb))*log10(d);
f2=300:0.1:1920;
Lb12=69.55+26.16*log10(f2)-13.82*log10(hb)-(3.2*(log10(11.75*hm)).^2-4.
97)+(44.9-6.55*log10(hb))*log10(d);
f=[f1f2];
Lb1=[Lb11Lb12];
figure(1);
hold on;
plot(f,Lb1,'r');
end
grid;
title('大城市');
xlabel('频率/MHz');
ylabel('损耗中值/dB');
0200400600800100012001400160018002000
90100
110
120
130
140150160
170
180
190
大城市
频率/MHz 损耗中值/d B 2.郊区clear;hb=200;hm=3;for d=[125103060100]f1=150:0.1:300;Lb11=69.55+26.16*log10(f1)-13.82*log10(hb)-(8.29*log10(1.54*hm).^2-1.1)+(44.9-6.55*log10(hb))*log10(d);Lb21=Lb11-2*(log10(f1/28)).^2-5.4;
f2=300:0.1:1920;Lb12=69.55+26.16*log10(f2)-13.82*log10(hb)-(3.2*(log10(11.75*hm)).^2-4.97)+(44.9-6.55*log10(hb))*log10(d);Lb22=Lb12-2*(log10(f2/28)).^2-5.4;f=[f1f2];Lb2=[Lb21Lb22];figure(1);hold on ;plot(f,Lb2,'b');end grid;title('郊区');xlabel('频率/MHz');ylabel('损耗中值/dB');
0200400600800100012001400160018002000
8090
100
110
120
130140
150
160
170
郊区
频率/MHz 损耗中值/d B 3.开阔区clear;hb=200;hm=3;for d=[125103060100]f1=150:0.1:300;
Lb11=69.55+26.16*log10(f1)-13.82*log10(hb)-(8.29*log10(1.54*hm).^2-1.1)+(44.9-6.55*log10(hb))*log10(d);Lb21=Lb11-2*(log10(f1/28)).^2-5.4;Lb31=Lb11-4.78*(log10(f1)).^2+18.33*log10(f1)-40.98;f2=300:0.1:1920;Lb12=69.55+26.16*log10(f2)-13.82*log10(hb)-(3.2*(log10(11.75*hm)).^2-4.97)+(44.9-6.55*log10(hb))*log10(d);Lb22=Lb12-2*(log10(f2/28)).^2-5.4;Lb32=Lb12-4.78*(log10(f2)).^2+18.33*log10(f2)-40.98;f=[f1f2];Lb3=[Lb31Lb32];figure(3);hold on ;plot(f,Lb3,'g');end grid;title('开阔区');xlabel('频率/MHz');ylabel('损耗中值/dB');
0200400600800100012001400160018002000
6070
80
90
100
110120
130
140
150
开阔区
频率/MHz 损耗中值/d B 4.clear;hb=200;hm=3;for d=[125103060100]f1=150:0.1:300;
Lb11=69.55+26.16*log10(f1)-13.82*log10(hb)-(8.29*log10(1.54*hm).^2-1.1)+(44.9-6.55*log10(hb))*log10(d);Lb21=Lb11-2*(log10(f1/28)).^2-5.4;Lb31=Lb11-4.78*(log10(f1)).^2+18.33*log10(f1)-40.98;f2=300:0.1:1920;Lb12=69.55+26.16*log10(f2)-13.82*log10(hb)-(3.2*(log10(11.75*hm)).^2-4.97)+(44.9-6.55*log10(hb))*log10(d);Lb22=Lb12-2*(log10(f2/28)).^2-5.4;Lb32=Lb12-4.78*(log10(f2)).^2+18.33*log10(f2)-40.98;f=[f1f2];Lb1=[Lb11Lb12];Lb2=[Lb21Lb22];Lb3=[Lb31Lb32];figure(1);hold on ;plot(f,Lb1,'r');plot(f,Lb2,'b');plot(f,Lb3,'g');end grid;title('大城市—红，郊区—蓝，开阔区—绿');xlabel('频率/MHz');