实验三_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=[1 2 5 10 30 60 100]
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=[f1 f2];
Lb1=[Lb11 Lb12];
figure(1);
hold on;
plot(f,Lb1,'r');
end
grid;
title('大城市' );
xlabel('频率/MHz');
ylabel('损耗中值/dB');
020040060080010001200
1400160018002000
90100
110
120
130
140150160
170
180
190
大城市
频率/MHz 损耗中值/d B
2.郊区
clear;
hb=200;
hm=3;
for d=[1 2 5 10 30 60 100]
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=[f1 f2];
Lb2=[Lb21 Lb22];
figure(1);
hold on ;
plot(f,Lb2,'b');
end
grid;
title('郊区' );
xlabel('频率/MHz');
ylabel('损耗中值/dB');
020040060080010001200
1400160018002000
8090
100
110
120
130140
150
160
170
频率/MHz 损耗中值/d B
3.开阔区
clear;
hb=200;
hm=3;
for d=[1 2 5 10 30 60 100]
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=[f1 f2];
Lb3=[Lb31 Lb32];
figure(3);
hold on ;
plot(f,Lb3,'g');
end
grid;
title('开阔区' );
xlabel('频率/MHz');
ylabel('损耗中值/dB');