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