[发明专利]一种FFT和IFFT逆序数表的并行处理方法

权利要求书

1.一种FFT和IFFT逆序数表的并行处理方法，其特征在于，包括以下步骤：

(1)预存储任意正整数N阶的逆序数索引表A，A的长度为n＝2^N；

(2)确定所需处理计算的二进制逆序数索引表B的阶数为M，B的长度m＝2^M；

(3)根据预存储的逆序数索引表A计算所需处理计算的二进制逆序数索引表B；

其是采用并行处理方式计算所需处理计算的二进制逆序数索引表B，对需并行计算的数据分段整理，并确定数据的位长及存储空间；具体过程为：

设D为一个无符号整形变量，用于递推计算中，对数值进行辅助修正，

D＝(1＜＜(M-N-1))

根据无符号整形变量D构造数组γ，

γ＝D|(D＜＜T)|(D＜＜2T)|…|(D＜＜((P-1)T))

P为一次并行计算能处理的数据量；P值的计算方式为：

W表示CPU位处理指令能支持的最长位宽，T表示逆序数索引表A与数组γ中数值最大的元素按二进制存储至少需要用到的比特位数；

T的计算方式是：

n表示逆序数索引表A中数值的个数；即逆序数索引表A为包含n个整数的有序序列；

对预存储的逆序数索引表A中的数据采用压缩存储法进行存储，

则预存储的逆序数索引表A中的n个数能够分成n′组表达，逆序数索引表B结果中的m个数能够分成m ′ 批来计算；

设定循环变量j＝n′；

则预存储逆序数索引表A中的n个数能够分成n′组表达，α_i表示第i组，X_i*n表示预存储逆序数索引表A中第i*n个数据；

A＝{α₀,α₁,…,α_n′-1},

α_i＝{X_i*n|(X_i*n+1＜＜T)|(X_i*n+2＜＜2T)|…|(X_(i+1)*n-1＜＜((P-1)T))}

B中的m个数能够分成m′批来表达，β_i表示第i组，Y_i*m表示预存储逆序数索引表B中第i*m个数据；

B＝{β₀,β₁,…,β_m′-1},

β_i＝{Y_i*m|(Y_i*m+1＜＜T)|(Y_i*m+2＜＜2T)|…|(Y_(i+1)*m-1＜＜((P-1)T))}

采用如下方式快速的求取到B；

k用来表示外层循环控制变量；i为内层循环控制变量。

2.根据权利要求1所述的FFT和IFFT逆序数表的并行处理方法，其特征在于，对预存储的逆序数索引表A中的数据采用压缩存储法进行存储的具体过程为：

(21)建立一个数组α_i，数组中每个元素的位宽均为T，且均为无符号整数；

(22)将P个元素按下标顺序依次致密写入内存中，元素之间不留空隙；

(23)把以α_i所在地址为首地址的内容看作位宽为T*P的无符号整数来取值。