说明：  function [Sx,f,alpha] = Sa(s,fs,M,N,read_f,del_f) 计算信号的循环谱 [Sx,f,alpha] = Sa(s,fs,M,N,read_f,del_f) 输入： s 序列 fs 采用频率 M 滑动频率窗长 N fft计算点数 read_f 可视频率（一般为载频的整数倍） del_f 可视频率间隔 输出： Sx 信号的循环谱 f 频率轴 alpha 循环频率轴
(Cyclosationary Spectrum)