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
Links
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 230000000306 recurrent effect Effects 0.000 abstract 1
- 239000000047 product Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 229920000535 Tan II Polymers 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods 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/49—Computations with a radix, other than binary, 8, 16 or decimal, e.g. ternary, negative or imaginary radices, mixed radix non-linear PCM
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03C—MODULATION
- H03C3/00—Angle modulation
- H03C3/38—Angle modulation by converting amplitude modulation to angle modulation
- H03C3/40—Angle 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods 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/4806—Computations with complex numbers
- G06F7/4812—Complex multiplication
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods 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/4806—Computations with complex numbers
- G06F7/4818—Computations with complex numbers using coordinate rotation digital computer [CORDIC]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods 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/544—Methods 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/5446—Methods 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)を利用した新規な
手段に関する。
複素数の乗算を行なう為に、少なくとも1つの座標回転
ディジタル計算機(CORDI C)を利用した新規な
手段に関する。
いろいろな形式の今日の電子装置では、相関の検出、離
散的なフーリエ変換等の様な比較的複雑な信号処理機能
が利用されている。基本の信号処理素子は乗算−累算セ
ルである。ベースバンド信号処理等を利用する超音波作
像の様なある装置では、完全に複素数の信号を乗算しな
ければならない。電子信号の同様な複素数乗算は、レー
ダ、ソナー等の様な他の多くの信号処理分野で、特にデ
ィジタル信号、更に特に2進形の信号の処理に見られる
。従って、何れもx+jy又はR2Oと云う形をした複
素数を表わす1つのディジタル信号の複素数乗算を実施
するディジタル信号手段を提供することが非常に望まし
い。
散的なフーリエ変換等の様な比較的複雑な信号処理機能
が利用されている。基本の信号処理素子は乗算−累算セ
ルである。ベースバンド信号処理等を利用する超音波作
像の様なある装置では、完全に複素数の信号を乗算しな
ければならない。電子信号の同様な複素数乗算は、レー
ダ、ソナー等の様な他の多くの信号処理分野で、特にデ
ィジタル信号、更に特に2進形の信号の処理に見られる
。従って、何れもx+jy又はR2Oと云う形をした複
素数を表わす1つのディジタル信号の複素数乗算を実施
するディジタル信号手段を提供することが非常に望まし
い。
従来、IREトランザクションズ・オン・エレクトロニ
ック・コンピュータズ誌EC−8、第330頁乃至第3
34頁(1959年)所載のJ。
ック・コンピュータズ誌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対の方程式で表わされる様
な関係を持つ。
」に最初に記載された様な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ビッ
トのシフトによって行なわれる。
とであり、2進数では、このn回の各々に対して1ビッ
トのシフトによって行なわれる。
この為、複素数乗算は、(必要とする場合の)倍率K(
θ)の乗算を別として、1組のシフトレジスタ及び加算
器を用いて実行することが出来る。
θ)の乗算を別として、1組のシフトレジスタ及び加算
器を用いて実行することが出来る。
この基本的なC0RD I C方式を利用して、1対の
複素数を乗算する新規な装置を提供することが非常に望
ましい。
複素数を乗算する新規な装置を提供することが非常に望
ましい。
発明の要約
この発明の現在好ましいと考えられる2つのディジタル
複索数B及びCを乗算するC0RD I C装置は、直
角座標の形(例えば、CR及びC+)で表わした一方の
数の実数データ部分及び虚数データ部分の各々を、極座
標の形で表わした他方の数(例えば、IB+、φ)の位
相角φだけ回転させるために、乗算器のない、再帰形又
はパイプライン逐次形の何れかのN段C0RD I C
回転手段を用いる。回転の後、CORDIC回転手段の
出力の実数又は虚数データ部分の各々に他方の数の大き
さのデータ|B|のスカラー乗算を行なう手段を設ける
。こうして計算された最終的なデータが、積の実数及び
虚数部分である。
複索数B及びCを乗算するC0RD I C装置は、直
角座標の形(例えば、CR及びC+)で表わした一方の
数の実数データ部分及び虚数データ部分の各々を、極座
標の形で表わした他方の数(例えば、IB+、φ)の位
相角φだけ回転させるために、乗算器のない、再帰形又
はパイプライン逐次形の何れかのN段C0RD I C
回転手段を用いる。回転の後、CORDIC回転手段の
出力の実数又は虚数データ部分の各々に他方の数の大き
さのデータ|B|のスカラー乗算を行なう手段を設ける
。こうして計算された最終的なデータが、積の実数及び
虚数部分である。
現在好ましいと考えられる別の複素数乗算CORDIC
装置では、1対の再帰形又はバイブライン逐次形位相回
転手段の各々が、第1及び第2の直角座標形式の複素数
の実数及び虚数部分に作用する。各々の複素数の位相は
、符号検出器で判定して、第1の数がゼロの位相角にな
るまで平等に回転し、これにより出力データの位相が第
1及び第2の複素数の位相角の和に等しくなる様にする
。
装置では、1対の再帰形又はバイブライン逐次形位相回
転手段の各々が、第1及び第2の直角座標形式の複素数
の実数及び虚数部分に作用する。各々の複素数の位相は
、符号検出器で判定して、第1の数がゼロの位相角にな
るまで平等に回転し、これにより出力データの位相が第
1及び第2の複素数の位相角の和に等しくなる様にする
。
必要によっては、全体の複素数の積を決定する為に、積
データのスカラー乗算を使うことが出来る。
データのスカラー乗算を使うことが出来る。
相関の統計等の計算の様な多くの場合、正確な位相情報
を累算することだけが必要であって、位相情報の精度が
保たれている限り、積の振幅項には比較的精度が要求さ
れない。使う場合、各々のスカラー乗算器はシフト及び
累算部分にすることが出来る。
を累算することだけが必要であって、位相情報の精度が
保たれている限り、積の振幅項には比較的精度が要求さ
れない。使う場合、各々のスカラー乗算器はシフト及び
累算部分にすることが出来る。
従って、この発明の目的は、1対のディジタル複素数を
乗算する新規なC0RD I C回路を提供することで
ある。
乗算する新規な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を設ける
場合、その部品は比較的精度の低いものであってよい。
られる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に発生する。
の成分に対する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に現れる。
データ乗算器手段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によって決定される。
力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の全
ての符号制御ビットを設定する。
わされる場合、符号反転及びマルチプレクサ手段は一層
簡単な構造に置換えることが出来る。この構造では、入
力数の各ビットが排他的オア(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つに現れる。
実数及び虚数の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に現れる。
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から夫々シフト制御パルスを受取る。
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に戻さ
れる。
各々のパルスに応答して、シフタ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’の精度で実
施することが出来る。
) 前に述べた様に、所望の計算精度に対して選ばれた繰返
し回数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にも現れる。
±90″の回転を実行し、符号を修正したl又はQデー
タを各々サイクルN−8クロック−パルスの第1のクロ
ック−パルス(即ち、クロック・パルス番号C−1)に
応答して、最初にクリアされた関連する累算器手段22
に個別にロードする。この第1パスの累算データが夫々
節txt及び11Qに現れ、夫々関連する加算器手段の
入力26−1b及び26−2bにも現れる。
符号選択手段14−1又は14−2で表わされる第1段
は、項α1の±90°の回転を表わす。
は、項α1の±90°の回転を表わす。
実数及び虚数軸の両方のチャンネルにある累算器22か
ら始まる第2段を(N−1)回の動作に対して再帰形で
利用し、その為(N−1)個のCLKパルスしか必要と
しない。相次ぐクロック・サイクルで、C−2,3,・
・・・・・8の時、各々の累算器の内容を他方の累算器
の内容と加算する。これは、反対のチャンネルの符号選
択手段14−3及び14−4とシフタ信号24−1又は
24−2に於ける符号の修正及び切捨ての後に行なう。
ら始まる第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個のクロック・サイクルとを必要とする。
最初の(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で表わされる。
った回転を行なう為に、最初のパルスと5個のクロック
・サイクルとを必要とする。5段の手順では、7.1゜
の精度を持つ回転の為に、最初のパルスと4個のクロッ
ク・サイクルとを必要とする。5用1zの入力データ速
度を利用すると、関連する最低クロツタ周波数Fは弐F
−(N−1)Dで表わされる。
こ\でDは入力データ速度である。こう云う周波数及び
データ速度が、今日側われている最も多い形式の半導体
集積回路で実現するのに見合ったものであることが認め
られよう。
データ速度が、今日側われている最も多い形式の半導体
集積回路で実現するのに見合ったものであることが認め
られよう。
加算器及び累算器は45°からO″までの回転に対処す
る位の奥行のビット密度を持っていなければならないし
、各段の回転に対する倍率をも考慮に入れなければなら
ない。この倍率はI065の値に漸近的に近付く。この
発明では、加算器及び累算器を入力データのビット密度
よりも更に2ビット奥行が深くなる様に設計すれば、両
方の因子に対処する十分な余裕が得られることが判った
。
る位の奥行のビット密度を持っていなければならないし
、各段の回転に対する倍率をも考慮に入れなければなら
ない。この倍率はI065の値に漸近的に近付く。この
発明では、加算器及び累算器を入力データのビット密度
よりも更に2ビット奥行が深くなる様に設計すれば、両
方の因子に対処する十分な余裕が得られることが判った
。
即ち、入力信号が7ビットのデータ・ワード(例えば、
P−7)であれば、(P+2)−9ビットの奥行の加算
器、累算器及びシフタを用いる。
P−7)であれば、(P+2)−9ビットの奥行の加算
器、累算器及びシフタを用いる。
この複素数乗算器10は、0. 9mmXI0 2mm
の面積を持つCMOSシリコン回路として集積されてい
るが、1例として、これはAを検出器の複素数出力デー
タ、Bを複素数基準データ入力及びCを複素数検出器デ
ータ入力として、次の式%式%) によって左右されるMタップのベースバンド相関検出器
に利用することが出来る。同様に、M個の点を持つ離散
的なフーリエ変換を左右する方程式4式% こ\でAは複素数順序aの複素数の離散的なフーリエ変
換であり、W、Mは1のM番目の複素数根である。従っ
て、相関の場合でも離散的なフーリエ変換の場合でも、
基本的な処理工程は複素数の乗算−累算であることが理
解されよう。相関の例では、複素数乗算した出力は次の
形である。
の面積を持つ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) この複素数乗算は次の様に書換えることが出来る。
で添字のRは実数成分、添字の1は虚数成分を表わし、 A、−BRCr−B、CR(8) この複素数乗算は次の様に書換えることが出来る。
AR= l B l (CRcosφ−CI sln
φ)AI−|B| (CRsinφ+Ch Cos
φ)基準uBの係数を実数及び虚数部分のデータではな
く、・大きさIB+及び位相φのデータとしてロードす
る場合、第1図の回路を乗算−累算セルに・対する複素
数乗算器として使うことが出来ることが理解されよう。
φ)AI−|B| (CRsinφ+Ch Cos
φ)基準uBの係数を実数及び虚数部分のデータではな
く、・大きさIB+及び位相φのデータとしてロードす
る場合、第1図の回路を乗算−累算セルに・対する複素
数乗算器として使うことが出来ることが理解されよう。
即ち、各々の乗算−累算セルに対するC0RDI Cが
同じ段数を持っているから、倍率K(φ)は全てのセル
に対して同じであり、個々の大きさの係数を変更するこ
とにより、又は相関器の最終的な出力に倍率をかけるこ
とにより、倍率をはっきりと考慮に入れることが出来る
。第1図の構造は、CORDICプロセツサ10を使う
ことによって、2つの乗算器及び2つの加算器を省略す
ることが出来るから、従来の複素数乗算セルよりも効率
がずっとよいことが理解されよう。
同じ段数を持っているから、倍率K(φ)は全てのセル
に対して同じであり、個々の大きさの係数を変更するこ
とにより、又は相関器の最終的な出力に倍率をかけるこ
とにより、倍率をはっきりと考慮に入れることが出来る
。第1図の構造は、CORDICプロセツサ10を使う
ことによって、2つの乗算器及び2つの加算器を省略す
ることが出来るから、従来の複素数乗算セルよりも効率
がずっとよいことが理解されよう。
複素数相関の例では、第1図の回路は、離散的なフーリ
エ変換(D F T)アルゴリズムに於ける複素数乗算
は大きさが1、即ちIWM+−1であるから、離散的な
フーリエ変換(D F T)の計算を効率よくする為に
、更に簡単にすることが出来る。
エ変換(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を使う。
の第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を使う。
び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の実施例で得られ
る結果よりも、殆んど常にずっと速い計算になる。
の除算機能では、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に加えるべき回転の符号を決定する。
形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の直角位相
データ・ワドは夫々次の様に表わされる。
及び第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)のデータで
ある。
φ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に略
ゼロの大きさの剰余を発生する。
すると、第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データのスカラー乗算を行なうことが出来る。
第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図の構造は、他のどの乗算器の構成よりも、複素数乗
算に対してずっと融通性のある方式を表わす。
器に比べて多数の利点がある。第1に、最も重要なこと
であるが、位相と振幅の精度を切離すことが出来る。振
幅の切捨てが最終的な答の位相の精度に影響しないし、
逆も真である。第2に、こう云う実施例の複素数乗算器
10’ (1対のC0RDIG装置11−1及び11
−2を用いるか又は1個の装置1I″を用いる)は、中
間段階の結果を生ずるが、これは用途によっては役に立
つことがある。例えば、これらの構造に対する一方の入
力が単に他方の入力の共役複素数である場合、こう云う
装置は同時に出力10′又は10′eに入力信号の振幅
を出し、出力10′ f又は10′fに入力信号のエネ
ルギを出す。従って、この実施例は大きさ及びエネルギ
の同時の検出器として使うことが出来る。第3図及び第
4図の構造は、他のどの乗算器の構成よりも、複素数乗
算に対してずっと融通性のある方式を表わす。
この発明の新規な複素数CORDIC乗算器の現在好ま
しいと考えられる幾つかの形式並びにその中で用いられ
る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:加算器
器の簡略ブロック図、第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の複素数の実数部分及び虚数部分の夫々一方を
、+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記載の装置。
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)
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)
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)
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 |
-
1988
- 1988-05-31 US US07/200,491 patent/US4896287A/en not_active Expired - Fee Related
-
1989
- 1989-03-23 CN CN89101615A patent/CN1017283B/zh not_active Expired
- 1989-05-10 FR FR8906124A patent/FR2632088B1/fr not_active Expired - Lifetime
- 1989-05-25 DE DE3917059A patent/DE3917059A1/de not_active Withdrawn
- 1989-05-30 KR KR1019890007529A patent/KR0146334B1/ko not_active IP Right Cessation
- 1989-05-30 GB GB8912385A patent/GB2220090B/en not_active Expired - Fee Related
- 1989-05-31 JP JP13629589A patent/JP3283504B2/ja not_active Expired - Fee Related
Cited By (2)
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 |