内容发布更新时间 : 2024/11/10 12:55:05星期一 下面是文章的全部内容请认真阅读。

姓名：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');