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