• 大小: 0M
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2021-06-17
  • 语言: 其他
  • 标签: 其他  

资源简介

energydetection.zip

资源截图

代码片段和文件信息

clear all;
clc;
n =10000;                     %采样点数
pfa = 0.01;                   %虚警概率
snr1=-25:1:0;               %信噪比
k  = length(snr1);
for i = 1:k                 %不同信噪比下仿真
    s = 0;                  %1个数初值设为0,每次循环都清零
    snr(i) = power(10snr1(i)/10);%db转化十进制
    for kk = 1:1         %每种信噪比下仿真次数 
        t = 1/n:1/n:1;              %离散点1000
        x = sin(pi*t);
        x = ( x > 0);                %大于0判为1,否则0,方波产生方法
        x  =2*x-1;                   %信源数据产生
        noise = randn(1length(t));  %均值0方差1的噪声
        pn = (std(noise)).^2;        %噪声功率
        a = sqrt(2*snr(i));       %信号幅值
        xx = a.*sin(2*pi*4*t);       %调制信号    
        y  = x.*xx;                  %bpsk调制信号
        ps = mean(abs(y).^2);        %信号平均功率
        yy = y+noise;                %输出信号
       
        r = (n+sqrt(2*n)*sqrt(2)*erfcinv(2*pfa));
        %r = gaminv(1-pafn/21)*2;
        sum1 = sum(abs(yy).^2);      %输出信号的累积功率
        if sum1>r
            s = s+1;                 %输出大于阈值及判为1,否则0,统计1的个数
        end
    end
 p(i) = s/100;             %检测概率:1的个数/n
pd(i) = 0.5*erfc((erfcinv(2*pfa)-snr(i)*sqrt(n/2))/sqrt(2+4*snr(i)));
%awgn channels
%pd(i) = marcumq(sqrt(snr(i)*2)sqrt(r)n/2);
%pd(i)=pfa+exp(1/((2*n)*(snr(i)^2))-sqrt(2)*erfcinv(2*pfa)/(sqrt(n)*snr(i)))
%rayleigh fading(pfa要很小《0.01)
end
plot(snr1pd‘o-‘);
hold on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1517  2007-10-23 16:21  energydetection.m

评论

共有 条评论