• 検索結果がありません。

12 重積分の変数変換

N/A
N/A
Protected

Academic year: 2021

シェア "12 重積分の変数変換 "

Copied!
2
0
0

読み込み中.... (全文を見る)

全文

(1)

82 (20130715) 12

12 重積分の変数変換

■線形変換と面積 R2の線形変換

LA:R2∋x7−→X=Ax∈R2 (A2次の正方行列)

を考える.行列A が正則,すなわちdetA̸= 0ならばLAは逆写像をもつ.

とくにLA 11の写像 (単射)である.行列Aが正則であるときLA 正則な線形変換とよぶ.

補題 12.1. 線形変換LA によるR2 の直線の像は直線または一点である.と

くにLAが正則ならば直線の像は直線になる.

証明.異なる2P,Q∈R2を結ぶ直線lの像を調べよう.P,Qの位置ベクトルを それぞれp,qとすると直線l

l={(1−t)p+tq|t∈R}

と表される.ここでLAの線形性から LA

((1−t)p+tq)

= (1−t)Ap+tAq なので,lLAによる像は

l={(1−t)˜p+t˜q|t∈R} ˜p=Ap, ˜q=Aq とかける.とくに−−→

OP= ˜p,−−→

OQ= ˜qとなる点P,Q をとると(1)P̸=Qのとき,

l P,Qを通る直線となる.(2)P=QのときlP 1点からなる集合である.

さらにdetA̸= 0なら写像LA11であるから(2)のケースは起こりえない.

補題12.2. 正則な線形変換LAによるR2の平行な2直線の像は平行な2 線である.

証明.平行な2直線の像は2つの直線であるが,これらが交わるとするとLA1 1であることに反する.

201379

12 (20130715) 83

補題 12.3. 直線 l 上の異なる2 P, Q をとっておく.直線 l にない2 R, S が直線 l の同じ側にあるための必要十分条件は,det(−→P R,−−→P Q) det(−→P S,−−→P Q)が同じ符号をもつことである.ここでR2 のベクトルは列ベク トルとみなし,det2つの2次列ベクトルを並べてできる行列の行列式を 表す.

証明.t(a, b) =−−→

P Qとおき,n=t(−b, a)とすると,(1) det(−−→

P Q,v) = (v,n)である.

ただし右辺はR2の内積を表す.(2)nは直線lに直交する零でないベクトルである.

直線l上にない点Rが,直線lnが指し示す側にあるための必要十分条件は−→P R nが鋭角をなすことである:(−→P R,n)>0.このことと(1)から結論が得られる.

補題 12.4. 線形変換LAによって,R2 の平行四辺形とその内部はR2の平 行四辺形とその内部,または線分に移る.とくにLA が正則ならば平行四辺 形の像は平行四辺形である.

証明.簡単のためLAが正則であるとし,平行四辺形P QRSの像を求める:p=−−→OP, q=−−→OQとすると,線分P Q{(1−t)p+tq|0≦t≦1}となるので,その像は線 P,Q となる.ただしP,Q はそれぞれLAによるP,Qの像.各辺に対して同 様のことを考えれば,平行四辺形の像が平行四辺形となることがわかる.さらに,平行 四辺形の内部は4つの辺を含む直線の一方の側の共通部分なので,補題12.3から結論 を得る(すこし端折った)

補題12.5. 平行四辺形P QRSの面積は|det(a,b)|である.ただしa=−−→P Q, b=−→P Rで,これらを2次の列べクトルとみなしている.

証明.ベクトルa,bのなす角をθとすると,求める面積は (12.1) |a| |b| |sinθ|=√

|a|2|b|2− |a|2|b|2cos2θ=√

|a|2|b|2−(a,b)2.

ただし (a,b) a, bの内積を表す.ここで a=t(a1, a2),b =t(b1, b2) とおいて (12.1)を計算すれば結論を得る.

補題12.6. 線形変換LA による平行四辺形D の像の面積は,|detA| |D| ある.ただし|D|D の面積である.

(2)

84 (20130715) 12 証明.平行四辺形D=P QRSの各頂点の位置ベクトルをp,q,r,sとし,

a=−−→

P Q=q−p, b=−→

P R=r−p

とおく.P,Q,RLAによる像をそれぞれP,Q,R と書くと,

−−−→PQ=Aq−Ap=A(q−p) =Aa, −−−→

PR=Ab であるから

|D|=|det(Aa, Ab)|=det(

A(a,b))=detA·det(a,b)=|detA| |D|.

■2変数の変数変換 R2 の領域上で定義されたC1-級写像 F:R2⊃(u, v)7−→F(u, v) =(

x(u, v), y(u, v))

∈R2 を考える.微分可能性から

F(a+h, b+k)

=F(a, b) +

(xu(a, b) xv(a, b) yu(a, b) yv(a, b)

) (h k )

+√

h2+k2ε(h, k)

|ε(h, k)| →0 as (h, k)→(0,0) と書ける.このt(h, k)の係数行列は,F の微分dF またはヤコビ行列(第6 節)である.このことから,(h, k)が十分小さいときは,近似式

Φ(h, k) :=F(a+h, b+k)−F(a, b)≑

(xu(a, b) xv(a, b) yu(a, b) yv(a, b)

) (h k )

が成り立つ.

記号. ヤコビ行列の行列式を

∂(x, y)

∂(u, v) = det

(xu xv

yu yv

)

と書き,ヤコビ行列式,Jacobianという.

12 (20130715) 85

以下,

∂(x, y)

∂(u, v) ̸= 0 が至るところ成立しているとする.

定理 12.7(重積分の変数変換;テキスト92ページ). 上の状況でxy平面上

のコンパクト集合Duv平面上のコンパクト集合EF によって11 に対応しているとき,

∫∫

D

f(x, y)dx dy=

∫∫

E

f(

x(u, v), y(u, v))∂(x, y)

∂(u, v) du dv が成り立つ.

問題

12

12-1 テキスト94ページ,問6,7; 95ページ,問8, 9.

12-2 テキスト97ページ問10.

12-3 テキスト100ページ,問11, 12; 101ページ,問13, 14 12-4 テキスト119ページ(章末問題)1, 2, 3, 4.

参照

関連したドキュメント

劣モジュラ解析 (Submodular Analysis) 劣モジュラ関数は,凸関数か? 凹関数か?... LP ニュートン法 ( の変種

WAV/AIFF ファイルから BR シリーズのデータへの変換(Import)において、サンプリング周波 数が 44.1kHz 以外の WAV ファイルが選択されました。.

[r]

LLVM から Haskell への変換は、各 LLVM 命令をそれと 同等な処理を行う Haskell のプログラムに変換することに より、実現される。

[r]

凡例及び面積 全体敷地 2,800㎡面積 土地の形質の変更をしよ うとする場所 1,050㎡面積 うち掘削を行う場所

各テーマ領域ではすべての変数につきできるだけ連続変量に表現してある。そのため

変更前変更後備考 (2) 浸水防護重点化範囲の境界における浸水対策 【検討方針】