xn=signal;;%信号cxn=xcorr(xn,'unbiased'); %计算序列的自相关函数CXk=fft(cxn,Nfft);%对自相关函数进行傅里叶变换Pxx=abs(CXk);index=0:round(Nfft/2-1);k=index*Fs/Nfft;%频谱序列subplot(3,1,3),plot(k,Pxx(index+1));%绘制功率谱曲线 首先确定一下数据位数
