• PUDN用户
    了解作者
  • matlab
    开发工具
  • 1KB
    文件大小
  • rar
    文件格式
  • 1
    收藏次数
  • 10 积分
    下载积分
  • 725
    下载次数
  • 2010-04-15 17:34
    上传日期
认知无线电中频谱感知技术的AM信号能量检测技术和周期平稳检测技术
chengxu.rar
  • energy.m
    595B
  • zhouqipinwen.m
    1.1KB
内容介绍
clc clear all close all dt=0.01; fs=1/dt; fm=1; fc=30; t=0:dt:30; mt=sqrt(2)*cos(2*pi*fm*t); %AM A=2; xn=(A+mt).*cos(2*pi*fc*t); noise=randn(size(xn)); N=256; cn=-20;%%% dB cn1=10^(-cn/10); po3=mean(xn.*conj(xn)); c3=po3/2*cn1; zxc=sqrt(c3).*randn(size(t));; y = xn + zxc a='hamming'; z='hamming'; if ~exist('L') L=400;%平滑的长度6 end lx=length(xn); if (length(y)~=lx) error('Time series x and y must be same length') end n=0:floor((lx-N)/L); ln=length(n);%滑动次数 a=(hamming(N))'; g=(hamming(ln))'; X=zeros(N,ln); Y=zeros(N,ln); Ts=1/N; for i=1:ln n_r=n(i)*L+(1:N); X(:,i)=fftshift(fft(a.*xn(n_r)))'; Y(:,i)=fftshift(conj(fft(a.*y(n_r))))'; end lnd2=floor(ln/2); lnd4=floor(ln/4)+1; ln3d4=lnd4+lnd2-2; for alpha=-N/2+1:N/2-1 for f=-N/2:N/2-1 f1=f+alpha; f2=f-alpha; if ((abs(f1)<N/2)&(abs(f2)<N/2)) f1=f1+N/2; f2=f2+N/2; fsh=fftshift(g.*(X(f1,:).*Y(f2,:))); fsh=sum(fsh)/ln; S(alpha+N/2,N/2+1+f)=fsh; end end end df=fs/N; alpha=(-N/2+1:N/2-1)*df*2; f=(-N/2:N/2-1)*df*2; R=abs(S); mesh(f,alpha,R)
评论
    相关推荐