CN101331479A - 循环快速傅里叶变换 - Google Patents

循环快速傅里叶变换 Download PDF

Info

Publication number
CN101331479A
CN101331479A CNA2006800472468A CN200680047246A CN101331479A CN 101331479 A CN101331479 A CN 101331479A CN A2006800472468 A CNA2006800472468 A CN A2006800472468A CN 200680047246 A CN200680047246 A CN 200680047246A CN 101331479 A CN101331479 A CN 101331479A
Authority
CN
China
Prior art keywords
output port
logic output
fft
exchange logic
input
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CNA2006800472468A
Other languages
English (en)
Other versions
CN101331479B (zh
Inventor
利夫·威廉松
吉姆·斯文森
安德斯·贝尔克曼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Clastres LLC
Telefonaktiebolaget LM Ericsson AB
Original Assignee
Telefonaktiebolaget LM Ericsson AB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Telefonaktiebolaget LM Ericsson AB filed Critical Telefonaktiebolaget LM Ericsson AB
Publication of CN101331479A publication Critical patent/CN101331479A/zh
Application granted granted Critical
Publication of CN101331479B publication Critical patent/CN101331479B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms
    • G06F17/142Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2649Demodulators
    • H04L27/265Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
    • H04L27/2651Modification of fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators for performance improvement
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2662Symbol synchronisation

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Discrete Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Signal Processing (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Complex Calculations (AREA)

Abstract

一种在N频点FFT中使用的DIF FFT级,其中N是偶数。该DIF FFT级包括交换逻辑,该交换逻辑接收第一输入样本x(v)和第二输入样本x(v+N/2),并选择性地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,或者另选地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,其中0≤v<N/2。该DIF FFT级还包括:求和单元,用于对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相加;求差单元,用于对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相减;以及旋转因子逻辑,其将该求差单元所供应的值乘以旋转因子WN (v+s)mod(N/2),其中s是表示N个输入样本的循环移位量的整数。

Description

循环快速傅里叶变换
技术领域
本发明总体上涉及快速傅里叶变换(FFT),更具体地涉及利用输入数据的循环移位来实现FFT或等价地利用每个输出数据的旋转来实现FFT。
背景技术
FFT在从低级电信信令到语音和图像处理的现代数字信号处理中有大量应用。在某些应用中,出于特定的目的,对变换数据组进行前移位或后旋转。
通过变换理论我们已知前移位和后旋转是等价操作。然而,在常规布置中,它们的实施彼此不同。前移位器由大小至少为变换大小的缓冲器和某种寻址逻辑构成。与此对照,后旋转器被实施为复数乘法器或例如使用坐标旋转数字计算机(CORDIC)算法的纯旋转器和某种旋转角度发生器(例如存储器中存储的表)。
结合FFT使用的前移位和/或后旋转的典型应用是利用正交频分复用(OFDM)来传送信息的通信系统。目前,OFDM被用在各种无线局域网(WLAN)标准(例如,IEEE 802.11a和WiMAX)和数字电视(DVB)中。OFDM还被认为是未来的通信标准,如超宽带(UWB)和对“3G”(第三代)蜂窝系统的增强。
前移位的一种用途是改善OFDM中的时间同步。图1中例示了带有循环前缀的单时域OFDM符号。在发射器端,数据被编码成复数,它们在时域中被逆傅里叶变换且附带上循环前缀。在接收器端,通过应用正向傅里叶变换来还原出数据。即使不是全部,大多数实施都依赖于针对时频变换的某种FFT算法。
接收器时间同步(即,找到放置FFT窗口的最佳位置)对于实现良好的接收器性能是至关重要的。绝不能将窗口放置为,使得多于一个符号(包括其循环前缀)被覆盖,因为这将导致符号间干扰(ISI),即,两个独立的符号混合成一个符号。
窗口的最佳位置是将它放置得“尽可能晚”,而不包括来自下一符号的循环前缀的数据。然而,过于迅速(aggressive)的方法可能会导致ISI,所以必须有一些安全裕度。只要信道的脉冲响应足够短,尽早放置包括循环前缀部分的窗口就可以消除ISI的风险。
FFT窗口的早放置对应于变换输入数据的循环移位(见图1)。位于发射器中最右侧位置处的数据块此时出现在接收器FFT窗口的左侧。对于在信道估计器中使用频率插值的系统而言,移位越大,插值器就越复杂。为了保持信道评估简单,应当通过前移位或后旋转使输入数据在FFT窗口中对准。
OFDM中的前移位和/或后旋转的另一用途是当通信信道具有较大的延迟扩展时。这种信道经历了频域中的旋转,这能够通过所提出的FFT方案来消除。
前移位和/或后旋转的益处并不能毫无代价地获得。对于硬件实施而言,时域中的循环移位需要较大的缓冲器,以及相应的延迟和能耗。相应的后旋转需要作用于所有变换输出数据的高精度转动抵消器(de-rotator),从而导致芯片面积以及能耗的增大。
现有解决方案存在以上介绍的问题。现有技术循环移位要么在时域实施要么在频域实施。尽管这两种方法导致了不同的实现,但是它们的相同之处在于,和单独的FFT相比,面积、延迟和能耗增大到了显著的程度。
对于时域解决方案而言,FFT的输入数据在变换之前被循环移位。为了进行这种移位,如果FFT不是使用流水线方法实施的,则待变换的所有数据必须都可用。因而,就需要大小为N个字的额外缓冲器。即使对FFT使用了流水线方法,取决于这种结构能够处理的循环移位的量,也可能必须提供缓冲器来存储多达N-1个字。
在频域方法中,移位对应于每个输出数据的旋转。在这种情况下不需要缓冲器,但是该技术仍存在需要高精度旋转器的相关问题。此外,从FFT输出的数据可能是以位反转的顺序产生的,这意味着针对每个连续输出样本的旋转角度将表现为或多或少的随机性。因而,可能导致需要大表来存储这些角度。
鉴于上述讨论,很明显,现有技术方法面临着面积、延迟和能耗的显著负担。因此,希望提供利用数据的前移位/后旋转来实现FFT的改善技术和装置。
发明内容
应当强调,本说明书中使用的“包括”一词用于指定所陈述特征、整数、步骤或部件的存在,但是这些词的使用并不排除存在或添加一个或更多个其他特征、整数、步骤、部件或它们的组合。
根据本发明的一个实施方式,上述和其他目的在用于N频点(bin)FFT中的频率抽取(DIF)FFT级中实现,其中N是偶数。DIF FFT级包括交换逻辑(swap logic),该交换逻辑接收第一输入样本x(v)和第二输入样本x(v+N/2),并选择性地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,或者另选地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,其中0≤v<N/2。DIF FFT级还包括:用于对第一交换逻辑输出端口和第二交换逻辑输出端口供应的值进行相加的求和单元;用于对第一交换逻辑输出端口和第二交换逻辑输出端口供应的值进行相减的求差单元;以及将求差单元供应的值乘以旋转因子WN (v+s)mod(N/2)的旋转因子逻辑,其中s是表示N个输入样本的循环移位量的整数。
另一方面,FFT处理器包括上述DIF FFT级和一种逻辑,该逻辑通过当(v+s)mod N<N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)mod N≥N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,来控制该交换逻辑。
在另选实施方式中,FFT处理器包括上述DIF FFT级和一种逻辑,该逻辑通过当(v+s)∧N/2=0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)∧N/2≠0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,来控制该交换逻辑,其中∧表示位逻辑AND运算。
附图说明
通过阅读以下结合附图的详细描述将理解本发明的这些目的和优点,附图中:
图1是带有循环前缀的单时域OFDM符号。
图2是用于根据两个大小为N/2的变换来确定大小为N的傅里叶变换的布置。
图3是具有4个蝶形级(butterfly stage)的N=16的DIF FFT的示意图,每个蝶形级都有8个蝶形。
图4(a)例示了针对折叠FFT实施的基-2蝶形级的结构。
图4(b)例示了移位工作模式下的折叠FFT实施。
图4(c)例示了计算工作模式下的折叠FFT实施。
图5示意性地例示了包括N/2个蝶形操作的布置中的两个蝶形操作。
图6示意性地例示了在进行了对输入数据索引的循环左移位之后对图5的蝶形操作的影响。
图7示意性地例示了在进行了对输入数据索引的循环右移位之后对图5的蝶形操作的影响。
图8是根据本发明的观点的示例性基-2DIF FFT级的示意图。
具体实施方式
下面将参照附图来描述本发明的各个特征,附图中,使用相同的标号来表示相同的部分。
现将结合多个示例性实施方式更为详细地描述本发明的各个方面。为了有助于理解本发明,按照计算机系统或能够执行编程指令的其他硬件中的元件要执行的动作序列来描述本发明的很多方面。应当意识到,在每个实施方式中,各种动作可以由专用电路(例如,互连以实现专用功能的离散逻辑门)、由一个或更多个处理器所执行的程序指令,或者由二者的组合来执行。此外,本发明还可以被认为是整体实施在任意形式的计算机可读载体中,如固态存储器、磁盘、光盘或载波(如射频、音频或光频载波),这些计算机可读载体包含了可以使处理器执行这里描述的技术的适当计算机指令组。因而,本发明的各个方面可以按照很多不同形式来实施,且所有这些形式都被理解为落在本发明的范围之内。对于本发明的各个方面中的每一个,任意这种形式的实施方式在此都可以称为“被配置成执行所述动作的逻辑”或者另选地称为“执行所述动作的逻辑”。
本发明的一个方面是在同一操作中实现移位/旋转和FFT,从而获得硬件、时间和能量的更有效利用。下面详细描述的该方法适用于可以使用移位FFT的所有技术领域,包括但不限于图像处理,无线局域网(WLAN)、超宽带(UWB)通信、回波消除等。该方法还适用于所有类型的FFT结构,而与流水线、基或并行化(parallelization)无关。
基-2FFT算法是利用分治(divide and conquer)方法从离散傅里叶变换得出的。该算法有两种基本版本:一种是通过时域抽取得到的,另一种是通过频域抽取得到的。为了便于理解本发明的各个方面,现在开始描述基-2频域抽取(DIF)FFT算法的推导。
首先,令x(n)为n=0...N-1的样本序列,其中N是偶数,且优选是2的幂。然后,将序列x(n)的离散傅里叶变换(DFT)表示为X(k),其中k=0...N-1。使用旋转因子符号
WN=e-J2π/N
来简化表达,从x(n)到X(k)的DFT可以写为:
DFT N { x ( n ) } = X ( k ) = Σ n = 0 N - 1 x ( n ) W N kn . - - - ( 1 )
然后针对偶数和奇数频率k分别来分析公式(1)。对于偶数频率,该表达式为:
X ( 2 k ) = Σ n = 0 N - 1 x ( n ) W N 2 kn = { W N 2 k = W N / 2 kn } = Σ n = 0 N - 1 x ( n ) W N / 2 kn .
和原始公式(1)相比,和的旋转因子对于n=0...N-1的轮数完成了两倍。因此,和被分为从0到N/2-1的两个半程和(halfrange sum),
X ( 2 k ) = Σ n = 0 N / 2 - 1 ( x ( n ) + x ( n + N / 2 ) ) W N / 2 kn
通过在0≤n<N/2的情况下定义u(n)=x(n)+x(n+N/2),该表达式比较容易理解。实际上,得出u(n)的大小为N/2的FFT为
X ( 2 k ) = Σ n = 0 N / 2 - 1 u ( n ) W N / 2 kn = DFT N / 2 { u ( n ) } .
现在将侧重点改为变换的奇数频率,并且应用相同的技术,给出:
X ( 2 k + 1 ) = Σ n = 0 N - 1 x ( n ) W N ( 2 k + 1 ) n = Σ n = 0 N - 1 x ( n ) W N 2 kn W N n = Σ n = 0 N - 1 ( W N n x ( n ) ) W N / 2 kn .
同样,与公式(1)相比,旋转因子在n从0到N-1变化时进行了两轮,并且与上述相似,和被分成从0到N/2-1的两个半程和,
X ( 2 k + 1 ) = Σ n = 0 N / 2 - 1 ( W N n x ( n ) + W N n + N / 2 x ( n + N / 2 ) ) W N / 2 kn ,
它等于
X ( 2 k + 1 ) = Σ n = 0 N / 2 - 1 ( W N n { x ( n ) - x ( n + N / 2 ) } ) W N / 2 kn .
为n=0至N/2-1引入新的辅助变量 v ( n ) = W N n ( x ( n ) - x ( n + N / 2 ) ) , 得出
X ( 2 k + 1 ) = Σ n = 0 N / 2 - 1 v ( n ) W N / 2 kn = DF T N / 2 { v ( n ) } .
结论是,信号x(n)的大小为N的FFT被分割成计算序列u(n)和v(n)的两个大小为N/2的FFT,其中u和v是x的简单函数
DFTN{x(n)}→DFTN/2{u(n)}且DFTN/2{v(n)},其中
u(n)=(x(n)+x(n+N/2)),
v ( n ) = W N n ( x ( n ) - x ( n + N / 2 ) ) .
图2中例示了用于根据两个大小为N/2的FFT来确定大小为N的FFT的布置。该示意图中根据x来生成u和v的布置因为图中的加减结构(包围在虚线框201中)的形状而一般被称为“FFT蝶形”。
只要被分割成两个半大小变换的变换是相等的,就可以重复应用这种分治方法。每次重复总会导致一组N/2个新蝶形。这种组被称为蝶形级。为说明这一点,图3示出了具有4个蝶形级的N=16的DIF FFT的示意图,每个蝶形级都有8个蝶形。注意输出怎样以位反转的顺序产生:假设输入顺序为0,1,2,3,...,13,14,15(以二进制表示为0000,0001,0010,0011,...,1101,1110,1111),则输出以顺序0,8,4,12,...,11,7,15(以二进制表示为0000,1000,0100,1100,...,1011,0111,1111)出现。
还注意到,在每次将一个变换独立地分成两个半大小的变换时,N不必是2的幂。然而,如果完整的变换的大小不是2的幂,则在某些级N将是偶数,但N/2不是,于是使用上述技术不能进一步分割该变换。
现在来描述前移位和后旋转的效果。如前所述,时间离散信号x(n)(n=0...N-1)的DFT被定义为:
X ( k ) = Σ n = 0 N - 1 x ( n ) W N kn ,
其中k=0...N-1,令xs(n)是信号x(n)向左移位了s个循环步骤的结果,则
xs(n)=x((n+s)mod N).
现在,xs(n)的傅里叶变换由下式定义:
X s ( k ) = Σ n = 0 N - 1 x s ( n ) W N kn = Σ n = 0 N - 1 x ( ( n + s ) mod N ) W N kn
= { l = n + s } = Σ l = s N - 1 x ( l ) W N k ( l - s ) + Σ l = 0 s - 1 x ( l ) W N k ( l - s )
= Σ l = 0 N - 1 x ( l ) W N k ( l - s ) = W N - ks Σ l = 0 N - 1 x ( l ) W N kl .
代入n=l,得到:
X s ( k ) = Σ n = 0 N - 1 x s ( n ) W N kn = W N - ks Σ n = 0 N - 1 x ( n ) W N kn
它是x(n)乘以旋转因子WN -ks的变换。因而,时域中的移位对应于频率中的旋转。对于逆傅里叶变换,情况相反。
FFT的硬件实施可以基于其数据流向图(data-flowgraph)的直接映射,如图3所示。然而,这种结构将需要支持硬件,在大多数情况下,这些支持硬件太快,因为它需要对于每次变换所有的输入数据都并行可用。换句话说,为了利用根据数据流向图直接设计出的FFT,数据必须以比实际传输速度快N倍的步调到达。如果数据速率较慢,则硬件将面临明显的等待输入的空闲时间,造成该实施不必要的芯片面积消耗。
对于很多数字信号处理应用而言,数据是以一个样本接着一个样本的串行方式到达的。例如,对于数字化语音信号或来自于天线的采样无线基带信号,情况就是这样。并行FFT实施对于这种应用明显太快。为了避免这种失配,数据流向图可以按这样的方式折叠:对于要实施的每个蝶形级仅需要一个蝶形处理器。
图4(a)中例示了针对这种折叠FFT实施的基-2蝶形级400的结构。大小为N的FFT由串联的log2(N)个这种级构成。基-2蝶形级400具有延迟线401和可选择的数据路径逻辑403。通过基-2蝶形级400的数据路径具有两种工作模式:移位模式(图4(b)中所示)和计算模式(图4(c)所示)。
在移位模式下,数据路径逻辑403被配置成从输入接收N/2个样本,并将它们馈入到延迟线401中。同时,数据路径逻辑403接收延迟线401的内容,并将其作为基-2蝶形级400的输出进行供应。因而,该模式下的基-2蝶形级400仅充当延迟线。
在计算模式下,数据路径逻辑403根据当前输入的样本和从延迟线移位输出的数据来计算蝶形操作。蝶形操作的一个输出被作为基-2蝶形级400的输出进行供应,而蝶形操作的另一输出被馈入到延迟线401的输入。
在本发明的一方面,对标准FFT处理器进行改动,使得它独立于作用于输入数据的循环移位量而计算出相同的结果。这种改动包括通过添加输入数据交换器以及向旋转因子索引添加偏移来改变蝶形处理器元件。包括某种简单控制逻辑以确定旋转因子偏移应该是多少以及应该何时交换输入数据。现在将对此进行更为详细的描述。
为了有助于理解本发明的各个方面,我们分析了针对N频点FFT的示例性单蝶形级。考虑N个输入样本x(0),...x(N-1)的索引0...N-1以及执行FFT所需的相应N/2个蝶形,图5中示意性示出了其蝶形0和1(“BF0”和“BF1”)。从FFT算法的推导明显得出,每个蝶形BFi都具有分离开N/2个输入样本的输入;即,BFi以x(i)和x(i+N/2)作为其输入。利用这种布置,能够适当地计算FFT算法,因为BFi的输出能够重写变量x(i)和x(i+N/2)。
在输入数据索引的循环左移位之后,情况如图6所示。如果希望如同没有发生循环左移位那样来计算FFT,则蝶形操作也必须移位。例如,之前接收与索引0和N/2相对应的输入的最左侧蝶形此时被馈入了与索引1和N/2+1相对应的输入。因此,该蝶形不是充当BF0,而是变为充当BF1
类似地,之前接收与索引N/2-1和N-1相对应的输入的最右侧蝶形此时被馈入了与索引N/2和0相对应的输入。因此,该蝶形级的功能不是充当BFN/2-1,而是变得与其输入很好地匹配。然而,可以观察到,因为样本x(0)从最左侧交换到了最右侧位置,所以这两个样本x(0)和x(N/2)的输入顺序与标准BF0蝶形相比是相反的。为了表示输入的这种反转,蝶形被表示为BF0’。
对于输入数据的循环右移位,情况如图7所描绘。可以看出,因为输入的变化,最左侧的蝶形可以表示为BF′N/2-1,而其右边的邻居变成BF0。剩余的蝶形相应地变化,最右侧的蝶形变成BFN/2-2
图8是根据本发明的方面的示例性基-2DIF FFT级800(此后,称为“级800”)的示意图。为了完成所需的循环FFT功能,级800在若干方面与常规蝶形级不同。首先,两个输入被供应到两个复用器801、803的相应一个。这使得能够根据蝶形级是处于常规BFx位置(没有交换必要)还是处于BFx’位置(必需进行交换以撤销输入的反转)来有条件地交换输入。
用于复用器801、803的控制逻辑(未示出)可以如下运行:假定两个输入x(v)和x(v+N/2),其中0≤v<N/2,如果(v+s)mod N≥N/2,或等价地(当N是2的幂时),如果(v+s)∧N/2≠0,则复用器801、803会交换输入,其中s是移位数,而∧表示位逻辑AND操作。否则,不进行交换。s的符号可用于控制移位的方向,但是为方便起见,可以将s限制在0≤s<N的范围内,因为一个方向(例如,左)上的s位的移位等于相反方向(例如,右)上N-s位的移位。
第一复用器801的输出被供应到求和单元805的第一输入,还被供应到求差单元807的第一输入。类似地,最后一个复用器803的输出被供应到求和单元805的第二输入,还被供应到求差单元807的第二输入。第二求差单元807的输出被供应到旋转因子逻辑809。
级800不同于常规蝶形级的另一方面在于旋转因子逻辑809中使用的旋转乘法系数偏移了因子s,从而其等于WN (v+s)mod(N-2)。如果N是2的幂,则取模运算变得没有成本,因为它导致二进制数系统中的线回(wrap-around)。
级800供应了两个输出:第一个输出x′(v)由求和单元805供应。第二个输出x′(v+N/2)由旋转因子逻辑809供应。
多个级800可以按照图3所示的方式互连以形成完整的DIF FFT单元。另选的是,可以将级800并入到如图4(a)所示的设计中(用作数据路径逻辑403的一部分)以生成如上所述的流水线基-2FFT结构。
级800所示例的蝶形级使得能够建立这样的FFT处理器(这里称为“循环FFT处理器”),该处理器直接计算前移位或等价的后旋转数据的变换而无需单独的前移位或后旋转电路。循环FFT处理器在很多应用中是有用的。这些应用包括但不限于:补偿OFDM接收器中的失位FFT窗口;以及抵消由延迟扩展较大的信道导致的旋转。
根据本发明的设计提供了优于常规技术的很多优点。例如,这些设计不会面临附加的面积、能量和时间增长,如果应用现有技术的前移位或后旋转方法来执行相同的操作则将会导致这种附加的面积、能量和时间增长。因为移位是在FFT处理器内部进行的,所以处理器的整体控制也被简化,且除了传达移位量的输入信号s以外,不需要添加附加的块或引线。对于系统设计者而言,附加的移位功能隐藏在FFT硬件中。
根据本发明的循环FFT在这样的情况下尤其有用:样本并不是以恰当的顺序x(0),x(1),x(2),...x(N-1)到达的,而是循环移位了某个量s,因而以x(s),x(s+1),x(s+2),...x(N-1),x(0),x(1),...x(s-1)的顺序到达。例如,这可能在OFDM接收器中存在不恰当放置的FFT窗口时发生。而且,可能由于通过延迟扩展较大的信道的传播而造成输入样本的不希望的循环移位。此处教示的循环FFT可用于解决这些问题,因为它能够如同输入样本没有任何移位那样生成相同的FFT输出。从概念上来讲,它就好像在计算之前FFT处理器已经在相反的方向使输入循环移位了量s,尽管如图所示,不需要实际的移位来完成这点。这些实施方式将包括用于确定输入样本的循环移位量和控制FFT处理器来使用与输入样本的循环移位量相同但是方向相反的s值的逻辑。
已经参照特定实施方式描述了本发明,然而,对于本领域技术人员而言,很明显,可以按照不同于上述实施方式的特定形式来实施本发明。所述实施方式仅是说明性的且不应被认为是任意形式的限制。本发明的范围由所附权利要求而不是前面的描述给出,且此处旨在涵盖落在权利要求范围内的所有的变型和等价形式。

Claims (12)

1、一种在N频点FFT中使用的频率抽取(DIF)快速傅里叶变换(FFT)级,其中N是整数,该DIF FFT级包括:
交换逻辑,其接收第一输入样本x(v)和第二输入样本x(v+N/2),并选择性地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,或者另选地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,其中0≤v<N/2;
求和单元,用于对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相加;
求差单元,用于对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相减;以及
旋转因子逻辑,其将该求差单元所供应的值乘以旋转因子WN (v+s)mod(N/2),其中s是表示N个输入样本的循环移位量的整数。
2、一种快速傅里叶变换(FFT)处理器,该FFT处理器包括:
权利要求1所述的DIF FFT级;以及
逻辑,其通过当(v+s)modN<N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)modN≥N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,来控制该交换逻辑。
3、一种正交频分复用(OFDM)接收器,该OFDM接收器包括:
权利要求2所述的FFT处理器;以及
逻辑,其确定这些输入样本的循环移位量并控制该FFT处理器来使用与这些输入样本的旋转量相等但是方向相反的s值。
4、根据权利要求3所述的OFDM接收器,其中这些输入样本的循环移位是由不恰当放置的FFT窗口所造成的。
5、根据权利要求3所述的OFDM接收器,其中这些输入样本的循环移位是由这些输入样本通过具有较大延迟扩展的信道而传播所造成的。
6、一种快速傅里叶变换(FFT)处理器,该FFT处理器包括:
权利要求1所述的DIF FFT级;以及
逻辑,其通过当(v+s)∧N/2=0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)∧N/2≠0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,来控制该交换逻辑,其中∧表示位逻辑AND运算。
7、一种执行在N频点傅里叶变换中使用的频率抽取(DIF)快速傅里叶变换(FFT)的方法,其中N是偶数,该DIF FFT方法包括以下步骤:
接收第一输入样本x(v)和第二输入样本x(v+N/2),并选择性地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,或者另选地在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,其中0≤v<N/2;
对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相加;
通过对第一交换逻辑输出端口和第二交换逻辑输出端口所供应的值进行相减而生成差值;以及
将该差值乘以旋转因子WN (v+s)mod(N/2),其中s是表示N个输入样本的循环移位量的整数。
8、根据权利要求7所述的方法,该方法包括以下步骤:
当(v+s)modN<N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)modN≥N/2时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本。
9、根据权利要求8所述的方法,该方法包括以下步骤:
确定这些输入样本的循环移位量,并将s设置为与这些输入样本的旋转量相等但方向相反的值。
10、根据权利要求9所述的方法,其中这些输入样本的循环移位是由不恰当放置的FFT窗口所造成的。
11、根据权利要求9所述的方法,其中这些输入样本的循环移位是由这些输入样本通过具有较大延迟扩展的信道而传播所造成的。
12、根据权利要求7所述的方法,该方法包括以下步骤:
当(v+s)∧N/2=0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第一输入样本和第二输入样本,而当(v+s)∧N/2≠0时在相应的第一交换逻辑输出端口和第二交换逻辑输出端口供应第二输入样本和第一输入样本,其中∧表示位逻辑AND运算。
CN2006800472468A 2005-12-14 2006-12-12 快速傅里叶变换级电路、处理器和方法及正交频分复用接收器 Expired - Fee Related CN101331479B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11/302,243 2005-12-14
US11/302,243 US7685220B2 (en) 2005-12-14 2005-12-14 Circular fast fourier transform
PCT/EP2006/011935 WO2007068438A2 (en) 2005-12-14 2006-12-12 Circular fast fourier transform

Publications (2)

Publication Number Publication Date
CN101331479A true CN101331479A (zh) 2008-12-24
CN101331479B CN101331479B (zh) 2012-01-11

Family

ID=38051780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006800472468A Expired - Fee Related CN101331479B (zh) 2005-12-14 2006-12-12 快速傅里叶变换级电路、处理器和方法及正交频分复用接收器

Country Status (8)

Country Link
US (1) US7685220B2 (zh)
EP (1) EP1960905A2 (zh)
JP (1) JP4903810B2 (zh)
KR (1) KR101229648B1 (zh)
CN (1) CN101331479B (zh)
BR (1) BRPI0619920A2 (zh)
TW (1) TW200741487A (zh)
WO (1) WO2007068438A2 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102152910A (zh) * 2011-04-27 2011-08-17 山东电力研究院 滤膜存放箱
CN102955760A (zh) * 2011-08-23 2013-03-06 上海华魏光纤传感技术有限公司 基于dif的基-2并行fft处理器及其处理方法
CN107913083A (zh) * 2016-10-06 2018-04-17 通用电气公司 用于超声复用的系统和方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005045519A1 (de) * 2005-09-23 2007-03-29 Newlogic Technologies Ag Verfahren und Vorrichtung zur FFT Berechnung
DE602006013964D1 (de) * 2006-07-18 2010-06-10 Ericsson Telefon Ab L M Faltung von Eingabedatenwerten für eine Transformation
CN101930426B (zh) * 2009-06-24 2015-08-05 华为技术有限公司 信号处理方法、数据处理方法及装置
CN102592601B (zh) 2011-01-10 2014-09-17 华为技术有限公司 一种信号处理方法和装置

Family Cites Families (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3638004A (en) * 1968-10-28 1972-01-25 Time Data Corp Fourier transform computer
US3767905A (en) * 1971-05-12 1973-10-23 Solartron Electronic Group Addressable memory fft processor with exponential term generation
US3927312A (en) * 1974-10-31 1975-12-16 Us Army Vector rotator
US4327367A (en) * 1978-04-24 1982-04-27 Electric Power Research Institute, Inc. Thyristor with even turn-on line potential and method with 1-micron to 5-mil wide alignment region band
US4231102A (en) * 1978-12-21 1980-10-28 Raytheon Company Cordic FFT processor
US4563750A (en) * 1983-03-04 1986-01-07 Clarke William L Fast Fourier transform apparatus with data timing schedule decoupling
FR2584213B1 (fr) * 1985-06-28 1994-03-11 Thomson Csf Dispositif de calcul d'une transformee de fourier discrete, glissante, et son application a un systeme radar.
FR2587819B1 (fr) * 1985-09-24 1989-10-06 Thomson Csf Dispositif de calcul d'une transformee de fourier discrete, glissante et non recursive, et son application a un systeme radar
FR2588680B1 (fr) * 1985-10-16 1989-08-25 Thomson Csf Dispositif de calcul d'une transformee de fourier discrete, et son application a la compression d'impulsion dans un systeme radar
US5034910A (en) * 1988-03-14 1991-07-23 E-Systems, Inc. Systolic fast Fourier transform method and apparatus
JPH05158967A (ja) 1991-12-05 1993-06-25 Nec Corp 計算機システム
US5371696A (en) * 1992-12-24 1994-12-06 Sundararajan; Duraisamy Computational structures for the fast Fourier transform analyzers
CN1117681A (zh) * 1994-03-07 1996-02-28 现代电子产业株式会社 定时相位探测器及其定时回复
JP3541526B2 (ja) * 1995-10-24 2004-07-14 日本ビクター株式会社 周波数分割多重信号生成方法及び復号方法
US5717620A (en) * 1995-10-24 1998-02-10 Airnet Communications Corporation Improved-accuracy fast-Fourier-transform butterfly circuit
JP3821316B2 (ja) * 1996-08-06 2006-09-13 ソニー株式会社 演算装置および方法
JP3797400B2 (ja) 1997-05-23 2006-07-19 ソニー株式会社 演算装置および方法
US6370188B1 (en) * 1999-03-31 2002-04-09 Texas Instruments Incorporated Phase and frequency offset compensation in a telecommunications receiver
FR2794921B1 (fr) * 1999-06-14 2001-09-14 St Microelectronics Sa Procede et dispositif de transformation de donnees reelles en symboles complexes, notamment pour la reception de porteuses modulees en phase et en amplitude et transmises sur une ligne telephonique
US6658441B1 (en) * 1999-08-02 2003-12-02 Seung Pil Kim Apparatus and method for recursive parallel and pipelined fast fourier transform
US6772181B1 (en) * 1999-10-29 2004-08-03 Pentomics, Inc. Apparatus and method for trigonometric interpolation
US6567042B2 (en) * 2000-12-04 2003-05-20 The United States Of America As Represented By The Secretary Of The Air Force Acquisition through circular correlation by partition for GPS C/A code and P(Y) code
US6633616B2 (en) * 2001-02-21 2003-10-14 Magis Networks, Inc. OFDM pilot tone tracking for wireless LAN
US7020218B2 (en) * 2001-06-18 2006-03-28 Arnesen David M Sliding-window transform with integrated windowing
US7020222B2 (en) * 2001-10-24 2006-03-28 Texas Instruments Incorporated Efficient method and system for offset phasor determination
US20030182342A1 (en) * 2002-03-25 2003-09-25 Murphy Charles Douglas Low-cost multiplication in small DFT butterflies
CN1180592C (zh) * 2002-04-01 2004-12-15 北京六合万通微电子技术有限公司 快速傅立叶变换和反变换的模拟数字混合电路及其在通信系统中的应用
CN100563226C (zh) * 2002-06-27 2009-11-25 三星电子株式会社 利用混合基数快速付里叶变换的调制设备
JP2004032568A (ja) * 2002-06-28 2004-01-29 Hitachi Kokusai Electric Inc 相関検出装置およびフーリエ変換装置
US7082451B2 (en) * 2002-09-09 2006-07-25 Freescale Semiconductor, Inc. Reconfigurable vector-FFT/IFFT, vector-multiplier/divider
US20050015420A1 (en) * 2003-07-18 2005-01-20 Gibb Sean G. Recoded radix-2 pipeline FFT processor
US8266196B2 (en) * 2005-03-11 2012-09-11 Qualcomm Incorporated Fast Fourier transform twiddle multiplication
US8229014B2 (en) * 2005-03-11 2012-07-24 Qualcomm Incorporated Fast fourier transform processing in an OFDM system
JP4698394B2 (ja) * 2005-11-25 2011-06-08 パナソニック株式会社 高速フーリエ変換回路

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102152910A (zh) * 2011-04-27 2011-08-17 山东电力研究院 滤膜存放箱
CN102152910B (zh) * 2011-04-27 2013-04-10 山东电力研究院 滤膜存放箱
CN102955760A (zh) * 2011-08-23 2013-03-06 上海华魏光纤传感技术有限公司 基于dif的基-2并行fft处理器及其处理方法
CN107913083A (zh) * 2016-10-06 2018-04-17 通用电气公司 用于超声复用的系统和方法

Also Published As

Publication number Publication date
CN101331479B (zh) 2012-01-11
KR20080078062A (ko) 2008-08-26
WO2007068438A2 (en) 2007-06-21
TW200741487A (en) 2007-11-01
EP1960905A2 (en) 2008-08-27
JP2009519637A (ja) 2009-05-14
WO2007068438A3 (en) 2008-06-05
US7685220B2 (en) 2010-03-23
BRPI0619920A2 (pt) 2011-10-25
US20070133389A1 (en) 2007-06-14
JP4903810B2 (ja) 2012-03-28
KR101229648B1 (ko) 2013-02-04

Similar Documents

Publication Publication Date Title
CN101331479B (zh) 快速傅里叶变换级电路、处理器和方法及正交频分复用接收器
US7693034B2 (en) Combined inverse fast fourier transform and guard interval processing for efficient implementation of OFDM based systems
CN103970718B (zh) 一种快速傅里叶变换实现装置及方法
CN105635022B (zh) 一种偏移正交多载波基带系统
CN103733516A (zh) 接收机、发射机以及用于数字多子频带处理的方法
KR20100014317A (ko) Ofdm 시스템용 보간 방법과 채널 추정 방법 및 장치
WO1999031603A1 (en) Computationally efficient analysis and synthesis of real signals using discrete fourier transforms and inverse discrete fourier transforms
US20110164490A1 (en) Fast fourier transform and inverse fast fourier transform (fft/ifft) operating core
CN101789918B (zh) 并行信道均衡方法
CN105095152A (zh) 一种可配置的128点fft装置
CN101937332A (zh) 基于基24算法的多路fft处理器中乘法器的复用方法
KR100720949B1 (ko) 직교 주파수 분할 다중화 시스템에서의 고속 푸리에 변환프로세서 및 그 변환 방법
Srinivasaiah et al. Low power and area efficient FFT architecture through decomposition technique
CN103873219A (zh) 一种prach基带信号的idft及子载波映射实现方法
CN103457895B (zh) 一种ofdm传输方法、装置和传输系统
CN103188192A (zh) 应用于视频传感器的基带处理装置
WO2014127663A1 (zh) 一种插值滤波方法及插值滤波器
CN102023963B (zh) 高速多模式时频域变换方法
US8977885B1 (en) Programmable logic device data rate booster for digital signal processing
Yuan et al. A 256-point dataflow scheduling 2× 2 MIMO FFT/IFFT processor for IEEE 802.16 WMAN
Zhang et al. A novel reconfigurable architecture for generic OFDM modulator based on FPGA
CN111291315A (zh) 一种数据处理方法、装置及设备
US6751642B1 (en) Process and device for transforming complex symbols into real data, especially for the transmission over a telephone line of phase-modulated and amplitude-modulated carriers
CN107454030B (zh) 一种电力线宽带载波半并行发射机及其实现方法
CN107979556A (zh) 信号处理方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: CLUSTER CO., LTD.

Free format text: FORMER OWNER: LM ERICSSON TELEPHONE CO., LTD.

Effective date: 20150420

Owner name: OPTIS WIRELESS TECHNOLOGY LLC

Free format text: FORMER OWNER: CLUSTER CO., LTD.

Effective date: 20150420

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20150420

Address after: Texas, USA

Patentee after: Telefonaktiebolaget LM Ericsson (publ)

Address before: Delaware

Patentee before: Clastres LLC

Effective date of registration: 20150420

Address after: Delaware

Patentee after: Clastres LLC

Address before: Stockholm

Patentee before: Telefonaktiebolaget LM Ericsson (publ)

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120111

Termination date: 20171212

CF01 Termination of patent right due to non-payment of annual fee