next up previous contents
Next: 4 QR 分解 Up: 3.2 Givens 変換 Previous: 計算の工夫 (1) を与えられた行列 にかける

計算の工夫 (2) 与えられた行列 $ A$ $ G$ で 相似変換 ($ G^T A G$ )

$\displaystyle B:= G^T A G, \quad G:=G(p,q,\theta)
$

$ p$ , $ q$ 行目, $ p$ , $ q$ 列目を除き $ A$ と変わらない。
  $\displaystyle b_{ij}$ $\displaystyle =$ $\displaystyle a_{ij}$   $\displaystyle \mbox{($i\ne p,q$ かつ $j\ne p,q$)}$$\displaystyle ,$
  $\displaystyle b_{pj}$ $\displaystyle =$ $\displaystyle a_{pj}\cos\theta-a_{qj}\sin\theta$   $\displaystyle \mbox{($j\ne p,q$)}$$\displaystyle ,$
  $\displaystyle b_{qj}$ $\displaystyle =$ $\displaystyle a_{pj}\sin\theta+a_{qj}\cos\theta$   $\displaystyle \mbox{($j\ne p,q$)}$$\displaystyle ,$
  $\displaystyle b_{ip}$ $\displaystyle =$ $\displaystyle a_{ip}\cos\theta-a_{iq}\sin\theta=b_{pi}$   $\displaystyle \mbox{($i\ne p,q$)}$$\displaystyle ,$
  $\displaystyle b_{iq}$ $\displaystyle =$ $\displaystyle a_{ip}\sin\theta+a_{iq}\cos\theta=b_{qi}$   $\displaystyle \mbox{($i\ne p,q$)}$$\displaystyle ,$
  $\displaystyle b_{pp}$ $\displaystyle =$ $\displaystyle a_{pp}\cos^2\theta-2a_{pq}\cos\theta\sin\theta+a_{qq}\sin^2\theta,$
  $\displaystyle b_{qq}$ $\displaystyle =$ $\displaystyle a_{pp}\sin^2\theta+2a_{pq}\cos\theta\sin\theta+a_{qq}\cos^2\theta
= a_{pp}+a_{qq}-b_{pp},$
  $\displaystyle b_{pq}$ $\displaystyle =$ $\displaystyle (a_{pp}-a_{qq})\cos\theta\sin\theta
+a_{pq}(\cos^2\theta-\sin^2\theta),$
  $\displaystyle b_{qp}$ $\displaystyle =$ $\displaystyle b_{pq}.$


next up previous contents
Next: 4 QR 分解 Up: 3.2 Givens 変換 Previous: 計算の工夫 (1) を与えられた行列 にかける
桂田 祐史
2015-12-22