ex1.m

循环自相关函数工具箱

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));