False alarm rate(错误报警率, Rayleigh分布）

function const_false_alarm_rate

% 生成 均值为10 的 Rayleigh 随机数
% 均值为 b*sqrt(pi/2)
% 方差为 （4-pi)*b^2/2
b = 10 / sqrt(pi/2);

% 生成 1,000,000 个随机数
ran = raylrnd(b, [1, 1e6]);

histogram(ran,60);
mean(ran)


% 阈值 db
Tdb = 5:13;
Tv = mean(ran) * 10.^(Tdb/20);

% number of false alarm
Nfa = zeros(size(Tv));
for k=1:length(Tv)
    Nfa(k) = length( find(ran > Tv(k)) );
end

far = Nfa / length(ran);

I = find(far < 1e-4,1);

Tdb(I)
end