Fs =5120; % Sampling frequency T1 = 1/Fs; % Sample time L = size(x,1); % Length of signal NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(x,NFFT)/L; Z=2*abs(Y(1:NFFT/2+1)); 非常好用。
研究过程中实现的,分享出来。 是一种多种分析图绘制合一的完整函数并对生成谱进行自动的peaks筛选输出。已经形成自定义函数,放到特定路径下即可直接调用,配有较详细的使用说明 style=1,plot amplitude value spectrumï¼›%style=2,plot power spectrum;%style=3,plot loglog PSD %style=4,pwelch with segements number 此绘图函数可以输入多种调用参数,如'FrequencyBan
离散频谱的比值校正法-离散频谱分析的一种新校正方法.pdf 离散频谱的比值校正法 在下帖子中讨论了“Matlab中FFT求正玹序列的振幅”https://www.ilovematlab.cn/thread-50611-1-1.html 我在帖子中指出了当“正弦信号的频率不与FFT后的某条谱线相合”,可用校正法来求正弦信号的频率。在这里提供一个比值校正法的程序,它已编写成一个函数,该方法的理论可参看以下附件。 function Z=Specorr [nx,mx]=size; if mx==1,
离散频谱的比值校正法-频谱分析的校正方法.pdf 离散频谱的比值校正法 在下帖子中讨论了“Matlab中FFT求正玹序列的振幅”https://www.ilovematlab.cn/thread-50611-1-1.html 我在帖子中指出了当“正弦信号的频率不与FFT后的某条谱线相合”,可用校正法来求正弦信号的频率。在这里提供一个比值校正法的程序,它已编写成一个函数,该方法的理论可参看以下附件。 function Z=Specorr [nx,mx]=size; if mx==1, x=x;