有一组数据，用FFT作频谱分析，并画出幅频特性、相频特性曲线

myinternet2008

示波器记录了一组数据，有2500个，想用matlab的FFT进行频谱分析，并画出幅频特性、相频特性曲线。
第一次接触matlab，希望各位老师同学帮忙。

ChaChing

养成好习惯, 先请动手搜一下!:loveliness:

myinternet2008

搜了，在各个论坛蹲两天了，楼上的帮帮我。。。

ChaChing

看看这里有没有你要问的问题－－信号处理专栏话题索引
http://forum.vibunion.com/forum/thread-49413-1-1.html

myinternet2008

楼上的，帮我看看我写的这程序对不对？
clear
x=xlsread('D:\sensor.xls');    %数据我放在excel里了，只取了2048个
N=length(x); fs=25000; t=(0:N-1)/fs;
y=fft(x); F=y(1:N/2+1); f=fs*(0:N/2)/N;
subplot(311); plot(t,x); subplot(312);plot(f,abs(F))
subplot(313);plot(f,phase(F)*180/pi)
谢谢改正！

做出来的图，总觉得它们的纵坐标有问题

[ 本帖最后由 ChaChing 于 2010-6-25 21:14 编辑 ]

ChaChing

有什么问题?

myinternet2008

第一图是用origin软件，进行快速傅里叶变化得到的图，第二图是用matlab进行FFT变换得到的图，两幅图的纵坐标不一样，是不是上面编写的matlab程序有问题？
还有一个问题就是：matlab的FFT能处理这个衰减振荡信号吗？

ChaChing

个人水平专业有限, 仅几点意见参考下!
1.仅取2048点是前2048点, 还是后2048点
2.幅值改为F=y(1:N/2+1)*2/N;
3.origin软件所用参数相同吗?
4.基本上, 个人以为fft不适宜处理这个衰减振荡信号, 但不同工具的结果应不至于相差太大

myinternet2008

取的是前2048个点，幅值改为F=y(1:N/2+1)/N会不会好点？
相频特性的纵坐标还是不对，不是在-180~180内变化，需要怎么在程序中修改？
把plot(f,phase(F)*180/pi)   改成  plot(f,angle(F)*180/pi) 好像对了，楼上的帮我分析下对吗？

[ 本帖最后由 ChaChing 于 2010-7-1 00:52 编辑 ]

ChaChing

个人知道的幅值是F=y(1:N/2+1)*2/N;
在-180~180内变化, 的确应该使用angle

个人水平专业有限, 看看高手们的帖
[原创]关于FFT结果为什么要除以N (5F)
http://forum.vibunion.com/forum/viewthread.php?tid=23665
FFT求频谱和相位的问题
http://forum.vibunion.com/forum/viewthread.php?tid=52269

myinternet2008

我把程序改成：
clear
x=xlsread('D:\sensor.xls');    %数据我放在excel里了，只取了前2048个
N=length(x); fs=25000; t=(0:N-1)/fs;
y=fft(x); F=y(1:N/2+1)*2/N; f=fs*(0:N/2)/N;
subplot(311); plot(t,x);
subplot(312); plot(f,abs(F))
subplot(313); plot(f,angle(F)*180/pi)
好像对了，谢谢楼上的ChaChing这几天的帮助，受益匪浅。

ChaChing

这个个人懂的不多, 没帮什么忙, 共同学习吧!

cui7659

LZ还在不，我也有一组时域衰减数据，要画出幅频曲线，你能把你的代码给我注释下吗，有些看不懂，谢谢

stao136

顶，支持，学习学习

stao136

，还不错