psd via fft and pwelch
2015-10-19 15:34
113 查看
%fft and pwelch方法求取功率谱 load x.mat Fs = 1; t = (0:1/Fs:1-1/Fs).'; Nx = length(x); % Window data w = hanning(Nx); xw = x.*w; % Calculate power nfft = Nx; X = fft(xw,nfft); mx = abs(X).^2; % Normalize by window power. Multiply by 2 (except DC & Nyquist) % to calculate one-sided spectrum. Divide by Fs to calculate % spectral density. mx = mx/(w'*w); NumUniquePts = nfft/2+1; mx = mx(1:NumUniquePts); mx(2:end-1) = mx(2:end-1)*2; Pxx1 = mx/Fs; Fx1 = (0:NumUniquePts-1)*Fs/nfft; [Pxx2,Fx2] = pwelch(x,[],[],[],Fs); plot(1./Fx1,Pxx1,1./Fx2,Pxx2,'r:'); legend('PSD via FFT','PSD v ia pwelch')
相关文章推荐
- 处理机调度
- C++编程注意事项
- iOS navigationBar导航栏底部与self.view的分界线的隐藏
- 中文乱码的原因及解决方法
- 滚轮事件的防冒泡、阻止默认行为的代码(效果是:只让当前div滚动,连当前文档都不滚动的效果)
- 论实践和认识的关系--来自生活体验和感悟
- android 线程初识
- 错误代码:1100 Table 't_depart_info' was not locked with LOCK TABLES的解决方法
- ZigBee使用Z-Stack协议栈时如何关闭ACK验证
- AIX下取前一天时间
- php模拟post提交数据
- Java基础知识学习(九)
- 浅谈3D Touch 在Xamarin.iOS上的应用
- 如何对map中的value进行排序
- Unity3d中BlinnPhong光照模型注解
- mongodb集群搭建-replica set模式
- Android应用程序发送广播(sendBroadcast)的过程分析
- EMV内核使用中的常见问题
- PMAC运动程序例程(二)
- 每天一个linux命令(19):find 命令概览