このページの翻訳は最新ではありません。ここをクリックして、英語の最新版を参照してください。

Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor

Q-less QR 分解を使用した無限の行数をもつ複素数値の行列のために、A'AX = B で X の値を計算する

• ライブラリ:

• Fixed-Point Designer HDL Support / Matrices and Linear Algebra / Linear System Solvers

• 

×

説明

Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor ブロックは Q-less QR 分解を使用して線形方程式系 A'AX = B を解きます。ここで A と B は複素数値の行列です。A は、ストリーミング データを表す無限大の tall 行列です。

正則化パラメーターが非ゼロの場合、Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor は A の行を因数分解する前に最初の上三角因子 R を λI_n に初期化します。ここで、λ は正則化パラメーター、I_n = eye(n) です。

端子

入力

A(i,:) — 行列 A の行

ベクトル

行列 A の行。ベクトルとして指定します。A は m 行 n 列の行列で、m ≥ 2 および m ≥ n です。B が single または double の場合、A は B と同じデータ型でなければなりません。A が固定小数点データ型の場合、A は符号付きで、2 進小数点スケーリングを使用し、B と同じ語長をもたなければなりません。勾配とバイアス表現は固定小数点データ型ではサポートされていません。

データ型: single | double | fixed point
複素数のサポート: あり

B — 行列 B

行列 | ベクトル

行列 B。ベクトルまたは行列として指定します。B は m 行 p 列の行列で、m ≥ 2 です。A が single または double の場合、B は A と同じデータ型でなければなりません。B が固定小数点データ型の場合、B は符号付きで、2 進小数点スケーリングを使用し、A と同じ語長をもたなければなりません。勾配とバイアス表現は固定小数点データ型ではサポートされていません。

データ型: single | double | fixed point

validInA — A 入力が有効であるかどうか

Boolean スカラー

A(i, ;) 入力が有効であるかどうか。boolean スカラーとして指定します。この制御信号は A(i,:) 入力端子からのデータが有効であるかどうかを示します。この値が 1 (true) で readyA の値が 1 (true) の場合、ブロックは A(i,:) 入力端子の値を取得します。この値が 0 (false) の場合、ブロックは入力サンプルを無視します。

true の validInA 信号を送信してから、readyA が false に設定されるまでに、多少の遅延が生じることがあります。すべてのデータが確実に処理されるように、別の true の validInA 信号を送信するときは、readyA が false に設定されるまで待たなければなりません。

データ型: Boolean

validInB — 入力 B が有効であるかどうか

Boolean スカラー

入力 B が有効であるかどうか。boolean スカラーとして指定します。この制御信号は B 入力端子からのデータが有効であるかどうかを示します。この値が 1 (true) で readyB の値が 1 (true) の場合、ブロックは B 入力端子の値を取得します。この値が 0 (false) の場合、ブロックは入力サンプルを無視します。

true の validInB 信号を送信してから、readyB が false に設定されるまでに、多少の遅延が生じることがあります。すべてのデータが確実に処理されるように、別の true の validInB 信号を送信するときは、readyB が false に設定されるまで待たなければなりません。

データ型: Boolean

restart — 内部状態をクリアするかどうか

Boolean スカラー

内部状態をクリアするかどうか。boolean スカラーとして指定します。この値が 1 (true) の場合、ブロックは現在の計算を停止し、すべての内部状態をオフにします。この値が 0 (false) で validInA と validInB の値が 1 (true) の場合、ブロックは新しいサブフレームを開始します。

データ型: Boolean

出力

X — 行列 X

行列 | ベクトル

行列 X。行列またはベクトルとして返されます。

データ型: single | double | fixed point

validOut — 出力データが有効であるかどうか

Boolean スカラー

出力データが有効であるかどうか。boolean スカラーとして返されます。この制御信号は出力端子 X のデータが有効であるかどうかを示します。この値が 1 (true) の場合、ブロックは X の行を正常に計算しています。この値が 0 (false) の場合、出力データは有効ではありません。

データ型: Boolean

readyA — 入力 A に対するブロックの準備が整っているかどうか

Boolean スカラー

入力 A に対するブロックの準備が整っているかどうか。boolean スカラーとして返されます。この制御信号は新しい入力データに対するブロックの準備が整っているかどうかを示します。この値が 1 (true) で validInA の値が 1 (true) の場合、ブロックは入力データを次のタイム ステップで受け入れます。この値が 0 (false) の場合、ブロックは入力データを次のタイム ステップで無視します。

true の validInA 信号を送信してから、readyA が false に設定されるまでに、多少の遅延が生じることがあります。すべてのデータが確実に処理されるように、別の true の validInA 信号を送信するときは、readyA が false に設定されるまで待たなければなりません。

データ型: Boolean

readyB — 入力 B に対するブロックの準備が整っているかどうか

Boolean スカラー

入力 B に対するブロックの準備が整っているかどうか。boolean スカラーとして返されます。この制御信号は新しい入力データに対するブロックの準備が整っているかどうかを示します。この値が 1 (true) で validInB の値が 1 (true) の場合、ブロックは入力データを次のタイム ステップで受け入れます。この値が 0 (false) の場合、ブロックは入力データを次のタイム ステップで無視します。

true の validInB 信号を送信してから、readyB が false に設定されるまでに、多少の遅延が生じることがあります。すべてのデータが確実に処理されるように、別の true の validInB 信号を送信するときは、readyB が false に設定されるまで待たなければなりません。

データ型: Boolean

パラメーター

行列 A の列数と行列 B の行数 — 行列 A の列数と行列 B の行数

4 (既定値) | 正の整数値のスカラー

行列 A の列数と行列 B の行数。正の整数値のスカラーとして指定します。

プログラムでの使用

ブロック パラメーター: n

型: 文字ベクトル

値: 正の整数値のスカラー

既定の設定: 4

行列 B の列数 — 行列 B の列数

1 (既定値) | 正の整数値のスカラー

行列 B の列数。正の整数値のスカラーとして指定します。

プログラムでの使用

ブロック パラメーター: p

型: 文字ベクトル

値: 正の整数値のスカラー

既定の設定: 1

忘却係数 — 行列の各行を因数分解した後に適用される忘却係数

0.99 (既定値) | 実数の正のスカラー

行列の各行を因数分解した後に適用される忘却係数。実数の正のスカラーとして指定します。A の各行に無限に入力されるたびに、出力が更新されます。

プログラムでの使用

ブロック パラメーター: forgettingFactor

型: 文字ベクトル

値: 正の整数値のスカラー

既定の設定: 0.99

正則化パラメーター — 正則化パラメーター

0 (既定値) | 非負の実数スカラー

正則化パラメーター。非負のスカラーとして指定します。正則化パラメーターを小さい正の値にすると、問題の調整を向上させ、推定値の分散を小さくすることができます。バイアスがありますが、推定値の分散が小さいと、多くの場合、最小二乗推定値と比べて平均二乗誤差が小さくなります。

プログラムでの使用

ブロック パラメーター: regularizationParameter

型: 文字ベクトル

値: 非負の実数スカラー

既定の設定: 0

出力データ型 — 出力行列 X のデータ型

fixdt(1,18,14) (既定値) | double | single | fixdt(1,16,0) | <data type expression>

出力行列 X のデータ型。fixdt(1,18,14)、double、single、fixdt(1,16,0)、またはユーザー指定のデータ型式として指定します。型は直接指定されるか、Simulink.NumericType などのデータ型オブジェクトとして表現されます。

プログラムでの使用

ブロック パラメーター: OutputType

型: 文字ベクトル

値: 'fixdt(1,18,14)' | 'double' | 'single' | 'fixdt(1,16,0)' | '<data type expression>'

既定の設定: 'fixdt(1,18,14)'

モデルの例

Implement Hardware-Efficient Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor

How to use the Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor block.

スクリプトを開く

Algorithms to Determine Fixed-Point Types for Complex Q-less QR Matrix Solve A'AX=B

Derivation of algorithms for determining fixed-point types for complex Q-less QR matrix solve.

ライブ スクリプトを開く

Determine Fixed-Point Types for Complex Q-less QR Matrix Solve A'AX=B

Use fixed.complexQlessQRFixedpointTypes to determine fixed-point types for computation of the complex least-squares matrix equation.

ライブ スクリプトを開く

固定小数点の HDL に最適化された最小分散無歪応答 (MVDR) ビームフォーマー

ハードウェア効率に優れた MVDR ビームフォーマーを実装します。

スクリプトを開く

Compute Forgetting Factor Required for Streaming Input Data

Use fixed.forgettingFactor and fixed.forgettingFactorInverse to compute forgetting factor.

ライブ スクリプトを開く

アルゴリズム

忘却係数を使用する Q-less QR 分解

Complex Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor ブロックは、次の再帰を実装し、忘却係数 α を使用して連続的にストリーミングする n 行 1 列の行ベクトル A(k,:) の上三角因子 R を計算します。これは行列 A が無限大の tall である場合と同様です。範囲 0 < α < 1 の忘却係数を使用すると、無制限に積分されなくなります。

$$\begin{array}{c}{R}_0=\mathrm{zeros}(n,n)\\ \left[\sim ,R_1\right]=\mathrm{qr}\left(\left[\begin{array}{c}{R}_0\\ A\left(1,:\right)\end{array}\right],0\right)\\ R_1=\alpha R_1\\ \left[\sim ,R_2\right]=\mathrm{qr}\left(\left[\begin{array}{c}{R}_1\\ A\left(2,:\right)\end{array}\right],0\right)\\ R_2=\alpha R_2\\ ⋮\\ \left[\sim ,R_k\right]=\mathrm{qr}\left(\left[\left[\begin{array}{c}{R}_{k-1}\\ A\left(k,:\right)\end{array}\right]\right],0\right)\\ R_k=\alpha R_k\\ ⋮\end{array}$$

忘却係数と Tikhonov 正則化を使用する Q-less QR 分解

k 番目の入力 A(k,:) を処理した後の出力 X_k は、次の反復を使用して計算されます。

$$\begin{array}{c}{R}_0=\lambda I_n\\ \left[~,R_1\right]=\mathrm{qr}\left(\left[\begin{array}{c}{R}_0\\ A\left(1,:\right)\end{array}\right],0\right)\\ R_1=\alpha R_1\\ X_1=R_1\backslash \left(R{\text{'}}_1\backslash B\right)\\ \left[~,R_2\right]=\mathrm{qr}\left(\left[\begin{array}{c}{R}_1\\ A\left(2,:\right)\end{array}\right],0\right)\\ R_2=\alpha R_2\\ X_2=R_2\backslash \left(R{\text{'}}_2\backslash B\right)\\ ⋮\\ \left[~,R_k\right]=\mathrm{qr}\left(\left[\begin{array}{c}{R}_{k-1}\\ A\left(k,:\right)\end{array}\right],0\right)\\ R_k=\alpha R_k\\ X_k=R_k\backslash \left(R{\text{'}}_k\backslash B\right)\\ ⋮\end{array}$$

これは A'_kA_kX = B を計算することと数学的に等価です。ここで、A_k は次のように定義されますが、実際にはこのブロックで A_k は作成されません。

$$A_k=\left[\begin{array}{c}{\alpha }^k\lambda I_n\\ \left[\begin{array}{cccc}{\alpha }^k& & & \\ & {\alpha }^{k-1}& & \\ & & \ddots & \\ & & & \alpha \end{array}\right]A\left(1:k,:\right)\end{array}\right]$$

前進代入と後退代入

上三角因子の準備が整うと、前進代入と後退代入が現在の入力 B で計算され、出力 X が生成されます。

$$X=R_k\backslash \left(R_k^{\text{'}}\backslash B\right)$$

実装方法の選択

部分的なシストリック実装はスペースの制約よりも計算速度を優先し、バーストによる実装は演算速度を犠牲にしてスペースの制約を優先します。次の表は、行列の分解と線形方程式系の求解に使用できる実装間のトレードオフを示しています。

実装	Ready	レイテンシ	領域	サンプル ブロックまたは例
シストリック	C	O(n)	O(mn2)	シストリック アレイでの CORDIC を使用したハードウェア効率に優れた QR 分解の実装
部分的なシストリック	C	O(m)	O(n2)	Real Partial-Systolic QR Decomposition Real Partial-Systolic Matrix Solve Using QR Decomposition
忘却係数を使用する部分的なシストリック	C	O(n)	O(n2)	固定小数点の HDL に最適化された最小分散無歪応答 (MVDR) ビームフォーマー
バースト	O(n)	O(mn2)	O(n)	Real Burst QR Decomposition Real Burst Matrix Solve Using QR Decomposition

ここで、C はデータの語長に比例する定数、m は行列 A の行数、n は行列 A の列数です。

拡張機能

C/C++ コード生成

Simulink® Coder™ を使用して C および C++ コードを生成します。

勾配とバイアス表現は固定小数点データ型ではサポートされていません。

HDL コード生成

HDL Coder™ を使用して FPGA 設計および ASIC 設計のための Verilog および VHDL のコードを生成します。

HDL Coder™ には、HDL の実装および合成されたロジックに影響する追加のコンフィギュレーション オプションがあります。

HDL アーキテクチャ

このブロックには単一の既定の HDL アーキテクチャがあります。

HDL ブロック プロパティ

一般
ConstrainedOutputPipeline	既存の遅延を設計内で移動することによって出力に配置するレジスタの数。分散型パイプラインではこれらのレジスタは再分散されません。既定の設定は 0 です。詳細については、ConstrainedOutputPipeline (HDL Coder)を参照してください。
InputPipeline	生成されたコードに挿入する入力パイプライン ステージ数。分散型パイプラインと制約付き出力パイプラインでは、これらのレジスタを移動できます。既定の設定は 0 です。詳細については、InputPipeline (HDL Coder)を参照してください。
OutputPipeline	生成されたコードに挿入する出力パイプライン ステージ数。分散型パイプラインと制約付き出力パイプラインでは、これらのレジスタを移動できます。既定の設定は 0 です。詳細については、OutputPipeline (HDL Coder)を参照してください。

制限

固定小数点データ型のみをサポートします。

固定小数点の変換
Fixed-Point Designer™ を使用して固定小数点システムの設計とシミュレーションを行います。

バージョン履歴

R2020b で導入

参考

ブロック

• Real Partial-Systolic Matrix Solve Using Q-less QR Decomposition with Forgetting Factor | Real Partial-Systolic Q-less QR Decomposition with Forgetting Factor | Complex Partial-Systolic Q-less QR Decomposition | Complex Burst Q-less QR Decomposition

関数

• fixed.qlessQRMatrixSolve