| 查看: 1541 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助傅里叶变换对不对？ 已有1人参与
|
||
|
我得到了一组数据，第一列是时间，单位是ps（1.0e-12s)，第二列是速度自相关函数。 0 1 10 0.04994 20 0.004652 30 0.01519 40 0.06579 50 0.01538 60 0.02599 70 0.0297 80 0.04439 90 0.02218 100 0.01489 110 0.0004905 120 0.02008 130 0.009628 140 0.02908 150 -0.001292 用matlab进行FFT变换： xlsread VACF.xls y=fft(ans) N=16 f=(0:N-1)*1.0E+11/N mag=abs(y) plot(fs,mag) 得到的图像感觉不对，请教各位哪里出了问题？谢谢~~  FFT.jpg |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助Matlab的fft2与ifft2
已经有7人回复
- 连续函数和不连续的函数的傅里叶变换的主要差别是什么？
已经有6人回复
- 怎么求传递函数
已经有8人回复
- origin傅里叶变换求助
已经有3人回复
- 帮忙看下我的化简对不对？
已经有7人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 请教Laplace变换问题
已经有8人回复
- 求助：谁用过origin傅里叶变换图（FFT)
已经有7人回复
- 对傅里叶变换的理解…
已经有22人回复
- 傅立叶变换的物理意义
已经有9人回复
- 怎么消除虚频，细节
已经有10人回复
- 求膜的弹性模量，用傅里叶变换请问这样做对吗？
已经有5人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- 求助：谁会傅里叶变换图
已经有8人回复
- 【求助】红外光谱解析问题？不过无法确定具体是哪一种，希望能得到帮助。
已经有4人回复
- 【求助】紧急求助：三角函数的傅里叶变换
已经有14人回复
- 【求助】如何做炒蚬子？？看看我的方法对不对啊。
已经有9人回复
- 【求助】原位傅里叶变换红外光谱仪
已经有5人回复
- 求助岛津傅里叶变换红外光谱仪FTIR-8400S型 的使用方法
已经有7人回复
3楼2014-05-05 15:54:54
2楼2014-05-05 15:51:46
4楼2014-05-05 16:55:37
|
我把程序改了，但感觉还是不对，请帮我看看哪里还有问题？谢谢 其中VACF.xls里只有一列16个数据。 1 0.01273 0.01196 0.005036 0.005826 -0.008611 0.003049 0.004471 -0.0008861 6.05E-03 0.006342 0.009512 0.01197 -0.005031 -0.01906 0.00776 0.02477 -0.02541 0.008003 0.05522 -0.08248 修改后的程序为： xlsread VACF.xls Fs=1.0E+14 N=16 N=pow2(nextpow2(N)) %ans=ans-mean(ans) y=fft(ans,N) f=(1:N)*Fs/N f=f(1:N/2) mag=abs(y) mag=mag(1:N/2) plot(f,mag) grid xlabel('Frequency (Hz)') ylabel('Mag') title('{\bf FFT}') saveas(gcf,'FFT','png') P=y.*conj(y)/N P=P(1:N/2) plot(f,P) % plot(f,10*log10(P)) grid xlabel('Frequency (Hz)') ylabel('Power(dB)') title('{\bf Periodogram}') saveas(gcf,'Periodogram','png')  FFT.png  Periodogram.png |
5楼2014-05-05 17:00:29
100