|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 wdhd 于 2016-9-14 10:33 编辑
因为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);
复制代码 |
|