% 功能：對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','歸一化直方圖')