c语言搜索算法伪代码,FFT算法伪代码
FFT简介(来源:百度百科):
快速傅立叶变换(FFT)是离散傅立叶(DFT)的快速算法,它是根据离散傅立叶变换的奇、偶、虚、实等特性,对离散傅立叶变换的算法进行改进获得的。它对傅立叶变换的理论并没有新的发现,但是对于在计算机系统或者说数字系统中应用离散傅立叶变换,可以说是进了一大步。
设x(n)为N项的复数序列,由DFT变换,任一X(m)的计算都需要N次复数乘法和N-1次复数加法,而一次复数乘法等于四次实数乘法和两次实数加法,一次复数加法等于两次实数加法,即使把一次复数乘法和一次复数加法定义成一次“运算”(四次实数乘法和四次实数加法),那么求出N项复数序列的X(m),即N点DFT变换大约就需要N2次运算。当N=1024点甚至更多的时候,需要N2=1048576次运算,在FFT中,利用WN的周期性和对称性,把一个N项序列(设N=2k,k为正整数),分为两个N/2项的子序列,每个N/2点DFT变换需要(N/2)2次运算,再用N次运算把两个N/2点的DFT变换组合成一个N点的DFT变换。这样变换以后,总的运算次数就变成N+2(N/2)2=N+N2/2。继续上面的例子,N=1024时,总的运算次数就变成了525312次,节省了大约50%的运算量。而如果我们将这种“一分为二”的思想不断进行下去,直到分成两两一组的DFT运算单元,那么N点的DFT变换就只需要Nlog2N次的运算,N在1024点时,运算量仅有10240次,是先前的直接算法的1%,点数越多,运算量的节约就越大,这就是FFT的优越性。
Iterative_FFT
参数 n ——a的元素个数
a[0,1,…,(n-1)]——系数
局部变量
wd——单位元的n次原根
w——当前点的估值多项式
y——变换结果
y = BitReverse(a)
for j = 1 to logn
d =
wd =
w = 1
for k = 0 to d/2-1
for m = k to n-1 step d
t = w×y[m + d/2]
x = y[m]
y[m] = x + t
y[m + d/2] = x – t
endfor
w = w×wd
endfor
endfor
return y
C语言实现:
void Iterative_FFT(double *ir,double *ii,double *or,double
*oi,int len)
{
int i,j,n,d,k,m;
double
t;
double wdr,wdi,wr,wi;
double wtr,wti;
double tr,ti;
double xr,xi;
n=pow2(len);
for(i=0;i
{
j = BitReverse(i,len);
or[j] = ir[i];
oi[j] = ii[i];
}
for(j=1;j<=len;j++)
{
d = pow2(j);
t = 2*pi/d;
wdr = cos(t); wdi = sin(t);
wr = 1; wi = 0;
for (k=0;k<=d/2-1;k++)
{
for (m=k;m<=n-1;m+=d)
{
tr = wr*or[m + d/2] - wi*oi[m + d/2];
ti = wr*oi[m + d/2] + wi*or[m + d/2];
xr = or[m]; xi = oi[m];
or[m] = xr + tr; oi[m] = xi + ti;
or[m+d/2] = xr - tr; oi[m+d/2] = xi -
ti;
}
wtr = wr*wdr - wi*wdi;
wti = wr*wdi + wi*wdr;
wr = wtr;
wi = wti;
}
}
}
FFT的串行递归实现伪代码:
Recursive_FFT
参数 n ——a的元素个数
a[0,1,…,(n-1)]——系数
局部变量
wn——单位元的n次原根
w——当前点的估值多项式
a0——偶数项系数
a1——奇数项系数
y——变换结果
y0——a0的FFT结果
y1——a1的FFT结果
if n = =1 then return a
else
wn =
w = 1
a0 = (a[0],a[2],…a[n-2])
a1 = (a[2],a[3],…a[n-1])
y0 = Recursive_FFT(a0,n/2)
y1 = Recursive_FFT(a1,n/2)
for k = 0 to n/2-1
y[k] = y0[k] + w×y1[k]
y[k +
n/2] = y0[k]-w×y1[k]
w =
w×wn
endfor
return y
endif
C语言实现:
void Recursive_FFT(double *ir,double *ii,double *or,double
*oi,int len)
{
int i,k;
int n;
double wnr,wni;
double wr,wi;
double wtr,wti;
double temp0_ir[MAX_N], temp0_ii[MAX_N], temp1_ir[MAX_N],
temp1_ii[MAX_N];
double temp0_or[MAX_N], temp0_oi[MAX_N], temp1_or[MAX_N],
temp1_oi[MAX_N];
memset(temp0_ir, 0 ,sizeof(temp0_ir));
memset(temp0_ii, 0 ,sizeof(temp0_ii));
memset(temp1_ir, 0 ,sizeof(temp1_ir));
memset(temp1_ii, 0 ,sizeof(temp1_ii));
memset(temp0_or, 0 ,sizeof(temp0_oi));
memset(temp0_or, 0 ,sizeof(temp0_oi));
n=pow2(len);
if(len == 0)
{
or = ir;
oi = ii;
}
else
{
wnr =
cos(2*pi/n);
wni = sin(2*pi/n);
//printf("wnr:%.2f wni:%.2f\n",wnr,wni);
wr = 1;
wi = 0;
for(i=0; i
{
temp0_ir[i/2] = ir[i];
temp0_ii[i/2] = ii[i];
temp1_ir[i/2] = ir[i+1];
temp1_ii[i/2] = ii[i+1];
}
//for(i=0; i<2; i++)
//{
// printf("%.2f %.2f %.2f
%.2f\n",temp0_ir[i],temp0_ii[i],temp1_ir[i],temp1_ii[i]);
//}
Recursive_FFT(temp0_ir, temp0_ii, temp0_or, temp0_oi, len-1);
Recursive_FFT(temp1_ir, temp1_ii, temp1_or, temp1_oi, len-1);
for(k=0; k<=n/2-1; k++)
{
or[k] = temp0_or[k] + (wr*temp1_or[k] - wi*temp1_oi[k]);
oi[k] = temp0_oi[k] + (wr*temp1_oi[k] + wi*temp1_or[k]);
or[k+n/2] = temp0_or[k] - (wr*temp1_or[k] - wi*temp1_oi[k]);
oi[k+n/2] = temp0_oi[k] - (wr*temp1_oi[k] + wi*temp1_or[k]);
wtr = wr*wnr - wi*wni;
wti = wr*wni + wi*wnr;
wr = wtr;
wi = wti;
}
}
}
c语言搜索算法伪代码,FFT算法伪代码相关推荐
- c代码实现 ifft运算_C语言系列之FFT算法实现
0x10 序言 长文预警,详细介绍FFT算法的编程原理和C实现,并在文章的最后附上了本文的所有源代码. 0x11 速览 1)FFT背后的数学原理 2)码位倒序 3)蝶形运算设计 4)利用复数FFT编写 ...
- c语言程序设计 算法,C语言程序设计第二章算法
<C语言程序设计第二章算法>由会员分享,可在线阅读,更多相关<C语言程序设计第二章算法(38页珍藏版)>请在人人文库网上搜索. 1.C程序设计,主讲人:袁丽,燕大里仁基础教学部 ...
- 用c语言实现的FFT
一.对FFT的介绍 1. FFT(Fast Fourier Transformation),即为快速傅里叶变换,是离散傅里叶变换的快速算法,它是根据离散傅里叶变换的奇.偶.虚.实等特性,对离散傅里叶变 ...
- 贝叶斯算法伪代码_字节跳动算法面经
写在前面 在某大厂做算法一年,本来安安稳稳还没想过跳槽,某天接到了猎头的电话说字节在招数据挖掘工程师,择日不如撞日那就去试一试呗. 一面 问贝叶斯网络的原理. 一些统计学的原理比如t-test, AU ...
- 算法描述怎么写伪代码java_伪代码描述算法
伪代码是自然语言和类编程语言组成的混合结构.它比自然语言更精确,描述算法很简洁:同时也可以很容易转换成计算机程序.下面就为大家介绍一下伪代码描述算法的介绍. 伪代码描述算法 一.算法描述是指对设计出的 ...
- Latex 编写算法伪代码,基于algorithmicx包的使用说明(人工翻译自CTAN)
目录 摘要 其他布局的package 简介 算法块 简单的一行 注释 标签和引用 分解较长的算法 同一文档中使用多布局 结构化语法 for语句块 while循环 repeat语句 if语句块 proc ...
- 算法描述怎么写伪代码java_算法描述---伪代码
算法描述 算法描述是指对设计出的算法,用一种方式进行详细的描述,以便与人交流.描述可以使用自然语言.伪代码,也可使用程序流程图,但描述的结果必须满足算法的五个特征. 使用自然语言描述算法显然很有吸引力 ...
- c代码实现 ifft运算_fft算法c语言_matlab fft算法_ifft c语言
FFT快速算法C程序_工学_高等教育_教育专区.电子信息工程综合课程设计报告书 DSP 课程设计 报告 题学 目: 院: FFT 快速算法 C 程序 计算机与信息工程学院 09 ... fft算法代码 ...
- 写出TREE-PREDECESSOR的伪代码(算法导论第三版12.2-3)
写出TREE-PREDECESSOR的伪代码(算法导论第三版12.2-3) TREE-PREDECESSOR(x)if x.left != NILreturn TREE-MAXIMUM(x.left) ...
- dsp实现快速傅里叶的C语言程序,DSP-快速傅立叶变换(FFT)算法实验
<DSP-快速傅立叶变换(FFT)算法实验>由会员分享,可在线阅读,更多相关<DSP-快速傅立叶变换(FFT)算法实验(10页珍藏版)>请在人人文库网上搜索. 1.中 南 大 ...
最新文章
- Java项目:设计管理系统(java+SSM+JSP+MYSQL+layui+Maven)
- 为啥程序员下班后只关显示器从不关电脑?看看各大网站的答案~
- boost::multiprecision模块cpp_complex相关的测试程序
- G6 图可视化引擎——核心概念——节点/边/Combo——内置节点——Rect
- java集合的遍历_Java集合之遍历
- Captaris Workflow开发系列课程介绍。
- bat 调用class文件_拯救写框架的程序员!用字节码替代反射,实现任意函数调用...
- 情侣在公交车上接吻遭乘客指责
- 《中国人史纲》读书笔记:第二章 神话时代 第三章 传说时代
- GetLastError错误码
- 高德sdk android加载3857,GCJ02-Correct
- RAID技术及应用总结
- 2T硬盘安装WIN7分区格式
- 基于FOC电路低次谐波抑制Simulink仿真
- STM32开发基础知识——定时器
- 柴胡加龙骨牡蛎汤去大黄合当归芍药散治顽固心悸案(李跃海)
- Flutter-解决Try catch出现异常:type ‘_TypeError‘ is not a subtype of type ‘Exception‘ in type cast
- 在移动硬盘里移动视频文件到移动硬盘 另外一个文件夹 显示正在计算_稳定可靠的数据之仓 柯达X200 SSD固态移动硬盘体验评测...
- 2015.3.30第一次博客测试
- som聚类用matlab实现,使用SOM对数据进行聚类