【技术实现步骤摘要】
本专利技术属于电数字数据处理
,具体地说是一种基于三叉树的分裂基FFT算法的优化方法。
技术介绍
分裂基FFT算法是一种快速计算离散傅里叶变换的算法,是由P.Duhamel和H.Hollmann发明的,可计算一维的实数和复数数据以及规模为2的整数次幂的离散傅里叶变换。和基2的FFT算法不同,分裂基FFT算法将计算大规模的傅里叶变换转换为计算三个小规模的子傅里叶变换过程,从而递归地将子傅里叶变换进行划分,最终完成FFT的计算过程。相对于基2算法,它拥有更少的计算量,能很好地利用在信号处理等领域。但是,现在分裂基FFT算法主要的实现方式是递归调用子过程,该方式在递归调用时带来大量的系统开销和过程调用的开销。
技术实现思路
本专利技术是为了克服现有技术存在的不足之处,提供一种基于三叉树的分裂基FFT算法的优化方法,以期能解决降低分裂基算法计算过程中的系统开销,从而能提升傅里叶变换计算过程中性能。本专利技术为解决技术问题采用如下技术方案:本专 ...
【技术保护点】
一种基于三叉树的分裂基FFT算法的优化方法,其特征是按如下步骤进行:步骤1、全局定义定义分裂基FFT算法的源向量为X[X0,X1,…,Xk,…,XN‑1];N表示源向量的长度;Xk表示FFT算法中的第k+1个复数数据;定义旋转因子数组为W[W0,W1,…,Wk,…,WN‑1];Wk表示第k+1个旋转因子;0≤k≤N‑1;定义全局变量为β,γ;定义链表L;定义所述三叉树中每个节点的计算信息包含:计算规模n、计算的起始位置pos、旋转因子在所述旋转因子数组W中的位置posOmega、父亲节点father、以及左子树left、中子树mid和右子树right;步骤2、使用迭代过程建 ...
【技术特征摘要】
1.一种基于三叉树的分裂基FFT算法的优化方法,其特征是按如下步骤进行:
步骤1、全局定义
定义分裂基FFT算法的源向量为X[X0,X1,…,Xk,…,XN-1];N表示源向量的长度;Xk表
示FFT算法中的第k+1个复数数据;
定义旋转因子数组为W[W0,W1,…,Wk,…,WN-1];Wk表示第k+1个旋转因子;0≤k≤N-1;
定义全局变量为β,γ;
定义链表L;
定义所述三叉树中每个节点的计算信息包含:计算规模n、计算的起始位置pos、旋转因
子在所述旋转因子数组W中的位置posOmega、父亲节点father、以及左子树left、中子树mid
和右子树right;
步骤2、使用迭代过程建立一棵三叉树Tree:
步骤2.1、定义在创建每个节点时的传入参数包含:
传入参数a为整数类型,表示当前计算规模;
传入参数b为整数类型,表示当前计算的起始位置;
传入参数c为整数类型,表示当前旋转因子在所述旋转因子数组W中的位置;
传入参数d为指针类型,表示指向当前处理节点的父亲节点;
步骤1.3、初始化所述传入参数和全局变量:a=N;b=0;c=1;β=1;γ=1;
步骤1.4、创建第β个节点Rootβ:
更新所述第β个节点Rootβ的计算规模nβ为a、计算的起始位置posβ为b、旋转因子在所
述旋转因子数组W中的位置posOmegaβ为c、父亲节点fatherβ为d;
步骤1.5、判断γ>β是否成立,如果成立,则跳入步骤2;否则,跳入步骤1.6;
步骤1.6、判断nγ=2或nγ=4是否成立,若成立,则将所述第γ个节点Rootγ的左子树leftγ、
中子树midγ和右子树rightγ均设置为空;否则跳入步骤1.7;
步骤1.7、令b=posγ、c=2×posOmegaγ、d=Rootγ;
步骤1.8、将β+1赋值给β后,创建第β个节点Rootβ:
更新所述第β个节点Rootβ的计算规模nβ为a、计算的起始位置posβ为b、旋转因子在所
述旋转因子数组W中的位置posOmegaβ为c、父亲节点fatherβ为d;
步骤1.9、令c=4×posOmegaγ、d=Rootγ;
步骤1.10、将β+1赋值给β后,创建第β个节点Rootβ:
更新所述第β个节点Rootβ的计算规模nβ为a、计算的起始位置posβ为b、旋转因子在所
述旋转因子数组W中的位置posOmegaβ为c、父亲节点fatherβ为d;
步骤1.11、令c=4×posOmegaγ、d=Rootγ;
步骤1.12、将β+1赋值给β后,创建第β个节点Rootβ:
更新所述第β个节点Rootβ的计算规模nβ为a、计算的起始位置posβ为b、旋转因子在所
述旋转因子数组W中的位置posOmegaβ为c、父亲节点fathe...
【专利技术属性】
技术研发人员:顾乃杰,张明,任开新,
申请(专利权)人:合肥康捷信息科技有限公司,
类型:发明
国别省市:安徽;34
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。