[size=150][b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br][/b][size=100]ガウスの定理、ストークスの定理の共通点は[br]中で打ち消しあうから、表面だけみても同じだという内容だけでない。[br]両方とも∇(ナブラ)というベクトルの微分を積分するという点もある。[br][br][/size][/size][b][size=150]そこで、まずナブラを思い出しながら広げてみよう。[br][br]<∇の思い起こし>[/size][/b][br][color=#0000ff]・[b]ヤングの定理[/b]が便利。U(x,y)が連続微分できるなら、2つの偏微分は入れ替え可能だ。(Uxy=Uyx)[br][/color]・fの勾配ベクトルgrad f=(∂f/∂x, ∂f/∂y) とか、(∂f/∂x, ∂f/∂y, ∂f/∂z) [br]のfの前の部分を演算子として取り出した(∂/∂x, ∂/∂y)とか(∂/∂x, ∂/∂y, ∂/∂z)を∇と名付けたね。[br][br]・基本の使い道は2つで、[br] ベクトル[b]A[/b]=(Ax,Ay,Az)に対して[br][b] 内積∇・A=div A(Aの 流れ出し、発散)[br] 外積∇×A=rot A(Aのうず、回転)[br]・rot(grad f) =0 [br][/b] (∂/∂x, ∂/∂y, ∂/∂z)×(∂f/∂x, ∂f/∂y, ∂f/∂z)[br] =( fzy-fyz, fxz- fzx, fxy-fyx) =(0,0,0)=[b]0[br][/b] 勾配ベクトルの回転は0[br][b]・div(rot A) =0[br][/b] (∂/∂x, ∂/∂y, ∂/∂z)・(∂/∂x, ∂/∂y, ∂/∂z)×(Ax,Ay,Az)[br] =(∂/∂x, ∂/∂y, ∂/∂z)・(∂Az/∂y-∂Ay/∂z, ∂Ax/∂z-∂Az/∂x, ∂Ay/∂x-∂Ax/∂y)[br] = (Az)yx-(Ay)zx+ (Ax)zy-(Az)xy+(Ay)xz-(Ax)yz =0[br] うずによる立体的な流れ出しは0[br][b]・∇[/b][sup]2[/sup][b]=∇・∇=(∂[/b][sup]2[/sup][b]/∂x[/b][sup]2[/sup][b], ∂[/b][sup]2[/sup][b]/∂y[/b][sup]2[/sup][b], ∂[/b][sup]2[/sup][b]/∂z[/b][sup]2[/sup][b])をラプラシアンという。[br][/b] 2回偏微分は上に凸なグラフなら負、下に凸なグラフなら正。 直線のグラフはゼロ。[br] 成分ごとの総和はx,y,zの増減で元に戻す力、粘りに対応している。[br][b]・∇・∇[sup][/sup]×A=0[br][/b] div (rot A)だから、0[br][b]・∇×∇f=0[br][/b] rot(grad f) だから、[b]0[br][color=#0000ff]・div(fA)= ∇・(fA) = (∇f)・A+ f(∇・A) スカラーとベクトルの積の∇微分公式[br][/color][/b] ∇・(fA)=(∂/∂x, ∂/∂y, ∂/∂z)・(fAx,fAy,fAz)[br] =∂/∂x (fAx)+ ∂/∂y(fAy) + ∂/∂z(fAz)[br] =∂f/∂x (Ax)+f ∂(Ax)/∂x + ∂f/∂y(Ay)+ f ∂(Ay)/∂y+ ∂f/∂z(Az)+ f∂(Az)/∂z[br] =∂f/∂x (Ax)+ ∂f/∂y(Ay)+ ∂f/∂z(Az)+f ( ∂(Ax)/∂x + ∂(Ay)/∂y+ ∂(Az)/∂z )[br]=(∂f/∂x, ∂f/∂y, ∂f/∂z)・(Ax+Ay+Az)+ f((∂/∂x, ∂/∂y, ∂/∂z)・(Ax,Ay,Az)= (∇f)・A+f(∇・A)
[b][size=150][br]<ガウスの定理>[/size][/b][br][color=#0000ff] ガウスの定理は、発散の積分について定理。これって、ところてんの定理?[br][/color] 閉曲面Sで囲まれた領域Vがベクトル場f=(f,g,h)のとき、[br] 微小体積の熱変化量[∂f/∂x +∂g/∂y+∂h/∂z]⊿x⊿y⊿z=div f を集めると、[br] 有限体積Vの熱変化量div f の体積分が出せる。[br][br][b][color=#0000ff] 領域V全体のfの総和(体積分)= 表面Sだけのfの内積の総和(面積分)[br][/color][/b] 数式にすると、[br][color=#0000ff] ∫[sub]V[/sub] div [b]f [/b]dV=∫[sub]S[/sub] [b]f[/b] ・d[b]S[/b] (ガウスの発散定理)[br][br][/color] 微小立体がつまっている状態では、隣接面での熱の出入りは相殺される。[br] だから、[color=#0000ff]集めたあとの全体の立体の表面だけで熱の出入りを求めることができるはず[/color]だね。[br][br]・ 数式で考えると [br] 積分の線形性、 f=∫∂f/∂x dx, g=∫∂g/∂y dy, h =∫∂h/∂z dz、[br] d[b]S[/b]=(dSx,dSy,dSz)=(dydz, dxdz, dxdy)[br] などから、[br] ∫[sub]v[/sub] div [b]f[/b] d[b]V[/b]=[br] ∫∫∫[sub]v[/sub] div [b]f[/b] dxdydz=∫∫∫[sub]v[/sub] (∂f/∂x +∂g/∂y+∂h/∂z ) dxdydz[br] =∫∫[sub]Sx [/sub](∫∂f/∂x dx)dydz+∫∫[sub]Sy [/sub](∫∂g/∂y dy)dxdz+∫∫[sub]Sz [/sub](∫∂h/∂z dz)dxdy[br] =∫∫[sub]Sx[/sub] f dydz+∫∫[sub]Sy[/sub] g dxdz+∫∫[sub]Sz[/sub] h dxdy[br] =∫[sub]S[/sub] f dSx+∫[sub]S[/sub] g dSy+∫[sub]S[/sub] h dSz[br] =∫[sub]S[/sub] (f dSx+g dSy+h dSz)[br] =∫[sub]S[/sub] [b]f[/b] ・d[b]S[/b] [br][color=#0000ff](例)[/color][br]閉曲面Sの領域Vで[br]全電荷をQ、電場をE、電荷密度をρ、真空の誘電率をε[sub]0[/sub]とすると、Sはベクトル場Eとなり、[br][b][color=#0000ff]∫[sub]S[/sub]E・dS=Q/ε[sub]0[/sub]=∫[sub]v[/sub] ρ dV /ε[sub]0[/sub] =∫[sub]v[/sub] 1/ε[sub]0[/sub] ρ dV[/color][/b] これを[b]ガウスの法則(積分形)[/b]という。[br]一方、ガウスの定理から∫[sub]S[/sub]E・dS=∫[sub]V [/sub]divE dVだから、[br][b]divE=1/ε0 ρ[/b]。これを[b]ガウスの法則(微分形)[/b]という。[br][color=#0000ff](例)[/color][br]半径aの球面S内の領域Vで一様に電荷ρがあるとき、中心からa以下のrの点の電場E(r)を求める。[br][b]ガウスの法則の積分形∫[sub]S[/sub]E・dS(面積分)=∫[sub]v[/sub] 1/ε[sub]0[/sub] ρ dV(体積分)[br][/b]半径rの体積4/3πa[sup]3[/sup]、半径rの表面積は4πr[sup]2 [/sup][br]体積分 ρ/ε[sub]0 [/sub]∫[sub]v[/sub] 1 dV =ρ/ε[sub]0 [/sub] 4/3πr[sup]3[/sup][br]面積分 ∫[sub]S[/sub]E(r)・dS=E(r) 4πr[sup]2 [/sup][br]半径r以下の領域では、上記の体積分と面積分が等しいので、[b]E(r)[/b]=ρ/ε[sub]0 [/sub] 4/3πr[sup]3[/sup]/4πr[sup]2[/sup]=[b]ρr /3ε[/b][sub][b]0[/b] [/sub] [br][color=#0000ff](例)[/color][br]半径1の球面Sはベクトル場[b]f[/b]=(ax, by,cz) の面積分∫[sub]S[/sub][b]f[/b]・d[b]S[/b]をガウスの定理で求める。[br]半径rの体積4/3πa[sup]3[/sup]、半径rの表面積は4πr[sup]2 [/sup][br]div f= ∇・f=∂(ax)/∂x +∂(by)/∂y+∂(cz)/∂z=a+b+c[br]面積分∫[sub]S[/sub][b]f[/b]・d[b]S[/b]=体積分∫[sub]v[/sub] div [b]f[/b] d[b]V[/b]=∫[sub]v[/sub] (a+b+c)d[b]V[/b]= (a+b+c )∫[sub]v[/sub] 1d[b]V[/b]=4/3π(a+b+c )
[size=150][b]<立体のグリーンの定理>[/b][br][/size][color=#0000ff]ガウスの法則は体積分を面積分に直せるという、[br]中抜き、空洞化の定理だった。[br]これを使うと、立体のグリーンの定理というものが証明できる。[br][/color][br] 閉曲面Sで囲まれた領域Vにスカラー場f(x,y,z)、g(x,y,z)があるとき、[br] ∂f/∂n, ∂g/∂nはSに対する外向き法線方向の方向微分係数としよう。[br] そのとき、つぎの「体積分=面積分」という関係が2つ成り立つことを「グリーンの定理」という。[br][b][color=#0000ff] (式1)∫[sub]v[/sub]{ f ∇[sup]2[/sup]g +∇f・∇g } dV = ∫[sub]S[/sub] f ∂g/∂n dS[br] (式2)∫[sub]v[/sub]{ f ∇[sup]2[/sup]g - g ∇[sup]2[/sup]f } dV = ∫[sub]S[/sub] (f ∂g/∂n- g ∂f/∂n)dS[br][/color][/b][color=#0000ff](理由)[br][/color] ガウスの定理∫[sub]V[/sub] div [b]f [/b]dV=∫[sub]S[/sub] [b]f[/b] ・d[b]S[/b]=∫[sub]S[/sub] [b]f [/b]・[b]n[/b] dS[br]・[b]f[/b]にf ∇gを代入すると[br] ∫[sub]V[/sub] div(f ∇g) dV = ∫[sub]S [/sub](f ∇g)・[b]n[/b] dS[br] 左辺の積分で、[b]積の∇微分公式[/b]から、div(f ∇g)=∇・(f∇g)=∇f・∇g+f∇[sup]2[/sup]g[br] 右辺の積分で、∇g・[b]n[/b] dS= ∂g/∂n [b]n・n[/b] dS=∂g/∂n dS [br] これから式1がでる。[br]・[b]f[/b]にg ∇fを代入すると[br] 式1のf,g入れかえ式ができる。[br] (式1R)∫[sub]v[/sub]{ g ∇[sup]2[/sup]f +∇g・∇f } dV = ∫[sub]S[/sub] g∂f/∂n dS[br] 式1Rー式1から、∫[sub]v[/sub]{ f ∇[sup]2[/sup]g - g ∇[sup]2[/sup]f } dV = ∫[sub]S[/sub] (f ∂g/∂n- g ∂f/∂n)dS[br] これが式2だね。[br]
[b][size=150]<平面のグリーンの定理>[/size][/b][br][color=#0000ff][b]閉曲線Cでかこかれた領域Sのスカラー場f,gに対して、[br] ∫∫[sub]S[/sub] (∂g/∂x-∂f/∂y) dxdy= ∳[sub]C[/sub] (f dx +g dy) [br][/b][/color]積分の線形性から、[br](式1)∫∫[sub]S[/sub] (-∂f/∂y) dxdy= ∳[sub]C[/sub] f dx [br](式2)∫∫[sub]S[/sub] (∂g/∂x) dxdy= ∳[sub]C[/sub] g dy [br]に分けて考えてよいね。[br]・式1[br]閉曲線Cでかこかれた領域Sのスカラー場fに対して、[br]面積分∫∫[sub]S[/sub] (-∂f/∂y) dxdyをすることは、[br]閉曲線Cをx値の下端a, 上端bの点A,Bで2分割して、[br]y座標の上側をC2(y=y2),下側をC1(y=y1)として反時計回りに線積分することと同じになる。[br]∫∫[sub]S[/sub] (-∂f/∂y) dxdy=-∫[sub]a[/sub][sup]b[/sup]∫[sub]y1[/sub][sup]y2[/sup](∂f/∂y dy) dx[br]=-∫[sub]a[/sub][sup]b[/sup][f] [sub]y1[/sub][sup]y2[/sup] dx = -∫[sub]a[/sub][sup]b[/sup][f(y2)-f(y1)] dx [br]= ∫[sub]b[/sub][sup]a[/sup]f(y2)dx+∫[sub]a[/sub][sup]b[/sup]f(y1)dx=∫[sub]C2[/sub]f(y2)dx+∫[sub]C1[/sub]f(y1)dx=∫[sub]C[/sub]f dx[br]上側C2がBA順の経路、下側C1がAB順の経路の積分だから、[br]AとBでつないで[b]閉曲線Cの1周の経路積分になる[/b]ことがわかるね。図をかくとよくわかるでしょう。[br]・式2[br]閉曲線Cでかこかれた領域Sのスカラー場gに対して、[br]面積分∫∫[sub]S[/sub] (∂g/∂y) dxdyをすることは、[br]閉曲線Cをy値の下端a, 上端bの点A,Bで2分割して、[br]x座標の右側をC2(x=x2),下側をC1(x=x1)として反時計回りに線積分することと同じになる。[br]∫∫[sub]S[/sub] (∂g/∂x) dxdy=∫[sub]a[/sub][sup]b[/sup]∫[sub]x1[/sub][sup]x2[/sup](∂g/∂x dx) dy[br]=∫[sub]a[/sub][sup]b[/sup][g] [sub]x1[/sub][sup]x2[/sup] dy = ∫[sub]a[/sub][sup]b[/sup][g(x2)-g(x1)] dy [br]= ∫[sub]a[/sub][sup]b[/sup]g(x2)dy+∫[sub]b[/sub][sup]a[/sup]g(x1)dy=∫[sub]C2[/sub]g(x2)dy+∫[sub]C1[/sub]g(x1)dy=∫[sub]C [/sub]gdy[br]右側C2がAB順の経路、下側C1がBA順の経路の積分だから、[br]AとBでつないで閉曲線Cの1周の経路積分になることがわかるね。[br][br]★閉曲線にへこみがあるときは、適当に凸になるように領域を分割する。[br]各領域で平面のグリーンの定理が成り立ち、分割線で切った境界線での積分方向が逆向きになるから、分割した積分を合体することで打ち消しあう。[br]だから、へこみのある領域でも平面のグリーンの定理は成り立つことがわかるね。
[size=150][b]<ストークスの定理>[/b][br][/size][color=#0000ff]ストークスの定理は、うず、回転の積分の定理。これって、メロンパンの定理?[br][/color] 閉曲線Cを境界線とする曲面Sがベクトル場[b]f[/b]=(f,g,h)のとき、[br] 微小領域dSでのうずであるrot [b]f[/b] =[b]∇×f[/b] を集めると、曲面S全体の合計がでる。[br] ガウスの発散の積分と同様に隣接する微小面どうしのうずは打ち消しあうので、[br] 境界線Cの線素ベクトルd[b]p[/b]=(dx,dy,dz)にそって合計してもよい。[br][br][b][color=#0000ff]曲面S全体のうず∇fの総和(面積分)= 境界Cだけのfの総和(線積分)[br][/color][/b][color=#0000ff] ∫[sub]S[/sub] (rot [b]f[/b])・ d[b]S[/b] =∫[sub]c[/sub] [b]f [/b]・d[b]p[/b] ストークスの定理。[/color][br][br]・数式で考える。[br] パラメータu,vが曲線C0でかこまれた領域S0を動くときに、[br] 空間内の曲面Sの境界線がCで、Sの方程式がr(X,Y,Z)=r(u,v)としよう。[br] 位置[b]r[/b]の2つの接ベクトルが∂[b]r[/b]/∂u,∂[b]r[/b]/∂vになるから、[br] 接平面がdS=|∂[b]r[/b]/∂u×∂[b]r[/b]/∂v|dudv, 単位法線ベクトル[b]n[/b]=∂[b]r[/b]/∂u×∂[b]r[/b]/∂v /|∂[b]r[/b]/∂u×∂[b]r[/b]/∂v|[br] ∇=(∂/∂x, ∂/∂y, ∂/∂z) , f=(f,g,h) から、[br] rot [b]f[/b]=∇×f = ((∂h/∂y-∂g/∂z), (∂f/∂z- ∂h/∂x) (∂g/∂x-∂f/∂y)) [br] ∂[b]r[/b]/∂u×∂[b]r[/b]/∂v=(Xu,Yu,Zu)×(Xv,Yv,Zv)=(YuZv-ZuYv, XvZu-XuZv,XuYv-YuXv)[br] [color=#0000ff](rot [b]f[/b])・d[b]S[/b][/color]= (rot [b]f[/b])・[b]n[/b] dS= (rot [b]f[/b])・∂[b]r[/b]/∂u×∂[b]r[/b]/∂v /|∂[b]r[/b]/∂u×∂[b]r[/b]/∂v| |∂[b]r[/b]/∂u×∂[b]r[/b]/∂v|dudv[br][color=#0000ff]=[/color][color=#0000ff](rot [b]f[/b])・∂[b]r[/b]/∂u×∂[b]r[/b]/∂v dudv[br][/color]= ((∂h/∂y-∂g/∂z), (∂f/∂z- ∂h/∂x) (∂g/∂x-∂f/∂y))・(YuZv-ZuYv, XvZu-XuZv,XuYv-YuXv) dudv[br]=[ (∂h/∂y-∂g/∂z)(YuZv-ZuYv)+ (∂f/∂z- ∂h/∂x) (XvZu-XuZv)+ (∂g/∂x-∂f/∂y)(XuYv-YuXv)]dudv[br]また、Uが偏微分可能な連続関数ならUxUy=UyUzだからUxUy-UyUx=0を数式にたしてよい。[br][br]・ここからが面積分の開始[br]∫[sub]S[/sub] (rot [b]f[/b])[b]・[/b]d[b]S[/b]=[br]∫[sub]S0[/sub][ (∂h/∂y-∂g/∂z)(YuZv-ZuYv)+ (∂f/∂z- ∂h/∂x) (XvZu-XuZv)+ (∂g/∂x-∂f/∂y)(XuYv-YuXv)]dudv[br]=∫[sub]S0[/sub][ (∂f/∂z (XvZu-XuZv) -∂f/∂y(XuYv-YuXv) [br] +∂g/∂x(XuYv-YuXv) -∂g/∂z)(YuZv-ZuYv)[br] +∂h/∂y(YuZv-ZuYv) - ∂h/∂x) (XvZu-XuZv)]dudv[br]=∫[sub]S0[/sub][ (∂f/∂z XvZu-∂f/∂zXuZv -∂f/∂yXuYv+∂f/∂yYuXv [br] +∂g/∂xXuYv-∂g/∂xYuXv -∂g/∂zYuZv+∂g/∂zZuYv[br] +∂h/∂yYuZv-∂h/∂yZuYv - ∂h/∂xXvZu+∂h/∂xXuZv)] dudv[br]=∫[sub]S0[/sub][ (Xv(∂f/∂z Zu+∂f/∂yYuXv) -Xu(∂f/∂zZv +∂f/∂yYv) +∂f/∂x(XuXv-XvXu) 0になる項を加える[br] +Yv(∂g/∂xXu+∂g/∂zZu)-Yu(∂g/∂xXv +∂g/∂zZv) +∂g/∂y(YuYv-YvYu)[br] +Zv(∂h/∂yYu+∂h/∂xXuZv)-Zu(∂h/∂yYv +∂h/∂xXv) +∂h/∂z(ZuZv-ZvZu) ]dudv[br]=∫[sub]S0[/sub][ (Xv(∂f/∂x Xu +∂f/∂yYu+∂f/∂z Zu) -Xu(∂f/∂x Xv+∂f/∂yYv+∂f/∂zZv) [br] +Yv(∂g/∂xXu+∂g/∂yYu+∂g/∂zZu)-Yu(∂g/∂xXv +∂g/∂yYv+∂g/∂zZv)[br] +Zv(∂h/∂xXuZv+∂h/∂yYu+∂h/∂zZu)-Zu(∂h/∂xXv+∂h/∂yYv +∂h/∂zZv) ]dudv[br][br]=∫[sub]S0[/sub][ (∂f/∂u Xv -∂f/∂v Xu) +(∂g/∂u Yv -∂g/∂v Yu)+(∂h/∂u Zv-∂h/∂v Zu) ]dud∂h/∂vv[br]=∫[sub]S0[/sub][ (∂/∂u (f Xv) -∂/∂v (f Xu)) +(∂/∂u(g Yv) -∂/∂v(g Yu))+(∂/∂u(h Zv)-∂/∂v(h Zu)) ]dud∂h/∂vv[br]=∫[sub]C0[/sub][ (f Xu)du +(f Xv) dv + (g Yu) du + (g Yv)dv +(h Zu)du +(h Zv)dv ]dudv グリーンの定理から[br]=∫[sub]C0[/sub][ (f Xu)du +(f Xv) dv] dudv +∫[sub]C0[/sub][ (g Yu) du + (g Yv)dv] dudv +∫[sub]C0[/sub][ (h Zu)du +(h Zv)dv ]dudv[br]=∫[sub]C [/sub](f dx + g dy + h dz) 全微分から[br]= ∫[sub]C [/sub][b]f[/b]・d[b]p[br][/b](例)[br]アンペールの法則[br]閉曲面Sで単位面積を流れる電流ベクトル(電流密度)i の面積分と[br]境界線Cでの磁場の強さHの線積分は等しい。[br][color=#0000ff] ∫[sub]S[/sub] [b]i [/b](x,y,z)・ d[b]S[/b] =∫[sub]c[/sub] [b]H [/b]・d[b]s[/b] 積分形。[/color][br] rot [b]H = i[/b] (x,y,z) 微分形[br](例)[br]z軸を中心とする半径aの円柱に電流が流れていて、電流密度が i0, 電流からrの距離の磁場の強さH(r)[br]半径rの境界線Cでの磁場の強さの線積分は、Cの内部の閉曲面の電流ベクトルの面積分に等しい。[br](円柱の外部)[br]半径aと半径rの間の電流は0なので、半径aの円Aで面積分すればよい。[br]面積分∫[sub]S[/sub] [b]i [/b](x,y,z)・ d[b]S[/b] =∫[sub]A[/sub] i0 dS=i0 ∫[sub]A[/sub] 1dS =i0 πa[sup]2[/sup] [br]線積分∫[sub]c[/sub] H(r)ds=H(r)∫[sub]c[/sub] 1ds = H(r)2πr [br]この2つが等しいので、 H(r)=i0 πa[sup]2[/sup] /2πr=[b]i0 a[sup]2[/sup]/2r[/b] [br](円柱の内部)[br]半径r以下はすべて電流が流れているので、半径rの円Rで面積分すればよい。[br]面積分∫[sub]S[/sub] [b]i [/b](x,y,z)・ d[b]S[/b] =∫[sub]R[/sub] i0 dS=i0 ∫[sub]R[/sub] 1dS =i0 πr[sup]2[/sup] [br]線積分∫[sub]c[/sub] H(r)ds=H(r)∫[sub]c[/sub] 1ds = H(r)2πr [br]この2つが等しいので、 H(r)=i0 πr[sup]2[/sup] /2πr=[b]i0 r/2[/b]