COST231-WI模型及MATLAB仿真

COST231-WI模型介绍

COST231-WI模型（COST231-Walfish-Ikegami）是基于大量现场测试和数据分析获取的Walfisch模型与Ikegami-Bertonin模型组合，适用于建筑物高度一致的平坦郊区、市区及建筑物密集和人口密集区域，小区覆盖半径可以小于1km。该模型考虑了天线多屏正向衍射损耗，提供更好的路径损耗预测。若基站和移动接收机之间的信号在视距（LOS）环境下传播，则LOS环境下的路径损耗公式为：

其中，为收发机的实际距离(km)；为发射信号频率(MHz)；为LOS环境下模型预测的路径损耗值(dB)。

图 1 城市及郊区典型 NLOS 传播下定义COST231-WI模型使用的相关参数示意图

若基站与移动接收机之间的信号在非视距（NLOS）环境下传播，其采用COST231-WI模型中的相关参数如图 1 所示。其中，为收发机之间的实际距离（km），表示基站发射天线的高度（m），表示周围建筑物的平均高度（m），表示移动接收机天线的高度(m)，表示周围建筑物的平均间隔(m)，表示街道的宽度(m)，表示基站发射天线与建筑物的高度差(m),表示建筑物与接收机天线的高度差(m),表示信号传播直接路径与道路方向的夹角（°）。该模型在NLOS环境下预测的路径损耗值为，其路径损耗公式为：

为自由空间的路径损耗值，其与距离、频率成正比关系：

为信号在传播过程中绕射与散射引起的衰落损耗值，与及成正比，与成反比：

在上式中，是街道方向衰落函数，取决于参数：

为信号沿屋顶传播时多屏绕射引起的衰落损耗值，与、成正比，与成反比：

 

其各个参数的计算公式如下：

 其中，

当周围建筑物屋顶高度比基站天线所在位置高时，产生的路径损耗增加值用表示：

表示对应基站天线与周围建筑物屋顶高度差导致的多屏绕射损耗值：

表示对应基站发出射频信号的多屏绕射损耗值：

 MATLAB仿真

在LOS和NLOS两种情况下，分别对800MHz、1500MHz、2000MHz三种频率的无线电信号的传播损耗进行仿真。

function Lp = cost231_WI(propagation_type,city_type,fc,d,ht,hm,hroof,alpha,w,s)
% propagation_type:LOS为0，NLOS为1
% city_type:0代表中小城市，1代表大都市
% fc:发射信号频率（MHz）
% d:收发机的实际距离（km）
% ht:基站天线的高度（m）
% hm:辐射源天线的高度（m）
% hroof:周围建筑物的平均高度（m）
% alpha:信号传播直接路径与道路方向的夹角（°）
% w:街道的宽度(m)
% s:周围建筑物的平均间隔(m)
if propagation_type == 0
    Lp = 42.6 + 26*log10(d) + 20*log10(fc);
    else
        dht = ht - hroof;
        dhm = hroof - hm;
        Lfs = 32.4 + 20*log10(d) + 20*log10(fc); %自由空间的路径损耗值
        if (0 <= alpha) && (alpha < 35)
            Lori = -10 + 0.354*alpha; %街道方向衰落函数
            elseif (35 <= alpha) && (alpha < 55)
                Lori = 2.5 + 0.075*(alpha - 35);
                else
                    Lori = 4 - 0.114*(alpha - 55);
        end
        Lrts = -16.9 - 10*log10(w) + 10*log10(fc) + 20*log10(dhm) + Lori;% 绕射与散射引起的衰落损耗值

        if dht > 0
            Lbhs = -18*log10(1 + dht);
            ka = 54;
            kd = 18;
        else
            Lbhs = 0;
            kd = 18 - 15*dht/hroof;
            if d >= 0.5
                ka = 54 - 0.8*dht;
            else
                ka = 54 - 1.6*dht*d;
            end
        end

        if city_type == 0
            kf = -4 + 0.7*(fc/925 - 1);
        else
            kf = -4 + 1.5*(f/925 - 1);
        end
        Lmsd = Lbhs + ka + kd*log10(d) + kf*log10(fc) - 9*log10(s);% 信号沿屋顶传播时多屏绕射引起的衰落损耗值
        Lp = Lfs + Lrts + Lmsd;
end

end

上面这段代码将COST231-WI模型封装成了一个函数，下面对其进行调用及结果显示。

clc;clear;close all;
propagation_type = [0,1];
city_type = 0;
fc = [800,1500,2000];
ht = 24;
hm = 1.5;
hroof = 21;
alpha = 45;
w = 15;
s = 30;

d = 0.01:0.1:5;
Lp1 = cost231_WI(propagation_type(1),city_type,fc(1),d,ht,hm,hroof,alpha,w,s);
Lp2 = cost231_WI(propagation_type(1),city_type,fc(2),d,ht,hm,hroof,alpha,w,s);
Lp3 = cost231_WI(propagation_type(1),city_type,fc(3),d,ht,hm,hroof,alpha,w,s);

Lp4 = cost231_WI(propagation_type(2),city_type,fc(1),d,ht,hm,hroof,alpha,w,s);
Lp5 = cost231_WI(propagation_type(2),city_type,fc(2),d,ht,hm,hroof,alpha,w,s);
Lp6 = cost231_WI(propagation_type(2),city_type,fc(3),d,ht,hm,hroof,alpha,w,s);

figure
plot(d,Lp1,'r--')
hold on;
plot(d,Lp2,'b--')
hold on;
plot(d,Lp3,'c--')
hold on;
plot(d,Lp4,'g-.')
hold on;
plot(d,Lp5,'m-.')
hold on;
plot(d,Lp6,'k-.')
xlabel('距离(km)');
ylabel('路径损耗(dB)');
title('cost231-WI模型路径损耗');
legend({'LOS   fc=800MHz','LOS   fc=1500MHz','LOS   fc=2000MHz','NLOS fc=800MHz','NLOS fc=1500MHz','NLOS fc=2000MHz'},Location="southeast");

运行结果如下：

可以看到，NLOS的路径损耗比LOS要大得多，而且在相同的条件下，路径损耗随着信号频率的增加而增加。