% 根據(jù)下述已知條件生成GBSBEM模型多徑分量參數(shù)，并繪制功率－延遲－角度分布圖，功率－延遲分布圖和功率－角度分布圖
% 已知條件為：
% 1、發(fā)射機、接收機間距 d0＝500m
% 2、路徑損耗指數(shù) n=2.5
% 3、參考距離和參考功率 dref＝100m，Pref=10
% 4、散射體反射損耗 Lr=6dB
% 5、歸一化最大路徑延遲 rm＝2
% 6、多徑數(shù) 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);   % 得到復(fù)數(shù)幅度 
end

Z = plot3Dfigure(r,a,phi,tao,L);    % 繪制功率－延遲－角度分布圖
plot2Dfigure(Z);    % 計算并繪制功率－延遲分布圖和功率－角度分布圖