| 查看: 1345 | 回复: 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 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有3人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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人回复
2楼2014-05-05 15:51:46
3楼2014-05-05 15:54:54
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
6楼2014-05-06 23:41:09