CN103870437A - 数字信号处理装置及其处理方法 - Google Patents
数字信号处理装置及其处理方法 Download PDFInfo
- Publication number
- CN103870437A CN103870437A CN201310056976.6A CN201310056976A CN103870437A CN 103870437 A CN103870437 A CN 103870437A CN 201310056976 A CN201310056976 A CN 201310056976A CN 103870437 A CN103870437 A CN 103870437A
- Authority
- CN
- China
- Prior art keywords
- rotation coefficient
- sub
- data
- butterfly
- signal processing
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- 238000007620 mathematical function Methods 0.000 claims description 10
- 230000005055 memory storage Effects 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (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)
- Discrete Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
一种数字信号处理装置及其处理方法。数字信号处理装置包括存储器、控制逻辑单元、蝴蝶运算器、选择器、第一转动系数产生器、第二转动系数产生器以及转动系数闩锁器。第一转动系数产生器与第二转动系数产生器分别提供第一子转动系数与第二子转动系数。将第一子转动系数与第二子转动系数通过蝴蝶运算器进行复数乘法运算以产生权重值(转动系数)。
Description
技术领域
本发明实施例是有关于一种信号处理技术,且特别是有关于一种数字信号处理装置及其处理方法。
背景技术
由于数字储存或计算机处理技术的流行,一般信号处理常用快速傅立叶转换(Fast Fourier Transform,FFT)来求得信号所对应的频谱,而离散傅立叶转换(Discrete Fourier Transform,DFT)提供了一个对连续傅立叶转换的离散逼近技术,并在许多领域中它是非常有用的。基于DFT演算法的快速傅立叶转换的应用范围相当广泛,例如雷达、无线通信、医疗成像、频谱分析和声学。
发明内容
有鉴于此,本发明实施例提出一种数字信号处理装置及其处理方法,通过傅立叶转换(Fourier Transform)的方式处理信号。
本发明实施例提出一种数字信号处理装置,其包括存储器、控制逻辑单元、蝴蝶运算器、选择器、第一转动系数产生器、第二转动系数产生器以及转动系数闩锁器。存储器具有数据输出端、数据输入端与地址端。控制逻辑单元电连接数据地址端。蝴蝶运算器电连接数据输出端、数据输入端与控制逻辑单元。选择器电连接蝴蝶运算器与控制逻辑单元。第一转动系数产生器与第二转动系数产生器电连接控制逻辑单元与选择器。第一转动系数产生器与第二转动系数产生器分别提供第一子转动系数与第二子转动系数。转动系数闩锁器电连接控制逻辑单元与蝴蝶运算器。转动系数闩锁器用以闩锁第一子转动系数与第二子转动系数通过蝴蝶运算器的复数乘法运算后的权重值。
在本发明的一实施例中,数字信号处理装置更包括第一数据闩锁器、第二数据闩锁器与数据多工器。第一数据闩锁器与第二数据闩锁器,电连接数据输出端、控制逻辑单元与蝴蝶运算器。数据多工器电连接蝴蝶运算器、转动系数闩锁器与控制逻辑单元。数据多工器根据控制逻辑单元的控制将运算结果写回至存储器中对应的位置。
在本发明的一实施例中,蝴蝶运算器包括复数乘法器、复数加法器与复数减法器。复数乘法器电连接选择器、转动系数闩锁器与数据多工器,复数乘法器进行前述的复数乘法运算并产生权重值。复数加法器电连接第一数据闩锁器、第二数据闩锁器与数据多工器。复数减法器电连接第一数据闩锁器、第二数据闩锁器与选择器。
在本发明的一实施例中,蝴蝶运算器执行基数为2的蝴蝶运算。
在本发明的一实施例中,若存储器储存N点数据且假设Q是表示对应N/2点所需的地址位元数,则时,任一权重值,其中对任意的k值,可以用一个Q-bit二进为表示法来表示 假设 且 ai与aj+R∈{1,0}、R为整数且0≤R≤(Q-1).则k可改写成权重值可改写成其中为第一子转动系数,为第二子转动系数。
在本发明的一实施例中,若存储器储存N点数据且假设Q是表示对应N/2点所需的地址位元数,则时,第一转动系数产生器与第二转动系数产生器的存储器空间分别储存个第一子转动系数与个第二子转动系数或是分别储存个第一子转动系数与个第二子转动系数,其中floor()为一数学函数而以无条件舍去方式来取得整数,而ceiling()为一数学函数而以无条件进位方式来取得整数。
本发明另一实施例提出一种数字信号处理方法,包括:提供储存N点数据的一存储器;提供第一子转动系数与第二子转动系数;将第一子转动系数与第二子转动系数通过蝴蝶运算器的复数乘法器进行复数乘法运算以产生权重值;将权重值应用至N点数据的快速傅立叶转换以产生多个经转置值;以及将这些经转置值写回至所述存储器。
为让本发明的上述特征和优点能更明显易懂,下文特举实施例,并配合所附图式作详细说明如下。
附图说明
下面的所附图式是本发明的说明书的一部分,绘示了本发明的示例实施例,所附图式与说明书的描述一起说明本发明的原理。
图1为快速傅立叶转换的蝴蝶运算示意图。
图2是说明16点快速傅立叶转换的信号处理流程图。
图3是说明一般单一处理单元的架构。
图4是说明一般基于蝴蝶运算的数字信号处理装置的架构图。
图5是依照本案实施例的数字信号处理装置的示意图。
图6是本案实施例的信号处理方法的流程图。
附图标号说明:
110:存储器 120:处理单元
122:蝴蝶运算器 124:控制逻辑单元
126:转动系数产生器 400、500:数字信号处理装置
510:存储器 520:蝴蝶运算器
522、524:数据闩锁器 526:复数加法器
528:复数减法器 530:复数乘法器
532:数据多工器 540:控制逻辑单元
550、560:转动系数产生器 570:选择器
580:转动系数闩锁器 addr:地址端
B1、B2:子转动系数 din:数据输入端
dout:数据输出端 W1:权重值(转动系数)
S01~S1、S601~S609:步骤
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
FFT是以DFT为基础来大量节省计算量。对于具有N点数据的输入向量x,进行离散傅立叶转换(Discrete Fourier Transform,DFT)之后为具有N点数据的输出向量X,离散傅立叶转换式子如等式(1)所示:
其中系数称为转动系数(twiddle factor)。
倘若根据上面的等式(1)来直接计算这些N点的DFT,直接演算所需的操作数是约N2个复数乘法和N*(N-1)个复数加法。显然,这是一个非常低效率且不切实际的方式。如果没有任何优化执行,DFT将因消耗巨大的计算资源而不能够应用于大部分的实际应用。
在下面的段落中,将以基数为2的分频(Decimation in Frequency,DIF)算法作为一个例子,显示如何显著降低DFT的计算。当N点数据的输入向量x为2的幂次方时,上述等式(1)的N点数据的输出向量可以分解成偶数和奇数部分如下,其中偶数频率的点可写成如等式(2),而奇数频率的点可写成如等式(3),
以及
上述等式(2)和等式(3)揭示频率输出X(k)的偶数索引值和奇数索引值可以通过较短的点数据来运算DFT。如图1所示,可以得到以“基数为2的分频”的快速傅立叶转换的蝴蝶运算(butterfly operation)示意图,其中转动系数可写成如等式(4)。
图2是说明16点分频的快速傅立叶转换的信号处理流程图。图3是说明一般单一处理单元的架构。请同时参阅图2与图3,信号处理流程图的左边是欲进行快速傅立叶转换运算的原始数据(在此以16点为例),右边则是完成快速傅立叶转换运算的目的数据。图中表示进行复数加法,表示进行复数减法,表示进行复数乘法。进行快速傅立叶转换运算的原始数据需事先存放于存储器110中。在进行第一步骤S0时,处理单元120依序自存储器110中读取数据对[x(0)、x(8)]、[x(1)、x(9)]、[x(2)、x(10)]、…、[x(7)、x(15)]以分别进行一次蝴蝶运算(含复数加法、复数减法与复数乘法运算),并且将每次运算结果写回至存储器110。等到第一步骤(step)S0完成时,处理单元120开始进行第二步骤S1,亦即依序自存储器110中读取第一步骤S0运算结果的数据对[x(0)、x(4)]、…、[x(3)、x(7)]、[x(8)、x(12)]、…、[x(11)、x(15)]以分别进行一次蝴蝶运算,并且将每次运算结果写回至存储器110。以此类推,以继续完成步骤S2与S3。因此即可以单一处理单元完成16点快速傅立叶转换运算。同样的方法可以类推至任意长度的快速傅立叶转换运算。
由于FFT发挥越来越重要的作用,许多实施FFT的研究提出了实现更快、更准确的DFT结果。无论什么样的实施方法所需的存储器与转动系数产生器(twiddle factor generator)的空间始终是评估实施过程中必须考虑的一个因素。图4是说明一般基于蝴蝶运算的数字信号处理装置的架构图。数字信号处理装置400包括存储器110及处理单元120,而处理单元120包括蝴蝶运算器(butterfly arithmetic unit)122、控制逻辑单元124以及转动系数产生器126。对于N点的FFT运算,存储器110储存N点数据,转动系数产生器126需储存N/2点的转动系数。例如,对于储存131072点数据的存储器110,转动系数产生器126需储存65536点的转动系数。由此可见,当FFT所运算的取样点增加时,存储器110与转动系数产生器126的尺寸也显著增加,同时转动系数产生器126的功耗也跟着增加。
其中ai,aj∈{1,0}、R为整数且0≤R≤(Q-1) (7)。
将等式(7)代入等式(6)的项次“k”中,得到下述的等式(8)和等式(9)。
从等式(7)、(8)和(9),任何给定的二进位数字k可分解成两部分,第一部份包含最低有效位元(least significant bit,LSB)的R个位元,而第二部份包含最高有效位元(most significant bit,MSB)的“Q-R”个位元,且对应的转动系数可通过两个子转动系数(sub-twiddle factor)进行复数乘法运算而得之。由于0≤Ar≤(2R-1)且0≤Ar'≤(2Q-R-1),所需存储器空间用来储存的和的点数分别为2R和2Q-R,所以转动系数的点数由原先的2Q个可以减少为“2R+2Q-R”个。
对于最佳的R值可如下列等式(10)进行微分运算。
其中floor()为数学函数而以无条件舍去方式来取得整数(round-downoperation)。在本发明另一实施例中,亦可令其中ceiling()为数学函数而以无条件进位方式来取得整数(round-up operation)。
基于等式(11),若N点的FFT数据且时,所需转动系数产生的点数由原先的2Q个减少为“”个。例如,对于8192点的FFT数据,原始转动系数的点数为4096个,但通过分解为两个子转动系数的方式,所需的点数为 个,因此对于整体的转动系数的点数减少32倍(4096/128=32),故可大幅减少所需的存储器空间。
基于将转动系数分解为子转动系数的运算架构,图5是依照本发明一实施例的数字信号处理装置500的示意图。请参阅图5。数字信号处理装置500包括存储器510、蝴蝶运算器520、控制逻辑单元540、选择器570、转动系数产生器550和560以及转动系数闩锁器580。
在此值得一提的是,若存储器510储存N点数据且时,根据等式(7)至等式(9),转动系数产生器550针对最低有效位元的R个位元可提供子转动系数B1,子转动系数而转动系数产生器560针对最高有效位元的“Q-R”个位元可提供子转动系数B2,子转动系数 其中 ai与aj{1,0}、R为整数且0≤R≤(Q-1)。
此外,根据等式(11),若存储器510储存N点数据且时,转动系数产生器550的存储器空间可预先储存子转动系数B1的点数为个,子转动系数另外,转动系数产生器560的存储器空间可预先储存子转动系数B2的点数为个,子转动系数其中floor()为数学函数而以无条件舍去方式来取得整数。
承上述,存储器510具有数据输入端din、数据输出端dout与地址端addr。控制逻辑单元540电连接数据地址端addr。蝴蝶运算器520电连接数据输出端dout、数据输入端din与控制逻辑单元540。选择器570电连接蝴蝶运算器520与控制逻辑单元540。两个转动系数产生器550和560电连接控制逻辑单元540与选择器570。转动系数闩锁器580电连接控制逻辑单元540与蝴蝶运算器520。
经由控制逻辑单元540的控制,当子转动系数B1与子转动系数B2通过蝴蝶运算器520的复数乘法运算后,转动系数闩锁器580可闩锁B1和B2经复数乘法运算后的权重值W1(转动系数)。倘若根据等式(9),则权重值W1的通式可表示为
另外,数字信号处理装置500还包括数据闩锁器522、数据闩锁器524与数据多工器532。蝴蝶运算器520可包括复数乘法器530、复数加法器526与复数减法器528。复数乘法器530电连接选择器570与转动系数闩锁器580。数据多工器532电连接复数加法器526、复数乘法器530、控制逻辑单元540与数据输入端din。复数加法器526电连接数据闩锁器522和524。复数减法器528电连接数据闩锁器522和524与选择器570。选择器570例如可为乘法输入选择器,依据控制信号将输入信号输出至复数乘法器530以进行后续的运算。数据闩锁器522与数据闩锁器524电连接数据输出端dout与控制逻辑单元540与蝴蝶运算器520。数据闩锁器522和/或数据闩锁器524通过数据输出端dout由存储器510读取数据,再输出至复数加法器526和/或复数减法器528进行运算。
此外,在一变化实施例中亦可不需设置数据闩锁器522和/或数据闩锁器524。举例来说,可通过设置两个数据输出端dout的方式,经由一控制逻辑单元(例如控制逻辑单元540或另行设置的控制逻辑单元)产生例如两个存储器存取地址,由存储器510中对应的位置读取数据并通过数据输出端dout输出至复数加法器526和/或复数减法器528。
控制逻辑单元540可用来进行如下的相关控制:(a)产生存储器510的存取地址,以从存储器510中对应的位置读出数据及/或将数据写回至存储器510中对应的位置;(b)使转动系数产生器550和560根据存取地址分别产生子转动系数B1和B2;(c)控制选择器570中的两种输入对,以对其中一种输入对以在复数乘法器530进行复数乘法运算,例如选择器570具有两种输入对,第一种输入对是由子转动系数B1与B2进行复数乘法运算得到权重值W1,第二种输入对是由复数减法器528的运算结果与权重值W1进行复数乘法运算;(d)对转动系数闩锁器580的权重值W1进行闩锁控制;(e)对数据多工器532的运算结果写回至存储器510;以及(f)对存储器510的数据进行闩锁控制。
在此值得一提的是,经由运算得到权重值W1、或是由复数减法器528的运算结果与权重值W1进行复数乘法运算可通过同一个复数乘法器530。由于复数乘法器530的电路复杂,而本实施例可通过共用一个复数乘法器530的方式可避免数字信号处理装置500的增加庞大电路面积。
此外,数字信号处理装置500中的蝴蝶运算器520可搭配数据闩锁器522和524进行如图2的流程,以执行基数为2的蝴蝶运算。数据闩锁器522和524自存储器510中读取数据对[x(0)、x(8)]、[x(1)、x(9)]、[x(2)、x(10)]、…、[x(7)、x(15)]以分别进行一次蝴蝶运算(含复数加法、复数减法与复数乘法运算),并且数据多工器532将每次运算结果写回至存储器510中对应的位置。等到第一步骤(step)S0完成时,数字信号处理装置500开始进行第二步骤S1,亦即依序自存储器510中读取第一步骤S0运算结果的数据对[x(0)、x(4)]、…、[x(3)、x(7)]、[x(8)、x(12)]、…、[x(11)、x(15)]以分进行一次蝴蝶运算,并且将每次运算结果写回至存储器510。以此类推,以继续完成步骤S2与S3。基此教示内容,本领域技术人员应当可基于同样的方法推演/类推至任意长度的快速傅立叶转换运算,故而在此并不再加以赘述之。
基于上述实施例所揭示的内容,可以汇整出一种通用的数字信号处理方法。更清楚来说,图6绘示为本案实施例的数字信号处理方法的流程图。请合并参阅图5和图6,本实施例的数字信号处理方法可以包括以下步骤。
如步骤S601所示,提供储存N点数据的存储器510。
如步骤S603所示,提供(第一)子转动系数B1与(第二)子转动系数B2,若时,转动系数产生器550的存储器空间可预先储存子转动系数B1的点数为个,子转动系数而转动系数产生器560的存储器空间可预先储存子转动系数B2的点数为个,子转动系数其中floor()为数学函数而以无条件舍去方式来取得整数。
如步骤S605所示,将子转动系数B1与子转动系数B2通过蝴蝶运算器520的复数乘法器530进行复数乘法运算以产生权重值W1。
如步骤S607所示,将权重值W1应用至N点数据的快速傅立叶转换以产生多个经转置值(亦即经由蝴蝶运算的每次运算结果)。
如步骤S609所示,将这些经转置值写回至存储器510。
综上所述,本案的数字信号处理装置500可使整体的转动系数的点数骤减,而欲求得的转动系数可通过将转动系数产生器550和560的输出做一次复数乘法运算而得之。显然地,本发明实施例因采用指数特性将转动系数拆解成两种转动系数,再以快速傅立叶转换架构中的蝴蝶运算器算出此转动系数而不用增加庞大的额外电路,故不会大幅增加电路面积,另一方面可以有效地减少储存所需用于蝴蝶运算的转动系数的数量。
虽然本发明已以实施例揭露如上,然其并非用以限定本发明,任何所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作些许更动与润饰,故本发明的保护范围当视权利要求范围所界定者为准。
Claims (10)
1.一种数字信号处理装置,其特征在于,包括:
一存储器,具有一数据输出端、一数据输入端与一地址端;
一控制逻辑单元,电连接所述数据地址端;
一蝴蝶运算器,电连接所述数据输出端、所述数据输入端与所述控制逻辑单元;
一选择器,电连接所述蝴蝶运算器与所述控制逻辑单元;
一第一转动系数产生器与一第二转动系数产生器,电连接所述控制逻辑单元与所述选择器,所述第一转动系数产生器与所述第二转动系数产生器分别提供一第一子转动系数与一第二子转动系数;以及
一转动系数闩锁器,电连接所述控制逻辑单元与所述蝴蝶运算器,所述转动系数闩锁器用以闩锁所述第一子转动系数与所述第二子转动系数通过所述蝴蝶运算器的复数乘法运算后的一权重值。
2.如权利要求1所述的数字信号处理装置,其特征在于,更包括:
一第一数据闩锁器与一第二数据闩锁器,电连接所述数据输出端、所述控制逻辑单元与所述蝴蝶运算器;以及
一数据多工器,电连接所述蝴蝶运算器、所述转动系数闩锁器与所述控制逻辑单元,所述数据多工器根据所述控制逻辑单元的控制将运算结果写回至所述存储器中对应的位置。
3.如权利要求2所述的数字信号处理装置,其特征在于,所述蝴蝶运算器包括:
一复数乘法器,电连接所述选择器、所述转动系数闩锁器与所述数据多工器,所述复数乘法器进行所述复数乘法运算并产生所述权重值;
一复数加法器,电连接所述第一数据闩锁器、所述第二数据闩锁器与所述数据多工器;以及
一复数减法器,电连接所述第一数据闩锁器、所述第二数据闩锁器与所述选择器。
4.如权利要求1所述的数字信号处理装置,其特征在于,所述蝴蝶运算器执行基数为2的蝴蝶运算。
7.一种数字信号处理方法,其特征在于,包括:
提供储存N点数据的一存储器;
提供一第一子转动系数与一第二子转动系数;
将所述第一子转动系数与所述第二子转动系数通过一蝴蝶运算器的一复数乘法器进行复数乘法运算以产生一权重值;
将权重值应用至N点数据的快速傅立叶转换以产生多个经转置值;以及
将所述经转置值写回至所述存储器。
8.如权利要求7所述的数字信号处理方法,其特征在于,所述蝴蝶运算器执行基数为2的蝴蝶运算。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW101146111A TWI472932B (zh) | 2012-12-07 | 2012-12-07 | 數位訊號處理裝置及其處理方法 |
TW101146111 | 2012-12-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103870437A true CN103870437A (zh) | 2014-06-18 |
Family
ID=50882179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310056976.6A Pending CN103870437A (zh) | 2012-12-07 | 2013-02-22 | 数字信号处理装置及其处理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20140164460A1 (zh) |
CN (1) | CN103870437A (zh) |
TW (1) | TWI472932B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693724A (zh) * | 2020-05-25 | 2020-09-22 | 五邑大学 | 车速侦测系统、车速侦测方法及计算机可读存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9977116B2 (en) * | 2015-10-05 | 2018-05-22 | Analog Devices, Inc. | Scaling fixed-point fast Fourier transforms in radar and sonar applications |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1169083A (zh) * | 1997-03-06 | 1997-12-31 | 华邦电子股份有限公司 | 可巡回执行的离散余弦转换及其逆转换集成电路处理器 |
US6081821A (en) * | 1993-08-05 | 2000-06-27 | The Mitre Corporation | Pipelined, high-precision fast fourier transform processor |
TW580623B (en) * | 2000-02-26 | 2004-03-21 | Walter E Pelton | Apparatus, methods, and computer program products for accurately determining the coefficients of a function |
TW200717261A (en) * | 2005-10-19 | 2007-05-01 | Sunplus Technology Co Ltd | Digital signal processing apparatus |
CN1959666A (zh) * | 2005-10-31 | 2007-05-09 | 凌阳科技股份有限公司 | 数字信号处理装置 |
CN101340414A (zh) * | 2007-07-06 | 2009-01-07 | 联发科技股份有限公司 | 可变长度的快速傅立叶转换系统及其执行方法 |
US8194532B1 (en) * | 2010-01-25 | 2012-06-05 | Xilinx, Inc. | Mixed radix discrete fourier transform |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6735610B1 (en) * | 1999-04-29 | 2004-05-11 | Walter E. Pelton | Apparatus, methods, and computer program products for determining the coefficients of a function with decreased latency |
US20040172435A1 (en) * | 2003-02-27 | 2004-09-02 | Texas Instruments Incorporated | Architecture and method for performing a fast fourier transform and OFDM reciever employing the same |
CN102768654A (zh) * | 2011-05-05 | 2012-11-07 | 中兴通讯股份有限公司 | 具有fft基2蝶运算处理能力的装置及其实现运算的方法 |
-
2012
- 2012-12-07 TW TW101146111A patent/TWI472932B/zh active
-
2013
- 2013-02-22 CN CN201310056976.6A patent/CN103870437A/zh active Pending
- 2013-02-26 US US13/776,724 patent/US20140164460A1/en not_active Abandoned
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6081821A (en) * | 1993-08-05 | 2000-06-27 | The Mitre Corporation | Pipelined, high-precision fast fourier transform processor |
CN1169083A (zh) * | 1997-03-06 | 1997-12-31 | 华邦电子股份有限公司 | 可巡回执行的离散余弦转换及其逆转换集成电路处理器 |
TW580623B (en) * | 2000-02-26 | 2004-03-21 | Walter E Pelton | Apparatus, methods, and computer program products for accurately determining the coefficients of a function |
TW200717261A (en) * | 2005-10-19 | 2007-05-01 | Sunplus Technology Co Ltd | Digital signal processing apparatus |
CN1959666A (zh) * | 2005-10-31 | 2007-05-09 | 凌阳科技股份有限公司 | 数字信号处理装置 |
CN101340414A (zh) * | 2007-07-06 | 2009-01-07 | 联发科技股份有限公司 | 可变长度的快速傅立叶转换系统及其执行方法 |
US8194532B1 (en) * | 2010-01-25 | 2012-06-05 | Xilinx, Inc. | Mixed radix discrete fourier transform |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693724A (zh) * | 2020-05-25 | 2020-09-22 | 五邑大学 | 车速侦测系统、车速侦测方法及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
TWI472932B (zh) | 2015-02-11 |
TW201423438A (zh) | 2014-06-16 |
US20140164460A1 (en) | 2014-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Blahut | Fast algorithms for signal processing | |
US3517173A (en) | Digital processor for performing fast fourier transforms | |
Garrido et al. | Hardware architectures for the fast Fourier transform | |
Park | Guaranteed-stable sliding DFT algorithm with minimal computational requirements | |
Nascimento et al. | A new solution to the hyperbolic tangent implementation in hardware: Polynomial modeling of the fractional exponential part | |
CN101937424A (zh) | 基于fpga实现高速fft处理的方法 | |
KR102376492B1 (ko) | 실수값을 입력으로 하는 고속푸리에 변환장치 및 방법 | |
CN103870437A (zh) | 数字信号处理装置及其处理方法 | |
Selianinau | Computationally efficient approach to implementation of the Chinese Remainder Theorem algorithm in minimally redundant Residue Number System | |
Mankar et al. | FPGA implementation of 16-point radix-4 complex FFT core using NEDA | |
CN115544438A (zh) | 数字通信系统中的旋转因子生成方法、装置和计算机设备 | |
Bansal et al. | Memory-efficient Radix-2 FFT processor using CORDIC algorithm | |
Patronik et al. | Design of Reverse Converters for a New Flexible RNS Five-Moduli Set {2^ k, 2^ n-1, 2^ n+ 1, 2^ n+ 1-1, 2^ n-1-1\} 2 k, 2 n-1, 2 n+ 1, 2 n+ 1-1, 2 n-1-1 (n Even) | |
Venkataramanan et al. | Estimation of frequency offset using warped discrete-Fourier transform | |
Banerjee et al. | A Novel Paradigm of CORDIC-Based FFT Architecture Framed on the Optimality of High-Radix Computation | |
Wang et al. | A novel fast modular multiplier architecture for 8,192-bit RSA cryposystem | |
More et al. | FPGA implementation of FFT processor using vedic algorithm | |
Hazarika et al. | Energy efficient VLSI architecture of real‐valued serial pipelined FFT | |
WO2018170400A1 (en) | Apparatus and methods of providing an efficient radix-r fast fourier transform | |
JP3872724B2 (ja) | 高速フーリエ変換のための回転因子表およびそれを用いた高速フーリエ変換装置 | |
TWI402695B (zh) | 分裂基數-2/8快速傅立葉轉換裝置及方法 | |
Mamatha et al. | Triple-matrix product-based 2D systolic implementation of discrete Fourier transform | |
Du Pont et al. | Hardware Acceleration of the Prime-Factor and Rader NTT for BGV Fully Homomorphic Encryption | |
Chmielowiec | Fast, parallel algorithm for multiplying polynomials with integer coefficients | |
Das et al. | Fast fourier transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140618 |