一定周期で連続的に入力される複素相関行列をそれぞれ上三角行列に変換する信号処理回路であって、少なくともN個の複素相関行列を記憶する記憶手段と、記憶された複素相関行列の行列要素を読み出す読出手段と、CORDICアルゴリズムをパイプライン型回路方式で実装したCORDIC計算回路手段と、前記読出手段と前記CORDIC計算回路手段を制御して、1つの複素相関行列に対して前記CORDIC計算回路手段を繰り返し使用することで上三角行列を計算するQR分解制御手段であって、N個の複素相関行列の上三角行列への変換をインターリーブ形式で並行に計算するQR分解制御手段と、を備える、信号処理回路。前記QR分解制御手段は、一つの複素相関行列を入力してから上三角行列を出力するまでの計算処理をN個の処理ステップに分割し、N個の複素相関行列に対して、同じ時点でそれぞれ異なる処理ステップをインターリーブ形式で並行計算する、請求項1に記載の信号処理回路。前記QR分解制御手段は、前記読出手段及び前記CORDIC計算回路手段を制御するためのマイクロコードを記憶したプログラムメモリと、プログラムカウンタに従って読み出されたマイクロコードを記憶するマイクロコードレジスタと、を有し、マイクロコードレジスタに読み出されたマイクロコードに従って、前記読出手段及び前記CORDIC計算回路手段を制御する、請求項1又は2に記載の信号処理回路。前記CORDIC計算回路手段は、ベクタリング計算処理と回転計算処理の一方または双方を実行するパイプライン型計算回路を含む、請求項1乃至3のいずれか1項に記載の信号処理回路。被検体から放出された超音波を受信し受信信号に変換する複数の超音波受信素子と、前記複数の超音波受信素子から取得される複数の受信信号の位相を揃える整相遅延回路と、前記整相遅延回路から取得される位相が揃えられた受信信号に対してアダプティブ信号処理を施すアダプティブ信号処理手段と、を備え、前記アダプティブ信号処理手段が、位相が揃えられた受信信号をそれぞれ複素信号化し、複数の複素信号を取得する複素信号化回路と、前記複数の複素信号の複素相関行列を計算する相関行列計算回路と、前記複素相関行列を上三角行列に変換する請求項1乃至4のいずれか1項に記載の信号処理回路と、を含む、超音波診断装置。上記の問題点を考慮し、本発明は、計算量が特に多いQR分解を高速に計算可能な信号処理回路を提供し、アダプティブ信号処理を使用した高性能な超音波診断装置を実現することを目的とする。上記の課題を解決するために、本発明に係る信号処理回路は、一定周期で連続的に入力される複素相関行列をそれぞれ上三角行列に変換する信号処理回路であって、少なくともN個の複素相関行列を記憶する記憶手段と、記憶された複素相関行列の行列要素を読み出す読出手段と、CORDICアルゴリズムをパイプライン型回路方式で実装したCORDIC計算回路手段と、前記読出手段と前記CORDIC計算回路手段を制御して、1つの複素相関行列に対して前記CORDIC計算回路手段を繰り返し使用することで上三角行列を計算するQR分解制御手段であって、N個の複素相関行列の上三角行列への変換をインターリーブ形式で並行に計算するQR分解制御手段と、を備える。また、本発明に係る超音波診断装置は、被検体から放出された超音波を受信し受信信号に変換する複数の超音波受信素子と、前記複数の超音波受信素子から取得される複数の受信信号の位相を揃える整相遅延回路と、前記整相遅延手段から取得される位相が揃えられた受信信号に対してアダプティブ信号処理を施すアダプティブ信号処理手段と、を備え、前記アダプティブ信号処理手段が、位相が揃えられた受信信号をそれぞれ複素信号化し、複数の複素信号を取得する複素信号化回路と、前記複数の複素信号の複素相関行列を計算する相関行列計算回路と、前記複素相関行列を上三角行列に変換する請求項1乃至5のいずれか1項に記載の信号処理回路と、を含む。本発明によれば、QR分解を高速に計算可能な小型かつ低コストの信号処理回路を提供できる。また、アダプティブ信号処理を適用した高性能な超音波診断装置を提供できる。アダプティブ信号処理は、超音波診断装置の分野でも性能向上のための技術として適用することができる。その目的は、複数の受信素子で受信したエコー信号から雑音波の影響を排除し、フォーカス点からの反射波信号のみをSN良く受信し、エコー画像の分解能とコントラストの向上を図ることである。しかし、超音波エコー信号に重畳する雑音波は、被検体内部の反射物質の分布によって、方向と強度が極めて短時間に急速に変化してしまうため、時間方向の分解能を劣化させないために、数十nsから数百nsという極めて短い周期でアダプティブ信号処理固有の計算を繰り返す必要がある。更に複素相関行列Aの大きさも、少なくとも6×6~10×10程度が望ましく、アダプティブ信号処理を超音波診断装置に適用するためには、前述の非特許文献1と比較して実質数十倍から数百倍の計算速度を実現しなければならないという問題がある。ここで、計算速度を高速化するには、アダプティブ信号処理の大部分を占めるQR分解を高速に行うことが肝要である。QR分解の方法は種々知られているが、本発明ではギブンス回転法を用いる。ギブンス回転法は、相関行列Aを上三角化行列に変換するにあたって、2次元ベクトルの回転計算を繰り返しながら、相関行列の下三角領域の行列要素を一つずつ順番に0要素へと変換していく方法である。また、このギブンス回転法の基本計算である2次元ベクトル回転アルゴリズムには、良く知られたCORDICアルゴリズムを使用する。CORDICアルゴリズムは、予め定められた回転角の列に従って2次元ベクトルを右または左に回転させながら所望回転角に収束させることを基本原理とするアルゴリズムである。このCORDICアルゴリズムは、三角関数や平方根などの複雑な計算を必要とせずに、回転計算を整数型の加減算のみの簡単な処理で実現できる優れた特徴がある。しかしその反面、逐次計算を行うため原理的に結果が出力されるまでの時間が長いという欠点がある。QR分解処理は、ギブンス回転法とこのCORDICアルゴリズムを組み合わせた計算回路で容易に実現できる。しかし、ギブンス回転法によって行列要素を一つずつ0要素へと変換していく過程においては、一つの計算の結果を次の計算に使用するという計算過程が多く存在するため、一つの計算の結果が出るまでの時間が長いCORDICアルゴリズムは、従来高速計算のためには使用しづらかった。このため、非特許文献1においてはギブンス回転法を四則演算のみで実行可能なSGR法に修正し、浮動小数点型の計算回路で直接計算して高速化を実現している。このようにすれば、相関行列1個当たりのQR分解計算は大幅に高速化できる。しかし、浮動小数点型の計算回路は複雑で大きな回路になることから、超音波用途に適合できるまでの大幅な高速化と小型化は難しかった。本発明は、簡単な計算回路で実現可能なCORDICアルゴリズムの特徴を活かしながら高速回路を実現するために、1CORDICアルゴリズムの逐次処理をパイプライン型回路で実装し、2連続して入力される複数の複素相関行列のQR分解をインターリーブ形式で並行計算する、ようにする。ここで、インターリーブ形式での並行計算とは、一般に複数の計算ステップからなる複雑な複合計算を一つのプロセッサで複数回計算する時、複数の複合計算の各計算ステップが相互に挟まりあう形でプロセッサをタイムシェアしながら計算されるデータ処理形式のことである。パイプライン形式の回路は結果が出るまでの時間は長くなるが、1クロックごとに異なる回転計算を次々と入力して実行できるので、1個の回転計算当たりの所要時間は最高1クロックまで短縮が可能である。従って、複数の複素相関行列に対する計算のうち計算可能な回転計算を優先し、インターリーブ形式で相互に挟まりあう形で実行させることにより、単位時間あたりの計算量を大幅に向上させることが出来る。この結果、計算結果が出力されるまでの時間は長いものの、入力行列1個あたりの計算時間は大幅に短縮することが可能になる。<超音波診断装置の全体構成>図1はCMP法を用いた本発明の実施形態に係る超音波診断装置のブロック図である。超音波探触子1は、被検体に接触させるように用いられ、被検体に向けて超音波ビームを送信および受信する。超音波探触子1は、印加される駆動信号に基づいて超音波ビームを送信すると共に、被検体から放出反射され伝播してきた超音波エコーを受信してアナログ信号である受信信号に変換し出力する複数の超音波送受信素子5で構成される。これらの超音波送受信素子は1次元または2次元状に配列されてトランスデューサアレイ素子アレイになっている。なお、超音波探触子1は超音波送信素子と超音波受信素子とから構成されるようにしても良い。送信信号処理回路2はCPUの指示により送信信号3を発生し、スイッチ回路4を介して超音波送受信素子群5を駆動する。これにより、超音波探触子1からパルス状の超音波ビーム6が送信される。超音波ビーム6の反射によって発生した超音波エコー波7は、再び超音波送受信素子群5によって電気信号に変換され、スイッチ回路4を通って整相遅延回路8に送られる。整相遅延回路8は、受信フォーカス処理、すなわち、同一点Pからのエコー信号の到着時刻が揃うように各受信信号の遅延時間を調整する回路である。具体的には、複数チャンネルの受信超音波信号に対して、所望のフォーカス遅延を与えた後に加算処理を行う。この整相遅延回路8までの構成は、アダプティブ信号処理を行わない通常の超音波診断装置と共通であるが、これ以降の処理がアダプティブ信号処理CMP法固有の処理となっている。複素信号化回路9は各受信信号にそれぞれ90°位相のずれた信号を虚数信号として付加し、複数の複素信号を取得する回路である。複素相関行列計算回路10は複素信号化された受信信号間の複素相関行列Aを一定時間周期に生成し、記憶回路11に記憶させる回路である。例えば、送信信号処理回路2からの送信信号3が超音波送受信素子群5を介して発信されると、発信後約100μsの間、そのエコー信号が連続的に受信され、50MHzクロックでサンプリングされたディジタル信号が整相遅延回路8に転送される。整相遅延回路によって遅延時間の調整された信号は、複素信号化回路9と複素相関行列計算回路10によって10クロックごと5MHzに複素相関行列に変換され、変換された相関行列の系列A1A2A3が記憶回路11に一時記憶される。複素相関行列計算回路10では、実際には入力信号の相関行列に対して時間平均化処理、空間平均化処理などの付加的な処理を行って出力行列を生成するが、これらの詳細は本発明と直接関係しないので省略する。QR分解回路12、後退代入回路13、内積計算回路14はアダプティブ信号処理の特徴的な処理回路であり、記憶回路11に次々と記憶される相関行列Aiのそれぞれに対して、既知の拘束ベクトルCを用いた計算項を計算する。ここで、拘束ベクトルCは所望の超音波信号の受信方向を指定する方向ベクトルであり、整相遅延処理の行われた入力信号に対しては、通常全ての要素を1とする縦ベクトルである。また、右肩の記号Hは複素共役転置行列であることを意味する。計算項の計算は、通常計算項を計算し、その結果と拘束ベクトルCとの内積計算を行うようにする。この計算項は逆行列を含むため、アダプティブ信号処理の計算量の多さの主な原因となっており、この計算を高速に実行する計算手段の実現が本実施形態の主眼となっている。LOG変換回路15は、その計算結果を用いてCMP法の拘束付最小電力値pをのように計算し、エコー画像の1画素信号としてCPUに転送する。すなわち、前述のように10クロックごとに複素相関行列Aが生成されると、10クロックに1個の画素値が計算される。CPUは転送されたエコー画像信号を編集して表示装置16に、エコー画像として表示する。図1の構成により、アダプティブ信号処理を用いた高性能超音波診断装置を実現することができる。しかし、この構成の装置でリアルタイム画像表示を行うためには、前述のQR分解回路12、後退代入回路13、内積計算回路14の三つの回路が、次々と生成される相関行列の繰り返し時間周期内にQR分解処理の計算を完了する必要がある。この三つの回路のうち特に計算量の多いのはQR分解回路12である。相関行列の繰り返し時間周期は、前述の実施例では50MHzクロックの10クロック分の時間に相当する200nsである。つまり、非特許文献1に開示された81μsに対しても数十倍の高速計算が必要である。本実施形態では、アダプティブ信号処理のうち特に計算量の多いQR分解処理を行うQR分解回路12を工夫して、従来例に対して数十倍から数百倍高速に計算可能とする。また、QR分解回路12を小型で低コストに作成可能とする。このような小型かつ低コストな高性能なQR分解回路12により、アダプティブ信号処理を適用した高性能な超音波診断装置の実現ができる。なお、上記の例では、アダプティブ信号処理としてCMP法を採用した構成を説明したが、MMSE法、MSN法などのその他の手法を採用する場合でも、やはり高速な計算処理が必要となる点は変わらない。そして、いずれの手法においても、QR分解計算が特に計算量の多い部分であるので、本実施形態におけるQR分解回路12を使って、MMSE法やMSE法などのアダプティブ信号処理を適用した高性能な超音波診断装置を実現できる。ギブンス回転アルゴリズム本実施形態におけるQR分解回路12の内容を具体的に説明するために、以下、ギブンス回転アルゴリズムによるQR分解処理について更に詳細に説明する。ギブンス回転法は、前述のごとく2次元ベクトルの回転処理を繰り返し実施して、行列要素をひとつずつ0化する計算方法である。すなわち、まず式に示すように相関行列Aの左から回転行列Qを掛け相関行列Aの第2行1列要素を0化する。このとき、式の回転行列Qによる乗算は、相関行列Aの第1行と第2行の縦方向に配列した二つの要素から構成されるn個の二次元ベクトルの回転計算を実行することに相当している。次に、式のように左から回転行列Qを掛け第3行1列要素を0化する。式の回転行列Qによる乗算は、相関行列Aの第1行と第3行の縦方向に配列した二つの要素から構成されるn個の二次元ベクトルの回転計算を実行することに相当している。同様の処理を第n行1列要素まで繰り返すことにより、式のように第1列要素の上三角行列化を完了することができる。第1列要素の上三角行列化が完了したら、式のように回転行列Qを掛け第2行2列要素を0化する。このような処理を次々と繰り返して、第2列要素の上三角化が完了したら、第3列、第4列を順次上三角化し、最終的に入力行列の上三角化を完了する。以上のようにギブンス回転によるQR分解は0化する要素の数だけ行列計算を繰り返しているが、その行列計算は全て2次元ベクトルの回転を基本としている。そこで、基本的な計算として式に示す2×2行列の上三角化計算を例として、ギブンス回転の計算手順を更に詳しく説明する。2×2行列の場合、ギブンス回転処理は、まず相関行列Aの第2行1列要素を0とするような回転行列Qの各係数を求め、次にその回転行列Qによる行列乗算計算を行う。もし相関行列が実数行列であれば、回転行列Qは、第1列要素からなるベクトルの絶対値と偏角を式のように計算し、式のように求められる。また、相関行列Aと回転行列Qとの乗算結果は、相関行列Aの第2列要素からなるベクトルの回転計算結果を式のようにおくと、式のようになる。ここでは、2次元ベクトルの絶対値と偏角を求める式の処理をVEC処理、式の回転計算をROT処理と称する。式から式に示した任意のn×n実数行列に対するギブンス回転法によるQR分解は、このVEC処理とROT処理を多数回繰返すことにより実行できる。相関行列Aが複素数の場合は、複素数平面内での回転処理が必要になるため、ギブンス回転法の計算手順はやや複雑になる。しかし、複素数の複素数平面内での回転計算や、複素数の絶対値と偏角を求める計算も2次元ベクトルの計算とまったく同様に計算できることから、詳細は省略するが、相関行列が複素数の場合であっても、ギブンス回転法はVEC処理とROT処理を多数回繰返すことで実行できることが知られている。アダプティブ信号処理の場合、使用される相関行列は通常複素数であるが、複素数であっても実数行列の場合と同様にVEC処理とROT処理の繰り返しで実行できるので、説明をわかりやすくするために、ここでは実数行列を仮定して説明することにする。CORDICアルゴリズムギブンス回転法のVEC処理、ROT処理は、良く知られたようにCORDICアルゴリズムを適用した回路で計算する方法が最も簡単である。そこで、次にCORDICアルゴリズムの概略を説明する。CORDICアルゴリズムは、回転行列をのように変形し、を満足する角度θjの回転計算を組み合わせて所望の回転計算を行う方法である。すなわち、トータルの回転角が所望の回転角に満たなければ正方向に角度θjだけ回転計算し、所望の回転角を超えていれば負方向に角度θjだけ回転する。図2は、2次元ベクトルV0を角度θだけ回転してベクトルV1とするROT処理の手順を、具体的に示したものである。まず、角度θが正なのでベクトルV0を角度θ1だけ正方向に回転する。次に回転角θ1がθより大きいので、回転結果を角度θ2だけ負方向に回転する。更に、トータルの回転角θ1-θ2が角度θより小さいので、回転結果を角度θ3だけ正方向に回転する。このようにベクトルV0の回転計算をj123の角度系列に従って実行すれば、トータルの回転角は最終的に2jのオーダーで所望の角度θに収束し、角度θだけ回転したベクトルV1が計算できる。式を満足する角度θjの回転計算は、tanθjの乗算が単にビットシフト処理となるので、cosθjの乗算を無視すれば、整数型の加算だけで実行できる。cosθjの乗算を省略すると、ベクトルの大きさにcosθjに相当する誤差が生じるが、この誤差は最後に纏めて乗算処理することによって容易に修正できるので問題は生じない。以下、簡単のために最後の大きさ補正のための乗算は無視して説明する。VEC処理は、ROT処理の変形で実現できる。図2は2次元ベクトルV0の偏角θと大きさRを計算するVEC処理の手順を、具体的に示したものである。VEC処理においては、ベクトルV0のy方向要素の符号を監視し、y要素が正ならベクトルを負方向に回転し、負なら正方向に回転することによって、ベクトルをx軸上に収れんさせる。ベクトルがx軸上に重なれば、その時のx要素の大きさがベクトルの大きさであり、トータルの回転角が偏角θに相当する。このように、CORDICアルゴリズムを用いれば、式、で示す複雑なVEC処理、ROT処理のどちらも基本的に簡単な加減算処理の繰り返しで計算できる。図3はCORDICアルゴリズムをパイプライン型回路方式で実装した時の実施例である。図において、PR0PR1PRmはm段からなるパイプライン回路を構成するパイプラインレジスタ、ROT1ROT2ROTmはそれぞれ式に示される固定回転角θ1θ2θmの回転計算を行う回転計算回路である。このパイプライン型計算回路においては、PR0に入力となるベクトルをセットすると、所定のパイプラインクロックに従って、順々に必要な回転計算がおこなわれて、mクロック後に回転結果が出力レジスタPRmから出力される。このとき、パイプライン段数mがある程度以上大きければ、回転角度の誤差は2mrad以下であることが保証されるので、パイプラインの段数mは必要な要求計算精度によって決定される。パイプライン型のCORDIC計算回路は、入力レジスタPR0に1クロックごとに異なるベクトルをセットすれば、それぞれの回転計算結果が所定のクロック後に次々と出力されるので、一つの回転計算当たりの計算時間が非常に短く実行できるという特徴がある。しかし、結果が出るまでの時間がパイプライン段数に相当するクロック数だけ遅れるという問題もある。本実施形態におけるQR分解回路図4は、以上に述べたギブンス回転法とCORDIC回路の特性に鑑みてなされたQR分解回路12の実施例である。図において、11は複数の相関行列を記憶する記憶回路、11a111a211a3は行列データ記憶領域、21は記憶された相関行列の所望の要素を読みだす読み出し制御回路、22は前述のパイプライン型のCORDIC回路、23はCODIC計算回路22の出力をいったん記憶して再びCORDIC計算回路22に入力することのできる中間記憶回路である。CORDIC計算回路22は相関行列の所望の要素と中間記憶回路23に記憶された中間データを入力とし、VEC処理またはROT処理を行い、再び結果を中間メモリ23に記憶する。QR分解制御回路25は、プログラムカウンタ26を用いてプログラムメモリ27からマイクロコードをマイクロコードレジスタ28に順次読みだし、読み出されたマイクロコードによって、1記憶回路11及び中間記憶回路23からのデータ読出し、2CORDIC計算回路22への入力設定、3計算結果の中間記憶回路23への書込み、など、QR分解計算手順を制御する。このとき、計算結果が最終計算結果であれば、計算結果をQR分解結果として後段の後退代入回路13に出力する。ギブンス回転法は前述のようにVEC処理とROT処理を繰り返すことで実行できるので、プログラムメモリ27にあらかじめQR分解のためのマイクロコード列を記憶しておけば、この構成により、QR分解を実行することができる。このようにQR分解はCORDIC計算回路22を繰り返し駆動することで実行できる。ただし、ギブンス回転法によるQR分解処理は一つのCORDIC計算の結果を次の計算に使用する場合が多いので、一つのCORDIC計算が完了するまで次の計算の待たされる状態が頻繁に発生し、計算速度を十分に高速化することができない。そこで本実施形態では、複数の行列に対するQR分解をインターリーブ形式で並行して計算させることにより、パイプライン型CORDIC計算回路22の稼働率を向上し、行列1個当たりの計算時間を向上させる。インターリーブ形式の並行計算は、QR分解制御回路25のプログラムメモリ27の内容を適切に設定することにより容易に実現可能である。図5は、本実施形態の並行計算の効果を説明した図である。今、単一行列に対するQR分解処理が二つのVEC処理Pv1Pv2と三つのROT処理Pr1Pr2Pr3からなっていて、1VEC処理、ROT処理の計算にはそれぞれ5クロックの時間が必要、2Pr1Pr2はPv1計算結果を、Pv2はPr2の結果をそれぞれ使用し、3Pr3はPv2の計算結果を使用する、と仮定する。図5はこの場合のQR分解処理のCORDIC計算回路22への入出力シーケンスを具体的に示したものである。すなわち、まず、VEC処理Pv1inが入力されると、5クロック後にその結果Pv1outが出力され、その結果を待ってROT処理Pr1in、Pr2inが入力される。5クロック後にROT処理Pr1、Pr2の結果Pr1outPr2outが出力されると、その結果を待ってVEC処理Pv2が入力される。このようにして、最後にROT処理Pr3の結果Pr3outが出力されて、単一行列に対するQR分解処理が完了する。この場合の計算時間は、図から推測できるように約25クロックである。一方、図5は6個の入力行列A1A2A6を3個ずつ、1回のVEC処理あるいはROT処理を単位としたインターリーブ形式で並行計算したときの動作シーケンスを示したものである。図では、1つの行列に対する計算時間を3つの時間区間に等分割し、1分割ずつ位相をずらせて並行計算するようにしている。図は各行列に対する計算のCORDIC計算回路入力のみを示しているが、この例ではCORDIC計算回路22への各入力はそれぞれ空き時間に挿入され、図のA1~A6合成結果のように合成される。この結果、CORDIC計算回路22の稼働率が向上し、3個の行列が1個当たり単一行列計算時の約13の9クロックで計算可能となる。図6は、並行計算する行列数Nを3とし、一定の時間周期で連続的に入力される行列A1A2A3をQR分解する動作シーケンスを、時間区間T1T2T3ごとに示した図である。ただし、図においてAiPjは行列Aiに対して区間処理Pjを行うことを意味し、区間処理Pjは一つの行列に対する全計算時間をN個の時間区間に等分割した時の第j番目の時間区間の計算処理を意味する。また、図の左の11a111a211a3は相関行列記憶回路11のN個の行列記憶領域を示す。すなわち、時間区間T1では、行列記憶領域11aに記憶された行列A1に対して区間処理P1を行う。時間区間T2では、行列記憶領域11aに記憶された行列A1に対する区間処理P2と、行列記憶領域11bに記憶された行列A2に対する区間処理P1とをインターリーブ形式で同時並行的に計算する。時間区間T3以降も同様であり、T3以降はN個の行列に対してそれぞれ異なる区間処理P1P2PNが割り当てられて、同時並行的に処理される。このようにN個の複素行列に対して、同じ時点でそれぞれ異なる処理ステップをインターリーブ形式で平行計算する。この結果、各時間区間でのトータルのCORDIC回路計算数は、各区間処理の計算数に大きなばらつきがあったとしても、開始および終了区間を除いて均一となり、並行計算が容易になるという利点がある。なお、同一時間区間でN個の区間処理P1P2PNを合成する時、複数の区間処理においてCORDIC回路入力を重複して実行するクロックがあり得るが、その場合には重複する区間処理を1クロックずつ遅延させて合成するようにする。このようにすれば、実際の各時間区間の幅は若干大きくなるものの区間処理の合成は可能であり、各時間区間の長さも均一にできる。超音波診断装置のように相関行列が一定周期で連続的に入力される場合には、入力周期を前記の時間区間に一致させることで、QR分解処理のリアルタイム処理が容易に可能になる。前述の実施例では、CORDIC回路の計算機能をVEC処理またはROT処理に切り替えて使用していたが、QR分解の総処理時間を高速化するために、CORDICの機能を複合して一つのCORDIC回路に合成することは技術的に容易である。例えば、CORDIC回路を、二つのベクトルに対する同一角度のROT処理を同時に実施するようにしたり、VEC処理とROT処理との双方を同時に実施するように構成することもできる。更に一つのROT処理の結果に対するROT処理を一つの回路で実行したり、その一方の回転処理を通常の乗算回路で実施したりすることも可能である。このような高速化や小型化のための変形をしたとしても、CORDICアルゴリズムをパイプライン回路方式で実現し、複数行列をインターリーブ方式で並行的に計算するようにすれば大幅な高速化は可能であり、本発明に含まれる。シミュレーション評価によれば、1個のQR分解回路で8個の行列を並行計算することにより、CORDICパイプライン回路の稼働率を約7倍に向上できる。その結果、本実施形態のQR分解回路12により、6×6複素行列に対するQR分解計算を120クロック以内に完了できる小型回路が実現可能である。回路の駆動クロックを200MHzとすれば1行列当たりの計算時間は600nsである。この計算時間は、計算時間が行列寸法の3乗に比例することを考慮すると、4×4複素行列を81μsで計算する非特許文献1の公知例に対して約40倍高速である。本実施形態のQR分解回路を3個並列に駆動するようにすれば、1行列当たりの計算時間は200nsとなり、前述の超音波診断装置に必要な5MHz周期の計算が可能となる。QR分解回路の中心となるCORDIC計算回路は、整数型の加減算回路を主体とする比較的小型の回路なので、3回路並列またはそれ以上の並列度の構成であっても容易に簡易LSIであるFPGA内部に実装が可能である。このように、本発明によればアダプティブ信号処理に必要な高速QR分解処理を小型低価格の回路で実現することができる。特に高速計算の必要な超音波診断装置へのアダプティブ信号処理の適用が容易になり、超音波診断装置の高性能化を実現する上で極めて効果が大きい。アダプティブ信号処理技術は、実施例で述べたCMP法以外にも、前述のように最小2乗誤差法、最大SNR法など多くの手法が知られている。これらの手法の多くは相関行列Aの逆行列に既知のベクトルyを掛けたなる計算を含み、この計算を実行するために実施例で述べたQR分解手段が必要となる。したがって、CMP法以外を採用する場合にも、上記のQR分解処理回路を利用することで、アダプティブ信号処理を高速に行える。また、超音波診断装置以外の無線通信、無線LANなどの装置であっても、アダプティブ信号処理の適用には同様な計算回路の実装が必要であり、本発明の信号処理回路がそれらの装置の小型高速化に有効であることは言うまでもない。