cost_231_walfisch_ikegami.m

移动通信中的信道仿真COST_231_Walfisch_Ikegami

function pld=COST_231_Walfisch_Ikegami(distance)

carry_f=848.97;%(MHz)
h_Base=50;%(m)4-50m
h_Mob=2;%(m)1-3m
h_Roof=12;
w_street=15;
ang=90;
dis_betwen_build=30;
width_str=15;
city=1;% 1 represent big city ;0 represent middle city or subburd
deltah_Base=h_Base-h_Roof;
deltah_Mob=h_Roof-h_Mob;
L0_no_d=32.4+20*log10(carry_f);
L_rts=-16.9-10*log10(width_str)+10*log10(carry_f)+20*log10(deltah_Mob);
if  0<ang<35
    L_ori=-10+0.354*ang;
else
    if  35<=ang<=65
        L_ori=2.5+0.075*(ang-35);
    else
        L_ori=4.0-0.114*(ang-55);
    end
end
L_rts=L_rts+L_ori;
if h_Base>h_Roof
    L_bsh=-18*log10(1+deltah_Base);
else
    L_bsh==0;
end
if h_Base>h_Roof
    k_d=18;
else
    k_d=18-15*deltah_Base/h_Roof;
end
if city==1
    k_f=-4+1.5*(carry_f/925-1);
else
    k_f=0.7*(carry_f/925)-1;
end  
Lmsd_no_d=L_bsh+k_f*log10(carry_f)-9*log10(dis_betwen_build);
    if h_Base>h_Roof
        k_a=54;
    else
        if distance>=0.500&h_Base<=h_Roof
            k_a=54-0.8*deltah_Base;
        else
            k_a=54-0.8*deltah_Base/0.5;
        end
    end   
    L_msd=Lmsd_no_d+k_a+k_d*log10(distance);
    L0=L0_no_d+20*log10(distance);
    L_city_bar=L0+L_rts+L_msd;
   pld=-L_city_bar;