rar，dan.zip_energy detection_failedenm_单门限检测_能量检测

%多节点基于能量检测的频谱感知概率仿真 close all clear all clc SNR_dB=0;%信噪比dB N=3; K=10;%采样点数 sigma=1;%噪声功率 lambdas=0:2:12;%阈值 L=length(lambdas); Num_Run=50000;%仿真次数 %初始化输出变量 Pd=zeros(1,L); Pf=zeros(1,L); for ii=1:L Es=10^(SNR_dB/10)*sigma/K; lambda=lambdas(ii); for jj=1:Num_Run %存在主用户信号 bits=rand(1,K)>0.5;%随机二进制比特序列 x=sqrt(Es)*(2*bits-1);%BPSK调制，且进行功率调节 result=0; for n=1:N z=randn(1,K);%噪声，均值为0，方差为1 y=raylrnd(1/sqrt(2),size(x)).*x+z;%接收信号 energy=sum(abs(y).^2);%采样点能量和 if(energy>lambda)%能量大于门限 Pd(ii)=Pd(ii)+1/Num_Run; break; end end %不存在主用户信号 for n=1:N y=randn(1,K);%噪声，均值为0，方差为1,只有噪声 energy=sum(abs(y).^2); if(energy>lambda) Pf(ii)=Pf(ii)+1/Num_Run; break; end end end end figure, subplot(2,1,1) plot(lambdas,Pd,'-r') xlabel('Threshold') ylabel('概率') ylim([0,1]) grid on legend('Pd仿真值') subplot(2,1,2) plot(lambdas,Pf,'-g') xlabel('Threshold') ylabel('概率') ylim([0,1]) grid on legend('Pf仿真值')