16,472
社区成员
发帖
与我相关
我的任务
分享求声音的频谱分析算法，就像winamp的柱状显示。
求声音的频谱分析算法，就像winamp的柱状显示。
大概又叫光谱分析吧，有函数代码最好，估计理论上应该是较复杂的,听说用到个什么傅立叶。
有熟悉的朋友帮帮忙吧。
大概又叫光谱分析吧，有函数代码最好，估计理论上应该是较复杂的,听说用到个什么傅立叶。
有熟悉的朋友帮帮忙吧。
...全文
请发表友善的回复…
发表回复
Kudeet 2005-01-08
- 打赏
- 举报
http://www.codeproject.com/audio/waveInFFT.asp
EnochShen 2005-01-08
- 打赏
- 举报
波形图可不可以?可以的话,发邮件给我,我有这么个例子 EnochShen@QQ.com
oyljerry 2005-01-07
- 打赏
- 举报
算法就是用FFT
SpeedChina 2005-01-07
- 打赏
- 举报
不是是否有用：
我这儿也有一个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);
}
}
我这儿也有一个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);
}
}
zhangqu_980371 2005-01-07
- 打赏
- 举报
希望那位高人有就给出来看看啊。我和楼主一样也很想知道。
chuanke 2005-01-07
- 打赏
- 举报
关注
oldboy1234 2005-01-07
- 打赏
- 举报
用傅立叶变换啊，求出波形的频率响应，winamp的柱状图的轮廓就是频响曲线，FFT可以用matlab算
saliors 2005-01-07
- 打赏
- 举报
光谱分析与频谱分析没有关系的啊。。。
msm 2005-01-07
- 打赏
- 举报
关注！学习！
