已有程序将时域信号(*.wav文件和pcm文件)使用傅立叶算法把时域信号变成
频域信号了,但怎么处理频域信号呢?
我要得到文件中的最大频率及该频率的出现的次数; 还有静音的多久,是单频还是双频?对傅里叶已经忘记了,希望做过这方面的同仁能告知,谢谢!附:傅里叶变换程序:回复人:ytandwy(云天) ( ) 信誉:100 2003-7-8 23:37:54 得分:50
?
/**************************************************************************
* 文件名:FreTrans.cpp
*
* 正交变换API函数库:
*
* FFT() - 快速付立叶变换
* IFFT() - 快速付立叶反变换
* DCT() - 离散余弦变换
* WALSH() - 沃尔什-哈达玛变换
*
* Fourier() - 图像的付立叶变换
* DIBDct() - 图像的离散余弦变换
* DIBWalsh() - 图像的沃尔什-哈达玛变换
*
*************************************************************************/#include "FreTrans.h"
#include "stdafx.h"
#include "DIBAPI.h"#include <math.h>
#include <direct.h>
#include <complex>
using namespace std;// 常数π
#define PI 3.1415926535/*************************************************************************
*
* 函数名称:
* FFT()
*
* 参数:
* complex<double> * TD - 指向时域数组的指针
* complex<double> * FD - 指向频域数组的指针
* r -2的幂数,即迭代次数
*
* 返回值:
* 无。
*
* 说明:
* 该函数用来实现快速付立叶变换。
*
************************************************************************/
VOID WINAPI FFT(complex<double> * TD, complex<double> * FD, int r)
{
// 付立叶变换点数
LONG count;// 循环变量
int i,j,k;// 中间变量
int bfsize,p;// 角度
double angle;complex<double> *W,*X1,*X2,*X;// 计算付立叶变换点数
count = 1 << r;// 分配运算所需存储器
W = new complex<double>[count / 2];
X1 = new complex<double>[count];
X2 = new complex<double>[count];// 计算加权系数
for(i = 0; i < count / 2; i++)
{
angle = -i * PI * 2 / count;
W[i] = complex<double> (cos(angle), sin(angle));
}// 将时域点写入X1
memcpy(X1, TD, sizeof(complex<double>) * count);// 采用蝶形算法进行快速付立叶变换
for(k = 0; k < r; k++)
{
for(j = 0; j < 1 << k; j++)
{
bfsize = 1 << (r-k);
for(i = 0; i < bfsize / 2; i++)
{
p = j * bfsize;
X2[i + p] = X1[i + p] + X1[i + p + bfsize / 2];
X2[i + p + bfsize / 2] = (X1[i + p] - X1[i + p + bfsize / 2]) * W[i * (1<<k)];
}
}
X = X1;
X1 = X2;
X2 = X;
}// 重新排序
for(j = 0; j < count; j++)
{
p = 0;
for(i = 0; i < r; i++)
{
if (j&(1<<i))
{
p+=1<<(r-i-1);
}
}
FD[j]=X1[p];
}// 释放内存
delete W;
delete X1;
delete X2;
}
* complex<double> * FD - 指向频域数组的指针
现在的问题是如何使用FD中的数据进行其他处理,如找出最大频率,次最大频率,静音时间,是单频还是双频等?

解决方案 »

  1.   


    :Ryan1005(从左向右走)代码是从网上拷贝的,可能就是您整理的。
    :qinxiaogang(qinxg) 
       俺确实不知道傅里叶变换为何物了,只是隐约知道在高等数学里学过,
    具体是什么早就忘记了, 你说的复数的模, 我也快忘记了,不过你一说我也回忆起了一些。。
       你有这方面的资料吗,能否让参考一下。
       
       [email protected]  or [email protected]
       
      

  2.   

    把时域里的东西转换到频域里主要是看它的频率分布,然后做一点滤波处理,比如低通滤波,带通滤波等。滤去你不想要的频率分量,保留你所需要的频率成分。这个和电路上的滤波器原理差不多。用软件处理时还可以自行设计许多滤波函数,应该说还是比较方便的。
    完成的频域里的处理以后再用IFFT变换,变换到时域里去,就可以看到处理的效果了。
      

  3.   

    可惜我不用FFT都没有能够从频域恢复到实域!