一、求解的问题
从一个坐标系(x,y)中的已知三点坐标
(
x
i
,
j
,
k
,
y
i
,
j
,
k
)
(x_{i,j,k},y_{i,j,k})
(xi,j,k,yi,j,k),映射到另一个坐标系(u,v)中的三点坐标
[
(
1
,
0
)
,
(
0
,
1
)
,
(
0
,
0
)
]
[(1,0),(0,1),(0,0)]
[(1,0),(0,1),(0,0)]的问题,求解其中的映射关系,讨论不同坐标下的一个微分问题。
简单来说,就是求解一个坐标变换矩阵
M
M
M,将坐标系(x,y)中的三点坐标矩阵
A
A
A 转换到坐标系(u,v)
B
B
B 中去:
B
=
M
×
A
A
=
[
x
i
x
j
x
k
y
i
y
j
y
k
1
1
1
]
,
B
=
[
1
0
0
0
1
0
1
1
1
]
M
=
B
⋅
A
−
1
=
[
x
i
−
x
k
x
j
−
x
k
x
k
y
i
−
y
k
y
j
−
y
k
y
k
0
0
1
]
x
=
(
x
i
−
x
k
)
u
+
(
x
j
−
x
k
)
v
+
x
k
y
=
(
y
i
−
y
k
)
u
+
(
y
j
−
y
k
)
v
+
x
k
\begin{aligned} B&=M\times A\\ A=\begin{bmatrix} x_i &x_j& x_k\\ y_i &y_j& y_k\\ 1 &1 &1 \end{bmatrix}&,B=\begin{bmatrix} 1&0&0\\ 0&1&0\\ 1&1&1 \end{bmatrix}\\ M&=B\cdot A^{-1}\\ &=\begin{bmatrix} x_i-x_k&x_j-x_k&x_k\\ y_i-y_k&y_j-y_k&y_k\\ 0&0&1 \end{bmatrix}\\ x&=(x_i-x_k)u+(x_j-x_k)v+x_k\\ y&=(y_i-y_k)u+(y_j-y_k)v+x_k \end{aligned}
BA=⎣⎡xiyi1xjyj1xkyk1⎦⎤Mxy=M×A,B=⎣⎡101011001⎦⎤=B⋅A−1=⎣⎡xi−xkyi−yk0xj−xkyj−yk0xkyk1⎦⎤=(xi−xk)u+(xj−xk)v+xk=(yi−yk)u+(yj−yk)v+xk
则它们之间的微分关系为:
∣
d
e
t
(
J
)
∣
=
∣
d
x
d
u
d
x
d
v
d
y
d
u
d
y
d
v
∣
=
∣
x
i
−
x
k
x
j
−
x
k
y
i
−
y
k
y
j
−
y
k
∣
\begin{aligned} \left|det(J)\right|=\left|\begin{matrix} \dfrac{dx}{du}&\dfrac{dx}{dv}\\ \dfrac{dy}{du}&\dfrac{dy}{dv}\\ \end{matrix}\right|=\left|\begin{matrix} x_i-x_k&x_j-x_k\\ y_i-y_k&y_j-y_k\\ \end{matrix}\right| \end{aligned}
∣det(J)∣=∣∣∣∣∣∣∣dudxdudydvdxdvdy∣∣∣∣∣∣∣=∣∣∣∣xi−xkyi−ykxj−xkyj−yk∣∣∣∣
他们之间的面积微元关系为:
d
x
×
d
y
=
(
d
x
d
u
d
u
+
d
x
d
v
d
v
)
×
(
d
y
d
u
d
u
+
d
y
d
v
d
v
)
=
d
x
d
v
d
v
×
d
y
d
u
d
u
+
d
x
d
u
d
u
×
d
y
d
v
d
v
=
(
d
x
d
u
d
y
d
v
−
d
x
d
v
d
y
d
u
)
d
u
×
d
v
=
∣
d
x
d
u
d
x
d
v
d
y
d
u
d
y
d
v
∣
d
u
×
d
v
=
∣
det
(
J
)
∣
d
u
×
d
v
\begin{aligned} d x \times d y&=\left(\frac{d x}{d u} d u+\frac{d x}{d v} d v\right) \times\left(\frac{d y}{d u} d u+\frac{d y}{d v} d v\right)\\ &=\frac{d x}{d v} d v \times \frac{d y}{d u} d u+\frac{d x}{d u} d u \times \frac{d y}{d v} d v\\ &=\left(\frac{d x}{d u} \frac{d y}{d v}-\frac{d x}{d v} \frac{d y}{d u}\right) d u \times d v\\ &=\left|\begin{array}{ll}\dfrac{d x}{d u} & \dfrac{d x}{d v} \\ \dfrac{d y}{d u} & \dfrac{d y}{d v}\end{array}\right| d u \times d v\\ &=|\operatorname{det}(J)| d u \times d v \end{aligned}
dx×dy=(dudxdu+dvdxdv)×(dudydu+dvdydv)=dvdxdv×dudydu+dudxdu×dvdydv=(dudxdvdy−dvdxdudy)du×dv=∣∣∣∣∣∣∣dudxdudydvdxdvdy∣∣∣∣∣∣∣du×dv=∣det(J)∣du×dv