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