Gamma cdf 函数返回 1 的数组

如何解决Gamma cdf 函数返回 1 的数组

clc
clear all
close all
dth_dB_vec = -5:5:15
epsilon = 5
v = 4
qp = 4
var = 10
Mp = 2^(qp)
Detection_per_vec = [];
for i = 1:length(dth_dB_vec)
    dthdB = dth_dB_vec(i);
    dth_Natural = 10^(dthdB/10);
    Detection_per = 0;
    Pfa = 1-(gamcdf(dth_Natural,Mp,(Mp/2*v*epsilon*var)));
    Detection_per = Pfa
    Detection_per_vec = [Detection_per_vec Detection_per];
end
figure(1)
semilogy(dth_dB_vec,Detection_per_vec,'--b')
hold on 
grid on
ylabel('Probability')
xlabel('Detection threshold [dB over noise level')

我正在尝试绘制 gamma CDF 函数的 FA 率，但是当结果向量只返回 1 时，我不知道我哪里出错了。想要的图随着 dth_dB_vec 的每次迭代而向下弯曲。所有参数都正确。这是原始的数学公式。 False Alarm formula。感谢您的帮助