JP2015503785A - Fft/dftの逆順ソーティングシステム、方法およびその演算システム - Google Patents

Fft/dftの逆順ソーティングシステム、方法およびその演算システム Download PDF

Info

Publication number
JP2015503785A
JP2015503785A JP2014549307A JP2014549307A JP2015503785A JP 2015503785 A JP2015503785 A JP 2015503785A JP 2014549307 A JP2014549307 A JP 2014549307A JP 2014549307 A JP2014549307 A JP 2014549307A JP 2015503785 A JP2015503785 A JP 2015503785A
Authority
JP
Japan
Prior art keywords
order
data
control signal
reversed
reverse
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
Application number
JP2014549307A
Other languages
English (en)
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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Publication of JP2015503785A publication Critical patent/JP2015503785A/ja
Pending legal-status Critical Current

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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • 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

Landscapes

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

Abstract

本発明は、高速フーリエ変換(FFT)/離散フーリエ変換(DFT)の逆順ソーティングシステムを開示する。逆順ソーティングシステムは、各ポイントの順番逆転すべきデータに対応する制御信号を取得し、順番逆転装置に制御信号を送信して、順番逆転装置が受信された制御信号に応じて各ポイントのデータを順番逆転することに用いられる信号制御装置と、受信された制御信号に応じて、順番逆転すべきデータを順番逆転することに用いられる順番逆転装置とを備える。また、本発明は、高速フーリエ変換(FFT)/離散フーリエ変換(DFT)の逆順ソーティングシステム方法と高速フーリエ変換(FFT)/離散フーリエ変換(DFT)の逆順ソーティング演算システムを開示する。本発明に開示されたシステム及び方法により、ベクトルプロセッサでFFT/DFTを計算する場合、FFT/DFTの逆順ソーティング時間が長すぎ、ベクトルプロセッサの利用効率が低い問題を解決し、且つFFT/DFTのベクトル計算効率を向上させることができる。

Description

本発明は、通信技術分野に関し、特に、高速フーリエ変換(Fast Fourier Transformation、FFT)/離散フーリエ変換(DFT)の逆順ソーティングシステム、方法およびその演算システムに関する。
FFT及びDFTは、信号処理、通信システム、特に、直交周波数分割多重技術(OFDM:Orthogonal Frequency Division Multiplexing)に関する長期発展型(Long term Evolution、LTE)、デジタルオーディオ放送(Digital Audio Broadcast、DAB)などの通信システムに広くに応用されている。従来のFFT/DFTアルゴリズムは、時間間引き(DIT)アルゴリズムや周波数間引き(DIF)アルゴリズムに関わらず、入力データまたは出力データに対して専門的な順番逆転操作を一回行わなければならず、且つDITアルゴリズムやDIFアルゴリズムに関わらず、ベクトルプロセッサで実現される場合に、ベクトル処理ユニットに一定の交差接続線を導入しなければならなく、そうでなければ、ベクトルプロセッサの効率を十分に発揮することができず、消費電力の増加をもたらす。
図1に示すように、16ポイントの正順番で入力されたFFTの時間間引きアルゴリズムのフローチャートであり、直接に4つのDIFバタフライベクトルプロセッサで実現される場合、最高の効率は、各段のFFTが2回のベクトルバタフライ演算により完了されるものである。図2から分かるように、最後の2段のFFTが交差接続線を導入しなければ、すべての4つのバタフライを十分に利用することができず、そうでなければ、ベクトルプロセッサの効率を十分に発揮することができない。特に最後の1段に毎回一つだけのバタフライを演算できるので、効率が大幅に低下し、消費電力が増加する。
一方、交差接続線の導入により演算の効率を向上させることができるが、計算ユニットが多い場合、接続線の数量が膨大になり、例えば、ベクトル計算ユニットが16個のDITユニットを備える場合、FFT演算は、最後の4段で交差接続線を導入しなければならず、交差接続線の総数が4000と多く、ここで、データのディジット(digit)幅が16とすると、大量の交差接続線により充放電の電流増加と消費電力の増加をもたらす。
また、上記の交差接続線の問題のほかに、標準的なFFT/DFTアルゴリズムは、バタフライ演算に加えて、一回の専門的なデータ順番逆転操作を必要とする。データの順番逆転は、シリアルまたはパラレルで行われることができ、シリアル方式において、ソフトウェアまたはハードウェアで各データの順番逆転後のアドレスを計算することができ、この方式の順番逆転は、大量の時間がかかり、ベクトルバタフライユニットが多い場合、順番逆転時間がバタフライ演算のための時間の数倍かもしれないので、ベクトルプロセッサの効率が低下する。
本発明の実施例は、FFT/DFT演算における逆順ソーティングを高速で並列に完了することができるFFT/DFTの逆順ソーティングシステム、方法およびその演算システムを提供する。
本発明の実施例に係るFFT/DFTの逆順ソーティングシステムは
各ポイントの順番逆転すべきデータに対応する制御信号を取得し、前記順番逆転装置に前記制御信号を送信することに用いられ、前記制御信号が特定の順番逆転ルールに基づいて設定され、前記特定の順番逆転ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番でm段のFFT/DFTの逆順ソーティングを行われるこを判断し、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにし、データに対応するアドレスディジット(digit)をc...cに表し、その第iディジットのcがその値に応じて上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされるものである信号制御装置と、
受信された前記制御信号に応じて前記順番逆転すべきデータを順番逆転することに用いられる順番逆転装置とを備える。
好ましくは、前記信号制御装置に制御信号検索テーブルが記憶されており、それに応じて、信号制御装置は、前記制御信号検索テーブルで前記順番逆転すべきデータに対応する制御信号を検索して送信することに用いられる。
また、前記システムは、制御信号検索テーブル記憶装置をさらに備える。
それに応じて、前記信号制御装置は、前記制御信号検索テーブル記憶装置から記憶された制御信号検索テーブルを取得し、前記制御信号検索テーブルで前記順番逆転すべきデータに対応する制御信号を検索して送信することに用いられる。
また、前記順番逆転装置は、
前記信号制御装置から送信された制御信号を受信して対応する交換ユニットに送信することに用いられる制御信号処理ユニットと、
各ポイントの順番逆転すべきデータを記憶することに用いられる記憶ユニットと、
前記制御信号に応じて記憶ユニットに記憶された順番逆転すべきデータを交換することに用いられる交換ユニットとをさらに含む。
好ましくは、前記交換ユニットは、多段の網状構造を構成し、段数がベクトルプロセッサによりアクセスされたデータの幅によって決定され、
各段の交換ユニットの数量が同じであり、
各段の前記交換ユニットの数量がベクトルプロセッサにおける演算ユニットの数量によって決定される。
好ましくは、各段の前記交換ユニットの数量がベクトルプロセッサにおける演算ユニットの数量によって決定されることは、
各段の交換ユニットの数量がベクトルプロセッサにおける演算ユニットの数量の半分であることである。
好ましくは、前記交換ユニットは、具体的に2つのマルチプレクサを含み、
各マルチプレクサの制御信号入力側は、2つのポイントの順番逆転すべきデータに対応する制御信号を入力することに用いられ、
各マルチプレクサの2つのデータ入力側は、それぞれ2つのポイントの順番逆転すべきデータに対応する記憶ユニットに接続され、入力された前記制御信号に応じて、入力された2つのポイントのデータを交換し、または直接通過させる。
本発明に係るFFT/DFT演算システムは、ベクトル演算装置と、ベクトルレジスタと、メモリとを備え、前記ベクトル演算装置が各段のFFT/DFT演算後のデータを前記ベクトルレジスタに送信してキャッシュし、前記ベクトルレジスタによりキャッシュされたデータが順番逆転されて前記メモリに送信される。前記システムは、上述した順番逆転システムをさらに備える。
前記順番逆転システムは、前記ベクトル演算装置により各段のFFT/DFT演算を行われた後、前記ベクトルレジスタにキャッシュされたデータを順番逆転して、前記メモリに送信することに用いられる。
前記メモリは、前記順番逆転システムにより順番逆転された各段のデータを記憶し、次の段の演算すべきデータを前記ベクトルレジスタに送信してキャッシュし、次の段のFFT/DFT演算を待っていることに用いられる。
本発明に係るFFT/DFTの逆順ソーティング方法は、
各ポイントの順番逆転すべきデータに対応する制御信号を取得し、前記制御信号に応じて、前記逆順すべきデータを順番逆転するようにソーティングするステップを含む。
ここで、前記制御信号は、特定の順番逆転ルールに基づいて設定され、前記特定の順番逆転ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番でm段のFFT/DFTの逆順ソーティングを行われることを確定し、順番逆転ルールが、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにして、データに対応するアドレスディジット(digit)をc...cに表し、その第iディジットのcがその値に応じて上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされるものである。
また、初めて順番逆転すべきデータを取得する場合、演算段を0に初期化し、
次の段の逆順データを受信するたびに、演算段に1を加える。
また、各段のデータを逆順ソーティングする場合、記憶された制御信号検索テーブルを検索することで、前記順番逆転すべきデータに対応する制御信号を検索し、前記制御信号に応じてこの段のデータを逆順ソーティングする。m回のバタフライ演算とm-1回のソーティングを経た後、FFT/DFT演算を完了させる。
本発明は、従来技術に比べて以下のようなメリットを持っている。
本発明の実施例に係るFFT/DFTの逆順ソーティングシステム、方法およびその演算システムによれば、FFT/DFTの逆順ソーティングは、ベクトルパイプラインデータチャネル中のハードウェアソーティングネットワークにおける交換ユニットにより行われ、各段のバタフライ演算後にパイプラインにおける順番逆転装置がデータをソーティングするように制御することにより、ベクトルプロセッサでFFT/DFTを計算する場合、FFT/DFTの逆順配列時間が長すぎ、ベクトルプロセッサの利用効率が低い問題を解決して、FFT/DFTのベクトル演算の効率を向上させることができる。
従来技術における16ポイントの正順番で入力されたFFTの時間間引きのフローチャートである。 図1の4つのDIFバタフライベクトルプロセッサの必要な交差線のレイアウト図である。 本発明の実施例に係る逆順ソーティングシステムの構造を示す図である。 本発明の実施例に係る演算システムの構造を示す図である。 本発明の実施例に係る方法のフローチャートである。 本発明の実施例1に係る順番逆転装置の構造を示す図である。 本発明の実施例1における交換ユニットの構造を示す図である。 本発明の実施例2に係る16ポイントのFFTの反復順番逆転ルールを示す図である。 本発明の実施例3に係る24ポイントのDFTの反復順番逆転ルールを示す図である。 本発明の実施例4に係る16ポイントのデータを段間にソーティングするフローチャートである。 本発明の実施例5に係る24ポイントのデータを段間にソーティングするフローチャートである。 本発明の実施例6における16個のDIFバタフライユニットによりソーティングするフローチャートである。 本発明の実施例6における第1段のデータチャネルの構造を示す図である。 本発明の実施例6における第2段のデータチャネルの構造を示す図である。 本発明の実施例6における第3段のデータチャネルの構造を示す図である。 本発明の実施例6における第4段のデータチャネルの構造を示す図である。
以下、各図面を参照して本発明の実施例に係る技術的スキームの主な実現原理、具体的な実施形態およびそれに対応する有益な効果を詳細に説明する。
図3は本発明の実施例に係る逆順ソーティングシステムの構造を示す図である。当該逆順ソーティングシステムは、信号制御装置31と、順番逆転装置32とを備える。
信号制御装置31は、各ポイントの順番逆転すべきデータに対応する制御信号を取得し、順番逆転装置32に前記制御信号を送信して、前記順番逆転装置32が受信された前記制御信号に応じて各ポイントのデータを、特定の順番逆転ルールに基づいて順番逆転することに用いられ、前記制御信号が前記特定の順番逆転ルールに基づいて設定され、前記特定の逆順ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番でm段のFFT/DFTの逆順ソーティングを行われることを判断し、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにして、データに対応するアドレスディジット(digit)をc...cに表し、その第iディジットのcがその値に応じて上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされる。
逆順装置32は、受信された前記制御信号に応じて前記順番逆転すべきデータを順番逆転することに用いられる。
前記信号制御装置に制御信号検索テーブルが記憶されている。
それに応じて、信号制御装置は、前記制御信号検索テーブルで前記順番が逆転すべきデータに対応する制御信号を検査して送信することに用いられる。
前記逆順システムは、制御信号検索テーブル記憶装置をさらに備える。
それに応じて、前記信号制御装置は、前記制御信号検索テーブル記憶装置から記憶された制御信号検索テーブルを取得し、前記制御信号検索テーブルで前記順番逆転すべきデータに対応する制御信号を検出して送信することに用いられる。
ここで、前記順番逆転装置は、制御信号処理ユニットと、記憶ユニットと、交換ユニットとをさらに備える。
前記制御信号処理ユニットは、前記信号制御装置から送信された制御信号を受信して交換ユニットに送信することに用いられる。
前記記憶ユニットは、各ポイントの順番逆転すべきデータを記憶することに用いられる。
前記交換ユニットは、前記制御信号に応じて記憶ユニットに記憶された順番逆転すべきデータを交換することに用いられる。
ここで、前記交換ユニットは、m段の網状構造を構成し、各段の交換ユニットの数量が同じである。
各段の交換ユニットの数量は、ベクトルプロセッサにおける演算ユニット数によって決定される。
すなわち、各段の交換ユニットの数量がベクトルプロセッサにおける演算ユニットの数量の半分である。
ここで、前記交換ユニットは、2つのマルチプレクサを含む。
各マルチプレクサの制御信号入力側は、2つのポイントの順番逆転すべきデータに対応する制御信号を入力することに用いられる。
各マルチプレクサの2つのデータ入力側は、それぞれ2つのポイントの順番逆転すべきデータに対応する記憶ユニットに接続され、入力された前記制御信号に応じて、入力された2つのポイントのデータを交換し、または直接通過させる。
図4は本発明の実施例に係る演算システムの構造を示す図である。前記演算システムは、ベクトル演算装置41と、ベクトルレジスタ42と、メモリ43とを備え、ベクトル演算装置41が各段のFFT/DFT演算後のデータをベクトルレジスタ42に送信してキャッシュし、ベクトルレジスタ42によりキャッシュされたデータが順番逆転されて前記メモリ43に送信され、前記演算システムは、順番逆転システム44をさらに備える。
順番逆転システム44は、前記ベクトル演算装置41により各段のFFT/DFTの演算を行われて前記ベクトルレジスタ42にキャッシュされたデータを順番逆転して、前記メモリ43に送信することに用いられる。
メモリ43は、各段の前記順番逆転システム44により順番逆転されたデータを記憶し、次の段の演算すべきデータを前記ベクトルレジスタ42に送信してキャッシュし、次の段FFT/DFTの演算を待っていることに用いられる。
図5は本発明の実施例に係る方法のフローチャートである。前記方法は、以下のステップを含む。
ステップ1:各ポイントの順番逆転すべきデータに対応する制御信号を取得する。
ステップ2:前記制御信号に応じて、前記順番逆転すべきデータを順番逆転する。
ここで、前記制御信号は、特定の順番逆転ルールに基づいて設定され、前記特定の逆順ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番でm段のFFT/DFTの逆順ソーティングを行行われることを確定し、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにして、データに対応するアドレスディジット(digit)をc...cに表し、その第iディジットのcがその値に応じて上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされるものである。
図6は本発明の実施例1による逆順装置の構造を示す図である。図6に示すように、各黒点は、1つの交換ユニットを表し、5段の交換ユニットを含み、各段の交換ユニットの数量が16個であり、32個のサブスクリプト装置、すなわち、記憶ユニットで順番逆転すべきデータを記憶することにより、各ポイントデータは、直接的又は交差して出力されることができる。具体的なチャネルは、交換ユニットにより受信された制御信号によって制御される。図7は交換ユニットの構造を示す図である。前記交換ユニットは、図7の矩形ブロック内に示される二つのマルチプレクサを含み、各マルチプレクサの制御信号入力側Ct1が二つのポイントの順番逆転すべきデータに対応する制御信号を入力することに用いられ、各マルチプレクサの二つのデータ入力側In0、In1がそれぞれ二つのポイントの順番逆転すべきデータに対応する記憶ユニットに接続され、入力された前記制御信号によって、入力された二つのポイントのデータを交換し、その出力側Out0、Out1を介して交換された二つのポイントのデータを出力する。
図8は本発明の実施例2に係る16ポイントのFFTの反復順番逆転ルールを示す図である。そのプロセスは、以下のとおりである。16個のデータが4個のアドレスディジット(digit)Cで表され、左から右に、データを一回ソーティングするたびに、4つのアドレスディジット(digit)が毎回1つのディジットを増やして最終の順番逆転結果に一致するようにする。第1回のソーティングにより、Cが順番逆転後のアドレス配列に一致し、そのプロセスがデータをC=0の上半分とC=1の下半分に分割し、その後上半分と下半分を交互にソーティングすることにより実現される。第2回のソーティングにより、Cが順番逆転後のアドレス配列に一致し、既に配列されたCが影響を受けないようにするために、隣接しているCをそれぞれ0と1であるデータに分割して二つずつ1つのデータグループと見なし、データグループを単位にしてデータをソーティングし、そのプロセスがCのソーティングプロセスと類似し、同様に、データをC=0の上半分とC=1の下半分に分割し、その後上半分と下半分のデータを、データグループを単位にして交互にソーティングすることにより実現される。第3回のソーティングにより、Cが順番逆転後のアドレス配列に一致し、既に配列されたCが影響を受けないことを保証するために、隣接しているCがそれぞれ00、01、10、11であるデータを1つのデータグループと見なし、データグループを単位にしてデータをソーティングし、同様に、データをC=0の上半分とC=1の下半分に分割し、その後上半分と下半分のデータを、データグループを単位にして交互にソーティングすることにより実現される。この場合、Cの配列順番と順番逆転後のCの配列順番は、自動的に一致し、順番逆転を完了する。
図9は本発明の実施例3に係る24ポイントのDFTの反復順番逆転ルールを示す図である。図9に示すように、ラディクス2、2、3、2の順番で計算された24ポイントのDFTの反復順番逆転プロセスは、以下のとおりである。そのアドレス混合ラディクスがC(C、CおよびCが2進数(binary)であり、Cが3進数(ternary)である)に表されると仮定した場合、第1回のソーティングにより、Cの順番が順番逆転後のアドレス配列に一致し、そのソーティングプロセスが上記の16ポイントのFFTと類似する。第2回のソーティングにより、Cの順番が順番逆転後のアドレス配列に一致し、そのソーティングプロセスは、上記の16ポイントのFFTと類似する。第3回のソーティングにより、Cの順番が順番逆転後のアドレス配列に一致し、そのソーティングプロセスが、ソーティングされたCへの影響を避けるように、隣接しているCがそれぞれ00、01、10、11であるデータを1つのデータグループと見なし、データグループを単位にしてデータをソーティングし、その後、データをCの値に応じてC=0、1、2の3つの部分に分割して、3つの部分のデータを、データグループを単位にして交互にソーティングすることである。この場合、Cの順番と順番逆転後のCの配列順番は、自動的に一致して、順番逆転を完了する。
図10は本発明の実施例4に係る16ポイントのデータを段間にソーティングするフローチャートである。図11は本発明の実施例5に係る24ポイント(ラディクス2、2、3、2)のデータを段間にソーティングするフローチャートである。図10と図11に示すように、各ポイントのデータが正順番で入力され、演算結果は、FFT/DFTの順番逆転すべき結果である。同時に、図10と図11から分かるように、各段の入力が要求に応じてソーティングされると、バタフライ演算のラディクスが同じであれば、異なる演算段の入力データが完全に同じ順番でベクトルユニットに送信されて演算されることができる。16ポイントのFFTを例にすると、各段で、最初の8個のデータと最後の8個のデータが対応してバタフライを構成して演算し、かつバタフライの入力データ間に等しい最大距離があり、バタフライのベクトル計算に非常に有利であり、且つベクトル計算ユニットにFFT/DFTのために余分なデータ接続線を導入する必要がないため、消費電力の低減にも有利である。順番逆転のために占用される時間を完全になくすために、上記のソーティングプロセスを、図4に示すようにデータをレジスタからメモリに記憶するパイプライン段に挿入することができる。ベクトル演算ユニットに16個のDITバタフライユニットがあり、かつFFT/DFTポイント数に少なくとも4つの因子2が含まれる場合、ソーティングネットワークにおいて図12に示すような4段のデータソーティングモードを実現する必要がある。ここで、ラディクス2演算を一番前に配列しなければならなく、対応するソーティングが図6に示すような5段のソーティングネットワークで実現されることができる。ここで、図6における各黒点は、一つの交換ユニットを表し、データは、直接通過または交差してもよい。交換ユニットは、図7に示すような回路の構造を採用してもよく、すなわち、2つのマルチプレクサから構成されることができる。当該ソーティングネットワークは、図12で要求された4段のソーティングを実現することができ、当該ソーティングネットワークにより最初の4段のソーティングを実現する場合のデータチャネルは、図13〜図16に示される。
また、より大きなポイント数のFFTのバタフライ演算、例えばLTEの必要な2048ポイントのFFTに対して、最初の4段の演算において毎回1024ポイントを隔てる連続した16個のデータをベクトル計算ユニットに入力してバタフライ演算を行い、次に結果をメモリに記憶する場合に淳番逆転システムは、それぞれ4つのソーティング方式でデータをソーティングし、4段後のFFT演算において、交互にソーティングされたデータグループの長さが16の倍数であるため、データを記憶する時にメモリへのデータの記憶アドレスを制御するだけで必要なソーティングを実現することができる。FFT/DFTにおける因子2の個数が4個未満の場合、例えば、LTEアップリンクDFTにおける因子2の個数が少なくとも2つであり、かつ因子2、3、5しかない場合、ラディクス2演算を先ず配列して計算する必要があり、ソーティングネットワークのソーティング機能を利用して演算結果をソーティングし、その後、1または2段の演算は、データがメモリに記憶される場合に長さ4、12または8のマスクとシフト制御を追加しなければならず、4、12または8つのバタフライ演算結果を記憶するたびに演算効率が少し低下するが、依然としてデータの記憶アドレスで制御してデータをソーティングするために占有された排他時間を避けることができる。
本発明に係るFFT/DFTのベクトル計算アーキテクチャとそれに対応するハードウェアソーティングネットワークでは、FFT/DFTの逆順を実行するため、ベクトル演算結果がメモリのパイプラインに記憶され、データの順番逆転のために占有された排他時間を削除し、且つベクトル計算ユニットに大量の交差接続線が導入されることを避けることができ、ハードウェアのオーバーヘッドを削減し、FFT/DFT演算の消費電力を削減することができる。このソーティングネットワークは、ベクトル演算に必要なその他のソーティング操作をサポートするため、一連のソーティング操作を実行する。
本分野の当業者は、本発明の実施例が方法、システム、またはコンピュータプログラム製品として提供されることができると理解すべきである。したがって、本発明は、完全なハードウェアの実施例、完全なソフトウェアの実施例、またはソフトウェアとハードウェアを組み合わせた実施例の携帯を適用することができる。また、本発明は、コンピュータ使用可能プログラムコードが含まれている1つ以上のコンピュータ使用可能な記憶媒体で実施されることができるコンピュータプログラム製品の形態を適用することができ、前記コンピュータ使用可能な記憶媒体は、ディスク記憶装置、CD-ROM、光記憶装置を含むがこれらに限定されない。
本発明は、本発明の実施例に係る方法、装置/システム、およびコンピュータプログラム製品のフローチャートおよび/またはブロック図を参照して説明する。コンピュータプログラムコマンドによって、フローチャートおよび/またはブロック図中の各プロセスおよび/またはブロック、およびフローチャートおよび/またはブロック図中のプロセスおよび/またはブロックの組み合わせを実現することができると理解すべきである。これらのコンピュータプログラムコマンドを、汎用コンピュータ、専用コンピュータ、組み込みプロセッサ、または他のプログラミング可能なデータ処理装置のプロセッサに提供して1つのマシンを生成することができ、それによってコンピュータまたは他のプログラム可能なデータ処理装置のプロセッサによって実行されるコマンドによってフローチャートにおける1つ以上のプロセスおよび/またはブロック図における1つ以上のブロック中に指定された機能を実現するためのデバイスを生成する。
これらのコンピュータプログラムコマンドは、コンピュータまたは他のプログラミング可能なデータ処理装置が特定の方式で動作するように案内することができるコンピュータ読み取り可能な記憶装置に記憶されることができ、それによってこのコンピュータ読み取り可能な記憶装置に記憶されたコマンドは、コマンドデバイスが含まれている製品を生成し、このコマンドデバイスがフローチャートにおける1つ以上のプロセスおよび/またはブロック図における1つ以上のブロック中に指定された機能を実現する。
これらのコンピュータプログラムコマンドは、コンピュータまたは他のプログラミング可能なデータ処理装置にロードされることができ、それによってコンピュータまたは他のプログラミング可能な装置で一連の操作手順を実行して、コンピュータで実現される処理を生成することにより、コンピュータまたは他のプログラミング可能な装置で実行されるコマンドによって、フローチャートにおける1つ以上のプロセスおよび/またはブロック図における1つ以上のブロック中に指定された機能を実現するための手順を提供される。
本発明の好ましい実施例を説明したが、本分野の当業者にとって一旦基本的な創造性の概念を知った場合、これらの実施例を別に変更および修正することができる。従って、添付された請求項は、好ましい実施例及び本発明の範囲内に含まれるすべての変更と修正が含まれていると解釈されることを意図している。。
明らかに、本分野の当業者にとって、本発明の精神と範囲を逸脱することなく、本発明に対して様々な変更及び変形を行うことができる。このようにして、本発明のこれらの変更及び変形が本発明の請求項およびその同等技術の範囲内に属すると、本発明は、これらの変更及び変形を含むことを意図している。

Claims (10)

  1. 各ポイントの順番逆転すべきデータに対応する制御信号を取得し、前記順番逆転装置に前記制御信号を送信することに用いられ、
    前記制御信号が特定の順番逆転ルールに基づいて設定され、前記特定の順番逆転ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番に応じてm段のFFT/DFTの逆順ソーティングを行われることを確定し、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにして、データに対応するアドレスディジットをc...cに表し、c...cの第iディジットのcがその値に応じて上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされるものである信号制御装置と、
    受信された前記制御信号に応じて、前記順番逆転すべきデータを順番逆転することに用いられる順番逆転装置とを備えることを特徴とする高速フーリエ変換(FFT)/離散フーリエ変換(DFT)の逆順ソーティングシステム。
  2. 前記信号制御装置に制御信号検索テーブルが記憶されており、
    それに応じて、前記信号制御装置は、前記制御信号検索テーブルで前記順番逆転すべきデータに対応する制御信号を検索して送信することに用いられることを特徴とする
    請求項1に記載のシステム。
  3. 制御信号検索テーブル記憶装置をさらに備え、
    それに応じて、前記信号制御装置は、前記制御信号検索テーブル記憶装置から記憶された制御信号検索テーブルを取得し、前記制御信号検索テーブルで前記順番逆転すべきデータに対応する制御信号を検索して送信することに用いられることを特徴とする
    請求項1記載のシステム。
  4. 前記信号制御装置から送信された制御信号を受信して交換ユニットに送信することに用いられる制御信号処理ユニットと、
    各ポイントの順番逆転すべきデータを記憶することに用いられる記憶ユニットと、
    前記交換ユニットは、前記制御信号に応じて記憶ユニットに記憶された順番逆転すべきデータを交換することに用いられる交換ユニットとをさらに備えることを特徴とする
    請求項1に記載のシステム。
  5. 前記交換ユニットは、m段の網状構造を構成し、各段の交換ユニットの数量が同じであり、
    各段の前記交換ユニットの数量がベクトルプロセッサにおける演算ユニット数によって決定されることを特徴とする
    請求項4に記載のシステム。
  6. 各段の交換ユニットの数量がベクトルプロセッサにおける演算ユニット数の半分であることであることを特徴とする
    請求項5に記載のシステム。
  7. 前記交換ユニットは、2つのマルチプレクサを含み、
    各マルチプレクサの制御信号入力側は、2つのポイントの順番逆転すべきデータに対応する制御信号を入力することに用いられ、
    各マルチプレクサの2つのデータ入力側は、それぞれ2つのポイントの順番逆転すべきデータに対応する記憶ユニットに接続され、入力された前記制御信号に応じて、入力された2つのポイントのデータを交換し、または直接通過させることを特徴とする
    請求項6に記載のシステム。
  8. ベクトル演算装置と、ベクトルレジスタと、メモリとを備え、前記ベクトル演算装置が各段のFFT/DFTの演算後のデータを前記ベクトルレジスタに送信してキャッシュし、前記ベクトルレジスタによりキャッシュされたデータが順番逆転されて前記メモリに送信されるFFT/DFTの演算システムであって、
    請求項1〜7に記載の逆順ソーティングシステムをさらに備え、
    前記逆順システムは、前記ベクトル演算装置により各段FFT/DFTの演算を行われて前記ベクトルレジスタにキャッシュされたデータを順番逆転して、前記メモリに送信することに用いられ、
    前記メモリは、前記逆順システムにより順番逆転された各段のデータを記憶し、次の段の演算すべきデータを前記ベクトルレジスタに送信してキャッシュし、次の段のFFT/DFTの演算を待っていることに用いられることを特徴とするシステム。
  9. FFT/DFTの逆順ソーティング方法であって、
    各ポイントの順番逆転すべきデータに対応する制御信号を取得し、前記制御信号に応じて、前記順番逆転すべきデータを順番逆転するステップとを含み、
    前記制御信号は、特定の順番逆転ルールに基づいて設定され、前記特定の逆順ルールが、順番逆転すべきデータのポイント数N=p...p...pに応じて、当該ポイント数がラディクスpm、...ラディクスpi、... ラディクスp2、ラディクスpの順番に応じてm段のFFT/DFTの逆順ソーティングを行われることを確定し、第m-i段のデータを順番逆転する場合、隣接しているp*pm-1*...*pi+1を一つのデータグループにして、データに対応するアドレスディジットをc...cに表し、c...cの第iディジットのcが、上から下にpi個の部分に分けられ、前記データグループを単位にして交互にソーティングされることを特徴とする方法。
  10. 各段のデータを逆順ソーティングする場合、記憶された制御信号検索テーブルを検索することにより、前記順番逆転すべきデータに対応する制御信号を検索し、前記制御信号に応じてこの段のデータを逆順ソーティングするステップをさらに含むことを特徴とする
    請求項9に記載の逆順ソーティング方法。
JP2014549307A 2011-12-27 2012-06-28 Fft/dftの逆順ソーティングシステム、方法およびその演算システム Pending JP2015503785A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201110444738.3 2011-12-27
CN201110444738.3A CN103186503B (zh) 2011-12-27 2011-12-27 一种fft/dft倒序排列系统与方法及其运算系统
PCT/CN2012/077800 WO2013097436A1 (zh) 2011-12-27 2012-06-28 一种fft/dft的倒序排列系统与方法及其运算系统

Publications (1)

Publication Number Publication Date
JP2015503785A true JP2015503785A (ja) 2015-02-02

Family

ID=48677674

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014549307A Pending JP2015503785A (ja) 2011-12-27 2012-06-28 Fft/dftの逆順ソーティングシステム、方法およびその演算システム

Country Status (5)

Country Link
EP (1) EP2800010A4 (ja)
JP (1) JP2015503785A (ja)
KR (1) KR101696987B1 (ja)
CN (1) CN103186503B (ja)
WO (1) WO2013097436A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10013253B2 (en) * 2014-12-23 2018-07-03 Intel Corporation Method and apparatus for performing a vector bit reversal
CN107948102A (zh) * 2016-10-12 2018-04-20 深圳市中兴微电子技术有限公司 一种矢量数据的并行排序网络及方法
CN109522674A (zh) * 2018-12-05 2019-03-26 天津大学 基2-2算法的快速傅里叶变换硬件设计方法
CN112667959B (zh) * 2020-12-30 2024-02-23 上海金卓科技有限公司 混合基dft处理方法、装置、设备及存储介质
CN114578093B (zh) * 2022-03-10 2023-08-18 中国计量科学研究院 一种基于混合基fft的激光多普勒测速仪测速方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06342449A (ja) * 1993-06-01 1994-12-13 Fujitsu Ltd データ並べ替え回路,及び、高速フーリエ変換回路
JP2011227559A (ja) * 2010-04-15 2011-11-10 Nec Corp データ並び替え回路に対する制御情報自動生成方法、データ並び替え回路、及び制御情報自動生成プログラム

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3731284A (en) * 1971-12-27 1973-05-01 Bell Telephone Labor Inc Method and apparatus for reordering data
US6609140B1 (en) * 1999-11-30 2003-08-19 Mercury Computer Systems, Inc. Methods and apparatus for fast fourier transforms
US7062523B1 (en) * 2000-08-01 2006-06-13 Analog Devices, Inc. Method for efficiently computing a fast fourier transform
KR100692997B1 (ko) * 2001-04-17 2007-03-12 삼성전자주식회사 패스트퓨리에변환 장치
CA2352683C (en) * 2001-07-09 2005-11-15 Tropic Networks Inc. Real-time method for bit-reversal of large size arrays
US20040059766A1 (en) * 2002-09-23 2004-03-25 Yeou-Min Yeh Pipelined low complexity FFT/IFFT processor
CN101083643A (zh) * 2006-05-31 2007-12-05 中国科学院微电子研究所 一种低存储器开销的混合基fft处理器及其方法
CN101404550B (zh) * 2008-07-08 2010-07-28 宁波中科集成电路设计中心有限公司 一种drm接收机中基于正交频分复用的解调器
CN102346728B (zh) * 2010-07-29 2016-02-10 中兴通讯股份有限公司 一种采用矢量处理器实现fft/dft倒序的方法和装置
CN102129419B (zh) * 2011-03-04 2016-02-03 中山大学 基于快速傅立叶变换的处理器

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06342449A (ja) * 1993-06-01 1994-12-13 Fujitsu Ltd データ並べ替え回路,及び、高速フーリエ変換回路
JP2011227559A (ja) * 2010-04-15 2011-11-10 Nec Corp データ並び替え回路に対する制御情報自動生成方法、データ並び替え回路、及び制御情報自動生成プログラム

Also Published As

Publication number Publication date
WO2013097436A1 (zh) 2013-07-04
KR101696987B1 (ko) 2017-01-16
EP2800010A1 (en) 2014-11-05
CN103186503A (zh) 2013-07-03
CN103186503B (zh) 2017-05-17
EP2800010A4 (en) 2015-08-05
KR20140116160A (ko) 2014-10-01

Similar Documents

Publication Publication Date Title
Kim et al. Hardware architecture of a number theoretic transform for a bootstrappable RNS-based homomorphic encryption scheme
CN105975251B (zh) 一种基于粗粒度可重构架构的des算法轮迭代系统及迭代方法
CN107957976B (zh) 一种计算方法及相关产品
CN102652315B (zh) 信息处理设备及其控制方法
JP2015503785A (ja) Fft/dftの逆順ソーティングシステム、方法およびその演算システム
CN111737638A (zh) 基于傅里叶变换的数据处理方法及相关装置
CN107957977B (zh) 一种计算方法及相关产品
CN104699624A (zh) 面向fft并行计算的无冲突存储访问方法
US20150301986A1 (en) Fast fourier transform circuit, fast fourier transform processing method, and program recording medium
JP4755610B2 (ja) 高速フーリエ変換装置
JP5601327B2 (ja) データ並べ替え回路、可変遅延回路、高速フーリエ変換回路、およびデータ並べ替え方法
CN104268124A (zh) 一种fft实现装置和方法
KR100836624B1 (ko) 가변 고속 푸리에 변환 장치 및 그 방법
CN108108189B (zh) 一种计算方法及相关产品
CN103493039B (zh) 数据处理方法、数据处理装置、接入设备和用户设备
US9268744B2 (en) Parallel bit reversal devices and methods
CN110766136B (zh) 一种稀疏矩阵与向量的压缩方法
CN109672524B (zh) 基于粗粒度可重构架构的sm3算法轮迭代系统及迭代方法
KR20140142927A (ko) 혼합 기수 파이프라인 fft 프로세서 및 이를 이용한 fft 프로세싱 방법
EP1076296A2 (en) Data storage for fast fourier transforms
CN108037908B (zh) 一种计算方法及相关产品
CN111368250B (zh) 基于傅里叶变换/逆变换的数据处理系统、方法及设备
CN102023963B (zh) 高速多模式时频域变换方法
CN109558638B (zh) Fft处理器
CN109753629B (zh) 多粒度并行fft计算装置

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150917

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160121

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20160520

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160920

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20160929

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20161014

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170809