求声音的频谱分析算法，就像winamp的柱状显示(2).

[AskLove7]

求声音的频谱分析算法，就像winamp的柱状显示。

大概又叫光谱分析吧，有函数代码最好，估计理论上应该是较复杂的,听说用到个什么傅立叶。
有熟悉的朋友帮帮忙吧。

由于较急，发了两个帖子。
另一个相同的帖子如下：
http://community.csdn.net/Expert/topic/3705/3705908.xml?temp=.683468

[oyljerry]

用快速傅立叶变缓算法

[SpeedChina]

下面的代码不知是否有用：

我这儿也有一个c版的源程序
[代码]FFT 快速傅立叶变换 (c语言)

算法参见徐世良《计算机常用算法》第二版 第294页

void KFFT(
BNU::vector<std::complex<double> > &x,
BNU::vector<std::complex<double> > &y,
int inverse)
{
int n = x.size();
int k = log((double)n)/log(2.);
int t ,m ,s , i, j, NV, r;
double angle;
std::complex<double> V,PODD;
for (t = 0 ; t <= n-1 ; t++) {
m = t;
s = 0;
for (i = 0 ; i <=k-1 ; i++) {
j = int(m/2);
s = 2*s + m - 2*j;
m = j;
}
if (&x!=&y) {
y[t] = x[s];
}else {
if (s>t) {
V=x[t];
y[t]=x[s];
y[s]=V;
}
}
}
for (t = 0 ; t <= n-2 ; t+=2) {
V = y[t];
y[t] = V + y[t+1];
y[t+1] = V - y[t+1];
}
m = n/2;
NV = 2;
for (r = k-2 ; r >=0 ; r-- ) {
m = m/2;
NV = 2*NV;
for (t = 0 ; t <= (m-1)*NV ; t=t+NV ){
for (j = 0 ; j <= (NV/2)-1 ; j++) {
angle = -(2.0*M_PI*inverse*j)/NV;
PODD = complex<double>(cos(angle),sin(angle))*y[t+j+(NV/2)];
y[t+j+(NV/2)] = y[t+j] - PODD;
y[t+j] = y[t+j] + PODD;
}
}
};
for (i = 0 ; i < n; i++) {
if (inverse==-1)
y = y/double(n);
}
}

[neilxp]

傅立叶变换不行吗

[EnochShen]

波形图可不可以?可以的话,发邮件给我,我有这么个例子 [e-mail]

[Kudeet]

http://www.codeproject.com/audio/waveInFFT.asp
http://www.codeguru.com/Cpp/G-M/multimedia/audio/article.php/c4709/

[hushuangyan74]

应该要用到ＤＦＴ，可以用ＤＣＴ来加快算法速度，不过没有具体搞过，帮不了忙！