DifferentialGeometry Lessons

はじめる前に
新機能一覧
Maple ワークシートの作成
Mapleワークシートを共有
Maple ウィンドウのカスタマイズ
Maple システムのカスタマイズ
コネクティビティ
Mathematics
数学
- テンソル解析
- パッケージ
- ファンクションアドバイザ
- ベキ級数
- ベクトル解析
- 不活性関数
- 代数
- 因数分解と方程式解法
- 基本数学
- 基本的な注意事項
- 変分法
- 変換
- 幾何
- 微分代数方程式
- 微分幾何学
  - Group Actions
  - JetCalculus
  - LieAlgebras
  - Tensor
  - ツール
  - ライブラリ
  - レッスンおよびチュートリアル
    - LieAlgebras
    - チュートリアル
    - 微分幾何学
    - チュートリアルの概要
    - レッスンの概要
  - 概要
  - &wedge
  - Annihilator
  - ApplyTransformation
  - ChangeFrame
  - ComplementaryBasis
  - ComposeTransformations
  - DeRhamHomotopy
  - DGbasis
  - DGDualBasis
  - DGsetup
  - DGzip
  - evalDG
  - ExteriorDerivative
  - Flow
  - FrameData
  - GetComponents
  - Hook
  - InfinitesimalTransformation
  - IntegrateForm
  - IntersectSubspaces
  - InverseTransformation
  - LieBracket
  - LieDerivative
  - Pullback
  - PullbackVector
  - Pushforward
  - RemoveFrame
  - Transformation
  - 参考文献
  - 変換
  - 環境設定
- 微分方程式
- 微積分
- 数
- 数値計算
- 数学関数
- 数論
- 最適化
- 特殊関数
- 統計
- 線形代数
- 群論
- 評価
- 論理
- 過去のパッケージ
- 金融
- 離散数学
- piecewise examples
- グラフ電卓
- 区分関数
物理
プログラミング
グラフィックス
Student パッケージ
Science and Engineering
科学・エンジニアリング
アプリケーションと例題
リファレンス
システム
MapleSim
MapleSim ツールボックス
MapleSim ヘルプ
Tasks
Toolboxes
エラーメッセージガイド
マニュアル
数学アプリ

Lesson 6: Transformations

Overview

In this lesson, you will learn to do the following:

–	Define a transformation between two manifolds.

–	Apply a transformation to a point.

–	Compose two transformations.

–	Find the inverse of a transformation.

–	Define the identity transformation.

–	Test for the equality of two transformations.

–	Compute the Jacobian of a transformation.

–	Find the flow of a vector field.

Transformation, ApplyTransformation

>	with(DifferentialGeometry):

>	DGsetup([x, y], M):

M >	DGsetup([u, v], N):

To define a transformation from the manifold M to the manifold N, use the Transformation command. The arguments are the domain frame name, the range frame name, and a list of equations defining the map, one equation for each coordinate in the range.

N >	phi := Transformation(M, N, [u = x^2 + y^2, v = 2xy]);

(2.1)

Like all objects in the DifferentialGeometry environment, a transformation has its own internal representation. This internal representation of a transformation includes the Jacobian of the map.

M >

op(phi);

(2.2)

The data contained in this internal representation can be obtained using the DGinfo command. For example:

M >	Tools:-DGinfo(phi,"DomainFrame");

(2.3)

M >	Tools:-DGinfo(phi, "RangeFrame");

(2.4)

M >	Tools:-DGinfo(phi, "JacobianMatrix");

(2.5)

To apply the transformation to a point, use the ApplyTransformation command -- the second argument is the coordinates of a point in the domain frame. Note that the active frame changes to the range space after the invocation of the ApplyTransformation command.

(2.6)

M >	ApplyTransformation(phi, [3, 4]);

(2.7)

N >	Tools:-DGinfo("CurrentFrame");

(2.8)

N >	ApplyTransformation(phi, [x = a, y = b]);

(2.9)

Note that a transformation can be defined from a coordinate system to itself.

N >	psi := Transformation(M, M,[x = 2x + 3y + 1, y = -x + 2*y + 3]);

(2.10)

Compose transformations

The composition of transformations is performed with the ComposeTransformations command.

M >	with(DifferentialGeometry):

M >	DGsetup([x, y], M1):

M1 >	DGsetup([u, v], M2):

M2 >	DGsetup([r, s], M3):

M3 >	DGsetup([t], M4);

(3.1)

M4 >	phi := Transformation(M1, M2, [u = 3*x, v = y/2]);

(3.2)

M1 >	psi := Transformation(M2, M3, [r = u^2, s = 1/v]);

(3.3)

M2 >	sigma := Transformation(M4, M1, [x = sin(t), y = cos(t)]);

(3.4)

M4 >	ComposeTransformations(psi, phi);

(3.5)

More than two transformations can be composed with a single call to ComposeTransformation.

M1 >	ComposeTransformations(psi, phi, sigma);

(3.6)

Find the inverse of a transformation

The DifferentialGeometry command InverseTransformation uses the Maple solve command to find the inverse of a transformation. All possible inverse functions are calculated. With one argument, the first solution is returned -- this may change from one Maple session to another.

With the second argument branch = "all",all solutions are returned. With the second argument branch = [pt1, pt2], where pt1 is in the range and point pt2 is in the domain, InverseTransformation returns the inverse transformation which sends pt1 to pt2.

M4 >	with(DifferentialGeometry):

M4 >	DGsetup([x, y], M1):

M1 >	DGsetup([u, v], M2):

Define a transformation psi and calculate its inverse.

M2 >	psi := Transformation(M1, M2, [u = x^2 + y^2, v = 2xy]);

(4.1)

M1 >	InverseTransformation(psi);

(4.2)

Set the Maple global variable _EnvExplicit to true to get explicit solutions.

M2 >	_EnvExplicit := true;

(4.3)

M2 >	InverseTransformation(psi);

[x = v/(2*((u+v)^(1/2)/2+(u-v)^(1/2)/2)), y = (u+v)^(1/2)/2+(u-v)^(1/2)/2]

(4.4)

M2 >	InverseTransformation(psi, branch = "all");

[[x = v/(2*((u+v)^(1/2)/2+(u-v)^(1/2)/2)), y = (u+v)^(1/2)/2+(u-v)^(1/2)/2], [x = v/(2*((u+v)^(1/2)/2-(u-v)^(1/2)/2)), y = (u+v)^(1/2)/2-(u-v)^(1/2)/2], [x = v/(2*(-(u+v)^(1/2)/2+(u-v)^(1/2)/2)), y = -(u+v)^(1/2)/2+(u-v)^(1/2)/2], [x = v/(2*(-(u+v)^(1/2)/2-(u-v)^(1/2)/2)), y = -(u+v)^(1/2)/2-(u-v)^(1/2)/2]]

(4.5)

M2 >	InverseTransformation(psi, branch = [[2, 2], [1, 1]]);

[x = v/(2*((u+v)^(1/2)/2+(u-v)^(1/2)/2)), y = (u+v)^(1/2)/2+(u-v)^(1/2)/2]

(4.6)

M1 >	InverseTransformation(psi, branch = [[2, 2], [-1, -1]]);

[x = v/(2*(-(u+v)^(1/2)/2+(u-v)^(1/2)/2)), y = -(u+v)^(1/2)/2+(u-v)^(1/2)/2]

(4.7)

Define the identity transformation. Test if two transformations are equal. DGequal

The commands IdentityTransformation and DGequal are found in the package Tools.

M1 >	with(DifferentialGeometry): with(Tools):

M1 >	DGsetup([x, y, z], E3):

E3 >	phi := IdentityTransformation();

(5.1)

E3 >	psi := Transformation(E3, E3, [x = x, y = y, z = z]);

(5.2)

E3 >	DGequal(phi, psi);

(5.3)

Find the Jacobian Matrix for a transformation

The command DGinfo/JacobianMatrix returns the Jacobian matrix of a transformation.

E3 >	with(DifferentialGeometry):

E3 >	DGsetup([x, y], M): DGsetup([u, v], N):

N >	psi := Transformation(M, N, [u = exp(x)sin(y), v = exp(x)cos(y)]):

M >	Tools:-DGinfo(psi, "JacobianMatrix");

(6.1)

Find the flow of a vector field. Find the infinitesimal transformations (vector fields) for a multi-parameter group of transformations.

The command Flow uses the Maple dsolve command to compute the flow of a vector field. The command InfinitesimalTransformation is the inverse to the command Flow, it calculates the vector field defined by a r-parameter group of transformations.

M >	with(DifferentialGeometry):

M >	DGsetup([x, y], E2);

(7.1)

E2 >	X := evalDG(-yD_x + xD_y);

(7.2)

E2 >	Rot := Flow(X, theta);

(7.3)

E2 >	Tx := Flow(D_x, a);

(7.4)

E2 >	Ty := Flow(D_y, b);

(7.5)

Composing these three transformations gives the group of Euclidean motions in the plane.

E2 >	psi := ComposeTransformations(Ty, Tx, Rot);

(7.6)

The command InfinitesimalTransformation is the inverse to the command Flow, it calculates the vector fields defined by an r-parameter group of transformations.

E2 >	InfinitesimalTransformation(psi, [theta, a, b]);

(7.7)

Exercises

Exercise 1

Check, by means of an example, the chain rule for the Jacobian of the composition of two maps.

Solution

E2 >	with(DifferentialGeometry):

E2 >	DGsetup([x, y], E2):

E2 >	DGsetup([u, v, w], E3):

E3 >	DGsetup([p, q], M):

M >	phi := Transformation(E2, E3,[u = xy, v = xz, w = y*z]);

(8.1.1.1)

E2 >	psi := Transformation(E3, M, [p = u/v, q = w/v]);

(8.1.1.2)

E3 >	sigma := ComposeTransformations(psi, phi);

(8.1.1.3)

E2 >	J1 := Tools:-DGinfo(phi, "JacobianMatrix");

(8.1.1.4)

E2 >	J2 := Tools:-DGinfo(psi, "JacobianMatrix");

(8.1.1.5)

To verify the chain rule, we have to evaluate this matrix at psi. In the next lesson (Lesson 6), we shall see how to do this with the Pushforward command, but for now we just use the subs command.

E2 >	J2psi := map2(subs, {u = xy, v = xz, w = y*z}, J2);

(8.1.1.6)

E2 >	J3 := Tools:-DGinfo(sigma, "JacobianMatrix");

(8.1.1.7)

E2 >	evalm(J3 - J2psi &* J1);

Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0})

(8.1.1.8)

Exercise 2

E2 >	with(DifferentialGeometry):

E2 >	DGsetup([x,y], M2):

M2 >	DGsetup([u,v], N2):

N2 >	DGsetup([x, y, z], M3):

M3 >	DGsetup([t], R):

[i] At what points do the following maps fail to be a local diffeomorphism?

(a)

R >	phi1 := Transformation(M2, N2, [u = x/(1 + x^2 + y^2), v = y/(1 + x^2 + y^2)]);

(8.2.1)

(b)

M2 >	phi2 := Transformation(M3, M3, [x = x + y + z, y = 1/2(x^2 + y^2 + z^2), z = 1/3(x^3 + y^3 + z^3)]);

(8.2.2)

[ii] At what points do the following maps fail to be local immersions?

(a)

M3 >	phi3 := Transformation(N2, M3, [x = u*v, y = u + v, z = u^2 + v^2]);

(8.2.3)

(b)

N2 >	phi4 := Transformation(R, M3, [x = t, y = 1/sqrt(2)t^2, z = 1/3t^3]);

(8.2.4)

Solution

[i](a)

R >	J := Tools:-DGinfo(phi1, "JacobianMatrix");

(8.2.1.1)

R >	d := LinearAlgebra:-Determinant(J);

d := -(x^2+y^2-1)/(1+x^2+y^2)^3

(8.2.1.2)

The map phi1 fails to be a local diffeomorphism at points on the unit circle x^2 + y^2 = 1.

[i](b)

M2 >	J := Tools:-DGinfo(phi2, "JacobianMatrix");

R >	d := factor(LinearAlgebra:-Determinant(J));

(8.2.1.3)

The map phi2 fails to be a local diffeomorphism at points on the lines x = y, x = z, or y = z.

[ii](a) The map phi2 fails to be a local immersion at points where the Jacobian has rank 1, that is, at points where the 2 x 2 sub-determinants of the Jacobian matrix vanish.

R >	J := Tools:-DGinfo(phi3, "JacobianMatrix");

(8.2.1.4)

R >	J1 := LinearAlgebra:-SubMatrix(J, [1, 2], [1, 2]);

(8.2.1.5)

R >	J2 := LinearAlgebra:-SubMatrix(J, [1, 3], [1, 2]);

(8.2.1.6)

R >	J3 := LinearAlgebra:-SubMatrix(J, [2, 3], [1, 2]);

(8.2.1.7)

R >	use LinearAlgebra in eq := Determinant(J1)^2 + Determinant(J2)^2 + Determinant(J3)^2 end;

(8.2.1.8)

R >	factor(%);

(8.2.1.9)

Thus phi3 fails to be a local immersion at points where u = v.

[ii](b)

R >	J := Tools:-DGinfo(phi4, "JacobianMatrix");

(8.2.1.10)

R >	factor((J[1, 1]^2 + J[2, 1]^2 + J[3, 1]^2));

(8.2.1.11)

The map phi4 is a local immersion at all points.

Exercise 3

Find the flow phi_t of each of the following vector fields and check that phi_t o phi_s = phi_(t + s) and (phi_t)^(-1) = phi_(-t).

R >	with(DifferentialGeometry):

R >	DGsetup([x, y], M):

M >	DGsetup([u, v, w], N):

[i]

N >	X1 := evalDG(xyD_x - y^2*D_y);

(8.3.1)

[ii]

M >	ChangeFrame(N):

N >	X2 := evalDG(vD_u - uD_v + w*D_w);

(8.3.2)

Solution

[i]

N >	phi := simplify(Flow(X, t));

(8.3.1.1)

E2 >	InfinitesimalTransformation(phi, [t]);

(8.3.1.2)

E2 >	phi1 := phi:

E2 >	phi2 := eval(phi, t = s);

(8.3.1.3)

E2 >	phi3 := eval(phi, t = t + s);

(8.3.1.4)

E2 >	phi4 := ComposeTransformations(phi1, phi2);

(8.3.1.5)

E2 >	phi4 := DGmap(1, collect, phi4, [x, y]);

(8.3.1.6)

Check that phi_(-t) is the inverse of phi_t.

E2 >	phi5 := eval(phi, t = -t);

(8.3.1.7)

E2 >	phi6 := ComposeTransformations(phi1, phi5);

(8.3.1.8)

E2 >	simplify(phi6);

(8.3.1.9)

[ii]

E2 >	psi := Flow(X2, t);

(8.3.1.10)

N >	InfinitesimalTransformation(psi, [t]);

(8.3.1.11)

N >	psi1 := psi:

N >	psi2 := eval(psi, t = s);

(8.3.1.12)

N >	psi3 := eval(psi, t = t + s);

(8.3.1.13)

N >	map(expand, psi3);

(8.3.1.14)

N >	psi4 := map(expand, ComposeTransformations(psi1, psi2));

psi4 := [u = v*sin(s)*cos(t)+v*cos(s)*sin(t)-u*sin(s)*sin(t)+u*cos(s)*cos(t), v = -v*sin(s)*sin(t)+v*cos(s)*cos(t)-u*sin(s)*cos(t)-u*cos(s)*sin(t), w = w*exp(t)*exp(s)]

(8.3.1.15)

Check that psi_(-t) is the inverse of psi_t.

N >	psi5 := eval(psi, t = -t);

(8.3.1.16)

N >	psi6 := ComposeTransformations(psi1, psi5);

psi6 := [u = (v*cos(t)+u*sin(t))*sin(t)+(-v*sin(t)+u*cos(t))*cos(t), v = (v*cos(t)+u*sin(t))*cos(t)-(-v*sin(t)+u*cos(t))*sin(t), w = w*exp(-t)*exp(t)]