频谱感知中，虚警概率与检测概率的关系曲线

pf=0.1:0.01:0.9;%虚警概率 
N=500;%采样点数
u=1;
s=0.0999;
snr=-10;
k  = length（pf）;
for i=1:k%不同pf下仿真
r（i）=0;%清零 
for kk=1:5000%仿真次数 
  t=0:0.001:0.5; 
  n=randn（1length（t））; 
  x=0.447*sin（2*pi*4*t）; 
  y=x+n; 
  sum1（i） = （1/N）*sum（abs（y）.^2）;      %输出信号的平均功率 
  g=u*（1+qfuncinv（pf（i））/sqrt（N））;%判决门限 
     if sum1（i）>g 
              r（i