FFT和功率谱估计

《FFT和功率谱估计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、FFT和功率谱估计1.用Fourier变换求取信号的功率谱---周期图法clf;Fs=1000;N=256;Nfft=256;%数据的长度和FFT所用的数据长度n=0:N-1;t=n/Fs;%采用的时间序列xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转化为dBf=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列subplot(2,1,1),plot(f,Pxx);%绘制功率谱曲线xlabel('频率/

2、Hz');ylabel('功率谱/dB');title('周期图N=256');gridon;Fs=1000;N=1024;Nfft=1024;%数据的长度和FFT所用的数据长度n=0:N-1;t=n/Fs;%采用的时间序列xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转化为dBf=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列subplot(2,1,2),plot(f,Pxx);%绘制功率谱曲线xl

3、abel('频率/Hz');ylabel('功率谱/dB');title('周期图N=256');gridon;2.用Fourier变换求取信号的功率谱---分段周期图法%思想：把信号分为重叠或不重叠的小段，对每小段信号序列进行功率谱估计，然后取平均值作为整个序列的功率谱clf;Fs=1000;N=1024;Nsec=256;%数据的长度和FFT所用的数据长度n=0:N-1;t=n/Fs;%采用的时间序列randn('state',0);xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);Pxx1=abs(fft(xn(1:256),Nsec).^2

4、)/Nsec;%第一段功率谱Pxx2=abs(fft(xn(257:512),Nsec).^2)/Nsec;%第二段功率谱Pxx3=abs(fft(xn(513:768),Nsec).^2)/Nsec;%第三段功率谱Pxx4=abs(fft(xn(769:1024),Nsec).^2)/Nsec;%第四段功率谱Pxx=10*log10(Pxx1+Pxx2+Pxx3+Pxx4/4);%Fourier振幅谱平方的平均值，并转化为dBf=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列subplot(2,1,1),plot(f(1:Nsec/2),Pxx(1:Nse

5、c/2));%绘制功率谱曲线xlabel('频率/Hz');ylabel('功率谱/dB');title('平均周期图（无重叠）N=4*256');gridon;%运用信号重叠分段估计功率谱Pxx1=abs(fft(xn(1:256),Nsec).^2)/Nsec;%第一段功率谱Pxx2=abs(fft(xn(129:384),Nsec).^2)/Nsec;%第二段功率谱Pxx3=abs(fft(xn(257:512),Nsec).^2)/Nsec;%第三段功率谱Pxx4=abs(fft(xn(385:640),Nsec).^2)/Nsec;%第四段功率谱Pxx5=abs(fft(xn(51

6、3:768),Nsec).^2)/Nsec;%第四段功率谱Pxx6=abs(fft(xn(641:896),Nsec).^2)/Nsec;%第四段功率谱Pxx7=abs(fft(xn(769:1024),Nsec).^2)/Nsec;%第四段功率谱Pxx=10*log10(Pxx1+Pxx2+Pxx3+Pxx4+Pxx5+Pxx6+Pxx7/7);%Fourier振幅谱平方的平均值，并转化为dBf=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列subplot(2,1,2),plot(f(1:Nsec/2),Pxx(1:Nsec/2));%绘制功率谱曲线xla

7、bel('频率/Hz');ylabel('功率谱/dB');title('平均周期图（重叠1/2）N=1024');gridon;1.用Fourier变换求取信号的功率谱---welch方法%思想：welch法采用信号重叠分段，加窗函数和FFT算法等计算一个信号序列的自功率谱（PSD）和两个信号序列的互功率谱（CSD），采用MATLAB自%带的函数psdclf;Fs=1000;N=1024;Nfft=256;