ex3.m

来自「循环自相关函数工具箱」· M 代码 · 共 18 行

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