检测概率、误报警概率（ROC曲线）

function prob_detect_and_prob_false_alarm

nlook = 100000;
TdB = 8:0.5:12;
Pd = zeros(1,length(TdB));
Pfa = zeros(1,length(TdB));
b = 2 / sqrt(pi/2);             % 均值为2的Rayleigh分布的参数
for k = 1:length(TdB)           %  遍历每个检测， dB
    Nd = 0;
    Nfa = 0;
    for look=1:nlook
        target = randn(1,100) + 7;          % 正态分布，mu=7,sigma=1
        system_noise = raylrnd(b,1,100);    % Rayleigh分布，均值 2
        signal = target + system_noise;      % 每个look， 在每个range gate里面采集的信号
        
        Td = mean(target) * 10^( TdB(k) / 20);
        Tfa = mean(system_noise) * 10^( TdB(k) / 20);
        Td = Tfa;
        %plot(signal);
        %pause;
        
        Id = find(target > Td, 1);
        Ifa = find(system_noise > Tfa, 1);
        
        if ~isempty(Id)
            Nd = Nd + 1;
        end
        if ~isempty(Ifa)
            N