powerspectrum.m

matlab平台下一组低维混沌随机序列产生的源程序

function PowerSpectrum(Type,gn)
%Chaotic mapping types
ChaosType={'Logistic','ImprovedLogistic','Tent','Chebyshev'};
SeqType={'L2GNS','IL2GNS','T2GNS','C2GNS'};


switch Type
case {'Logistic'}
    funHandle=@Logistic;
    titleHandle=SeqType(1);
    tHandle=ChaosType(1);
case {'ImprovedLogistic'}
    funHandle=@ImprovedLogistic;
    titleHandle=SeqType(2);
    tHandle=ChaosType(2);
case {'Tent'}
    funHandle=@Tent;
    titleHandle=SeqType(3);
    tHandle=ChaosType(3);
case {'Chebyshev'}
    funHandle=@Chebyshev;
    titleHandle=SeqType(4);
    tHandle=ChaosType(4);
otherwise
    disp('>>???Unkown Chaotic mapping!');
end

close all hidden;
figure
periodogram(gn,[],2048,length(gn));
axis([0,2500,-70,-25]);
title('Power Spectrum of Gauss Noise');
% title(['Power Spectrum of ',char(titleHandle)]);
xlabel('Frequency');
ylabel('Amplitude of Power Spectrum');

figure
plot(gn);grid on;
title('Time Curve of Gauss Noise');
% title(['Time Curve of ',char(titleHandle)]);
xlabel('Time');
ylabel('Amplitude');