http://blog.sciencenet.cn/blog-2777453-1000182.html
2017-05-06 15:46
232 查看
%--------------------------------------------------------------------------
%功能: 谱分析代码
% 何成文 2016年9月1日
%--------------------------------------------------------------------------
clc;clear;close all;
load('WPM10000.txt');
data=WPM10000;
N=length(data);
figure;plot([1:N],data,'r');title('原始数据');grid on;xlabel('T/s');
%--------------------------------------------------------------------------
F=fft(data,N);
R=abs(F)/N;
Angle=angle(F);
n=[0:1:N-1];
Hz=n/(N*1);
figure
subplot(1,2,1);plot(Hz(1:N/2),R(1:N/2),'r');title('傅里叶幅度谱');xlabel('X/Hz');ylabel('幅度');grid on;
subplot(1,2,2);plot(Hz(1:N/2),Angle(1:N/2),'b');title('傅里叶相位谱');xlabel('X/Hz');ylabel('相位');grid on;
%--------------------------------------------------------------------------
Y=fft(data,N); %傅里叶变换!
Y1=(abs(Y)).^2/N; %Y1为周期图谱计算公式!
n=[0:1:N-1];
Hz=n/(N*1); %1是抽样时间间隔,本文默认为1!
figure;plot(Hz(1:N/2),Y1(1:N/2),'b'); %表现形式1
xlabel('X/Hz');ylabel('幅度');title('WPM噪声信号的傅里叶周期谱图');grid on;
figure;plot(log(Hz(1:N/2)),log(Y1(1:N/2)),'r'); %表现形式2-取对数!
xlabel('log(f/Fs)');ylabel('log(PSD)');title('WPM噪声信号的傅里叶周期谱图');grid on;
%--------------------------------------------------------------------------
来源:http://blog.sciencenet.cn/blog-2777453-1000182.html
%功能: 谱分析代码
% 何成文 2016年9月1日
%--------------------------------------------------------------------------
clc;clear;close all;
load('WPM10000.txt');
data=WPM10000;
N=length(data);
figure;plot([1:N],data,'r');title('原始数据');grid on;xlabel('T/s');
%--------------------------------------------------------------------------
F=fft(data,N);
R=abs(F)/N;
Angle=angle(F);
n=[0:1:N-1];
Hz=n/(N*1);
figure
subplot(1,2,1);plot(Hz(1:N/2),R(1:N/2),'r');title('傅里叶幅度谱');xlabel('X/Hz');ylabel('幅度');grid on;
subplot(1,2,2);plot(Hz(1:N/2),Angle(1:N/2),'b');title('傅里叶相位谱');xlabel('X/Hz');ylabel('相位');grid on;
%--------------------------------------------------------------------------
Y=fft(data,N); %傅里叶变换!
Y1=(abs(Y)).^2/N; %Y1为周期图谱计算公式!
n=[0:1:N-1];
Hz=n/(N*1); %1是抽样时间间隔,本文默认为1!
figure;plot(Hz(1:N/2),Y1(1:N/2),'b'); %表现形式1
xlabel('X/Hz');ylabel('幅度');title('WPM噪声信号的傅里叶周期谱图');grid on;
figure;plot(log(Hz(1:N/2)),log(Y1(1:N/2)),'r'); %表现形式2-取对数!
xlabel('log(f/Fs)');ylabel('log(PSD)');title('WPM噪声信号的傅里叶周期谱图');grid on;
%--------------------------------------------------------------------------
来源:http://blog.sciencenet.cn/blog-2777453-1000182.html
相关文章推荐
- http://blog.sciencenet.cn/blog-630522-491356.html
- http://blog.sciencenet.cn/blog-557238-489823.html
- swat主流域文件(file.cio)参数详解——引自http://blog.sciencenet.cn/blog-922140-710636.html
- 本文引用地址:http://blog.sciencenet.cn/blog-816706-709171.html 此文来自科学网靳远博客,转载请注明出处。
- 矩阵按列归一化(转载http://blog.sciencenet.cn/blog-86695-785606.html)
- http://blog.sciencenet.cn/blog-362400-629730.html
- 微分流形与黎曼几何学习笔记(转自http://blog.sciencenet.cn/home.php?mod=space&uid=81613&do=blog&id=333317)
- ASP.Net设计时需要考虑的性能优化问题 转载自http://blog.sina.com.cn/s/blog_3d7bed650100055p.html
- 转自http://blog.sciencenet.cn/u/施一公 《学生如何提高专业英文阅读能力 精选》
- 整理牛人看文献的方法 http://bbs.sciencenet.cn/thread-523775-1-1.html
- 闭包(Closures)在各种语言中的例子 http://blog.csdn.net/ruby_cn/archive/2004/11/23/192588.aspx
- Eclipse Rcp系列 http://www.blogjava.net/dreamstone/archive/2007/02/08/98706.html
- [ASP.NET] Session 一些比较详细的知识(转自:http://blog.csdn.net/zhoufoxcn/archive/2006/11/08/1373685.aspx)
- java网络编程记一下 http://www.jspcn.net/htmlnews/B20070807102121.html
- http://www.blogjava.net/beansoft/archive/2007/03/09/102812.html
- 停止更新.全部转到http://blog.fedora.net.cn
- asp.net中验证码的实现 (转自http://blog.csdn.net/zhoufoxcn/archive/2007/07/19/1698615.aspx)
- Tooltip效果 http://cceer.xmu.edu.cn/blog/post/tooltip.html
- Webdings字体图案 http://www.blogjava.net/majianan/archive/2007/05/04/88726.html
- Eclipse_Wtp_HTML(http://hongjunli.blogjava.net/)