关于能量检测BPSK 的matlab仿真代码

代码说明：

NR=-24:2:-6; q=10.^(SNR/10); b1=["r","g","b","k"]; b2=[":",":",":",":"]; b3=[".","+",">","o"]; n=[512,1024,1024*2,1024*5]; pf=0.1; count=0; for i=1:length(n)      r=n(i)+sqrt(2*n(i)).*erfinv(1-pf)*2;     for p=1:length(q)     for k=1:10000        m=randint(1,n(i),2); %m=rand(1,n(i));        s=dmodce(m,1,1,"psk",2);        noise=random("Normal",0,1,1,n(i));        real=q(p).*s+noise;        T=sum(abs(real.^2));        if(T>r) count=count+1;        end     end     pd(p)=count/10000;     count=0;     end      plot(SNR,pd,[b1(i) b2(i) b3(i)])      hold on end legen