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

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=0:2:300-2; 
figure(2) 
plot(ff(1:150),y(1:150));