doapdf.m

来自「GBSBEM(几何单反射椭圆模型)可以用来仿真功率－延迟－角度（pda-图）、功」· M 代码 · 共 28 行

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

% 对于多径归一化延迟r=1.6，生成100000个DOA
y = rand(1,100000);
r = 1.6;
a = DOA(y,r,10000);

% 绘成直方图
[N_sample,x] = hist(a,200);  % 得到直方图参数
subplot(2,1,1)
bar(x*180/pi,N_sample,1)
ylabel('样本数')
xlabel('DOA（度）')

% 由直方图得到pdf，并与理论pdf进行比较
theorypdf = (r^2 - 1)^(3/2) * (r^2 - 2*r*cos(x) + 1) ./ (pi * (2*r^2 - 1) * (r - cos(x)).^3) ;     % 根据公式计算theoryh(x)
del_x = x(3) - x(2);
p_hist = N_sample/100000/del_x;    % 由直方图得到pdf

subplot(2,1,2)
plot(x*180/pi,theorypdf,'k',x*180/pi, p_hist,'+k')  % 比较
axis([-180 180 0 max(p_hist)+0.1])
grid on
ylabel('概率密度函数')
xlabel('DOA（度）')
legend('理论pdf','归一化直方图')