www.pudn.com > cyclostationary_toolbox.rar > ex1.m, change:2007-10-13,size:349b


clear; 
clc; 
t=0:1/600:1023/600; 
A=1;B=1;fb=7;fa=100; 
x=A*(1+B*cos(2*pi*fb*t)).*cos(2*pi*fa*t); 
figure(1) 
plot(t,x); 
axis([0 1 -2 2]); 
alpha=2*pi/300; 
max_tau=0; 
R=cyclic_autocorrelation(x,alpha,max_tau); 
y=sqrt(R.*conj(R)); 
% ff=600*(0:1023)/512; 
ff=0:1:300-1; 
figure(2) 
% plot(ff(1:150)/(2*pi),y(1:150)); 
plot(ff(1:150)*2,y(1:150));