toapdf.m

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

% 功能：对TOA的理论pdf与仿真pdf进行比较

% 将最大路径延迟rm设为2，生成10000个TOA
rm = 2;
r = TOA(10000,rm);

% 绘成直方图
[N_sample,x] = hist(r,200);  % 得到直方图参数
subplot(2,1,1)
bar(x,N_sample,1)
ylabel('样本数')
xlabel('TOA（归一化）')

% 由直方图得到pdf，并与理论pdf进行比较
beta = rm * sqrt(rm ^ 2 - 1);
theorypdf = (2 * x.^2 -1) ./ (beta * sqrt(x.^2 -1));     % 根据公式计算theoryh(x)
del_x = x(3) - x(2);
p_hist = N_sample/10000/del_x;    % 由直方图得到pdf
subplot(2,1,2)
plot(x,theorypdf,'k',x, p_hist,'+k')  % 比较
grid on
ylabel('概率密度函数')
xlabel('TOA（归一化）')
legend('理论pdf','归一化直方图')