大尺度衰落下的路径损耗问题（2）

        上次我们主要提到了自由空间下的路径损耗问题，但在实际通信中可能难以达到这种条件。这时就产生了一种新的信道模型——Okumura模型。

        20世纪60年代，奥村（Okumura）等研究人员在东京近郊进行了广泛的频率测量，涵盖多种基站和移动台天线高度，并考察了各种不规则地形和环境条件下的信号强度。他们绘制了一系列曲线图，展示了不同频率下场强与距离之间的关系，其中基站天线的高度作为曲线的参数。这些曲线图包括了在不同地形和环境条件下的结果，比如开阔地与市区的平均场强与距离的关系、市区平均场强与频率的关系，以及市区与郊区之间的差异。此外，他们提供了郊区修正因子的曲线、信号强度随基站天线高度变化的曲线，以及移动台天线高度对信号强度的相互关系曲线等。同时，还提供了不同地形的修正数据。

        为了简化计算机预测，研究人员引入了Hata模型，这是基于奥村先生大量测试数据拟合得到的公式，即Okumura-Hata模型。这种方法类似于科学研究中的做法，通过数学模型和公式来总结和简化复杂的实验数据，使得信号传播特性可以更便捷地进行计算和预测。hata模型将Okumura模型扩展到各种传播环境，实际上这种模型也是当今最常用的路径损耗模型。设置载波频率Hz，基站距离大约为1-1000m，发射机高度以及接收机高度分别为30和2m。对比出三种环境下的路径损耗模型。《MIMO-OFDM无线通信技术及MATLAB实现》中的程序1.4中对于legend的语法有一点错误，已做修改，下面是书中的源码以及图片。

        程序1.4 对比不同的路径损耗模型

clear, clf
fc=1.5e9;
htx=30;%发射天线高度
hrx=2;%发射天线高度
distance=[1:2:31].^2; %距离
y_urban=PL_Hata(fc,distance,htx,hrx,'urban');%市区
y_suburban=PL_Hata(fc,distance,htx,hrx,'suburban');%郊区
y_open=PL_Hata(fc,distance,htx,hrx,'open');%开阔地
semilogx(distance,y_urban,'b-s', distance,y_suburban,'r-o', distance,y_open,'k-^')
grid on, axis([1 1000 40 110]),
title(['Hata PL model, f_c=',num2str(fc/1e6),'MHz'])
xlabel('Distance[m]'), ylabel('Path loss[dB]')
legend('urban','suburban','open area')

function PL=PL_Hata(fc,d,htx,hrx,Etype)
% Hata Model
% Input
%       fc    : carrier frequency[Hz]
%       d     : between base station and mobile station[m]
%       htx   : height of transmitter[m]
%       hrx   : height of receiver[m]
%       Etype : Environment Type('urban','suburban','open')
% output
%       PL    : path loss[dB]
if nargin<5, Etype = 'URBAN'; end
fc=fc/(1e6);
if fc>=150&&fc<=200, C_Rx = 8.29*(log10(1.54*hrx))^2 - 1.1;
 elseif fc>200, C_Rx = 3.2*(log10(11.75*hrx))^2 - 4.97; 
 else   C_Rx = 0.8+(1.1*log10(fc)-0.7)*hrx-1.56*log10(fc);
end
PL = 69.55 +26.16*log10(fc) -13.82*log10(htx) -C_Rx ...
     +(44.9-6.55*log10(htx))*log10(d/1000);
EType = upper(Etype);
if EType(1)=='S',  PL = PL -2*(log10(fc/28))^2 -5.4;
  elseif EType(1)=='O' 
    PL=PL+(18.33-4.78*log10(fc))*log10(fc)-40.97;
end

        对于"PL_Hata"函数展示了Okumura-Hata模型的具体实现，具体的公式可以参照原书，函数已经很清晰的显示出书上的公式。当不输入环境类型时，函数将默认为城市模型。下面是具体的仿真图。

注:原书代码都为老版本代码，对于legend函数添加图例，增加了修改字体功能，如果后面直接加入字号可能会导致无法运行代码，后续代码也有同样的问题，不再一一赘述，大家注意一下就行，后续还会增加一些对于一些函数的改动与新增会单独做一个总结和博客，敬请期待。