qr

シンボリック行列の QR 分解

構文

R = qr(A)

[Q,R] = qr(A)

[Q,R,P] = qr(A)

[___] = qr(A,"econ")

[Q,R,P] = qr(___,outputForm)

[___] = qr(A,0)

[C,R] = qr(A,B)

[C,R,P] = qr(A,B)

[___] = qr(A,B,"econ")

[C,R,P] = qr(___,outputForm)

[___] = qr(A,B,0)

___ = qr(___,"real")

説明

R = qr(A) は、QR 分解 A = Q*R の R 因子を返します。ここで、A は m 行 n 列のシンボリック行列、R は m 行 n 列の上三角シンボリック行列、Q は m 行 m 列のユニタリ シンボリック行列です。

[Q,R] = qr(A) は、A = Q*R となるように上三角行列 R とユニタリ行列 Q を返します。

[Q,R,P] = qr(A) は、A*P = Q*R となる置換行列 P も返します。A のすべての要素が浮動小数点数で近似できる場合、この構文は abs(diag(R)) が降順になるように置換行列 P を選択します。そうでない場合は P = eye(n) を返します。

[___] = qr(A,"econ") は、前述の出力引数の組み合わせのいずれかを使用してエコノミーサイズの分解を返します。出力のサイズは、m 行 n 列の行列 A のサイズに応じて、次のように異なります。

• m > n の場合、qr は Q の最初の n 列と R の最初の n 行のみを計算する。

• m ≤ n の場合、エコノミーサイズの分解は、"econ" オプションを使用しない場合の QR 分解と同じになる。

既定では、"econ" を使用すると、qr は置換情報を行列 P として返します。

[Q,R,P] = qr(___,outputForm) は、置換情報 P を行列またはベクトルのどちらで返すかを指定します。たとえば、outputForm が "vector" の場合、A(:,P) = Q*R になります。outputForm の既定値は、A*P = Q*R となる "matrix" です。

[___] = qr(A,0) は、qr(A,"econ","vector") と等価です。この構文は推奨されていません。代わりに "econ" オプションを使用してください。

[C,R] = qr(A,B) は、C = Q'*B および A = Q*R となるように上三角行列 R と行列 C を返します。ここでは、A および B は同じ行数でなければなりません。

C と R は、行列方程式 A*X = B の解を X = R\C として表現します。

[C,R,P] = qr(A,B) は、A*P = Q*R となる置換行列 P も返します。A のすべての要素が浮動小数点数で近似できる場合、この構文は abs(diag(R)) が降順になるように置換行列 P を選択します。そうでない場合は P = eye(n) を返します。

C、R、および P は、行列方程式 A*X = B の解を X = P*(R\C) として表します。

[___] = qr(A,B,"econ") は、前述の出力引数の組み合わせのいずれかを使用してエコノミーサイズの分解を返します。出力のサイズは、m 行 n 列の行列 A のサイズに応じて、次のように異なります。

• m > n の場合、qr は C と R の最初の n 行のみを計算する。

• m ≤ n の場合、エコノミーサイズの分解は、"econ" オプションを使用しない場合の QR 分解と同じになる。

既定では、"econ" を使用すると、qr は置換情報を行列 P として返します。

[C,R,P] = qr(___,outputForm) は、置換情報 P を行列またはベクトルのどちらで返すかを指定します。たとえば、outputForm が "vector" の場合、行列方程式 A*X = B の解は X(P,:) = R\C になります。outputForm の既定値は "matrix" で、A*X = B の解が X = P*(R\C) となります。

[___] = qr(A,B,0) は、qr(A,B,"econ","vector") と等価です。この構文は推奨されていません。代わりに "econ" オプションを使用してください。

___ = qr(___,"real") は、入力引数と中間結果が実数であることを仮定するため、abs および conj への呼び出しを無効にします。このオプションを使用すると、qr はすべてのシンボリック変数が実数を表すと仮定します。この構文では、すべての数値引数が実数でなければなりません。

"real" を使用して、結果に複素共役が含まれないようにします。

例

QR 分解の上三角因子

4 行 4 列のウィルキンソンの固有値テスト行列の QR 分解の上三角 R 因子を計算します。

4 行 4 列のウィルキンソンの固有値テスト行列を作成します。

A = sym(wilkinson(4))

A = 
$$\left(\begin{array}{cccc}\frac{3}{2}& 1& 0& 0\\ 1& \frac{1}{2}& 1& 0\\ 0& 1& \frac{1}{2}& 1\\ 0& 0& 1& \frac{3}{2}\end{array}\right)$$

Q 因子を返さずに、QR 分解の R 因子を返します。

R = qr(A)

R = 
$$\left(\begin{array}{cccc}\frac{\sqrt{13}}{2}& \frac{4\hspace{0.17em}\sqrt{13}}{13}& \frac{2\hspace{0.17em}\sqrt{13}}{13}& 0\\ 0& \frac{\sqrt{13}\hspace{0.17em}\sqrt{53}}{26}& \frac{10\hspace{0.17em}\sqrt{13}\hspace{0.17em}\sqrt{53}}{689}& \frac{2\hspace{0.17em}\sqrt{13}\hspace{0.17em}\sqrt{53}}{53}\\ 0& 0& \frac{\sqrt{53}\hspace{0.17em}\sqrt{381}}{106}& \frac{172\hspace{0.17em}\sqrt{53}\hspace{0.17em}\sqrt{381}}{20193}\\ 0& 0& 0& \frac{35\hspace{0.17em}\sqrt{381}}{762}\end{array}\right)$$

パスカル行列の QR 分解

3 行 3 列のパスカル行列の QR 分解を計算します。

3 行 3 列のパスカル行列を作成します。

A = sym(pascal(3))

A = 
$$\left(\begin{array}{ccc}1& 1& 1\\ 1& 2& 3\\ 1& 3& 6\end{array}\right)$$

A の QR 分解を表現する Q 行列と R 行列を計算します。

[Q,R] = qr(A)

Q = 
$$\left(\begin{array}{ccc}\frac{\sqrt{3}}{3}& -\frac{\sqrt{2}}{2}& \frac{\sqrt{6}}{6}\\ \frac{\sqrt{3}}{3}& 0& -\frac{\sqrt{6}}{3}\\ \frac{\sqrt{3}}{3}& \frac{\sqrt{2}}{2}& \frac{\sqrt{6}}{6}\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}\sqrt{3}& 2\hspace{0.17em}\sqrt{3}& \frac{10\hspace{0.17em}\sqrt{3}}{3}\\ 0& \sqrt{2}& \frac{5\hspace{0.17em}\sqrt{2}}{2}\\ 0& 0& \frac{\sqrt{6}}{6}\end{array}\right)$$

isAlways を使用して、A = Q*R であることを検証します。

TF = isAlways(A == Q*R)

TF = 3x3 logical array

   1   1   1
   1   1   1
   1   1   1

置換された QR 分解

置換を使用して、浮動小数点行列の QR 分解の数値安定性を向上させます。関数 qr は、置換情報を行列またはベクトルとして返します。

可変精度の演算を行うための有効小数桁数を 10 に設定します。3 行 3 列のシンボリック ヒルベルト行列を浮動小数点数で近似します。

previoussetting = digits(10);
A = vpa(hilb(3))

A = 
$$\left(\begin{array}{ccc}1.0& 0.5& 0.3333333333\\ 0.5& 0.3333333333& 0.25\\ 0.3333333333& 0.25& 0.2\end{array}\right)$$

まず、A の QR 分解を置換なしで計算します。

[Q,R] = qr(A)

Q = 
$$\left(\begin{array}{ccc}0.8571428571& -0.5016049166& 0.1170411472\\ 0.4285714286& 0.5684855721& -0.7022468832\\ 0.2857142857& 0.6520863915& 0.7022468832\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}1.166666667& 0.6428571429& 0.45\\ 0& 0.1017143303& 0.1053370325\\ 0& 0& 0.003901371573\end{array}\right)$$

A と Q*R の差を計算します。丸め誤差のため、計算された行列 Q と R は、厳密に等式 A = Q*R を満たしません。

E = A - Q*R

E = 
$$\left(\begin{array}{ccc}-3.469446952e-18& -6.938893904e-18& -4.33680869e-18\\ 0& -1.734723476e-18& -8.67361738e-19\\ 0& -8.67361738e-19& -4.33680869e-19\end{array}\right)$$

QR 分解の数値安定性を上げるには、3 つの出力引数を使用する構文を指定することで置換を使用します。シンボリック変数、式または関数を含まない行列に対し、次の構文ではピボットがトリガーされ、R に返される行列の abs(diag(R)) が降順になります。

[Q,R,P] = qr(A)

Q = 
$$\left(\begin{array}{ccc}0.8571428571& -0.4969293466& -0.1355261854\\ 0.4285714286& 0.5421047417& 0.7228063223\\ 0.2857142857& 0.6776309272& -0.6776309272\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}1.166666667& 0.45& 0.6428571429\\ 0& 0.1054092553& 0.1016446391\\ 0& 0& 0.003764616262\end{array}\right)$$

P = 3×3

     1     0     0
     0     0     1
     0     1     0

等式 A*P = Q*R を確認します。置換を使用した QR 分解では、丸め誤差が小さくなります。

E = A*P - Q*R

E = 
$$\left(\begin{array}{ccc}-3.469446952e-18& -4.33680869e-18& -6.938893904e-18\\ 0& -8.67361738e-19& -1.734723476e-18\\ 0& -4.33680869e-19& -1.734723476e-18\end{array}\right)$$

次に、"vector" 引数を使用して置換情報をベクトルとして返します。

[Q,R,p] = qr(A,"vector")

Q = 
$$\left(\begin{array}{ccc}0.8571428571& -0.4969293466& -0.1355261854\\ 0.4285714286& 0.5421047417& 0.7228063223\\ 0.2857142857& 0.6776309272& -0.6776309272\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}1.166666667& 0.45& 0.6428571429\\ 0& 0.1054092553& 0.1016446391\\ 0& 0& 0.003764616262\end{array}\right)$$

p = 1×3

     1     3     2

A(:,p) = Q*R であることを検証します。

E = A(:,p) - Q*R

E = 
$$\left(\begin{array}{ccc}-3.469446952e-18& -4.33680869e-18& -6.938893904e-18\\ 0& -8.67361738e-19& -1.734723476e-18\\ 0& -4.33680869e-19& -1.734723476e-18\end{array}\right)$$

厳密なシンボリック計算では、丸め誤差は発生しません。

A = sym(hilb(3));
[Q,R] = qr(A);
E = A - Q*R

E = 
$$\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right)$$

10 進数の有効桁数を以前の設定に戻します。

digits(previoussetting)

QR 分解を使用した行列方程式の求解

qr を使用して、行列形式の方程式系を解くことができます。

方程式系 A*X = b を解きます。ここで、A は 5 行 5 列のシンボリック行列で、b は 5 行 1 列のシンボリック ベクトルです。

A = sym(invhilb(5))

A = 
$$\left(\begin{array}{ccccc}25& -300& 1050& -1400& 630\\ -300& 4800& -18900& 26880& -12600\\ 1050& -18900& 79380& -117600& 56700\\ -1400& 26880& -117600& 179200& -88200\\ 630& -12600& 56700& -88200& 44100\end{array}\right)$$

b = sym([1:5]')

b = 
$$\left(\begin{array}{c}1\\ 2\\ 3\\ 4\\ 5\end{array}\right)$$

C = Q'*b かつ A = Q*R となる行列 C と行列 R を求めます。

[C,R] = qr(A,b);

解 X を計算します。

X = R\C

X = 
$$\left(\begin{array}{c}5\\ \frac{71}{20}\\ \frac{197}{70}\\ \frac{657}{280}\\ \frac{1271}{630}\end{array}\right)$$

isAlways を使用して、X が系 A*X = b の解であることを検証します。

tf = isAlways(A*X == b)

tf = 5x1 logical array

   1
   1
   1
   1
   1

置換された QR 分解を使用した行列方程式の求解

浮動小数点数をもつ方程式系を解く際に、置換行列または置換ベクトルをもつ QR 分解を使用します。

方程式系 A*X = b を解きます。ここで、A は 3 行 3 列の可変精度のシンボリック行列で、b は 3 行 1 列の可変精度のシンボリック ベクトルです。C = Q'*b かつ A = Q*R となる行列 C と行列 R を求めます。

previoussetting = digits(10);
A = vpa([2 -3 -1; 1 1 -1; 0 1 -1]);
b = vpa([2; 0; -1]);
[C,R,P] = qr(A,b)

C = 
$$\left(\begin{array}{c}-2.110579412\\ -0.2132007164\\ 0.7071067812\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}3.31662479& 0.3015113446& -1.507556723\\ 0& 1.705605731& -1.492405014\\ 0& 0& 0.7071067812\end{array}\right)$$

P = 3×3

     0     0     1
     1     0     0
     0     1     0

解 X を計算します。

X = P*(R\C)

X = 
$$\left(\begin{array}{c}1.0\\ -0.25\\ 0.75\end{array}\right)$$

または、置換情報をベクトルとして返します。

[C,R,p] = qr(A,b,"vector")

C = 
$$\left(\begin{array}{c}-2.110579412\\ -0.2132007164\\ 0.7071067812\end{array}\right)$$

R = 
$$\left(\begin{array}{ccc}3.31662479& 0.3015113446& -1.507556723\\ 0& 1.705605731& -1.492405014\\ 0& 0& 0.7071067812\end{array}\right)$$

p = 1×3

     2     3     1

置換ベクトル p を使用して解 X を計算します。

X(p,:) = R\C

X = 
$$\left(\begin{array}{c}1.0\\ -0.25\\ 0.75\end{array}\right)$$

10 進数の有効桁数を以前の設定に戻します。

digits(previoussetting)

エコノミーサイズの QR 分解

"econ" オプションを使用して、エコノミーサイズの QR 分解を計算します。

4 行 4 列のパスカル行列の最初の 2 列の行列を作成します。

A = sym(pascal(4));
A = A(:,1:2)

A = 
$$\left(\begin{array}{cc}1& 1\\ 1& 2\\ 1& 3\\ 1& 4\end{array}\right)$$

この行列の QR 分解を計算します。

[Q,R] = qr(A)

Q = 
$$\begin{array}{l}\left(\begin{array}{cccc}\frac{1}{2}& -{\sigma }_1& \frac{\sqrt{3}\hspace{0.17em}\sqrt{10}}{10}& 0\\ \frac{1}{2}& -\frac{\sqrt{5}}{10}& -\frac{2\hspace{0.17em}\sqrt{3}\hspace{0.17em}\sqrt{10}}{15}& \frac{\sqrt{6}}{6}\\ \frac{1}{2}& \frac{\sqrt{5}}{10}& -\frac{\sqrt{3}\hspace{0.17em}\sqrt{10}}{30}& -\frac{\sqrt{6}}{3}\\ \frac{1}{2}& {\sigma }_1& \frac{\sqrt{3}\hspace{0.17em}\sqrt{10}}{15}& \frac{\sqrt{6}}{6}\end{array}\right)\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\sigma }_1=\frac{3\hspace{0.17em}\sqrt{5}}{10}\end{array}$$

R = 
$$\left(\begin{array}{cc}2& 5\\ 0& \sqrt{5}\\ 0& 0\\ 0& 0\end{array}\right)$$

次に、この行列のエコノミーサイズの QR 分解を計算します。行数が A の列数より多いため、qr は Q の最初の 2 列と R の最初の 2 行のみを計算します。

[Q,R] = qr(A,"econ")

Q = 
$$\left(\begin{array}{cc}\frac{1}{2}& -\frac{3\hspace{0.17em}\sqrt{5}}{10}\\ \frac{1}{2}& -\frac{\sqrt{5}}{10}\\ \frac{1}{2}& \frac{\sqrt{5}}{10}\\ \frac{1}{2}& \frac{3\hspace{0.17em}\sqrt{5}}{10}\end{array}\right)$$

R = 
$$\left(\begin{array}{cc}2& 5\\ 0& \sqrt{5}\end{array}\right)$$

複素共役の回避

"real" オプションを使用すると、QR 分解の結果で複素共役を回避できます。

シンボリック変数を要素として含む行列を作成します。

syms x
A = [1 2; 3 x]

A = 
$$\left(\begin{array}{cc}1& 2\\ 3& x\end{array}\right)$$

この行列の QR 分解を計算します。既定の設定では、qr は x が複素数を表現すると仮定するので、結果には関数 abs を使用する式が含まれます。

[Q,R] = qr(A)

Q = 
$$\begin{array}{l}\left(\begin{array}{cc}\frac{\sqrt{10}}{10}& -\frac{\frac{3\hspace{0.17em}x}{10}-\frac{9}{5}}{{\sigma }_1}\\ \frac{3\hspace{0.17em}\sqrt{10}}{10}& \frac{\frac{x}{10}-\frac{3}{5}}{{\sigma }_1}\end{array}\right)\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\sigma }_1=\sqrt{{\left|\frac{x}{10}-\frac{3}{5}\right|}^2+{\left|\frac{3\hspace{0.17em}x}{10}-\frac{9}{5}\right|}^2}\end{array}$$

R = 
$$\left(\begin{array}{cc}\sqrt{10}& \frac{\sqrt{10}\hspace{0.17em}\left(3\hspace{0.17em}x+2\right)}{10}\\ 0& \sqrt{{\left|\frac{x}{10}-\frac{3}{5}\right|}^2+{\left|\frac{3\hspace{0.17em}x}{10}-\frac{9}{5}\right|}^2}\end{array}\right)$$

"real" オプションを使用すると、qr はすべてのシンボリック変数が実数を表現すると仮定するので、より短い結果を返すことができます。

[Q,R] = qr(A,"real")

Q = 
$$\left(\begin{array}{cc}\frac{\sqrt{10}}{10}& -\frac{\frac{3\hspace{0.17em}x}{10}-\frac{9}{5}}{\sqrt{\frac{x^2}{10}-\frac{6\hspace{0.17em}x}{5}+\frac{18}{5}}}\\ \frac{3\hspace{0.17em}\sqrt{10}}{10}& \frac{\frac{x}{10}-\frac{3}{5}}{\sqrt{\frac{x^2}{10}-\frac{6\hspace{0.17em}x}{5}+\frac{18}{5}}}\end{array}\right)$$

R = 
$$\left(\begin{array}{cc}\sqrt{10}& \frac{\sqrt{10}\hspace{0.17em}\left(3\hspace{0.17em}x+2\right)}{10}\\ 0& \sqrt{\frac{x^2}{10}-\frac{6\hspace{0.17em}x}{5}+\frac{18}{5}}\end{array}\right)$$

入力引数

A — シンボリック入力

m 行 n 列のシンボリック行列。

シンボリック入力。m 行 n 列のシンボリック行列として指定します。

B — シンボリック入力

シンボリック ベクトル | シンボリック行列

シンボリック入力。シンボリック ベクトルまたはシンボリック行列として指定します。B は A と同じ行数でなければなりません。

outputForm — 置換出力の形状

"matrix" (既定値) | "vector"

置換出力の形状。"matrix" または "vector" として指定します。このオプションは、置換出力 P が置換行列またはベクトルのいずれとして返されるかを制御します。このオプションを使用するには、qr に 3 つの出力引数を指定しなければなりません。

• outputForm が "matrix" の場合、P は A*P = Q*R を満たす置換行列。

• outputForm が "vector" の場合、P は A(:,P) = Q*R を満たす置換ベクトル。

出力引数

R — QR 分解の R 因子

上三角シンボリック行列

QR 分解の R 因子。m 行 n 列の行列 A に対して A = Q*R を満たす上三角シンボリック行列として返されます。

Q — QR 分解の Q 因子

ユニタリ シンボリック行列

QR 分解の Q 因子。m 行 n 列の行列 A に対して A = Q*R を満たすユニタリ シンボリック行列として返されます。

• 完全な分解の場合、qr(A) は Q を m 行 m 列のユニタリ行列として返す。

• m > n の方形行列 A の場合、エコノミーサイズの分解 qr(A,"econ") は、Q の最初の n 列と R の最初の n 行のみを計算する。

P — 置換情報

行列 | ベクトル

置換情報。倍精度値の行列またはベクトルとして返されます。P の形状は、outputForm の値によって異なります。

C — 線形系因子

シンボリック行列

線形系因子。C = Q'*B および X = R\C となる、行列方程式 A*X = B の解を表すシンボリック行列として返されます。置換出力 P を指定した場合、解は、outputForm が "matrix" の場合に X = P*(R\C)、outputForm が "vector" の場合に X(P,:) = R\C となります。

詳細

行列の QR 分解

QR 分解 (QR decomposition/factorization) は、m 行 n 列の行列 A を A = Q*R として表現するものです。完全な分解の場合、Q は m 行 m 列のユニタリ行列で、R は m 行 n 列の上三角行列になります。A の成分が実数の場合、Q は直交行列です。

ヒント

• 上三角行列 R は、条件 R = chol(A'*A) を満たします。

• 引数 "econ" と 0 は、返される行列の形状のみに影響します。

• シンボリック オブジェクトではない数値行列 (sym、syms または vpa により作成されていない行列) で qr を呼び出すと、MATLAB^® 関数 qr が起動します。

• 多くのシンボリック変数が含まれる行列計算は低速になる可能性があります。計算速度を向上させるには、特定の値を変数に代入することでシンボリック変数の数を減らします。

バージョン履歴

R2014a で導入

R2023a: エコノミーサイズの分解は置換情報を行列として返す

シンボリック入力のエコノミーサイズの分解用の構文 [Q,R,P] = qr(A,"econ") および [C,R,P] = qr(A,B,"econ") は、置換情報を行列 P として返すようになりました。以前のリリースでは、これらの構文は置換情報をベクトル P として返しました。

• R2023a 以降では、構文 qr(__,"econ") と構文 qr(__,"econ","matrix") は等価です。これらはいずれも置換情報を行列 P として返します。さらに、構文 qr(__,0) と qr(__,"econ","vector") は等価です。これらはいずれも置換情報をベクトル P として返します。

• R2023a より前では、構文 qr(__,"econ")、qr(__,0)、および qr(__,"econ","vector") は等価でした。これらはすべて、置換情報をベクトル P として返しました。構文 qr(__,"econ","matrix") は、以前はエラーを返しました。

参考

chol | eig | lu | svd