专利名称：一种dft/idft的快速计算方法和装置的制作方法
技术领域：
本发明涉及，特别涉及一种离散傅立叶变换(DFT)/逆离散傅立叶变换(IDFT)的 快速计算方法和装置。
背景技术：
为了在移动通信领域保持领先地位，把握新一轮技术浪潮，3GPP组织在2004年底 启动了 3G长期演进系统(LTE)的研究与标准化工作。LTE系统为了避免多载波系统带来的 高峰均比(PAPR)，同时又能利用多载波系统的优势，上行系统采用单载波的DFT-OFDMA接 入方式。图l给出一个LTE上行发射系统框图。 如图l所示，在上行系统中每次分配给用户的资源数(子载波数)为12的整数倍， 因此不能直接利用FFT/IFFT完成子载波的映射和解映射。同时，直接计算DFT/IDFT带来 的计算复杂度难以承受，使得整个上行系统的实现复杂度很高。另外，在OFDM系统中采取 导频辅助信道估计时，需要采用DFT/IDFT进行插值或基于域变换的去噪处理，也同样需要 面对直接计算DFT/IDFT带来的计算复杂度问题，从而导致去噪处理时延很大。总之，在需 要进行DFT/IDFT处理的系统中，由于直接计算DFT/IDFT的计算复杂度很高，使得整个系统 的实现复杂度增大，时延较大，处理速度较低。
发明内容
有鉴于此，本发明提供一种DFT/IDFT的快速计算方法和装置，能够降低DFT/IDFT
的计算复杂度，降低处理时延，从而大大提高DFT/IDFT所在系统的处理速度。 为实现上述目的，本发明采用如下的技术方案 —种离散傅立叶变换/逆离散傅立叶变换的快速计算方法，包括 将待变换的比特序列补零，使补零后的序列长度符合FFT/IFFT处理器的处理范
围； 利用FFT/IFFT处理器对补零后的序列进行FFT/IFFT计算，并按照由补零前和补 零后的序列长度确定的采样位置对FFT/IFFT计算结果序列进行时域/频域重采样；
将重采样结果作为DFT/IDFT的结果序列输出。
较佳地，所述将待变换的比特序列补零为 预先确定FFT/IFFT处理器的处理点数，在所述处理点数中，为DFT/IDFT的变换点 数选择对应的FFT/IFFT的处理点数，建立并保存二者间的对应关系； 在进行补零操作时，根据待变换的比特序列长度，选择对应的处理点数，将待变换 的比特序列补零到比特序列长度等于所选择的处理点数。 较佳地，当所述FFT/IFFT处理器的处理点数为2048时，建立并保存的所述对应关 系为 若补零前的序列长度为 120，则补零后的序列长度为512 ;
若补零前的序列长度为132 300，则补零后的序列长度为1024 ;
若补零前的序列长度为312 600，则补零后的序列长度为2048 ; 若补零前的序列长度为612 900，则补零后的序列长度为2048 ; 若补零前的序列长度为312 1200，则补零后的序列长度为2048。 较佳地，当所述FFT/IFFT处理器的处理点数为4096时，建立并保存的所述对应关
系为 若补零前的序列长度为12 120，则补零后的序列长度为512 ; 若补零前的序列长度为132 300，则补零后的序列长度为1024 若补零前的序列长度为312 600，则补零后的序列长度为2048 若补零前的序列长度为612 900，则补零后的序列长度为4096 若补零前的序列长度为312 1200，则补零后的序列长度为4096。 较佳地，当进行DFT时，在待变换的序列尾部补零； 当进行IDFT时，在待变换的序列中间补零；若待变换序列长度N为奇数时，则在所
述待变换的序列中索引值为'L^/2」和索引值为^/2」+l的符号间补零，或者，在所述待 变换的序列中索引值为「^V/2，P索引值为「iV/2，+l的符号间补零。 较佳地，由补零前和补零后的序列长度确定采样位置的方式为index(i)=
i*NFFT/IFFT/NDFT/IDFT，其中，NFFT/IFFT为补零后的序列长度，NDFT/IDFT为补零前的序列长度，
index(i)为FFT/IFFT变换后的采样位置索引值，i为重采样结果的索引值。 较佳地，当所述index (i)为整数时，所述进行时域/频域重采样为xDFT/IDFT (i)=
xfft/ifft (index—int (i)) ， i = 0，1，2，. ，NDFT/IFD1~l，xDFT/IDFT (i)为输出的DFT/IDFT结果序列
中的第i个值，xFFT/IFFT(i)为FFT/IFFT计算结果序列中的第i个值。 较佳地，当所述index (i)不是整数时，所述进行时域/频域重采样为 将index (i)表示为整数部分index」nt (i)和小数部分index_frac (i)，并根据计
算精度截取所述小数部分index_frac(i)中的前M位index_frac_M(i); 当截取的前M位indeX_fraC_M(i)中所有比特均为0时，xDFT/IDFT (i) = xFFT/
丽(index—int(i));当截取的前M位index—frac—M(i)中所有比特均为1时，xDFT/IDFT(i)=
Xfft/ifft (index_int(i)+l); 当截取的前M位indeX_fraC_M(i)中的比特不全为0、且不全为1时，根据预先确 定的插值方式，对x附/肿t(indexjnt(i))禾P xFFT/IFFT (index_int (i)+1)进行小数插值得到
xdft/idft (i)。 较佳地，所述预先确定的插值方式为零阶保持，则所述对~1/,(111(^_11^(1)) 和xfft/ott(indexjnt(i)+l)进行小数插值得到xDFT/IDFT (i)为当截取的前M位index_frac_ M(i)中的最高位为0时，xDFT/IDFT(i) = xFFT/IFFT(index_int(i));当截取的前M位index— frac—M(i)中的最高位为1时，xDFT/IDFT(i) = xFFT/IFFT (index_int (i)+1);
所述预先确定的插值方式为一阶线性插值，则所述对xFFT/IFFT(indeX_int(i))和
xfft/ifft
(index_int(i)+l)进行小数插值得至Ll xDFT/IDFT (i)为
xDFT/IDFT(i) = xFFT/IFFT (index_int (i)) * (l-index_frac_M (i))+xFFT/IFFT (index— int(i)+l)*index_frac_M(i); 所述预先确定的插值方式为二阶拉格朗日插值，则所述对~1/1附(111(1^_11^(1))禾口 xFFT/IFFT(index_int(i)+l)进畔亍小数插值f寻至Ll xDFT/IDFT (i)为:xDFT/IDFT (i) = xFFT/IFFT (index— int (i) _1)承index—frac—M (i)承(index—frac—M (i) _1) /2+xFFT/IFFT (index_int (i))承(index— frac_M (i) +1) * (l_index_f rac_M (i)) +xFFT/IFFT (index_int(i)+l)* (index_f rac_ M (i) +1) *index_frac_M (i) /2 。 较佳地，计算的(l-index_frac_M(i))方式为，对index_frac_M(i)中的比特逐位 进行取反操作。
—种离散傅立叶变换/逆离散傅立叶变换的快速计算装置，包括 补零器，用于将输入的待变换比特序列补零，使补零后的序列长度符合FFT/IFFT
处理器的处理范围； FFT/IFFT处理器，用于对补零后的序列进行FFT/IFFT计算； 重采样器，用于按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT 计算结果序列进行时域/频域重采样，并将重采样结果作为DFT/IDFT的结果序列输出。
由上述技术方案可见，本发明中，将待变换的比特序列补零，使补零后的序列长度 符合FFT/IFFT处理器的处理范围；利用FFT/IFFT处理器对补零后的序列进行FFT/IFFT计 算，并按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算结果序列进行 时域/频域重采样；将重采样结果作为DFT/IDFT的结果序列输出。通过上述方式，本发明利 用计算速度较快的FFT/IFFT代替直接进行DFT/IDFT，并对FFT/IFFT处理结果进行重采用 以得到DFT/IDFT处理结果。由于FFT/IFFT的处理速度很快，因此能够大大提高DFT/IDFT 的计算速度，降低DFT/IDFT的计算复杂度，降低处理时延，从而大大提高DFT/IDFT所在系 统的处理速度。
图1为LTE上行发射系统的结构示意图。 图2为本发明提供的DFT快速计算方法的流程图。 图3为上述本发明的方法和装置在LTE上行发射基带系统中的一个具体应用实 例。
具体实施例方式
为使本发明的目的、技术手段和优点更加清楚明白，以下结合附图对本发明做进 一步详细说明。 本发明的基本思想是利用FFT的处理代替DFT的处理，利用IFFT的处理代替 IDFT的处理，从而提高运算速度，降低处理复杂度。 图2为本发明提供的DFT快速计算方法的流程图。如图2所示，该方法包括
步骤201，将待变换的比特序列补零，使补零后的序列长度符合FFT/IFFT处理器 的处理范围。 本发明中，利用FFT的处理代替DFT的处理，利用IFFT的处理代替IDFT的处理， 以下简写为利用FFT/IFFT处理代替DFT/IDFT处理。 如前所述，本发明中利用FFT/IFFT代替DFT/IDFT，从而降低计算复杂度，提高计 算速度。而在FFT/IFFT处理过程中，可处理的序列点数为2的整数次幂。因此，本步骤中将需要DFT/IDFT的序列补零到某个2的整数次幂，然后在后续步骤中对FFT/IFFT的计算 结果再分别进行频域和时域的重采样得到DFT/IDFT的计算结果。
具体的补零方式可以为 预先确定FFT/IFFT处理器的处理点数，在所述处理点数中，为DFT/IDFT的变换点 数选择对应的FFT/IFFT的处理点数，建立并保存二者间的对应关系； 在进行补零操作时，根据待变换的比特序列长度，选择对应的处理点数，将待变换 的比特序列补零到比特序列长度等于所选择的处理点数。 其中，在FFT/IFFT处理中，FFT/IFFT处理器均存在最大的处理点数M，该处理器可 以处理的序列长度为小于等于M的任意2的整数次幂。因此，在建立所述对应关系时，任 意一种DFT/IDFT的变换点数可以对应大于N、且小于等于M的任意2的整数次幂，其中，N 为DFT/IDFT的变换点数。并且，在这些可选的点数中，FFT/IFFT处理点数越高变换结果的 精度越高，但相应所需的处理时间也要增长。因此，在建立所述对应关系时，可以根据实际 应用的需要，在处理时间和变换精度两方面寻求均衡。例如，本发明给出在LTE系统中应用的DFT/IDFT变换点数和FFT/IFFT处理点数
的对应关系，如表l所示。
DFT/IDFT点数12 120132 300312 600612 900912 1200
FFT/IFFT点数(2048)5121024204820482048
FFT/IFFT点数(4096)5121024204840964096 表1 表1中给出了 FFT/IFFT最大处理点数分别为2048点和4096点时，DFT/IDFT变 换点数和FFT/IFFT处理点数的对应关系。在进行补零操作时，可以根据表1确定补零后的 序列长度。例如，若DFT/IDFT点数(即分配给用户的子载波数)为660，则对于最大处理点 数为2048的FFT/IFFT处理器，补零到序列长度为2048，对于最大处理点数为4096的FFT/ IFFT处理器，补零到序列长度为4096。 这里仅是举个简单的例子，在具体应用中，还可以根据系统中业务的时延要求和 数据传输精度，进行具体设置。 对于DFT和IDFT补零的位置稍有不同，进行DFT时，在待变换序列的尾部补零；进 行IDFT时，在待变换序列的中间补零。当在待变换序列中间补零时，若待变换序列长度N
为奇数时，则在所述待变换的序列中索引值为L^/2」和索引值为Lw/2」+i的符号间补零，
或者，在所述待变换的序列中索引值为「iV/2"l和索引值为pV/2， + l的符号间补零。例如，
设待变换序列长度为71，则可以在索引值为35和36的符号之间进行补零，或者，也可以在 索引值为36和37的符号之间进行补零。
如上，即完成了具体的补零操作。 步骤202，利用FFT/IFFT处理器对补零后的序列进行FFT/IFFT计算。
本步骤的操作与现有实现方式相同，这里就不再赘述。
7
步骤203，按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算结 果序列进行时域/频域重采样并输出。 假设补零前的序列长度为NDFT/IDFT，补零后的序列长度为NFFT/IFFT，则经过补零后的 序列经过步骤102进行FFT/IFFT处理后的序列长度仍为NFFT/IFFT。设FFT/IFFT处理的计算
结果为 XFFT/IFFT
(k) ， k = 0， 1 ， 2，. . . ， NFFT/IFFT-1 ， NDFT/IDFT点序列的DFT/IDFT结果为
XDFT/IDFT (n)，
n = 0， 1，2，. . . ，NDFT/IDFT-1。则需要在本步骤中，根据FFT/IFFT处理的计算结果xFFT/IFFT(k)， 进行重采样，得到序列xDFT/IDFT (n) ， n = 0， 1， 2，…，NDFT/IDFT-1，作为DFT/IDFT的处理结果进 行输出。对于FFT处理，则进行频域重采样，对于IFFT处理，则进行时域重采样。
假定利用i表示重采样后输出序列的索引值，index (i) = i*NFFT/IFFT/NDFT/IDFT为对 xFFT/IFFT(k)进行重采样的位置，也即序列xFFT/IFFT(k)的索引惶，从而将Xfft加ft(index(i))赋
值给 XDFT/IDFT 、1 )输出。 若重采样位置index (i)为整数，则可以直接令xDFT/IDFT (i) = xFFT/IFFT (index (i))输 出即可。若重采样位置index(i)为非整数，即进行非整数倍采样时，则需要通过插值方式
确定 XDFT/IDFT (i)。 具体地，可以将index(i)表示为整数部分index_Int(i)和小数部分index_ frac(i)。将index (i)用无符号二进制数表示，其中index—int(i)包含ll比特(采用2048 点FFT/IFFT处理器)，或者12bit (采用4096点FFT/IFFT处理器)。假设index_frac (i) 部分包含的比特数为fare个比特，即index—frac(i)的无符号二进制数表示表示形式为 bA. bfrac—!，其中b。为MSB, bfrac—丄为LSB。 为减小计算量，对index (i)的小数部分index_frac (i)进行简化处理。具体地，
首先，根据计算精度的要求截取小数部分indeX_fraC(i)中的前M位index— fracM(i)，精度要求越高，M值越大，处理时间越长，反之，精度要求越低，M值越小，处理时 间越短。因此，这里截取小数部分位数时也需要根据计算精度和处理时间进行均衡选择。具 体实现时，可以对x^/，T(k)的相邻两个样点之间的间隔进行8等分或者16，32等分。8等 分时截取index_frac (i)的前3比特，即b。b山2。 16等分时候截取index_frac (i)的前4比 特，即b。b山2b3。 32等分时候截取index_frac (i)的前5比特，即KbibAslv用index_frac_ M(i)表示上述三种对index—frac (i)的截取情况，M为3时表示截取index—frac (i)的前 3比特。M为4时表示截取index_frac (i)的前4比特。M为5时表示截取index_frac (i) 的前5比特。 接下来，当截取的前M位indeX_fraC_M(i)中所有比特均为0时，表示重采样位置 很接近indexjnt(i)，因此令x。FT/n)FT(i) = xFFT/IFFT (indexjnt (i))输出；当截取的前M位 index_fraC_M(i)中所有比特均为1时，表示重采样位置很接近index」nt (i)+1，因此令 xDFT/mFT(i) = xFFT/IFFT(index_int(i)+l)输出; 当截取的前M位indeX_fraC_M(i)中的比特不全为0、且不全为1时，根据预先确 定的插值方式，对x附/肿T(indexjnt(i))禾P xFFT/IFFT (indexjnt (i)+1)进行小数插值得到 xDFT/IDFT(i)。具体插值方式可以采样现有的任意插值方式，不同的插值方式计算复杂度和精 度均有所不同，通常高阶插值方式计算精度高，但处理复杂度也相对较高。本实施例中仅以 三种插值方式为例进行说明，这三种插值方式分别为零阶保持，一阶线性插值，二阶拉格朗 日插值。
若预先确定的插值方式为零阶保持插值，则本步骤中进行非整数倍重采样时的具 体插值操作为当index_frac_M(i)中的MSB比特b。为零时，
xDFT/iDFT(i) = xFFT/IFFT (index—
int (i)) ， i = 1，2， ， NDFT/IFDT—1 ;当b0为查时，xDFT/IDFT(i) = xFFT/IFFT(index_int (i)+1) ， i
—1，2， ， Ndft/ifdT_1。 若预先确定的插值方式为一阶线性插值，则本步骤中进行非整数倍重采样时的具 体插值操作为 xDFT/IDFT(i) = xFFT/IFFT (index_int (i)) * (l-index_frac_M (i))+xFFT/IFFT (index— int(i)+l)hndex—frac^M其中，(l-index_frac_M(i))的运算可—/rac —M(/) 代替，以减少计算量，Z^fec —/rac — M(0表示对indeX_frc_M(i)中的比特逐位进行取反操 作。 若预先确定的插值方式为二阶拉格朗日插值，则本步骤中进行非整数倍重采样时 的具体插值操作为 xDFT/IDFT(i) = xFFT/IFFT (index_int (i)-1) *index_frac_M (i) * (index_frac_ M(i) _1) /2+xFFT/IFFT(index_int (i)) * (index_frac_M (i) +1) * (l_index_frac_M (i)) +xFFT/ IFFT (index_int (i) +1) * (index_frac_M (i) +1) *index_frac_M (i) /2，与一阶线性插值的计算 类似，(l-index_fraC_M(i))的运算可以利用Z^fcc — /rac — M(/)代替，以减少计算量。
通过上述方式即可以通过插值算法确定非整数倍采样时XD^n^(i)的取值。按照 上述方式，逐个确定重采样位置index(i)，并按照所确定的位置进行重采样得到序列xd^ IDFT(n)并输出，即为DFT/IDFT的处理结果。
至此，本发明的快速计算方法流程结束。 本发明还提供了一种DFT/IDFT的快速计算装置，用于需要计算DFT/IDFT的系统 中。该装置具体包括补零器，用于将输入的待变换比特序列补零，使补零后的序列长度符 合FFT/IFFT处理器的处理范围；FFT/IFFT处理器，用于对补零后的序列进行FFT/IFFT计 算；重采样器，用于按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算结 果序列进行时域/频域重采样，并将重采样结果作为DFT/IDFT的结果序列输出。上述装置 可以用于实施上述图2所示的本发明的方法。 图3为上述本发明的方法和装置在LTE上行发射基带系统中的一个具体应用实 例。如图3所示，在上行发送端，将数据进行编码调制后，根据系统分配的子载波数目(在 LTE中为12的整数倍)进行分组。然后对分组后的序列按照上述步骤201中的方式进行补 零操作，得到FFT处理器可以处理的序列长度。补零后的序列经过FFT处理器后得到频域 序列，然后对这个频域序列按照上述步骤203中的方式进行重采样得到与分组后的发射数 据相同长度的序列。最后对这个序列进行子载波映射、IFFT处理和插入循环前缀后就完成 LTE上行基带信号的处理。 在上述应用实例中，上行发送端需要进行DFT处理，因此利用FFT处理器对补零后 序列进行处理。若在进行IDFT处理的系统中，则需要利用IFFT处理器对补零后序列进行 处理。 由上述本发明的具体实现过程可见，本发明通过FFT/IFFT处理器达到对DFT/ IDFT的快速计算，避免了对DFT/IDFT直接计算带来的巨大的计算量。通过对时域/频域重
采样点位置的不同量化方式以及采样不同的小数插值方法可以在计算复杂度和计算精度
9之间找到一个最优平衡点。在一阶线性插值算法和二阶拉格朗日插值算法中通过对小数部 分的取反操作达到减少计算量的目的。 本发明可以广泛地应用于3GPP LTE、MT-Advanced、WiMAX等移动通信系统以及其 它需要计算DFT/IDFT的数字信号处理系统中，以降低DFT/IDFT处理的计算复杂度，提高系 统处理速度。同时，在如3GPP LTE等原本就包括FFT/IFFT处理器的系统中，实现了 FFT/ IFFT处理器的分时复用，更有效地提高了资源利用率。 以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发 明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围 之内。
10
权利要求
一种离散傅立叶变换/逆离散傅立叶变换的快速计算方法，其特征在于，该方法包括将待变换的比特序列补零，使补零后的序列长度符合FFT/IFFT处理器的处理范围；利用FFT/IFFT处理器对补零后的序列进行FFT/IFFT计算，并按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算结果序列进行时域/频域重采样；将重采样结果作为DFT/IDFT的结果序列输出。
2. 根据权利要求1所述的方法，其特征在于，所述将待变换的比特序列补零为 预先确定FFT/IFFT处理器的处理点数，在所述处理点数中，为DFT/IDFT的变换点数选择对应的FFT/IFFT的处理点数，建立并保存二者间的对应关系；在进行补零操作时，根据待变换的比特序列长度，选择对应的处理点数，将待变换的比 特序列补零到比特序列长度等于所选择的处理点数。
3. 根据权利要求2所述的方法，其特征在于，当所述FFT/IFFT处理器的处理点数为 2048时，建立并保存的所述对应关系为若补零前的序列长度为12 120，则补零后的序列长度为512 ; 若补零前的序列长度为132 300，则补零后的序列长度为1024 ; 若补零前的序列长度为312 600，则补零后的序列长度为2048 ; 若补零前的序列长度为612 900，则补零后的序列长度为2048 ; 若补零前的序列长度为312 1200，则补零后的序列长度为2048。
4. 根据权利要求2所述的方法，其特征在于，当所述FFT/IFFT处理器的处理点数为 4096时，建立并保存的所述对应关系为若补零前的序列长度为12 120，则补零后的序列长度为512 ; 若补零前的序列长度为132 300，则补零后的序列长度为1024 ; 若补零前的序列长度为312 600，则补零后的序列长度为2048 ; 若补零前的序列长度为612 900，则补零后的序列长度为4096 ; 若补零前的序列长度为312 1200，则补零后的序列长度为4096。
5. 根据权利要求1到4中任一所述的方法，其特征在于，当进行DFT时，在待变换的序 列尾部补零；当进行IDFT时，在待变换的序列中间补零；若待变换序列长度N为奇数时，则在所述待 变换的序列中索引值为L^/2」和索引值为^V/2」+l的符号间补零，或者，在所述待变换 的序列中索引值为「 vV/2"l和索引值为「 7V/21 +1的符号间补零。
6. 根据权利要求1所述的方法，其特征在于，由补零前和补零后的序列长度确定采样位置的方式为<formula>formula see original document page 2</formula>，其中，NFFT/丽为补零后的序列长度，N。FT/丽为补零前的序列长度，index(i)为FFT/IFFT变换后的采样位置索引值，i为重采样结果的索 引值。
7. 根据权利要求6所述的方法，其特征在于，当所述index(i)为整数时，所述进行时 域/频域重采样为<formula>formula see original document page 2</formula>为输出的DFT/IDFT结果序列中的第i个值，xFFT/IFFT (i)为FFT/IFFT计算结果序列中 的第i个值。
8. 根据权利要求6或7所述的方法，其特征在于，当所述index(i)不是整数时，所述进 行时域/频域重采样为将index (i)表示为整数部分index」nt (i)和小数部分index—frac (i)，并根据计算精 度截取所述小数部分index_frac(i)中的前M位index_frac_M(i);当截取的前M位index_frac_M(i)中所有比特均为0时，xDFT/IDFT(i) = xFFT/IFFT(index_ int (i));当截取的前M位index_fraC_M(i)中所有比特均为1时，xDFT/IDFT (i) = xFFT/ IFFT (index_int (i)+1);当截取的前M位index_fraC_M(i)中的比特不全为0、且不全为1时，根据预先确定的 插值方式，对xFFT/IFFT(index_int(i))禾口 xFFT/IFFT (index_int (i)+1)进行小数插值得至Ll xDFT/idft ( i )。
9. 根据权利要求8所述的方法，其特征在于，所述预先确定的插值方式为零阶保持，则所述对xFFT/IFFT (indexjnt (i))和xFFT/ IFFT(index_int(i)+l)进行小数插值得到xDFT/IDFT (i)为当截取的前M位index_frac_M(i) 中的最高位为0时，x。fw諮(i) =xFFT/IFFT(index_int(i));当截取的前M位index_frac_M(i) 中的最高位为1时，xDFT/IDFT(i) = xFFT/IFFT(index_int(i)+l);所述预先确定的插值方式为一阶线性插值，则所述对xFFT/IFFT(index_int(i))和xFFT/ IFFT(index_int(i)+l)进行小数插值得到xDFT/IDFT (i)为xdft/idft (i) =xfft/ifft (index_int (i))氺(1_index_f rac_M (i)) +xFFT/IFFT (index— int(i)+l)*index_frac_M(i);所述预先确定的插值方式为二阶拉格朗日插值，则所述对xFFT/IFFT(index_int (i))和xfft/ifft(indexjnt(i)+l)进行小数插值得到x。ft/丽(i)为xDFT/IDFT(i) =xFFT/IFFT (index— int (i) _1)承index—frac—M (i)承(index—frac—M (i) _1) /2+xFFT/IFFT (index_int (i))承(index— frac_M (i) +1) * (l_index_f rac_M (i)) +xFFT/IFFT (index_int(i)+l)* (index_f rac_ M (i) +1) *index_frac_M (i) /2 。
10. 根据权利要求9所述的方法，其特征在于，计算的(l-index_frac_M(i))方式为，对 index_frac_M (i)中的比特逐位进行取反操作。
11. 一种离散傅立叶变换/逆离散傅立叶变换的快速计算装置，其特征在于，该装置包括补零器，用于将输入的待变换比特序列补零，使补零后的序列长度符合FFT/IFFT处理 器的处理范围；FFT/IFFT处理器，用于对补零后的序列进行FFT/IFFT计算；重采样器，用于按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算 结果序列进行时域/频域重采样，并将重采样结果作为DFT/IDFT的结果序列输出。
全文摘要
本发明公开了一种离散傅立叶变换/逆离散傅立叶变换的快速计算方法，包括将待变换的比特序列补零，使补零后的序列长度符合FFT/IFFT处理器的处理范围；利用FFT/IFFT处理器对补零后的序列进行FFT/IFFT计算，并按照由补零前和补零后的序列长度确定的采样位置对FFT/IFFT计算结果序列进行时域/频域重采样；将重采样结果作为DFT/IDFT的结果序列输出。本发明还公开了一种离散傅立叶变换/逆离散傅立叶变换的快速计算装置。应用本发明，能够降低DFT/IDFT的计算复杂度，提高计算速度。
文档编号G06F17/14GK101751375SQ200810239910
公开日2010年6月23日 申请日期2008年12月12日 优先权日2008年12月12日
发明者胡炜, 陆涛, 雷旭 申请人:普天信息技术研究院有限公司