ex4.m

循环自相关函数工具箱

clear;
clc;
fc=12;fb=7;fa=100;fs=600;N=1024;
t=0:1/fs:1024/fs;
x=(1+cos(2*pi*fb*t)).*cos((2*pi*fa*t)+cos(2*pi*fc*t));
figure(1)
plot(t,x);
axis([0 1 -2 2]);
xlabel('\fontsize{9}\fontname{标楷体}时间/s');ylabel('\fontsize{9}\fontname{标楷体}幅值/um');
x1=x.*x;
y=fft(x1,N);
Pyy=sqrt(y.*conj(y))/N;
figure(2)
ff=fs*(0:N/2)/N;
plot(ff(1:N/2+1),Pyy(1:N/2+1));
xlabel('\fontsize{9}\fontname{标楷体}循环频率 /Hz');ylabel('\fontsize{9}{R_{x}}^{\alpha}(0) /um^2');