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


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; 
T=300; 
max_tau=0; 
R=cyclic_autocorrelation_fast(x,T,max_tau); 
y=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,y);