实验一 移动信道传播模型的估算（Hata-Okumura模型）

一、实验目的

1、熟悉MATLAB程序编辑与开发环境，掌握使用该软件解决移动通信问题的具体方法；

2、掌握无线电波传播损耗计算的基本知识；

3、掌握运用Hata-Okumura模型预测的方法实现大城市、中小城市、郊区、开阔农村地带的无线传播损耗预测系统；

4、通过无线传播损耗预测仿真加深对Hata-Okumura无线传播损耗预测模型的理解。

二、实验内容

1、学习和理解Hata-Okumura无线传播损耗预测模型的原理及实现方法；

2、使用MATLAB软件产生范围300 ~ 1500MHz的若干工作频率点，以及范围1 ~ 10km的若干转播距离；

3、利用MATLAB软件仿真Hata-Okumura模型，分别实现大城市、中小城市、郊区、开阔农村地带的无线传播损耗预测系统，要求基站发射机的有效天线高度为200m，移动台接收机的有线天线高度为3m；

4、分别在大城市、中小城市和郊区三种环境情况下，对无线传播损耗预测值进行仿真和分析，在多组工作频率的条件下，画出传播距离-无线传播损耗的关系曲线图并进行性能分析；在多组传播距离的条件下，画出工作频率-无线传播损耗的关系曲线图并进行性能分析。

clc;clear;
%作者：YLS
%Hata-Okumura模型 f:中心频率; hb:基站天线有效高度; hm:移动台天线有效高度,单位:m; d:通信距离,单位:km; a_hm:移动台天线高度校正;
hm=3;i=0;hb=200;d1=1:0.5:10;
figure()
subplot(221);
for f1=300:300:1500
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld1=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;i=i+1;
    if i==1
        plot(d1,Ld1,'kx-');hold on;
    elseif i==2
        plot(d1,Ld1,'ks-');hold on;
    elseif i==3
        plot(d1,Ld1,'kO-');hold on;
    elseif i==4
        plot(d1,Ld1,'k>-');hold on;
    else
        plot(d1,Ld1,'k*-');hold on;
    end
end
legend('300MHz','600MHz','900MHz','1200MHz','1500MHz');xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');title('（大城市）传输损耗随传播距离的变化曲线');hold off;i=0;
subplot(222);
for f1=300:300:1500
    a_hm=(1.1*log10(f1)-0.7)*hm-(1.56*log10(f1)-0.8);Lx1=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;i=i+1;
    if i==1
        plot(d1,Lx1,'kx-');hold on;
    elseif i==2
        plot(d1,Lx1,'ks-');hold on;
    elseif i==3
        plot(d1,Lx1,'kO-');hold on;
    elseif i==4
        plot(d1,Lx1,'k>-');hold on;
    else
        plot(d1,Lx1,'k*-');hold on;
    end
end
legend('300MHz','600MHz','900MHz','1200MHz','1500MHz');xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');title('（中小城市）传输损耗随传播距离的变化曲线');hold off;i=0;
subplot(223);
for f1=300:300:1500
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;Lj1=Ld-2*(log10(f1/28))^2-5.4;i=i+1;
    if i==1
        plot(d1,Lj1,'kx-');hold on;
    elseif i==2
        plot(d1,Lj1,'ks-');hold on;
    elseif i==3
        plot(d1,Lj1,'kO-');hold on;
    elseif i==4
        plot(d1,Lj1,'k>-');hold on;
    else
        plot(d1,Lj1,'k*-');hold on;
    end
end
legend('300MHz','600MHz','900MHz','1200MHz','1500MHz');xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');title('（郊区）传输损耗随传播距离的变化曲线');hold off;i=0;
subplot(224);
for f1=300:300:1500
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f1)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d1)-a_hm;Ln1=Ld-4.78*(log10(f1))^2-18.33*log10(f1)-40.98;i=i+1;
    if i==1
        plot(d1,Ln1,'kx-');hold on;
    elseif i==2
        plot(d1,Ln1,'ks-');hold on;
    elseif i==3
        plot(d1,Ln1,'kO-');hold on;
    elseif i==4
        plot(d1,Ln1,'k>-');hold on;
    else
        plot(d1,Ln1,'k*-');hold on;
    end
end
legend('300MHz','600MHz','900MHz','1200MHz','1500MHz');xlabel('传播距离(单位km)');ylabel('传输损耗L(单位dB)');title('（开阔农村地带）传输损耗随传播距离的变化曲线');hold off;
%图二
hm=3;hb=200;i=0;f2=300:100:1500;
figure(2)
subplot(221);
for d2=1:3:10
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld2=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;i=i+1;
    if i==1
        plot(f2,Ld2,'kx-');hold on;
    elseif i==2
        plot(f2,Ld2,'ks-');hold on;
    elseif i==3
        plot(f2,Ld2,'kO-');hold on;
    else
        plot(f2,Ld2,'k*-');hold on;
    end
end
legend('1km','4km','7km','10km');xlabel('工作频率(单位Hz）');ylabel('传输损耗L(单位dB)');title('（大城市）传输损耗随工作频率的变化曲线');hold off;i=0;
subplot(222);
for d2=1:3:10
    a_hm=(1.1*log10(f2)-0.7)*hm-(1.56*log10(f2)-0.8);Lx2=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;i=i+1;
    if i==1
        plot(f2,Lx2,'kx-');hold on;
    elseif i==2
        plot(f2,Lx2,'ks-');hold on;
    elseif i==3
        plot(f2,Lx2,'kO-');hold on;
    else
        plot(f2,Lx2,'k*-');hold on;
    end
end
legend('1km','4km','7km','10km');xlabel('工作频率(单位Hz）');ylabel('传输损耗L(单位dB)');title('（中小城市）传输损耗随工作频率的变化曲线');hold off;i=0;
subplot(223);
for d2=1:3:10
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;Lj2=Ld-2*(log10(f2/28)).^2-5.4;i=i+1;
    if i==1
        plot(f2,Lj2,'kx-');hold on;
    elseif i==2
        plot(f2,Lj2,'ks-');hold on;
    elseif i==3
        plot(f2,Lj2,'kO-');hold on;
    else
        plot(f2,Lj2,'k*-');hold on;
    end
end
legend('1km','4km','7km','10km');xlabel('工作频率(单位Hz）');ylabel('传输损耗L(单位dB)');title('（郊区）传输损耗随工作频率的变化曲线');hold off;i=0;
subplot(224);
for d2=1:3:10
    a_hm=3.2*(log10(11.75*hm))^2-4.97;Ld=69.55+26.16*log10(f2)-13.82*log10(hb)+(44.9-6.55*log10(hb))*log10(d2)-a_hm;Ln2=Ld-4.78*(log10(f2)).^2+18.33*log10(f2)-40.98;i=i+1;
    if i==1
        plot(f2,Ln2,'kx-');hold on;
    elseif i==2
        plot(f2,Ln2,'ks-');hold on;
    elseif i==3
        plot(f2,Ln2,'kO-');hold on;
    else
        plot(f2,Ln2,'k*-');hold on;
    end
end
legend('1km','4km','7km','10km');xlabel('工作频率(单位Hz）');ylabel('传输损耗L(单位dB)');title('（开阔农村地带）传输损耗随工作频率的变化曲线');hold off;

运行结果如下：

图1.传播距离—无线传播损耗的关系

图2.工作频率—无线传播损耗的关系