fft_size=length(t);
F1=fft(s3,fft_size);
pxx=abs(F1).^2/fft_size;
P=10*log(pxx(1:fft_size/2));
f=2*pi*fs*(0:round(fft_size/2-1))/fft_size;
subplot(4,4,5);
plot(f,P);我一般是这样求的
----------------------------------------
这要看你的对象和目的.
1.如果你的对象是周期信号,想得到周期信号的参数,那么幅值谱就可以了.此时的步骤是:
对得到的复数取模,然后除以处理的数据长度,乘以2,得到单边的幅值谱.
2.如果你的对象是平稳随机过程,那么你能够定量是功率谱密度(功率在频带上"连续"分布),那么步骤是:对得到的复数取模,然后除以数据的长度,再平方,然后乘以2,再除以频率分辨率,得到的就是单边功率谱;!!!!!!!!!!必须重复N次后平均!!!!!!!!!!!!!!
-----------------------------------------
文章评论(0条评论)
登录后参与讨论