cooperative energy majoritypd 认知无线电协作频谱感知

认知无线电协作频谱感知，能量检测，硬判决，Voting模式

clear all;
close all;
clc
tic;
k=1000;
fb=1*k*k;%基带码速率
fc=10*k*k;%载波频率
fs=64*k*k;%采样频率
N=25;%基带码源数
Q=fs/fb;%一个基带码源的采样点数
NN=N*Q;%总点数
K=100000;%%%%%统计次数
T=zeros（1K）;%%%%%信息熵的比值
% Nband=8;
NNN=5;%%%%%%%%%%%%总共5个用户
% delta=0.1;
% Threshold1=1.1251;
% Threshold2=Threshold1-delta;
% Threshold3=Threshold1+delta;
Threshold=1.0310;
Ns=768;

Energya=zeros（NNNK）;
Energyb=zeros（NNNK）;
Energyc=zeros（NNNK）;
Energy=zeros（NNNK）;
Dfirst=zeros（NNNK）;
Dsecond=zeros（NNNK）;
numfirsta=zeros（1K）;%%%%%%%%%%%%%%%%%%%%第一次检测有用户的数目
numfirstb=zeros（1K）;%%%%%%%%%%%%%%%%%%%%第一次检测无用户的数目
numsecond=zeros（1K）;%%%%%%%%%%%%%%%%%%%%第二次检测有用户的数目
Dfinal=zeros（1K）; 
numseconddetection=0;%%%%%%%%%%%%%%%%%%需要进行第二次检测的次数
for kk=1:1:K%%%%%%%%%%%%%%%%%%%%%%%%%%产生信号
    msend=sign（randn（1N））;%产生随机分布的基带码元
    x=zeros（1NN）;
    t=1:1:Q;
    carrier=cos（2*pi*fc*（t-1）/fs）;
    for i=1:1:N
        x（（i-1）*Q+1:i*Q）=carrier*msend（i）;
    end
    % figure（1）;
    % plot（x）;


    SNRset=-12;
    
    noise=randn（NNNlength（x））;
    NOISE=sum（noise（1:）.^2）/length（x）;
    D=sqrt（2*10^（SNRset/10））;
    x=x*D;
    S=sum（x.^2）/length（x）;
    
    SNRreal=10*log10（S/NOISE）;

    xnoise=zeros（NNNNN）;
    H=1;%%%%%定义信号的有无
    for j=1:1:NNN
        xnoise（j:）=H*x+noise（j:）;
    end
    
   
    for j=1:1:NNN%%%%%%%%%%%%%%%%%%%%第一次检测
        Energya（jkk）=mean（xnoise（jNs+1:2*Ns）.^2）;
        if Energya（jkk）>=Threshold
            Dfirst（jkk）=1;
            numfirsta（kk）=numfirsta（kk）+1;
        elseif Energya（jkk）            Dfirst（jkk）=0;
            numfirstb（kk）=numfirstb（kk）+1;
        end
    end
    
    if numfirsta（kk）>2%%%%%%%%%%%%第一次便可以判决有用户
        Dfinal（kk）=1;
    else
        Dfinal（kk）=0;
    end
        

end

Pd=sum（Dfinal）/K




figure（3）;
plot（Energy（1:）‘k‘）;
xlabel（‘Simulation times‘）;ylabel（‘Energy‘）;





toc;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2058  2012-02-10 16:57  cooperative_energy_thre_5users_majoritypd.m

----------- ---------  ---------- -----  ----

                 2058                    1