|
楼主 |
发表于 2009-3-13 12:48
|
显示全部楼层
本帖最后由 wdhd 于 2016-6-30 09:20 编辑
谢谢大大的回覆,所以用STFT就可以找到共振範围噜?因为psd没概念,有找到论坛上的程序,但有些参数不知如何套用自己的资料想请教一下,Matlab代码如下:
- clear;
- Fs=15000; %我的取样率
- n=0:1/Fs:33; %我资料笔数有495000个,495000/15000=33
- nfft=1024; %这个不知要设多少
- window=boxcar(100); %矩形窗
- window1=hamming(100); %海明窗
- window2=blackman(100); %blackman窗,有叁种窗口,只要选一种就可以了吗?
- noverlap=20; %数据无重叠 ,如果想要重叠要设多少?
- range='half'; %频率间隔为[0 Fs/2],只计算一半的频率
- [Pxx,f]=pwelch(xn,window,noverlap,nfft,Fs,range);
- [Pxx1,f]=pwelch(xn,window1,noverlap,nfft,Fs,range);
- [Pxx2,f]=pwelch(xn,window2,noverlap,nfft,Fs,range);
- plot_Pxx=10*log10(Pxx);
- plot_Pxx1=10*log10(Pxx1);
- plot_Pxx2=10*log10(Pxx2);
- figure(1)
- plot(f,plot_Pxx);
- pause;
- figure(2)
- plot(f,plot_Pxx1);
- pause;
- figure(3)
- plot(f,plot_Pxx2);
复制代码 |
|