JPH0229821A - Cordic複素数乗算器 - Google Patents

Cordic複素数乗算器

Info

Publication number
JPH0229821A
JPH0229821A JP1136295A JP13629589A JPH0229821A JP H0229821 A JPH0229821 A JP H0229821A JP 1136295 A JP1136295 A JP 1136295A JP 13629589 A JP13629589 A JP 13629589A JP H0229821 A JPH0229821 A JP H0229821A
Authority
JP
Japan
Prior art keywords
digital data
data word
output
sign
word
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
JP1136295A
Other languages
English (en)
Other versions
JP3283504B2 (ja
Inventor
Matthew O'donnell
マシュー・オドネル
William E Engeler
ウィリアム・アーネスト・エンゲラー
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of JPH0229821A publication Critical patent/JPH0229821A/ja
Application granted granted Critical
Publication of JP3283504B2 publication Critical patent/JP3283504B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/49Computations with a radix, other than binary, 8, 16 or decimal, e.g. ternary, negative or imaginary radices, mixed radix non-linear PCM
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03CMODULATION
    • H03C3/00Angle modulation
    • H03C3/38Angle modulation by converting amplitude modulation to angle modulation
    • H03C3/40Angle modulation by converting amplitude modulation to angle modulation using two signal paths the outputs of which have a predetermined phase difference and at least one output being amplitude-modulated
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/4806Computations with complex numbers
    • G06F7/4812Complex multiplication
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/4806Computations with complex numbers
    • G06F7/4818Computations with complex numbers using coordinate rotation digital computer [CORDIC]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/544Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices for evaluating functions by calculation
    • G06F7/5446Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices for evaluating functions by calculation using crossaddition algorithms, e.g. CORDIC

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Nonlinear Science (AREA)
  • Complex Calculations (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 発明の背景 この発明は算術処理手段、更に具体的に云えば、1対の
複素数の乗算を行なう為に、少なくとも1つの座標回転
ディジタル計算機(CORDI C)を利用した新規な
手段に関する。
いろいろな形式の今日の電子装置では、相関の検出、離
散的なフーリエ変換等の様な比較的複雑な信号処理機能
が利用されている。基本の信号処理素子は乗算−累算セ
ルである。ベースバンド信号処理等を利用する超音波作
像の様なある装置では、完全に複素数の信号を乗算しな
ければならない。電子信号の同様な複素数乗算は、レー
ダ、ソナー等の様な他の多くの信号処理分野で、特にデ
ィジタル信号、更に特に2進形の信号の処理に見られる
。従って、何れもx+jy又はR2Oと云う形をした複
素数を表わす1つのディジタル信号の複素数乗算を実施
するディジタル信号手段を提供することが非常に望まし
い。
従来、IREトランザクションズ・オン・エレクトロニ
ック・コンピュータズ誌EC−8、第330頁乃至第3
34頁(1959年)所載のJ。
E、ボルダ−の論文rCORDI C三角関数計算方式
」に最初に記載された様なC0RD I C装置及び方
式を使うことが述べられている。これは、角度θにわた
る回転は、その何れの回転も特別の1組の角度αの内の
1つにわたる様な幾つかの回転の和として表わすことが
出来ると云う計算方式%式% こ\でξ1−+1又は−1である。C1−90’と定義
すると −1−〇 α  −tan  (2)   n−0,1,2・・・
n+2 即ち、全体の角度が、複数個(n個)の角度α1の全部
を使って順次近似され、角度の一層細かい夫々の近似が
直角座標の結果X  及びyn+1をn+1 生じ、これらは(次に粗い近似に対する)直角座標の値
X。及びy。に対して次の1対の方程式で表わされる様
な関係を持つ。
Xn+1−K(θ)(xn十ξ+Yn/2)(3a) yo+1−K(θ)(”/n−ξ+xn/2)(3b) ニーでK(θ)はCOS (θ)に等しい倍率である。
各々の係数2−0は事実上除数2の除算をn回行なうこ
とであり、2進数では、このn回の各々に対して1ビッ
トのシフトによって行なわれる。
この為、複素数乗算は、(必要とする場合の)倍率K(
θ)の乗算を別として、1組のシフトレジスタ及び加算
器を用いて実行することが出来る。
この基本的なC0RD I C方式を利用して、1対の
複素数を乗算する新規な装置を提供することが非常に望
ましい。
発明の要約 この発明の現在好ましいと考えられる2つのディジタル
複索数B及びCを乗算するC0RD I C装置は、直
角座標の形(例えば、CR及びC+)で表わした一方の
数の実数データ部分及び虚数データ部分の各々を、極座
標の形で表わした他方の数(例えば、IB+、φ)の位
相角φだけ回転させるために、乗算器のない、再帰形又
はパイプライン逐次形の何れかのN段C0RD I C
回転手段を用いる。回転の後、CORDIC回転手段の
出力の実数又は虚数データ部分の各々に他方の数の大き
さのデータ|B|のスカラー乗算を行なう手段を設ける
。こうして計算された最終的なデータが、積の実数及び
虚数部分である。
現在好ましいと考えられる別の複素数乗算CORDIC
装置では、1対の再帰形又はバイブライン逐次形位相回
転手段の各々が、第1及び第2の直角座標形式の複素数
の実数及び虚数部分に作用する。各々の複素数の位相は
、符号検出器で判定して、第1の数がゼロの位相角にな
るまで平等に回転し、これにより出力データの位相が第
1及び第2の複素数の位相角の和に等しくなる様にする
必要によっては、全体の複素数の積を決定する為に、積
データのスカラー乗算を使うことが出来る。
相関の統計等の計算の様な多くの場合、正確な位相情報
を累算することだけが必要であって、位相情報の精度が
保たれている限り、積の振幅項には比較的精度が要求さ
れない。使う場合、各々のスカラー乗算器はシフト及び
累算部分にすることが出来る。
従って、この発明の目的は、1対のディジタル複素数を
乗算する新規なC0RD I C回路を提供することで
ある。
この発明の上記並びにその他の目的は、以下図面につい
て詳しく説明する所を読めば、明らかになろう。
発明の詳細な説明 最初に第1図について説明すると、現在好ましいと考え
られる1実施例の複素数ディジタルCORDIC乗算器
10では、C0RDI Cディジタル乗算器手段11と
1対のスカラー・ディジタル乗算手段12とを利用して
いる。第1の複索数Cが直角座標の形で導入され、Pビ
ットの実数軸成分CRデータ・ワードが入力10aから
入って、CORDIC乗算器の第1の入力又はI入力1
1aに接続され、Pビットの虚数軸成分crが第2の入
力10bに加えられ、乗算器の第2の入力又はQ入力1
1bに接続される。使う場合、1対のスカラー乗算手段
12の各々は、第1のデータ入力12−1a又は12−
2aが、CORDIC手段の虚数軸データQ′出力11
f又は実数軸データ■′出力11gから夫々データを受
取る。各々の手段12は、ディジタル・データ計算の分
野で周知の様に、シフトレジスタ及び累算器からなる1
組であってよい。ある計算タスクでは、振幅の精度が低
下しても差支えないことがあるから、手段12を設ける
場合、その部品は比較的精度の低いものであってよい。
他方の複素数Bが極座標の形で導入され、大きさIB+
の成分に対するPビットのデータ・ワードが入力10c
に加えられ、両方のスカラー・データ乗算手段の第2の
入力12−1b及び12”−2bに接続され、角度φ成
分のデータ・ワードが入力10dに加えられる゛。別の
入力10eに周期的なりロックCLK信号が加えられる
。Cベクトル入力のスカラー乗算した実数及び虚数成分
が、夫々第1及び第2のスカラー乗算手段の出力12−
1c及び12−2cにディジタル・データとして得られ
る。複素数乗算器の出力は、慢素数積Aの夫々直角座標
の形をした実数軸及び虚数軸の項AR及びArを別々に
夫々の出力端子10g及び10fに発生する。
この発明の別の一面として、CORDICデイジタル・
データ乗算器手段11は、第1図に示す様な再帰形乗算
器であるか、又は第2図について説明する様なバイブラ
イン(逐次)形乗算器である。実数軸Iデータ・ワード
又は虚数軸Qデータ・ワードの各々が、手段11に入力
されると、夫々第1又は第2の符号選択手段14−1又
は14−2の入力14−1a又は14−2aに現れる。
入力データの大きさではなく、符号は、夫々符号選択水
力14−1b又は14−2bの符号選択2進制御信号の
状態に応じて、影響を受けないが又は反転される(実効
的に−1の乗算が行なわれる)。最初の回転が+90″
であるか一90°であるかを決定する符号選択ディジタ
ル・データが、符号手段14−1又は14−2の出力1
4−1c又は14−2cに夫々現れる。符号選択手段は
、符号選択手段14−1で示す様に、入力14−1aを
符号ビット反転器16を介して2人カマルチブレクサM
UX手段18の第1の入力18aに接続することによっ
て構成することが出来る。第2のMUX手段の入力18
bが入力14−1aからのもとの符号ディジタル・デー
タ・ワードを受取る。MUX手段の選択入力18cにあ
る符号選択2進信号の状態に応じて、出力18d(並び
に手段の出力14−1c)かもとの極性の入力18b又
は反転した極性の入力18aに接続される。符号選択信
号の2進状態が符号制御手段20によって決定される。
入力のディジタル・データ・ワードが2の補数として表
わされる場合、符号反転及びマルチプレクサ手段は一層
簡単な構造に置換えることが出来る。この構造では、入
力数の各ビットが排他的オア(XOR)ゲートの一方の
入力に接続され、他方の入力がMUX選択入力に接続さ
れる。各々のXORゲートの出力ビットが、順次加算器
チェーンの関連する逐次的な入力に接続され、このチェ
ーンはXOR出力に1を加算する桁上げビットを持って
いる。従って、MUX符号選択入力信号が入力数か又は
その入力の否定の何れかを選択する。この比較的低速の
再帰動作を行なうCORDIC乗算器11では、複数個
(図示の場合は6個)の符号選択信号が必要であり、各
々異なる信号が、節11cから符号制御手段の入力20
gに供給された符号制御データ・ワードに応答して、符
号制御手段の出力(今の場合は出力20a乃至20f)
の内の異なる1つに現れる。このSビットのデータ・ワ
ードは、角度φによって一意的に設定され、それに応じ
た2進パターンを持つデータ・ワードの1組の符号制御
ビットの中をクロックによって前進することに応答して
、プログラム可能な論理配列(PLA)手段の様な論理
手段21の出力21aに供給することが有利である。即
ち、PLA手段の1つの入力21bが、入力lie及び
装置の入力10eからCLKパルスを受取り、PLA手
段の2番目の入力21Cが(CORD I C手段の入
力lid及び装置の入力10dを介して)第2の複素数
の位相角φ情報を受取る。再帰形実施例の動作は、特に
位相入力10dでデータが安定した時、n回の動作の各
々に対して別々のCLKパルスが発生し、この動作を加
算して、所望の積の成分を発生する。この為、各々のク
ロック中パルスにより、次の最も小さいC0RDI C
角度αi=tanが、合計の回転角度θに対して正又は
負の何れの寄与になるかり決定され、こうして入力20
gの符号制御ワードにより、出力20a乃至2Ofの全
ての符号制御ビットを設定する。
第1段(+/−90°選択段)の出力に出る符号選択の
実数及び虚数のPビットのデータ・ワドが、第1及び第
2の符号選択手段の出力14−1c又は14−2cに現
れ、第1の累算器ACCUM手段22−1の入力22−
1a又は第2のACCUM手段22−2の入力22−2
aの所で、(n−1)再帰段に供給される。この入力デ
ータが夫々第2の入力22−1b又は22−2bの(P
+2) ビットのデータ・ワードと加算され、第1又は
第2の累算器の出力22−1c又は222Cに(P+2
)ビットのディジタル・データ・ワードを発生する。こ
の出力データ・ワードがデータ節111又はIIQの夫
々1つに現れる。
節111からのデータ・ワードが、第3の符号選択手段
14−3の第1の入力14−3aに現れる。
この符号選択手段は、符号制御手段20の第3の出力2
0cから、符号制御信号を入力14−3bに受取る。節
11Qのデータ・ワードが第4の符号選択手段14−4
の第1の入力14−4aに現れる。この符号選択手段は
、符号制御手段20の第4の出力20dからの符号制御
信号を受取る符号選択入力14−4bを持っている。「
実数」チャンネルからの符号選択して累算したデータ・
ワードが出力14−3cに現れ、第1のシフタ手段24
−1の入力24−1aに結合され、「虚数」チャンネル
からの符号選択して累算したデータ・ワードか出力14
−4cに現れ、第2のシック手段24−2の入力24−
2aに結合される。各々のシフタ手段24は例えば入力
24−1b又は24−2bの様なシフト制御入力を持っ
ていて、これが符号制御手段20の第5又は第6の出力
20e又は20fから夫々シフト制御パルスを受取る。
シフト制御入力24−1b又は24−2bの一方に出る
各々のパルスに応答して、シフタ24−1又は24−2
にあるデータ・ワードが1つの2進位置だけ右に回転し
又はシフトし、1ビットだけシフトしたデータが第1の
シフタ手段の出力24−1c又は第2のバーレル・シフ
タ手段の出力24−20に現れ、夫々第2のディジタル
加算器手段26−2又は第1のディジタル加算器手段2
6−1の第1の入力26−2a又は26−1aに夫々交
差結合される。これらの加算器手段の別の入力2B−2
b又は2B−1bが、夫々節11Q又は111からのデ
ィジタル・データ・ワードを受取る。第1の加算器手段
の出力26−1cの加算データのワードは、(P+2)
ビットのデータ・ワードQ′であって、C0RD I 
C乗算器手段の第1の出力11−fに供給され、その後
乗算器の出力10fに供給されるが、第2の加算器手段
の出力26−2cのディジタル・データφワードは、別
の(P+2)ビットの信号1′であって、CORDIC
の手段の出力11g及び乗算器手段の出力10gに供給
される。Q′データ信号が累算器の入力22−1bに戻
され l / データが累算器の入力22−2bに戻さ
れる。
動作について説明すると、式(3a)及び(3bは次の
様に書き直すことが出来る。
I′曙K(θ)(1+ξ+Q/2) (4a) Q’  −K  (θ)  (Q−ξ、l/2)(4b
) 前に述べた様に、所望の計算精度に対して選ばれた繰返
し回数nに関係する一定値であるから、共通の倍数K(
θ)swcosθは無視することが出来る。入力の回転
角度φは、任意の角度であるが、これが最初に式(3a
)及び(3b)を充たす1組の回転角度αiに分解され
る。各々の回転で倍率は異なるが、倍率の大きさは回転
の符号には無関係であり、従って、各々の段で符号が異
なるが、同じ大きさの回転を使う一定数の回転に対して
は、全体的な倍率は、合計の回転の角度に無関係であり
、倍率を無視してもよいし、或いは何回もの回転順序の
終りに加えてもよい。これは、複素散積の位相が重要で
ある様な用途で重要なことである(第2のff1Bの大
きさが1であって、スカラー乗算器12を省略すること
が出来、この為、乗算器を全く必要としない場合、特に
重要である)。例として、n−8段の回転は、下記の表
1の符号選択情報を利用して、±0. 6’の精度で実
施することが出来る。
入力のCR及びCIデータ・ワードの符号を修正して、
±90″の回転を実行し、符号を修正したl又はQデー
タを各々サイクルN−8クロック−パルスの第1のクロ
ック−パルス(即ち、クロック・パルス番号C−1)に
応答して、最初にクリアされた関連する累算器手段22
に個別にロードする。この第1パスの累算データが夫々
節txt及び11Qに現れ、夫々関連する加算器手段の
入力26−1b及び26−2bにも現れる。
符号選択手段14−1又は14−2で表わされる第1段
は、項α1の±90°の回転を表わす。
実数及び虚数軸の両方のチャンネルにある累算器22か
ら始まる第2段を(N−1)回の動作に対して再帰形で
利用し、その為(N−1)個のCLKパルスしか必要と
しない。相次ぐクロック・サイクルで、C−2,3,・
・・・・・8の時、各々の累算器の内容を他方の累算器
の内容と加算する。これは、反対のチャンネルの符号選
択手段14−3及び14−4とシフタ信号24−1又は
24−2に於ける符号の修正及び切捨ての後に行なう。
この為、N段の回転には、最初のクロック・サイクルの
最初の(90°)パルスと、その他の(N−1)個のク
ロック・サイクルとを0.9°の精度を持つ回転を行な
う為に必要とする。同様に、N−7段の回転では、I0
 F3°の精度を持つ回転を実施する為に、最初のパル
スと追加の6個のクロック・サイクルとを必要とする。
これに対してN−6段の回転では、3.6°の精度を持
った回転を行なう為に、最初のパルスと5個のクロック
・サイクルとを必要とする。5段の手順では、7.1゜
の精度を持つ回転の為に、最初のパルスと4個のクロッ
ク・サイクルとを必要とする。5用1zの入力データ速
度を利用すると、関連する最低クロツタ周波数Fは弐F
−(N−1)Dで表わされる。
こ\でDは入力データ速度である。こう云う周波数及び
データ速度が、今日側われている最も多い形式の半導体
集積回路で実現するのに見合ったものであることが認め
られよう。
加算器及び累算器は45°からO″までの回転に対処す
る位の奥行のビット密度を持っていなければならないし
、各段の回転に対する倍率をも考慮に入れなければなら
ない。この倍率はI065の値に漸近的に近付く。この
発明では、加算器及び累算器を入力データのビット密度
よりも更に2ビット奥行が深くなる様に設計すれば、両
方の因子に対処する十分な余裕が得られることが判った
即ち、入力信号が7ビットのデータ・ワード(例えば、
P−7)であれば、(P+2)−9ビットの奥行の加算
器、累算器及びシフタを用いる。
この複素数乗算器10は、0. 9mmXI0 2mm
の面積を持つCMOSシリコン回路として集積されてい
るが、1例として、これはAを検出器の複素数出力デー
タ、Bを複素数基準データ入力及びCを複素数検出器デ
ータ入力として、次の式%式%) によって左右されるMタップのベースバンド相関検出器
に利用することが出来る。同様に、M個の点を持つ離散
的なフーリエ変換を左右する方程式4式% こ\でAは複素数順序aの複素数の離散的なフーリエ変
換であり、W、Mは1のM番目の複素数根である。従っ
て、相関の場合でも離散的なフーリエ変換の場合でも、
基本的な処理工程は複素数の乗算−累算であることが理
解されよう。相関の例では、複素数乗算した出力は次の
形である。
AR−BRCR+BI CI       (7)こ\
で添字のRは実数成分、添字の1は虚数成分を表わし、 A、−BRCr−B、CR(8) この複素数乗算は次の様に書換えることが出来る。
AR= l B l  (CRcosφ−CI sln
φ)AI−|B|  (CRsinφ+Ch  Cos
φ)基準uBの係数を実数及び虚数部分のデータではな
く、・大きさIB+及び位相φのデータとしてロードす
る場合、第1図の回路を乗算−累算セルに・対する複素
数乗算器として使うことが出来ることが理解されよう。
即ち、各々の乗算−累算セルに対するC0RDI Cが
同じ段数を持っているから、倍率K(φ)は全てのセル
に対して同じであり、個々の大きさの係数を変更するこ
とにより、又は相関器の最終的な出力に倍率をかけるこ
とにより、倍率をはっきりと考慮に入れることが出来る
。第1図の構造は、CORDICプロセツサ10を使う
ことによって、2つの乗算器及び2つの加算器を省略す
ることが出来るから、従来の複素数乗算セルよりも効率
がずっとよいことが理解されよう。
複素数相関の例では、第1図の回路は、離散的なフーリ
エ変換(D F T)アルゴリズムに於ける複素数乗算
は大きさが1、即ちIWM+−1であるから、離散的な
フーリエ変換(D F T)の計算を効率よくする為に
、更に簡単にすることが出来る。
即ち、両方の乗算器12を省略し、非常に簡単にした回
路が得られる。
計算を更に速くする為、C0RDI C乗算器の再帰形
の第2の部分又は位相回転手段11は、逐次形又はパイ
プライン形のアーキテクチュアに置換えることが出来る
。現在好ましいと考えられる1つのパイプライン形CO
RDIC乗算器11′が第2図に示されている。Pビッ
トの実数軸lデータ・ワードが入力11′ aに加えら
れ、Pビットの虚数軸Qデータ番ワードが入力11′ 
bに加えられる。■又はQデータ・ワードの両方が、±
90@の回転を行なう第1の部分で、符号選択手段14
−1又は14−2の内の関連する一方の作用を別々に受
ける。第2の部分で、複数個(N−1個)の同一の段3
0を用いる。図面に示したN−5の実施例では、4段3
0a乃至30dを使う。
各段30は夫々同一の実数軸及び虚数軸部分30−1及
び30−2を有する。各々の部分で、入力31a又は3
1bのデータ・ワードがシフト手段32a又は32bの
入力と、加算器手段34a又は34bの入力とに結合さ
れる。ビット・シフタ手段32の出力が符号選択手段3
6a又は36bの入力に結合される。各々の符号選択手
段は、インバータ(−1)手段37とマルチプレクサM
UX手段38とで構成される。関連する1つの入力40
a、40b、 ・−−−−−40a”、40b”にある
符号選択制御信号が、符号制御手段41の関連する出力
41b乃至41iに出る。符号選択手段が符号制御入力
11′ CにSビット幅の符号制御データ・ワードを受
取る。このデータ・ワードがMUX手段の出力42 a
、  42 b、 ・・・−・42 a”又は42b1
に於ける各々の信号の2進状態を決定する。出力42a
又は42bが、同じ段の反対側のチャンネル部分に対す
る加算器手段34の関連する第2の入力44b又は44
aに交差結合される。例えば、第1段30aの実数軸部
分の加算器手段に対する第2の入力44aが、虚数軸チ
ャンネルの出力42bに交差結合され、虚数軸チャンネ
ルの加算器手段の第2の入力44bが実数軸チャンネル
の出力42aに接続されると云う様になる。1≦に≦(
N−1)として、各々に番目の段30kにあるシフタ手
段32kが、前の(k−1)番目の段30 (k−1)
にあるシフタ手段よりも1つ多くのビットをシフトする
。k番目の段30には(k−1)ビットのシフト手段3
2kを使う。
従って、第1段のシフタ手段32a、32bは、除数1
の除算機能では、0ビ°ツトだけシフトし、これは通抜
けの接続によって実効的に置換えることが出来る。即ち
、利用しない。第2段30bにあるビット−シフタ32
a′及び32b′が更に1つ余分のビットで除算し、従
ってa■1であり、除数2の除算機能が行なわれる。同
様に、第3段30cでは、シフトは、除数4の除算機能
では、b−2ビットであり、第4段30dのシフタ手段
32a′は、除数8の除算機能に対し、c−3ビットだ
けシフトする。このパイプライン形CORDICアーキ
テクチュアは第1段のパルスしか必要とせず、実質的に
合計N段を通る論理回路の遅延によって設定される速度
を持つ。これは再帰形C0RDI Cの実施例で得られ
る結果よりも、殆んど常にずっと速い計算になる。
次に第3図について説明すると、再帰又はパイプライン
形C0RD I C回転装置11又は11′の一方又は
両方は、2つの複素数の積を発生する別の実施例の乗算
装置10′に利用することが出来る。第1の複素数N1
をIVl 1、φ1−1φ1 Able     (第3a図)として表わし、第2の
複素数N2を1V21、φ2=lA2e′φ2 (第3
b図)として表わせば、(1+十iQ+ )   (I
2 +iQ2 )の積PはP−A+  l IA:  
lel (φ1+φ2)である。第1の複素数は、ベク
トルV1として、その位相角φまたけ回転させて、その
虚数部分を除き、実数の剰余たけで終わる様にすること
、即ち、11= l V+  l = l A+  l
及び(h’−Of、:すルコとが出来ることが理解され
よう。そうする時、位相角φ1を解析して、その結果書
られる回転した第2のベクトルV2′が1A2181 
(φ1+φ2)に等しくなり、この時ff1I+’ と
して利用し得るAllとスカラー積を求めさえすれば、
最終的な複素数の積Pが得られる様に、第2の複素数ベ
クトル■2に加えるべき回転の符号を決定する。
この完全に複素数の乗算器10′が、第1の複素数N1
及び第2の複素数N2の両方を、夫々実数又は同相部分
11又はI2並びに夫々の虚数又は直角位相部分Q1又
はQ2を表わす入力データとして受取る。従って、第1
の複素数N1が同相部分のデータ・ワード11として第
1の入力10′aに現れると共に、直角位相部分のデー
タ・ワードQ1として別の入力10′ bに現れる。第
2の複素数N2は、実数部分のデータ・ワードI2が入
力10′ cに現れ、虚数部分のデータ・ワードQ2が
入力10′ dに現れる。第1の複素数の実数部分■1
及び虚数部分Q1のデータ・ワードが、第1のCORD
IC回転手段11−1に入力され、第2の複素数の実数
部分及び虚数部分12.Q2のデータ・ワードが第2の
C0RD I C回転手段11−2に入力される。両方
の回転装置の回転角度は少なくとも部分的には、回転制
御入力111c又は1l−2cの信号(この入力は第1
図及び第2図の回転装置の符号制御入力11cに対応す
る)によって制御される。節11−IQ (これは第1
図の節11Q又は第2図の節11Q−1に対応する)の
信号が、符号制御手段50の入力50aに印加され、節
11−IQのデータを、第1の回転手段11−1のゼロ
の残留位相角を表わす一定データ・パターンと比較する
ことにより、各々の出力50b、50c (従って、そ
れに接続された回転制御入力1l−1c及び1l−2c
)の2進信号の状態を決定する。即ち、第1の回転手段
の位相角φ1が、近似的にOoの残留位相角まで相次い
でデクレメントされ、これに対して第2の回転手段の合
計位相角φTを同じ回転角度だけインクレメントして、
φ1Σ08の時、φTさ(φ1+φ2)になり、所望の
積Pの位相角になる様にする。節11−IQのデータが
、数N1が直角位相成分を持つこと、即ち角度φ1が0
9ではないことを示す時に、回転が開始される。入力5
0aのデータに応答して、入力1l−1c及び1l−2
cの信号の論理状態を決定して、それに対して入力され
たデータ・ワードの位相から、次のnに対する次の増分
角度θ’ −tan−’  (1/2”)を両方の回転
装置によって加算又は減算させる。従って、両方の回転
装置に於ける回転の符号が、第1の回転手段11−1に
ある剰余“Qlの符号によって決定され、0°に向って
順次強制的に近似する。実際には、第1の回転装置11
−1は、一連のC0RD I C回転を通じて、第1の
複素数データ・ワードN1を、第1の複素数の大きさを
表わす実数A1に変換する。C0RD I C回転が所
定の段数を完了すると、出力50b、50cは変化しな
くなり、出力50dに、そしてその後乗算器の出力10
′ rに泣相回転完了READY信号が発生される。こ
の時、第2の回転手段の12L 出力1l−2dの同相
データ・ワード及びQ2/ 出力1l−2eの直角位相
データ・ワドは夫々次の様に表わされる。
12 −RelV2’ l −Re (lA21e1(
φ1+φ2))−1A2 1cos  (φT)   
           (lla)Q2’ −1m1V
2’ l−Im(lA21ei(φ1+φ2))−lA
21sin(φ丁)                
  (llb)第1の回転装置の1出力1l−1dから
のIA+lデータφワードが、この時節10′ eを介
して、第1及び第2のスカラー乗算手段52゜54の入
力52a、54aに現れる。!2′データ・ワードが第
1の乗算手段52の第2の入力52bに結合され、Q2
’ データ・ワードが第2の乗算手段54の第2の入力
に結合される。第1の乗算手段の出力52cのI′デー
タ・ワードが複素数CORDIC乗算器の出力10′ 
fに現れ、これがRe (P)=lA+  l lA2
 l cos(φT)のデータであり、第2の乗算手段
の出力54cのQ′データ・ワードが複素数C0RD 
IC乗算器の出力50′ gに現れ、これが1m(P)
−IA+  I lA21sin (φT)のデータで
ある。
次に第4図(第4A及び4B図よりなる)について説明
すると、第3図の回転装置11−1及び11−2と符号
制御手段50は、計算速度の速いパイプライン形乗算器
集成体11′で実現することが出来る。第1のN段逐次
形C0RDI C乗算器60aが第1段6O−1aを持
ち、これは入力11′aから符号選択手段6l−1a(
インバータ62及びMUX手段64で構成される)の所
で、実数軸11データを受取ると共に、入力11bから
別の符号選択手段6l−1bの所で虚数軸Q1データを
受取り、符号選択ビット(この第1段に対するものだけ
)がQ1人カデータから取出される。第1の乗算器60
aは、この後略同一の複数個(N−1)の段6O−2a
乃至6O−naを持ち、その各々は、■及びQ部分に対
して、夫々符号選択手段6l−2a乃至5l−na又は
6l−2b乃至6l−nbを含む同一の部分と、除数(
2” ) (7)除算手段66−1乃至66− (n 
−1)と(こ\で1≦■≦Nであり、■は段の番号であ
る)、加算器手段68−1乃至68−(n−1)とを持
っている。第2のN段の逐次形CORDIC乗算器60
bは、第1段6O−1bと(N−1)個の後続の同一の
段6O−2b乃至6〇−nbとの同じ構成を有する。(
第1の乗算器60aの各段のMUX符号選択制御入力6
5−1乃至65−nの)符号ビットが、複数個(N個)
の論理インバータ70−1乃至70−nの内の関連する
1つによって反転され、この為、第1のCOR6Ic乗
算器60aの任意の段のQ出力の符号ビットが、両方の
乗算器60a及び60bの次段の回転の符号を決定する
。従って、符号ビットが2つのパイプライン形CORD
IC構造60a及び60bの間で反転され、この為、第
1の乗算器60aが第1の入力データ(1+及びQ+)
を出力I′及びQ′に於ける(φ1+φ2)に等しい位
相φTへ回転し、第2の乗算器60bが第2の入力デー
タをゼロの位相に回転し、実数出力11′eにIAl 
1の大きさを発生する(そして剰余出力11′ rに略
ゼロの大きさの剰余を発生する。
このQ出力は二重回転装置の誤差の目安である)。
I″及びQ″データには、その乗算を必要とする場合、
第3図の実施例の様に、手段52及び54によってIA
l 1データのスカラー乗算を行なうことが出来る。
これらの複素数乗算器構造の各々は、従来の複素数乗算
器に比べて多数の利点がある。第1に、最も重要なこと
であるが、位相と振幅の精度を切離すことが出来る。振
幅の切捨てが最終的な答の位相の精度に影響しないし、
逆も真である。第2に、こう云う実施例の複素数乗算器
10’  (1対のC0RDIG装置11−1及び11
−2を用いるか又は1個の装置1I″を用いる)は、中
間段階の結果を生ずるが、これは用途によっては役に立
つことがある。例えば、これらの構造に対する一方の入
力が単に他方の入力の共役複素数である場合、こう云う
装置は同時に出力10′又は10′eに入力信号の振幅
を出し、出力10′ f又は10′fに入力信号のエネ
ルギを出す。従って、この実施例は大きさ及びエネルギ
の同時の検出器として使うことが出来る。第3図及び第
4図の構造は、他のどの乗算器の構成よりも、複素数乗
算に対してずっと融通性のある方式を表わす。
この発明の新規な複素数CORDIC乗算器の現在好ま
しいと考えられる幾つかの形式並びにその中で用いられ
るCORDIC回転装置の形式を例として説明したが、
当業者には、種々の変更が考えられよう。従って、この
発明は特許請求の範囲によって限定されるものであって
、こ\で説明した好ましい実施例の説明によって同等制
約されないことを承知されたい。
【図面の簡単な説明】
第1図はこの発明の第1の形式の複素数ディジタル乗算
器の簡略ブロック図、第2図は第1図の複素数乗算器に
示した再帰形C0RD I C手段の代りに使われる逐
次形又はパイプライン形CORDIC手段の現在好まし
いと考えられる別の実施例の簡略ブロック図、第3図は
この発明の別の形式の複素数ディジタル乗算器のブロッ
ク図、第3a図及び第3b図は第3図の複素数ディジタ
ル乗算器の動作を説明するのに役立つ考えを示すベクト
ル図、第4A及び4B図は両者が合わさって、第3図の
装置に於ける複素数乗算を行なう、現在好ましいと考え
られる実施例の逐次形又はパイプライン形手段を示す簡
略ブロック図である。 [主な符号の説明] 1or、10g:出力 11a、llb:入力 14−I014−2;符号選択手段 14−3.14−4:符号選択手段(再帰部分)22−
I022−27累算器 24−I024−2ニジフタ 26−I026−2:加算器

Claims (1)

  1. 【特許請求の範囲】 1、第1の複素数の実数部分及び虚数部分の各々を表わ
    すディジタル・データ・ワードを受取る手段と、 前記第1の複素数の実数部分及び虚数部分の夫々一方を
    、+90°及び−90°の内の選ばれた一方の第1の角
    度増分α_1だけ回転させて、ディジタル・データ・ワ
    ードI及びQの内の関連する一方を形成する手段を含む
    第1の部分と、 該第1の部分から受取ったI及びQディジタル・データ
    ・ワードの各々を、Nを2より大きい正の整数として、
    夫々角度α_1より小さいが、次の増分角度α_i_+
    _1よりも大きな、2≦i≦Nに対する正及び負の増分
    角度α_iの選ばれた一方だけ再帰的に回転修正する手
    段を含む再帰部分と、振幅/回転角度の形で表わした第
    2の複素数の回転角度φ部分を表わすディジタル・デー
    タ・ワードを受取り、前記第1の部分に於ける角度α_
    1並びに前記再帰部分に於ける増分角度α_iの全ての
    符号を制御して、第2の複素数の回転角度φを近似する
    手段と、 N個の全ての角度にわたる回転修正の後、第1部及び第
    2部を、虚数部分の回転出力ディジタル・データ・ワー
    ドA_I及び実数部分の回転出力ディジタル・データ・
    ワードA_Rとして夫々供給する手段とを有する再帰C
    ORDIC回転装置。 2、前記第1の角度増分α_1だけ回転する手段が、前
    記符号を制御する手段からの第1及び第2の符号制御信
    号の夫々1つに応答して、前記第1の複素数の実数部分
    及び虚数部分の夫々に対して同じ符号及び反転した符号
    の一方を選択する手段を含む請求項1記載の再帰COR
    DIC回転装置。 3、各々の増分角度が、n=i−2として、α_i=t
    an^−^1(2^−^n)である請求項2記載の再帰
    CORDIC回転装置。 4、iが8未満である請求項3記載の再帰CORDIC
    回転装置。 5、再帰部分の回転する手段が、前記第1の部分からの
    1ディジタル・データ・ワード及び逐次的にシフトさせ
    た(N−1)個の第1のディジタル・データ・ワードの
    内の夫々逐次的な1つを累算する第1の手段と、前記第
    1の部分からのQディジタル・データ・ワード及び逐次
    的にシフトさせた(N−1)個の第2のディジタル・デ
    ータ・ワードの逐次的な1つを累算する第2の手段と、
    前記符号を制御する手段からの第1の符号選択手段符号
    制御信号に応答して、その時前記第1の累算手段から供
    給されるディジタル・データ・ワードの符号を選択的に
    反転する第1の符号選択手段と、前記符号を制御する手
    段からの第2の符号選択手段符号制御信号に応答して、
    その時前記第2の累算手段から供給されるディジタル・
    データ・ワードの符号を選択的に反転する第2の符号選
    択手段と、前記符号を制御する手段から第1のシフト手
    段制御信号が発生する度に、前記第1の符号選択手段か
    らのディジタル・データ・ワードのビットを予定の方向
    に1ビットだけシフトする第1の手段と、前記符号を制
    御する手段からの第2のシフト手段制御信号が発生する
    度に、前記第2の符号選択手段からのディジタル・デー
    タ・ワードのビットを予定の方向に1ビットだけシフト
    する第2の手段と、前記第1の累算手段及び前記第2の
    シフト手段からのディジタル・データ・ワードを加算し
    て、前記第1の累算手段に対する逐次的にシフトさせた
    (N−1)個の第1のディジタル・データ・ワードの逐
    次的な1つを供給すると共に、A_I出力ディジタル・
    データ・ワードとして供給する第1の加算器手段と、前
    記第2の累算手段及び前記第1のシフト手段からのディ
    ジタル・データ・ワードを加算して、前記第2の累算手
    段に対する逐次的にシフトさせた(N−1)個の第2の
    ディジタル・データ・ワードの逐次的な1つを供給する
    と共に、A_R出力ディジタル・データ・ワードとして
    供給する第2の加算器手段とを有し、前記符号を制御す
    る手段は別の出力を持ち、受取った回転角度φデータに
    応答して、(N−1)回の再帰の各々に対し、第1及び
    第2の符号制御手段及び第1及び第2のシフト手段の符
    号制御信号が全て前記別の出力に供給される様にした請
    求項3記載の再帰CORDIC回転装置。 6、第1の複素数の実数部分及び虚数部分の各々を表わ
    すディジタル・データ・ワードを受取る手段と、 第1及び第2の複素数の実数及び虚数部分の各々1つを
    、+90°及び−90°の内の選ばれた一方の第1の角
    度増分α_1だけ別々に回転させる段手段を含む第1の
    部分と、 Nを2より大きい正の整数として、複数個(N−1個)
    の段手段を含む第2の部分であって、各々の段手段は1
    対の入力の実数及び虚数部分からなるディジタル・デー
    タ・ワードを、2≦i≦Nとして、正及び負の増分角度
    α_iの内の選ばれた一方だけ回転修正し、各々のi番
    目の段手段の増分角度α_iは角度α_1より小さいが
    、次の段手段の増分角度α_i_+_1より大きく、第
    1の段手段に対する入力ディジタル・データ・ワードは
    前記第1の部分から受取り、他の任意の段手段に対する
    入力ディジタル・データ・ワードは直前の段手段から受
    取る様になっている当該第2の部分と、振幅/回転角度
    の形で表わした第2の複素数の回転角度φ部分を表わす
    ディジタル・データ・ワードを受取って、前記第1の部
    分に於ける角度α_1並びに前記複数個の段手段に於け
    る全ての増分角度α_iの符号を制御して、第2の複素
    数の回転角度φを近似する手段と、 N個の角度全部にわたる回転修正の後、第1部及び第2
    部を、虚数部分の回転出力ディジタル・データ・ワード
    及び実数部分の回転出力ディジタル・データ・ワードと
    して夫々供給する手段とを有する逐次CORDIC回転
    装置。 7、前記第1の角度増分α_1だけ回転する手段が、前
    記符号制御手段からの第1の増分符号制御信号に応答し
    て、第1及び第2の複素数の実数及び虚数部分の選ばれ
    た一方に対して同じ符号及び反転した符号の内の一方を
    選択する手段を含む請求項6記載の逐次CORDIC回
    転装置。 8、各々の増分角度が、n=i−2として、α_i=t
    an^−^1(2^−^n)である請求項7記載の逐次
    CORDIC回転装置。 9、iが8未満である請求項8記載の逐次CORDIC
    回転装置。 10、各々の段手段が、該段手段に入力される実数部分
    のディジタル・データ・ワード及び虚数部分のディジタ
    ル・データ・ワードの両方を受取る手段と、前記実数部
    分の入力ディジタル・データ・ワードのビットを予定の
    方向にnビットだけシフトする第1の手段と、虚数部分
    の入力ディジタル・データ・ワードのビットを前記予定
    の方向にnビットだけシフトする第2の手段と、前記符
    号を制御する手段からの第1の符号選択手段制御信号に
    応答して、前記第1のシフト手段からのシフト済みディ
    ジタル・データ・ワードの符号を選択的に反転する第1
    の符号選択手段と、前記符号を制御する手段からの第2
    の符号選択手段制御信号に応答して、第2のシフト手段
    からのシフト済みディジタル・データ・ワードの符号を
    選択的に反転する第2の符号選択手段と、実数部分の入
    力ディジタル・データ・ワード及び前記第2の符号選択
    手段からのディジタル・データ・ワードを加算して当該
    段手段からの実数部分の出力ディジタル・データ・ワー
    ドを発生する第1の加算器手段と、虚数部分の入力ディ
    ジタル・データ・ワード及び前記第1の符号選択手段か
    らのディジタル・データ・ワードを加算して、当該段手
    段からの虚数部分の出力ディジタル・データ・ワードを
    発生する第2の加算器手段とを有し、前記符号を制御す
    る手段は、(N−1)個の段手段全部に対する全ての第
    1及び第2の符号制御手段制御信号が、受取った回転角
    度φデータに応答して供給される付加的な出力を持って
    いる請求項8記載の逐次CORDIC回転装置。 11、実数部分のディジタル・データ・ワードCR及び
    虚数部分のディジタル・データ・ワードC_Iとして表
    わされた第1の複素数Cと、大きさ|B|のディジタル
    ・データ・ワード及び位相角φのディジタル・データ・
    ワードとして表わされた第2の複素数Bとの積である出
    力ディジタル・データを発生する装置に於て、 位相角φのディジタル・データ・ワードに応答して、入
    力の実数部分のディジタル・データ・ワードI及び入力
    の虚数部分のディジタル・データ・ワードQを夫々回転
    修正して、出力ディジタル・データ・ワードQ′及び出
    力ディジタル・データ・ワードI′を夫々求めるCOR
    DIC手段と、少なくとも前記第1の複素数の実数部分
    のディジタル・データ・ワードC_R及び虚数ディジタ
    ル・データ・ワードC_Iを受取って、夫々入力I及び
    Qディジタル・データ・ワードを供給する入力手段と、 夫々出力Q′、及びI′ディジタル・データ・ワードを
    受取って、出力ディジタル・データ・ワードの虚数部分
    A_I及び出力ディジタル・データ・ワードの実数部分
    A_Rを夫々発生する出力手段とを有する装置。 12、CORDIC手段が再帰CORDIC回転装置で
    ある請求項11記載の装置。 13、前記出力手段が、第2の複素数の大きさ|B|の
    ディジタル・データ・ワードを受取る手段と、実数部分
    のディジタル・データ・ワードI′及び虚数部分のディ
    ジタル・データ・ワードQ’の夫々1つに大きさ|B|
    のディジタル・データ・ワードのスカラー積を求めて、
    実数部分A_R及び虚数部分A_Iのディジタル・デー
    タ・ワードを求める手段とを有する請求項12記載の装
    置。 14、CORDIC手段が逐次CORDIC回転装置で
    ある請求項11記載の装置。 15、前記出力手段が更に、第2の複素数の大きさ|B
    |のディジタル・データ・ワードを受取る手段と、実数
    部分のディジタル・データ・ワードI′及び虚数部分の
    ディジタル・データ・ワードQ′の夫々1つに大きさ|
    B|のディジタル・データ・ワードのスカラー積を求め
    て、夫々実数部分A_R及び虚数部分A_Iのディジタ
    ル・データ・ワードを求める手段とを有する請求項14
    記載の装置。 16、実数部分のディジタル・データ・ワードI_1及
    び虚数部分のディジタル・データ・ワードQ_1として
    表わされた第1の複素数と、実数部分のディジタル・デ
    ータ・ワードI_2及び虚数部分のディジタル・データ
    ・ワードQ_2として表わされた第2の複素数との積で
    ある出力ディジタル・データを発生する装置に於て、 N個の期間の各々i番目の期間に、別のディジタル・デ
    ータ・ビットの状態に応答して、入力の実数部分のディ
    ジタル・データ・ワードI及び入力の虚数部分のディジ
    タル・データ・ワードQを、α_1を±90°の選ばれ
    た一方、そしてNを2より大として、2≦i≦Nに対し
    、α_iをtan^−^1(2^−^n)(但しn=i
    −2)として、角度α_iだけ回転修正して、出力ディ
    ジタル・データ・ワードQ′及び出力ディジタル・デー
    タ・ワードI′を夫々求める第1及び第2のCORDI
    C手段を有し、第1のCORDIC手段は夫々のデータ
    ・ワードI_1及びQ_1を夫々実数部分及び虚数部分
    のデータ・ワードとして受取って、少なくとも1つの出
    力データ・ワードI_0及び有効な出力角度が略ゼロで
    ある時に符号の状態を変える信号を発生し、第2のCO
    RDIC手段は夫々のデータ・ワードI_2及びQ_2
    を実数部分及び虚数部分のデータ・ワードとして受取っ
    て、夫々実数部分及び虚数部分の出力データ・ワードI
    ″及びQ″を発生し、 更に、各々のi番目の期間に、符号状態信号が符号を変
    えるまで、両方のCORDIC手段によって出力ディジ
    タル・データ・ワードを回転修正させる様に選ばれたデ
    ータ状態を前記別のディジタル・データ・ビットに持た
    せる手段と、夫々出力I″及びQ″ディジタル・データ
    ・ワードを受取って、出力ディジタル・データの虚数部
    分I′及び出力ディジタル・データ・ワードの実数部分
    Q′を夫々発生する出力手段とを有する装置。 17、各々のCORDIC手段が再帰CORDIC回転
    装置である請求項16記載の装置。 18、前記出力手段が更に、第1のCORDIC手段の
    出力データ・ワードI_0を受取る手段と、実数部分の
    出力I″データ・ワード及び虚数部分の出力Q″ディジ
    タル・データ・ワードにI_0データ・ワードを乗じた
    スカラー積を求めて、夫々実数部分I′及び虚数部分Q
    ′のディジタル・データ・ワードを求める手段とを有す
    る請求項17記載の装置。 19、実数部分のディジタル・データ・ワードI_1及
    び虚数部分のディジタル・データ・ワードQ_1として
    表わされた第1の複素数と、実数部分のディジタル・デ
    ータ・ワードI_2及び虚数部分のディジタル・データ
    ・ワードQ_2として表わされた第2の複素数との積で
    ある出力ディジタル・データを発生する装置に於て、 α_1を±90°の選ばれた一方とし、Nを2より大と
    して、2≦i≦Nに対し、n=i−2として、α_iが
    tan^−^1(2^−^n)である様な角度α_iだ
    け、入力の実数部分のディジタル・データ・ワードI及
    び入力の虚数部分のディジタル・データ・ワードQを回
    転修正して、出力ディジタル・データ・ワードQ′及び
    出力ディジタル・データ・ワードI′を求める第1及び
    第2のCORDIC手段を有し、ここで各々の角度の符
    号は関連する符号ビットの論理状態に応答し、第1の手
    段の各々の角度に対する符号ビットは、その角度に対す
    る計算段の1対のデータ入力の一方にその時存在するデ
    ータ・ビットであり、 更に、入力ビットの論理状態を反転する複数個(N個)
    の手段を有し、各々i番目の反転する手段は、関連する
    i番目の角度に対する第1の手段に対する段から符号ビ
    ットを受取って、第2の手段の同じi番目の段に対して
    別の符号ビットを供給し、 前記第2の手段は夫々のデータ・ワードI_2及びQ_
    2を夫々実数部分及び虚数部分のデータ・ワードとして
    受取って、第2の複素数の大きさ|A_2|である少な
    くとも1つの出力ワードを発生し、 更に、夫々出力I′及びQ′ディジタル・データ・ワー
    ドを受取って、出力ディジタル・データの虚数部分Q″
    及び出力ディジタル・データ・ワードの実数部分I″を
    夫々発生する出力手段を有する装置。 20、前記出力手段が更に、第2のCORDIC手段の
    出力データ・ワード|A_2|を受取る手段と、実数部
    分の出力I″ディジタル・データ・ワード及び虚数部分
    の出力Q″ディジタル・データ・ワードの夫々1つと|
    A_2|データ・ワードとのスカラー積を求めて、夫々
    実数部分I′及び虚数部分Q′のディジタル・データ・
    ワードを求める手段とを有する請求項19記載の装置。
JP13629589A 1988-05-31 1989-05-31 Cordic複素数乗算器 Expired - Fee Related JP3283504B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US07/200,491 US4896287A (en) 1988-05-31 1988-05-31 Cordic complex multiplier
US200,491 1988-05-31

Publications (2)

Publication Number Publication Date
JPH0229821A true JPH0229821A (ja) 1990-01-31
JP3283504B2 JP3283504B2 (ja) 2002-05-20

Family

ID=22741948

Family Applications (1)

Application Number Title Priority Date Filing Date
JP13629589A Expired - Fee Related JP3283504B2 (ja) 1988-05-31 1989-05-31 Cordic複素数乗算器

Country Status (7)

Country Link
US (1) US4896287A (ja)
JP (1) JP3283504B2 (ja)
KR (1) KR0146334B1 (ja)
CN (1) CN1017283B (ja)
DE (1) DE3917059A1 (ja)
FR (1) FR2632088B1 (ja)
GB (1) GB2220090B (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6922451B1 (en) 1999-04-15 2005-07-26 Nec Corporation Frequency shifting circuit and method
JP2012147873A (ja) * 2011-01-18 2012-08-09 Canon Inc 信号処理回路及び超音波診断装置

Families Citing this family (70)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4983970A (en) * 1990-03-28 1991-01-08 General Electric Company Method and apparatus for digital phased array imaging
EP0453641B1 (de) * 1990-04-20 1997-03-12 Siemens Aktiengesellschaft CORDIC-Prozessor für Vektordrehungen in Carry-Save-Architektur
US5115492A (en) * 1990-12-14 1992-05-19 General Electric Company Digital correlators incorporating analog neural network structures operated on a bit-sliced basis
US5167008A (en) * 1990-12-14 1992-11-24 General Electric Company Digital circuitry for approximating sigmoidal response in a neural network layer
US5121364A (en) * 1991-08-07 1992-06-09 General Electric Company Time frequency control filter for an ultrasonic imaging system
US5142649A (en) * 1991-08-07 1992-08-25 General Electric Company Ultrasonic imaging system with multiple, dynamically focused transmit beams
DE4126953C2 (de) * 1991-08-14 1995-04-13 Fraunhofer Ges Forschung Schaltungsanordnung zur Durchführung des CORDIC-Algorithmus
US5235982A (en) * 1991-09-30 1993-08-17 General Electric Company Dynamic transmit focusing of a steered ultrasonic beam
US5291892A (en) * 1991-11-04 1994-03-08 General Electric Company Ultrasonic flow imaging
US5269307A (en) * 1992-01-31 1993-12-14 Tetrad Corporation Medical ultrasonic imaging system with dynamic focusing
JP2818345B2 (ja) * 1992-12-25 1998-10-30 株式会社東芝 ディジタル正弦波発生回路
US5349525A (en) * 1993-01-08 1994-09-20 General Electric Company Color flow imaging system utilizing a frequency domain wall filter
US5349524A (en) * 1993-01-08 1994-09-20 General Electric Company Color flow imaging system utilizing a time domain adaptive wall filter
US5465222A (en) * 1994-02-14 1995-11-07 Tektronix, Inc. Barrel shifter or multiply/divide IC structure
US5473654A (en) * 1994-06-24 1995-12-05 General Electric Company Backprojection for x-ray CT system
US5549111A (en) * 1994-08-05 1996-08-27 Acuson Corporation Method and apparatus for adjustable frequency scanning in ultrasound imaging
US5675554A (en) * 1994-08-05 1997-10-07 Acuson Corporation Method and apparatus for transmit beamformer
US5928152A (en) * 1994-08-05 1999-07-27 Acuson Corporation Method and apparatus for a baseband processor of a receive beamformer system
US5685308A (en) * 1994-08-05 1997-11-11 Acuson Corporation Method and apparatus for receive beamformer system
US5555534A (en) * 1994-08-05 1996-09-10 Acuson Corporation Method and apparatus for doppler receive beamformer system
US5793701A (en) * 1995-04-07 1998-08-11 Acuson Corporation Method and apparatus for coherent image formation
US6029116A (en) * 1994-08-05 2000-02-22 Acuson Corporation Method and apparatus for a baseband processor of a receive beamformer system
US5581517A (en) * 1994-08-05 1996-12-03 Acuson Corporation Method and apparatus for focus control of transmit and receive beamformer systems
AU3361095A (en) 1994-08-05 1996-03-04 Acuson Corporation Method and apparatus for transmit beamformer system
DE4442959C2 (de) * 1994-12-02 2001-02-08 Sican Gmbh Monolithisch integrierbare Schaltungsanordnung zur komplexen Multiplikation serieller Datenströme
US5668749A (en) * 1995-05-04 1997-09-16 Motorola, Inc. Circuit for performing arithmetic operations in a demodulator
DE69619963T2 (de) * 1995-08-30 2002-11-21 Koninkl Philips Electronics Nv Signalprozessor mit reduzierter komplexität und empfänger mit einem derartigen signalprozessor
US5822967A (en) * 1996-06-11 1998-10-20 Vermeer Manufacturing Co. Baler with swing arm bale wrapper
US5802111A (en) * 1997-02-24 1998-09-01 Motorola, Inc. Complex constellation point multiplier
US6192089B1 (en) 1998-08-07 2001-02-20 Motorola, Inc. Electronic circuit and method for automatic frequency control
JP3201364B2 (ja) * 1998-11-20 2001-08-20 日本電気株式会社 角度演算回路
US6123671A (en) * 1998-12-31 2000-09-26 General Electric Company Method and apparatus for distributed, agile calculation of beamforming time delays and apodization values
US6349317B1 (en) * 1999-03-13 2002-02-19 Vitit Kantabutra Efficient radix-4 CORDIC vector rotators and computers of sine and cosine functions
AU3608199A (en) * 1999-04-23 2000-11-10 Nokia Networks Oy Qam modulator
US6560536B1 (en) 1999-07-12 2003-05-06 Eagle-Eye, Inc. System and method for rapid telepositioning
US8255149B2 (en) 1999-07-12 2012-08-28 Skybitz, Inc. System and method for dual-mode location determination
US6480788B2 (en) 1999-07-12 2002-11-12 Eagle-Eye, Inc. System and method for fast acquisition reporting using communication satellite range measurement
US20040143392A1 (en) * 1999-07-12 2004-07-22 Skybitz, Inc. System and method for fast acquisition reporting using communication satellite range measurement
AU2001283181A1 (en) 2000-08-09 2002-02-18 Skybitz, Inc System and method for fast code phase and carrier frequency acquisition in gps receiver
AU2001283182A1 (en) * 2000-08-09 2002-02-18 Skybitz, Inc. Frequency translator using a cordic phase rotator
AU2001289045A1 (en) * 2000-09-08 2002-03-22 Avaz Networks Hardware function generator support in a dsp
US6452961B1 (en) 2000-09-12 2002-09-17 Interstate Electronics Corporation Massively paralleled sequential test algorithm
US6466958B1 (en) 2000-09-12 2002-10-15 Interstate Electronics Corporation, A Division Of L3 Communications Corporation Parallel frequency searching in an acquisition correlator
JP2004509358A (ja) * 2000-09-18 2004-03-25 スカイビッツ,インコーポレイテッド Gps受信機におけるコード位相ならびにキャリア周波数の高速捕捉システム並びに方法
JP3870105B2 (ja) * 2002-02-22 2007-01-17 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 逆投影方法およびx線ct装置
US7082451B2 (en) * 2002-09-09 2006-07-25 Freescale Semiconductor, Inc. Reconfigurable vector-FFT/IFFT, vector-multiplier/divider
US7889787B2 (en) * 2003-08-04 2011-02-15 Supertex, Inc. Ultrasound transmit beamformer integrated circuit and method
DE102005045519A1 (de) * 2005-09-23 2007-03-29 Newlogic Technologies Ag Verfahren und Vorrichtung zur FFT Berechnung
EP2541432A1 (en) 2005-10-07 2013-01-02 Altera Corporation Data input and output in systolic array processors
CA2641597C (en) * 2006-02-06 2013-11-26 Masco Corporation Of Indiana Power sprayer
US9547327B2 (en) 2006-11-09 2017-01-17 Alan N. Willson, Jr. Excess-fours processing in direct digital synthesizer implementations
US8131793B2 (en) * 2006-11-09 2012-03-06 Pentomics, Inc. Efficient angle rotator configured for dynamic adjustment
US9244483B1 (en) 2006-11-09 2016-01-26 Pentomics, Inc. Excess-fours processing in direct digital synthesizer implementations
US9268529B2 (en) 2006-11-09 2016-02-23 Pentomics, Inc. Efficient angle rotator configured for dynamic adjustment
US8484278B2 (en) * 2007-05-11 2013-07-09 Synopsys, Inc. Digital architecture for DFT/IDFT hardware
WO2009040957A2 (en) * 2007-09-26 2009-04-02 Nec Corporation Cordic-based fft and ifft apparatus and method
EP2278714B1 (en) * 2009-07-02 2015-09-16 Nxp B.V. Power stage
US9363068B2 (en) * 2010-08-03 2016-06-07 Intel Corporation Vector processor having instruction set with sliding window non-linear convolutional function
US9157940B2 (en) * 2011-02-09 2015-10-13 Smart Energy Instruments, Inc. Power measurement device
CN102799565A (zh) * 2011-05-26 2012-11-28 联芯科技有限公司 坐标旋转数字计算的改进方法及其装置
US8706794B1 (en) * 2011-08-23 2014-04-22 Gregory K. Fleizach No-multiply digital signal processing method
US8867592B2 (en) 2012-05-09 2014-10-21 Nxp B.V. Capacitive isolated voltage domains
US9489342B2 (en) * 2012-12-24 2016-11-08 Intel Corporation Systems, methods, and computer program products for performing mathematical operations
GB2513882A (en) * 2013-05-08 2014-11-12 Nordic Semiconductor Asa Digital Radios
JP6335661B2 (ja) 2014-06-04 2018-05-30 キヤノン株式会社 演算装置およびその制御方法
KR102268110B1 (ko) * 2014-08-05 2021-06-22 삼성전자주식회사 데이터를 변조하는 방법 및 장치 및 기록 매체
CN110178321B (zh) * 2017-01-17 2021-05-18 华为技术有限公司 信号发射方法及装置、发射机、信号传输系统
US11385336B2 (en) * 2018-07-31 2022-07-12 Maxim Integrated Products, Inc. Time of flight sensors and sensing methods
US11764940B2 (en) 2019-01-10 2023-09-19 Duality Technologies, Inc. Secure search of secret data in a semi-trusted environment using homomorphic encryption
CN113711157A (zh) 2019-02-20 2021-11-26 优创半导体科技有限公司 用于使用坐标旋转数字计算机(cordic)对浮点三角函数进行硬件高效自适应计算的装置和方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3868680A (en) * 1974-02-04 1975-02-25 Rockwell International Corp Analog-to-digital converter apparatus
US3976869A (en) * 1974-09-27 1976-08-24 The Singer Company Solid state resolver coordinate converter unit
US3926367A (en) * 1974-09-27 1975-12-16 Us Navy Complex filters, convolvers, and multipliers
US3927312A (en) * 1974-10-31 1975-12-16 Us Army Vector rotator
US4231102A (en) * 1978-12-21 1980-10-28 Raytheon Company Cordic FFT processor
US4354249A (en) * 1980-03-24 1982-10-12 Motorola Inc. Processing unit for multiplying two mathematical quantities including at least one complex multiplier
US4344151A (en) * 1980-04-21 1982-08-10 Rockwell International Corporation ROM-Based complex multiplier useful for FFT butterfly arithmetic unit
JPS60501186A (ja) * 1983-02-02 1985-07-25 ジ−イ−シ− アビニヨニクス リミテツド 2進デジタルプロセッサ
DE3312796A1 (de) * 1983-04-09 1984-10-11 Licentia Patent-Verwaltungs-Gmbh, 6000 Frankfurt Digitaler oszillator zur erzeugung komplexer signale
US4680727A (en) * 1984-09-24 1987-07-14 Rockwell International Corporation Complex multiplier for binary two's complement numbers
US4769779A (en) * 1985-12-16 1988-09-06 Texas Instruments Incorporated Systolic complex multiplier

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6922451B1 (en) 1999-04-15 2005-07-26 Nec Corporation Frequency shifting circuit and method
JP2012147873A (ja) * 2011-01-18 2012-08-09 Canon Inc 信号処理回路及び超音波診断装置

Also Published As

Publication number Publication date
US4896287A (en) 1990-01-23
FR2632088B1 (fr) 1994-08-26
CN1017283B (zh) 1992-07-01
JP3283504B2 (ja) 2002-05-20
KR0146334B1 (ko) 1998-09-15
CN1045879A (zh) 1990-10-03
DE3917059A1 (de) 1989-12-07
GB2220090A (en) 1989-12-28
GB2220090B (en) 1992-12-09
FR2632088A1 (fr) 1989-12-01
KR890017608A (ko) 1989-12-16
GB8912385D0 (en) 1989-07-12

Similar Documents

Publication Publication Date Title
JPH0229821A (ja) Cordic複素数乗算器
JP4980834B2 (ja) 算術プロセッサ
US4231102A (en) Cordic FFT processor
US20060282489A1 (en) Hardware function generator support in a DSP
US5457804A (en) Accumulating multiplication circuit executing a double-precision multiplication at a high speed
US4937775A (en) Apparatus for the cross-correlation of a pair of complex sampled signals
US5737253A (en) Method and apparatus for direct digital frequency synthesizer
EP0477011B1 (en) Processor element for calculating accumulation of data, processing unit, and processor
EP0441121A2 (en) Arithmetic operation apparatus for elementary function
US4945505A (en) Cordic apparatus and method for approximating the magnitude and phase of a complex number
CA1258537A (en) Apparatus and method for approximating the magnitude of a complex number
JP3003467B2 (ja) 演算装置
US4841552A (en) Digital phase shifter
US8452830B2 (en) Programmable CORDIC processor with stage re-use
US6011448A (en) Method and apparatus for frequency modulation synthesis
US6681235B2 (en) Frequency synthesizer and Gaussian noise generator using the same
GB2135087A (en) Vector rotator/accumulator
JP3684314B2 (ja) 複素乗算器および複素相関器
US5412588A (en) Digital sine-wave generating circuit
US20220374201A1 (en) Digital Signal Processing Device
JPH0535773A (ja) ベクトル除算方式とその装置
US5889691A (en) Apparatus and method for a multiplier unit with high component utilization
JP2864598B2 (ja) ディジタル演算回路
SU1119006A1 (ru) Устройство дл делени чисел
SU783791A1 (ru) Устройство дл умножени многочленов

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090301

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees