h****y
发帖数: 61 | 1 谢谢～～
由于偶的程序中就只要一次FFT，所花时间也不多，最后偶决定还是用最土的办法，直
接算三重积分得到fourier transform coeffients. 对比用FFT得到的结果，发现两者
还是怎么都对不上，除了常数项外。
对于一个有inversion symmetry的周期性体系的势，做fourier transform 后，起
fourier transform coeffients的虚部应该为零，可 FFT后给出的结果则不是这样的。
偶怀疑到底FFT后直接给出的是不是就是fourier transform coeffients啊？这个FFT变
换后各个fourier transform coeffients究竟是怎么布局的？
一维的FFT情况已经清楚了，我的理解是三维的FFT不就在每一维度上先后做一次一维的
FFT吗？我的程序是用Fortran写的，这也是为什么最初选择用FFTPACK,而没有用FFTW的
原因。对3D的情况，FFTW中可以容易的找出如n1=-2，n2= 2, n3=0 对应的系数吗？ |
|
r****y
发帖数: 1437 | 2 【 以下文字转载自 Science 讨论区 】
【 原文由 rossby 所发表 】
个人体会，也许对上面的那个问题有用。用过几次fft,发现一不小心还真爱出
错，估计是我水平太不照了。
1。 for real function f(x), ifft(fft(f(x))) != f(x), only
real(ifft(fft(f(x)))) = f(x), 计算中的round off errors很容易
产生很小的虚数项，~1e-14 to 1e-12.
1。 for f(x), ifft(fft(f(x)))= f(x), 这个没问题。但是要是你想通过
fft做微分，就得小心。比如用matlab, suppose你的数据是128个点，
y = fft(f, 128)
所得到的y还是128个值，前65个对应wavenumber 0, 1, 64, 后面的对应
-63, -62, ..., -1. 你要是以为后面的对应65, 66, 127，做微分的就会出错， |
|
r****y
发帖数: 1437 | 3 个人体会，也许对上面的那个问题有用。用过几次fft,发现一不小心还真爱出
错，估计是我水平太不照了。
1。 for real function f(x), ifft(fft(f(x))) != f(x), only
real(ifft(fft(f(x)))) = f(x), 计算中的round off errors很容易
产生很小的虚数项，~1e-14 to 1e-12.
1。 for f(x), ifft(fft(f(x)))= f(x), 这个没问题。但是要是你想通过
fft做微分，就得小心。比如用matlab, suppose你的数据是128个点，
y = fft(f, 128)
所得到的y还是128个值，前65个对应wavenumber 0, 1, 64, 后面的对应
-63, -62, ..., -1. 你要是以为后面的对应65, 66, 127，做微分的就会出错，
因为 d(exp(-ikx))/dx= -ik exp(-ikx), 这个k你就错了 |
|
j********i
发帖数: 2412 | 4 投诉人：jinxiaomei
投诉理由：跨版封人
投诉目标：解封道歉版三下台
我在biology版跟IVYtony言语不合，这厮恼羞成怒，居然在overseas以人身攻击为由封
我及我的另一个id xiaomai 14天。我以及xiaomai已经很长时间没去overseas灌水了。
他不用IVYtony身份，借了版主fft的id封我。fft已经自杀变成新人（证据如下），将密码告诉了IVYtony.
fft和IVYtony IP相同(证据如下). 既然fft无意做版主，那就退位让贤吧。至于IVYtony这样的版三，站务应该立即解除他的版主职务。谢谢！
证据：
biology版讨论在此：http://www.mitbbs.com/article_t/Biology/28417842.html
发信人: deliver (自动发信系统), 信区:
标 题: fft 封 jinxiaomei 在 Overseas
发信站: BBS 未名空间站自动发信系统 (Mon Nov 12 01:07:21 2007)
【此篇文章是由自动发信系统所张贴】
由于 jinxiaomei 在 Overseas |
|
|
m*******s
发帖数: 3142 | 6 【 以下文字转载自 Mathematics 讨论区 】
发信人: manifolds (流形), 信区: Mathematics
标 题: 有没有区间[0，pi]上的FFT算法？
发信站: BBS 未名空间站 (Sat Jun 25 23:17:54 2011, 美东)
现在有一个数值积分问题，需要得到如下积分的虚部
int_0^\pi dx f(cosx)*sinx*e^{imx}/(i*sinx-a)
其中m是自然数，a是非0实数.
我的做法是把被积函数的虚部和实部分离出来，分别积分，可以用FFT的现成程序完成。
但是我觉得这个做法比较笨，我觉得应该可以只用一次类似FFT的计算，然后直接取结
果的虚部。
我所了解的FFT应该是对[0，2pi]这个区间积分，问题是我这个被积函数非基非偶，不
知道该如何是好。
请问大家有没有什么建议? 谢谢！ |
|
t*******r
发帖数: 22634 | 7 这个我在跟人讨论的时候，正好做了个 FFT Spectrogram / Spectrum 的音准分析。
FFT Spectrogram / Spectrum 的图形见附件。
（插一句，所有这些图形都是用免费开源软件 Audacity 生成的，如何生成，我
会在后面跟帖说明）。
如何解读这个图，需要稍稍解释一下。
（1）首先是区分人声和伴奏：
每一个单独的乐器和人声，在 spectrogram 上，基音和各阶上方泛音的图形，
会形成独特重复的 pattern。观察这些 pattern 就可以区分人声和伴奏
（其实可以精确到每个声部里单独的乐器）。
在这个图里，我选中 1:25 - 1:26 之间，在唱 “never be” 的那个 “be” 字。
观察 spectrogram 可知 “be” 字的基音大约是在 300 Hz 左右。
（2）然后是 plot spectrum 测量基音能量峰值的频率：
对选中的那段 plot FFT spectrum。在 spectrum 窗口，把鼠标移动到 300 Hz
附近，spectrum 窗口显示能量顶峰是在 297 Hz，音高是属于 D4 的范围内... 阅读全帖 |
|
t*******r
发帖数: 22634 | 8 现在就是如何解释这个 +3 Hz / +20 cent 的 FFT 数据的问题。。。
俗话说的好，不怕不识货，就怕货比货。。。Paul Anka 前辈，这次就
让您老当小白鼠了。。。俺用 FFT 测了一下 Paul Anka 的。。。图形
贴在附件里。
Paul Anka 的那个 “never be” 的 “be” 音唱的比较直，不那么 vibrato。
（总体上他的风格 vibrato 比较少）。但也有麻烦，他唱的那个音是用
略略向上的 pitch bend 的。（插一句，pitch bend 不属于走音，是
vocal nuance 的一种）。。。所以这个如何是好嘛。
没关系，很多 pitch bend 还是有一定的准则的。附件图里的那种微微向上的
pitch bend，一般而言是以结尾点的音高为准。。。所以俺就选了结尾点最高
的的一小段，plot spectrum 一测，能量峰值是 288 Hz。。。相对 294 Hz
的标准音高，数据上是唱低了 6 Hz。
我擦。。。人唱低 6 Hz （FFT 能量峰值数据而言）都发唱片了，都赚美刀
了。。。俺唱高 3 Hz 算个球嘛。。。hi... 阅读全帖 |
|
e***e
发帖数: 3872 | 9 【 以下文字转载自 Graphics 讨论区 】
发信人: etude (小铃铛), 信区: Graphics
标 题: CUDA FFT
发信站: BBS 未名空间站 (Sat Apr 14 11:50:57 2007), 转信
比较了一下CUDA FFT Lib和Intel Math Kernel的性能，计算2D FFT，single
precision float complex to complex, 硬件分别是8800GTS和P4 3.0G(w/ HT)
显存640M和内存1G，下表中时间单位为毫秒：
size： 256 * 256 512 * 512 1024 *1024 2048 * 2048
CUDA 9.093843 31.65965 96.12374 603.8209
9.508812 31.58298 96.13728 609.9658
8.999843 31.65233 9 |
|
d*******2
发帖数: 340 | 10 碰到一个fft的奇怪问题.一个384x384的矩阵,fft后有很多噪音.对矩阵加零变成
512x512的矩阵后,fft反而没有噪音了.按理说加零揭示出频谱的更多细节,应该有更多
噪音才对啊.
程序如下,麻烦哪位前辈帮忙看一下
f384=abs(fftshift(fft2(fftshift(matrix384))));
GN=512;
N=384;
matrix512=zerso(GN);
matrix512(GN/2-N/2+1:GN/2+N/2,GN/2-N/2+1:GN/2+N/2)=matrix384;
f512=abs(fftshift(fft2(fftshift(matrix512))));
f512比f384平滑得多?!
先谢了! |
|
q********w
发帖数: 6 | 11 Yes! It's really a problem to measure the accurate amplitude
of a SIN signal if the sampling window is not exactly
the mutiple of period of the SIN,because FFT is just a
sample of FT.to solve the problem,one of the effective way
is to get N sets of the samples,say,A,B,C....,then average
the power of these samples.
Equation:((abs(FFT(A))^2+(abs(FFT(B))^2+...)/N
and of course another way,increase the sampling frequency
as high as possible. |
|
r*******g
发帖数: 1335 | 12 Hello
我想请教一个FFT的问题，我也是EE的，但是现在做的东西比较偏向cs，做的是network
traffic分
析。
给定一个时间序列，有无办法求它的各个频段上的能量，我想通过能量来classify
network traffic.
这是个很open的问题，我自己也不清楚能量为何物，是不是就是CFT变换后频域上波形
的面积。比如说我
可以给定截止频率，想求某两个频率之间的能量。
FFT的结果不知道可否来做这个事情，我记得时间采样则频域上周期扩展，频域上采样
则时间域周期扩展，
所以FFT的几个值好像是频域上采样啊，这个和我要求的能量有什么关系呢？
谢谢了。 |
|
a****l
发帖数: 8211 | 13 1. I don't think Fourier Transform is designed for analyzing network traffic
. They are simply not relevant.
2. If you really want to do it, I think you mean you want to get the power
spectrum, or Power Spectrum Density. PSD is not part of FFT, but you can use
FFT to calculate traditional PSD. There are lots of other definitions of
PSD, which shall be calculated differently.
3. Talking about "power of signal", it is not "the area under fft". Do some
basic math of a sine signal, you can see why.... 阅读全帖 |
|
m*******s
发帖数: 3142 | 14 现在有一个数值积分问题，需要得到如下积分的虚部
int_0^\pi dx f(cosx)*sinx*e^{imx}/(i*sinx-a)
其中m是自然数，a是非0实数.
我的做法是把被积函数的虚部和实部分离出来，分别积分，可以用FFT的现成程序完成。
但是我觉得这个做法比较笨，我觉得应该可以只用一次类似FFT的计算，然后直接取结
果的虚部。
我所了解的FFT应该是对[0，2pi]这个区间积分，问题是我这个被积函数非基非偶，不
知道该如何是好。
请问大家有没有什么建议? 谢谢！ |
|
e***e
发帖数: 3872 | 15 【 以下文字转载自 Graphics 讨论区 】
发信人: etude (小铃铛), 信区: Graphics
标 题: CUDA FFT
发信站: BBS 未名空间站 (Sat Apr 14 11:50:57 2007), 转信
比较了一下CUDA FFT Lib和Intel Math Kernel的性能，计算2D FFT，single
precision float complex to complex, 硬件分别是8800GTS和P4 3.0G(w/ HT)
显存640M和内存1G，下表中时间单位为毫秒：
size： 256 * 256 512 * 512 1024 *1024 2048 * 2048
CUDA 9.093843 31.65965 96.12374 603.8209
9.508812 31.58298 96.13728 609.9658
8.999843 31.65233 9 |
|
t*******r
发帖数: 22634 | 16 最近讨论踊跃，俺就试着结合提笑肌和 FFT 试唱一下
you raise me up 一小段，反正就是有高音的那小段，
其他省了。没学，直接跟唱的，音准节奏很烂。后期
2:1 压缩以及加混响。
感想：提笑肌容易声音有表情，但是也容易紧张。FFT
有助于放松，但是缺少表情。要两头都能顾及，目前俺
搞不定。 |
|
t*******r
发帖数: 22634 | 17 这个坑爹的 FFT spectrogram，很多情况还是不如耳朵好使。
不过耳朵的缺点是受心理影响比较大，不客观。FFT spectrogram
当个提醒、以及客观比较，还是很不错的。 |
|
t*******r
发帖数: 22634 | 18 这个我在跟人讨论的时候，正好做了个 FFT Spectrogram / Spectrum 的音准分析。
FFT Spectrogram / Spectrum 的图形见附件。
（插一句，所有这些图形都是用免费开源软件 Audacity 生成的，如何生成，我
会在后面跟帖说明）。
如何解读这个图，需要稍稍解释一下。
首先是区分人声和伴奏。每一个单独的乐器和人声，在 spectrogram 上，
基音和各阶上方泛音的图形，会形成独特重复的 pattern。观察这些
pattern 就可以区分人声和伴奏（其实可以精确到每个声部里单独的乐器）。
在这个图里，我选中 Paul Anka 在唱 “never be” 的那个 “be” 字。
观察 spectrogram 可知 “be” 字的基音大约是在 300 Hz 左右。
Paul Anka 的那个 “never be” 的 “be” 音唱的比较直，不那么 vibrato。
（总体上他的风格 vibrato 比较少）。但也有麻烦，他唱的那个音是用
略略向上的 pitch bend 的。（插一句，pitch bend 不属于走音，是
vocal nu... 阅读全帖 |
|
S**I
发帖数: 15689 | 19 这个本质上没法并行啊：后一个FFT依赖于前一个FFT的结果...... |
|
a*********t
发帖数: 261 | 20 AMD X2 560 unlock 到 4cores后，prime95 small FFT 可以运行一个晚上没有问题，
但large FFT 在一个小时左右就会从新启动机子，求建议。多谢。 |
|
s*******k
发帖数: 71 | 21
Mathematical routines such as FFT shouldn't be machine dependent unless it
uses hardware accelerations (for example, Intel MMX instructions), so what
you have for SGI should be able to compile and run on any machine. If not,
try www.fftw.org, they have a free FFT library which works for almost all
platforms (well, maybe not PalmPC!)
SH |
|
c******o
发帖数: 8 | 22 I need to install a Fast Fourier Transform subroutine in dec
alpha machine. comlib.sgimath has FFT subroutine, but it is
for sgi machine. Can anyone recommend any web site to
download a FFT subroutine for alpha machine? Thank you. |
|
s*****k
发帖数: 3 | 23 i am trying to write code for fft myself. of cause i think it is ok if it can
run. someone can tell me some references on fft algorithm.
thanks! |
|
a***n
发帖数: 328 | 24 对于有n个点的信号，fft的时候fft size该取多少有什么基本原则么？
如果2^m〉n〉2^m-1，用哪个好？
另外，用matlab里的pwelch分析信号的时候，默认做法是把信号分成重叠50%的八段信
号。但是如果信号本身比较短的话，是不是应该/可以少分几段，有没什么原则可以
follow？
外行提问，请多指教！谢谢！ |
|
a***n
发帖数: 328 | 25 非常感谢！
分的段数少不正好增加输入每个fft的数据量么？这样应该是增加精度啊，不是么？
如果程序并不要求fft size是2^x，且不考虑程序运算效率，选用n是不是可以提供更准
确的spectrum? |
|
s******e
发帖数: 841 | 26 我有如下50个data point，我想用FFT找出data中的频率是多少。用了Matlab的fft函数
，但是magnitude最大的总是第一个数，那对应的频率应该是什么呢，还是我弄错了，
谢谢 |
|
d*****l
发帖数: 8441 | 27 怪事。你的50个数据点是实数还是复数点啊？
八成是复数的吧。不然的话，如果是实数，做FFT得到的频谱幅度是对称的啊，可是
你的图上看不出对称啊。若是实数矢量的话，50点FFT说实话只能得到25个bin，
除非内插，或用其他基于模型的谱估计方法。 |
|
f***r
发帖数: 44 | 28 用fftshift(fft(....))就得到你平时看的频谱了
变换回去用ifft(ifftshift(....))
ifft(ifftshift(fftshift(fft(anything))))==anything |
|
j********i
发帖数: 2412 | 29 请站务处理一下这个，我认为fft已经自杀戒网，所以变成新人，密码告诉了IVYtony.
IVYtony借fft跨版封人，企图混淆是非，情节十分恶劣。 |
|
m****s
发帖数: 18160 | 30 ☆─────────────────────────────────────☆
TrueCabbage (@圆白菜@) 于 (Wed Jan 24 11:53:19 2007) 提到:
发信人: agostic (灌水有害健康), 信区: Complain
标 题: 弹劾Overseas的版主：fft
发信站: BBS 未名空间站 (Tue Aug 1 17:14:00 2006)
投诉人：agostic
－－－－打抱不平而已！！
投诉对象：fft
投诉标题：不整理版务，不分青红皂白删除所有帖子（几千个一起删除）
投诉目标：弹劾板主
投诉理由：
1. 理由一：不顾班规，不给警告之前就删除帖子
证据：无论什么帖子都已经被删除
2. 理由二：对于违反班规者，没有警告，也没有处理，采用删除所有帖子的办法
证据二：帖子都消失了。很多求助贴，意见贴，信息贴，帮助贴，和平贴，中性贴，等等都被删除！！！
投诉日期：08/01/2006
这么懒的版主，只会一刀切。把所有的讨论帖子，无论好坏，按帖子号，不按内容全 |
|
e***e
发帖数: 3872 | 31 比较了一下CUDA FFT Lib和Intel Math Kernel的性能，计算2D FFT，single
precision float complex to complex, 硬件分别是8800GTS和P4 3.0G(w/ HT)
显存640M和内存1G，下表中时间单位为毫秒：
size： 256 * 256 512 * 512 1024 *1024 2048 * 2048
CUDA 9.093843 31.65965 96.12374 603.8209
9.508812 31.58298 96.13728 609.9658
8.999843 31.65233 96.37738 604.2612
INTEL 307.4331 1228.842 5498.506 22838.89
294.6351 1229 |
|
|
t*******r
发帖数: 22634 | 33 当然，这个 “捷克式轻机枪 FFT 测音准” 嘛，虽然火力是猛烈，但是毕竟
轻机枪有点小重是不是，还有点小费子弹是不是。。。尼玛这次累坏俺了。。。
所以俺主要是忽悠大伙儿，俺自己一般不测，凑合听听就完事了。。。当然，
如果遇上网络对钢，那捎上 88 毫米火炮也不算重是不是？。。。sorry，
跑题了。。。:-P
下面言归正传，如何用免费且开源的 Audacity 生成这些
spectrogram / spectrum：
（1）首先在 edit -> preference -> spectrogram 里设参数。具体建议的
参数见附图。
（2）然后在 waveform 左边的名字下拉菜单里，把 “Waveform” 改成
“Spectrogram log(f)”。也见附图。
（3）如果要 plot spectrum，选中一段后，在菜单里选
“Analyze -> Plot Spectrum”，然后用鼠标在 spectrum 上移动，
自动会显示鼠标附近的峰值。
先到这里，大家有意见可以畅所欲言、自由发话嘛。。。 |
|
t*******r
发帖数: 22634 | 34 当然，这个 “捷克式轻机枪 FFT 测音准” 嘛，虽然火力是猛烈，但是毕竟
轻机枪有点小重是不是，还有点小费子弹是不是。。。尼玛这次累坏俺了。。。
所以俺主要是忽悠大伙儿，俺自己一般不测，凑合听听就完事了。。。当然，
如果遇上网络对钢，那捎上 88 毫米火炮也不算重是不是？。。。sorry，
跑题了。。。:-P
下面言归正传，如何用免费且开源的 Audacity 生成这些
spectrogram / spectrum：
（1）首先在 edit -> preference -> spectrogram 里设参数。具体建议的
参数见附图。
（2）然后在 waveform 左边的名字下拉菜单里，把 “Waveform” 改成
“Spectrogram log(f)”。也见附图。
（3）如果要 plot spectrum，选中一段后，在菜单里选
“Analyze -> Plot Spectrum”，然后用鼠标在 spectrum 上移动，
自动会显示鼠标附近的峰值。
先到这里，大家有意见可以畅所欲言、自由发话嘛。。。 |
|
S**I
发帖数: 15689 | 35 应该可以吧，看这个API：
cufftResult
cufftPlanMany(cufftHandle *plan, int rank, int *n, int *inembed,
int istride, int idist, int *onembed, int ostride,
int odist, cufftType type, int batch);
Creates a FFT plan configuration of dimension rank, with sizes specified in
the array n. The batch input parameter tells CUFFT how many transforms to
configure. With this function, batched plans of 1, 2, or 3 dimensions may be
created.
... cufftPlanMany() allows use of more complicated data layouts... 阅读全帖 |
|
|
b**u
发帖数: 17 | 37 已知二个函数f,g,f有紧支集，g属于L2[0,1],如何用FFT算法求卷积F(x)在[0。1]的值。
怎么把F(x)拓展成周期函数？
3x |
|
t*****o
发帖数: 74 | 38 why bother to write it by youself?
there are many very good standard fft routines
very fast and stable
already developed tens year by supercomputer center |
|
s*****k
发帖数: 3 | 39 I ever tried two different formulas of Fourier transform. One worked well, but
the other didn't. To know them more clearly, I try to write the code myself.
I don't expect that it will be nice to work on every machine or in any code,
but I hope it works well when I test it.
From www.google.com, it seems not easy to find some useful algorithm except
the private codes. Here somebody can show me the FFT algorithm.
Thanks!
it |
|
|
r****y
发帖数: 1437 | 41
nod nod, strongly recommend fftw.
I heard the professor who taugh me signal processing said that,
at 1960s, the guys who invented fft can actually patent it. But they did not.
If they did, they should be much richer than Bill Gates now.
|
|
h****y
发帖数: 61 | 42 版上有人用过那个FFTPACK（by Hugh C. Pumphrey）？
这里面FFT变换后得到的fourier transform的系数究竟是怎么对应的？
看了doc 文件，并和解析结果比较，对1D 的情况好像是这样的： 第一个数为fourier
transform 的常数项,即n=0的项；第2个数对应n=1项的实部，第3个数对应n=1项的虚部
，第3个数对应n=2的实部，第四个数对应n=2的虚部 ...
可是对于2D、3D的情况死活对应不了 :-( |
|
D********g
发帖数: 533 | 43 MatLab里面应该有吧。。。 我记得是有2维FFT的。。。 |
|
H*****n
发帖数: 95 | 44 大家能推荐一下好的C 的 FFT 和FFT2 吗？
网上查到很多版本，不知道哪个好，希望高手指点。
谢谢 |
|
z*****n
发帖数: 7639 | 45 the first FFT result is the DC component. It is simply the sum
of all the x(n)'s. |
|
n****n
发帖数: 8 | 46 除了直流问题，你还要解决频谱泄漏问题：
如果你的FFT不满足相干采样（coherent sampling）出来的频谱会有泄漏，你的目标频
率分量会落在2个bin的中间，需要用tabei-ueda算法进行估计，而且是无偏估计 |
|
d*****l
发帖数: 8441 | 47 也不说清楚是原始数据图，还是FFT后的图。或许是数据吧。 |
|
|
r*******g
发帖数: 1335 | 49
hi，能不能具体的讲讲。
这个window filter是时域还是频域上的，fft的各点的值是怎么和能量对应上的。我觉
得一个f对应的
值最多是omiga=2pi f那一个omiga处对应的值，但是我想求的是一段连续的omiga对应
的能量。
谢谢了。 |
|
T****i
发帖数: 2190 | 50 DM有inverse FFT的功能，但是没有做过直接把EDP直接变成image的
有专门做从衍射回复图像的，而且似乎不是很简单，UIUC的左建民在做这个
Science 300, 1419 (2003) |
|
