直線の線形変換

【Pythonで学ぶ線形代数学講座(18)】直線の線形変換

直線の方程式と線形変換

下図のように、平面上に 2 点 $P,Q$ があって、それぞれの位置ベクトルを $\mathit{p},\mathit{q}$ とします。
 

このとき、$P$ を始点、$Q$ を終点とするベクトルは $\mathit{q}-\mathit{p}$ なので、$P$ と $Q$ を結ぶ直線上の位置ベクトルは、任意のスカラー $t$ によって、
 $$\mathit{r}=\mathit{p}+t(\mathit{q}-\mathit{p})$$
と表せます ($\mathit{q}-\mathit{p}$ の延長線上に $\mathit{r}$ があります)。この式を少し整理すると
 $$\begin{array}{}\text{(1)}& \mathit{r}=(1-t)\mathit{p}+t\mathit{q}\end{array}$$
となります。これが ベクトルで表された直線の方程式 です。

直線の線形変換

位置ベクトル $\mathit{p}$ と $\mathit{q}$ が線形変換 $f$ によって、それぞれ $f(\mathit{p})$ と $f(\mathit{q})$ に写され、かつ $\mathit{p}\ne \mathit{q}$ であるとき、方程式 (1) で表される直線 $\mathit{r}$ を $f$ で変換すると、
 $$\begin{array}{}\text{(2)}& f(\mathit{r})=(1-t)f(\mathit{p})+tf(\mathit{q})\end{array}$$
となり、$f(\mathit{p})$ と $f(\mathit{q})$ を結ぶ直線に写されることがわかります。
 

例として方程式
 $$\begin{array}{}\text{(3)}& 5x-3y=2\end{array}$$
で表される直線を行列 $A=\left[\begin{array}{cc}1& -1\\ 3& 2\end{array}\right]$ で変換してみましょう。2通りの方法で解いてみます。

[1] 2点を変換する方法

まず直線上の 2 点を見つけます。
$(x,y)=(1,1)$ が (3) を満たすことは明らかです。適当な $x$ を (3) に放り込んで、もう１つの点を見つけることもできますが、$y$ が分数になると面倒なので、ちょっとしたテクニックを使って整数解を得ておきましょう。(3) を
 $$3y=5x-2$$
のように書き直して、右辺を強引に $5$ で括って
 $$3y=5(x-1)+3$$
とします。$3$ を左辺に移項すれば、
 $$\begin{array}{}\text{(4)}& 3(y-1)=5(x-1)\end{array}$$
となります。$k$ を任意の実数として、
 $$\begin{array}{}\text{(5)}& y-1=5k,x-1=3k\end{array}$$
で表される点はすべて (4) を満たします。$k=1$ とすれば $(x,y)=(4,6)$ が得られます。$(1,1)$ と $(4,6)$ を行列 $A$ によって変換すると
 $$\begin{array}{rl}& \left[\begin{array}{cc}1& -1\\ 3& 2\end{array}\right]\left[\begin{array}{c}1\\ 1\end{array}\right]=\left[\begin{array}{c}0\\ 5\end{array}\right]\\ & \left[\begin{array}{cc}1& -1\\ 3& 2\end{array}\right]\left[\begin{array}{c}4\\ 6\end{array}\right]=\left[\begin{array}{c}-2\\ 24\end{array}\right]\end{array}$$
が得られます。$(0,5)$ と $(-2,24)$ を結ぶ直線の傾き $m$ を求めると、
 $$\begin{array}{}\text{(6)}& m=\frac{5-24}{0-(-2)}=-\frac{19}{2}\end{array}$$
したがって、求める直線の方程式は
 $$y=-\frac{19}{2}x+5$$
となります。分母を払って整理すると
 $$\begin{array}{}\text{(7)}& 19x+2y=10\end{array}$$
が得られます。

[2] 逆行列を用いる方法

直線 $5x-3y=2$ の上にある点を $(x,y)$ とし、行列 $A=\left[\begin{array}{cc}1& -1\\ 3& 2\end{array}\right]$ によって点 $(X,Y)$ に写されたとします。
 $$\begin{array}{}\text{(8)}& A\left[\begin{array}{c}x\\ y\end{array}\right]=\left[\begin{array}{c}X\\ Y\end{array}\right]\end{array}$$
$A$ の逆行列を計算すると
 $$\begin{array}{}\text{(9)}& A^{-1}=\frac{1}{5}\left[\begin{array}{cc}2& 1\\ -3& 1\end{array}\right]\end{array}$$
となるので、これを (8) の両辺に左側から掛けて、
 $$\begin{array}{}\text{(10)}& \left[\begin{array}{c}x\\ y\end{array}\right]=\frac{1}{5}\left[\begin{array}{cc}2& 1\\ -3& 1\end{array}\right]\left[\begin{array}{c}X\\ Y\end{array}\right]=\frac{1}{5}\left[\begin{array}{c}2X+Y\\ -3X+Y\end{array}\right]\end{array}$$
$(x,y)$ を元の方程式 $5x-3y=2$ に代入して整理すると、
 $$\begin{array}{}\text{(11)}& 19X+2Y=10\end{array}$$
が得られます。

Pythonによる直線の線形変換

Python で直線を変換してみましょう。
以下のコードを実行するためには、コードライブラリ の coordinate() が必要です。

# python_linear_transformation_line

# In[1]

import numpy as np
import matplotlib.pyplot as plt

# FigureとAxes
fig, ax = plt.subplots(1, 1, figsize=(6, 6))

# 座標を設定
coordinate(ax, [-6, 6], [-6, 6])

# 直線5x-3y=2上の点
x = np.linspace(-6, 6, 65)
y = (5*x - 2)/3

# xとyを縦軸方向に連結
p = np.vstack((x, y))

# 変換行列aを定義
a = np.array([[1, -1],
              [3, 2]])

# 直線5x-3y=2を行列aで変換
ap = np.dot(a, p)

# 直線を表示
ax.plot(x, y, color = "blue", label = "5x-3y=2")
ax.plot(ap[0], ap[1], color = "red", label = "19x+2y=10")
ax.legend()