gbsbem.m

GBSBEM(几何单反射椭圆模型)可以用来仿真功率－延迟－角度（pda-图）、功率延迟图、时间－角度的联合统计、波达方向的边缘特征和窄带衰落包络。适用于低层系统

% 根据下述已知条件生成GBSBEM模型多径分量参数，并绘制功率－延迟－角度分布图，功率－延迟分布图和功率－角度分布图
% 已知条件为：
% 1、发射机、接收机间距 d0＝500m
% 2、路径损耗指数 n=2.5
% 3、参考距离和参考功率 dref＝100m，Pref=10
% 4、散射体反射损耗 Lr=6dB
% 5、归一化最大路径延迟 rm＝2
% 6、多径数 L

clear all
d0 = 500;
n = 2.5;
Pref = 10;
dref = 100;
Lr = 4.5;
rm = 2;
L = 5;

% 计算LOS分量
[P0,tao] = LOS(Pref,dref,d0,n);     

for i = 1:L-1
    r(i) = TOA(1,rm) * tao * 1000000;      % 生成一个归一化多径延时
    y = rand(1);    % 生成一个服从U(0,1)的随机变量
    phi(i) = DOA(y,r(i),1000);  % 用y生成到达角（DOA）——phi
    theta = 2 * atan(sin(phi(i)) / ((r(i)+1)/(r(i)-1) - cos(phi(i))));  % 用DOA计算出离开角（DOD）——theta    
    Pi = P0 - 10 * n * log10(r(i)) - Lr;   % 计算其他多径分量功率    
    gamma = rand(1) * 2 * pi;   % 生成服从U(0,2pi)的随机相位
    a(i) = 10 ^ ((Pi - P0) / 20) * exp(j * gamma * pi /180);   % 得到复数幅度 
end

Z = plot3Dfigure(r,a,phi,tao,L);    % 绘制功率－延迟－角度分布图
plot2Dfigure(Z);    % 计算并绘制功率－延迟分布图和功率－角度分布图