姓名：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');
大城市190180170160损耗中值/dB15014013012011010090020040060080010001200频率/MHz1400160018002000
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');
郊区170160150140损耗中值/dB1301201101009080020040060080010001200频率/MHz1400160018002000
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');