fft-1.1（fft - 1.1）.docVIP

fft-1.1（fft - 1.1） #包括 iom128。” #包括内联函数。” / ********************************************************************* 快速福利叶变换C程序包 函数简介：此程序包是通用的快速傅里叶变换C语言函数，移植性强，以下部分不依 赖硬件此程序包采用联合体的形式表示一个复数，输入为自然顺序的复。 数（输入实数是可令复数虚部为0），输出为经过FFT变换的自然顺序的 复数此程序包可在初始化时调用create_sin_tab()函数创建正弦函数表， 以后的可采用查表法计算耗时较多的罪和COS运算，加快可计算速度 使用说明：使用此函数只需更改宏定义fft_n的值即可实现点数的改变，fft_n的 应该为2的N次方，不满足此条件时应在后面补0。若使用查表法计算罪值和 因为值，应在调用FFT函数前调用create_sin_tab()函数创建正弦表 函数调用：FFT（S）； 作者：吉帅虎 时间：2010-2-20 版本：ver1.1 参考文献： ********************************************************************** / #包括数学。” #定义fft_n 128 / /定义福利叶变换的点数 #定义PI 3.1415926535897932384626433832795028841971 /定义圆周率值 结构compx {浮房，图像；}；//定义一个复数结构 结构compx的[ fft_n ]；// FFT输入和输出：从的[ 0 ]开始存放，根据大小自己定义 浮sin_tab [ fft_n / 2 ]；//定义正弦表的存放空间 / ******************************************************************* 函数原型：结构compx EE（struct结构compx compx B1，B2） 函数功能：对两个复数进行乘法运算 输入参数：两个以联合体定义的复数A、B 输出参数：一和B的乘积，以联合体的形式输出 ******************************************************************* / 结构compx EE（struct compx，结构compx B） { 结构compx C； c.real = A. * * b.imag b.real-a.imag； c.imag = A. * b.imag + a.imag * b.real； 返回（C）； } / ****************************************************************** 函数原型：无效create_sin_tab（float * sin_t） 函数功能：创建一个正弦采样表，采样点数与福利叶变换点数相同 输入参数：* sin_t存放正弦表的数组指针 输出参数：无 ****************************************************************** / 无效create_sin_tab（float * sin_t） { int i； 为（i = 0；i fft_n／2；i++） sin_t [我] = sin（2 * pi *我/ fft_n）； } / ****************************************************************** 函数原型：无效sin_tab（浮PI） 函数功能：采用查表的方法计算一个数的正弦值 Input parameters: Pi to calculate the sine value of radian value, range 0--2*PI, does not meet the need to change Output parameters: input value PI sine value ******************************************************************/ Float sin_tab (float PI) { Int n; Float a; N= (int) (pi*FFT_N/2/PI); If (n=0nFFT_N/2) A=SIN_TAB[n]; Else if (n=FFT_N/2nFFT_N) A=-SIN_TAB[n-FFT_N/2]; Return a; } /*************************************