[size=150][b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][/b][/size][size=150][size=100][br]前回まで「[b]慣性系[/b]」の相対性原理、つまり[b]特殊[/b]な相対性理論を見てきたね。[br]ポイントは光速が絶対で一定であることだった。[br][br]相対性理論を計量まで深めるためには、時空のゆがみや変換を学ぶ必要があることに気づいたね。[br]そこで、時空を数式で表現するために、今回はテンソルについて学ぼう。[/size][/size]

[b]＜テンソルは複数添え字（方向）をもつ数の入れ物＞[br][/b]テンソルはデータが入れ物の型でベクトル・行列の1種ともいえます。[br]配列やリストが入れにすると作れます。[br]テンソルは自由度からみると、スカラー、ベクトルの次にくるもので、行列を広げたものです。[br]実際に、[b]方向(軸）の数、階数（ランク）[/b]は[br]スカラーは０、ベクトルは１、テンソルは2以上です。[br]テンソルは、広い意味ではスカラー、ベクトルも含み、狭い意味では[b]行列以上のものに限ります[/b]。[br]ベクトルと同様に、成分は座標系に連動して変わるけれどもテンソルそのものは変わらない。[br][size=150][color=#0000ff]・テンソルの型は[b]shape[/b]で表示できます。[br][/color][color=#0000ff]・テンソルの軸は[b]transpose[/b]で,行列と同じように転置できます。[br][/color][/size][br][color=#0000ff](例）[/color][br]Python[br]#テンソルの型============[br]import numpy as np[br]scalar_x=np.array(2)[br]vector_x=np.array([1,2]) [br]matrix_x=np.array([[1,2], [3,4]]) [br]tensor_x=np.array([[[1,2], [3,4]], [[2,4], [6,8]]])[br][b]print(scalar_x.shape,",",vector_x.shape,",",matrix_x.shape,",",tensor_x.shape)[br][/b]#================[br][OUT}[br][b]() , (2,) , (2, 2) , (2, 2, 2)[br][/b][br][color=#0000ff](例）[br][/color]1から24までの数の配列PをP.shape=(2,3,4)とすることで、[br]2×3×4のテンソルに整形しましょう。2段に区切り、3本に区切り、1本に4個ならべます。[br]順にz軸、y軸、x軸とすることで、ループをｚ、ｙを固定してｘだけ動かすと[br]ｘ軸にそって、１，２，３，４と数がならんだベクトルが1本できます。[br]それをｙ軸を動かしてから、ｘ軸にそって、５，６，７，８と数がならんだベクトルが1本できます。[br]ｚを動かすと、上の段にも数のベクトルを2本ならべられます。[br][b]・Q=P.transpose(0,2,1)とすることで、ｚ段数は変えずに行と列を転置できますね。[br]・Q=P.transpose(2,1,0)なら、ｙ[/b][b]たて数を変えずに、ｚ段とｘよこの入れ替えができます。[br][/b][br]Python[br]#テンソルの転置============[br]import matplotlib.pyplot as plt[br]import numpy as np[br]P=np.array(range(1,25))[br]P.shape=(2,3,4)[br]#段数をかえずに、たて・よこの入れ替え[br]Q=P.transpose(0,2,1)[br]Q[br]fig, ax = plt.subplots(1, 2, figsize=(10, 4), subplot_kw={'projection': '3d'})[br][br]def drawTensor(P,n):[br] print(P)[br] z_size, y_size, x_size = P.shape[0], P.shape[1],P.shape[2][br] voxels = np.ones((x_size, y_size, z_size), dtype=bool)[br] ax[n].voxels(voxels, edgecolor='k', facecolors='green', alpha=0.1)[br] #テンソルの内容をその位置にかく[br] for pz in range(z_size):[br] for py in range(y_size):[br] for px in range(x_size):[br] contents = P[pz,py,px][br] #print(f"{pz},{py},{px};contents")[br] ax[n].text(px, py, pz+0.5, f'{contents}', fontsize=15, verticalalignment="bottom")[br] ax[n].set_xlabel('x')[br] ax[n].set_ylabel('y')[br] ax[n].set_zlabel('z')[br] ax[n].set_aspect('equal')[br]drawTensor(P,0)[br]drawTensor(Q,1)[br][br]plt.show()[br]#====================================[br][OUT][br][[[ 1 2 3 4][br] [ 5 6 7 8][br] [ 9 10 11 12]][br][br] [[13 14 15 16][br] [17 18 19 20][br] [21 22 23 24]]][br][[[ 1 5 9][br] [ 2 6 10][br] [ 3 7 11][br] [ 4 8 12]][br][br] [[13 17 21][br] [14 18 22][br] [15 19 23][br] [16 20 24]]][br]

[b]＜成分ごとの和差積商もある＞[br][/b]テンソルはデータが入れ物の型でベクトル・行列の1種だから、[br]同じ位置のどうしのたし算、引き算、かけ算、わり算というものもできます。[br]A+B[br]A－B[br]A*B[br]A/B[br]これはベクトルの和と差、行列の和と差といっしょで、[br]同じ位置の成分どうしの演算なので、[br]テンソルの形状が同じものどうしの計算になり、[br]とくにテンソルらしい複雑さは何もありません。[br]なので、省略します。[br][br][b]＜テンソル積＞[/b][br]テンソル積は、２つのテンソルで共通する方向（軸）を決めてその内積をとるものです。[br][b]共通軸を指定しないとき[/b]、[br]次元が q のテンソル A と[br]次元が r のテンソル B があるとき、[br]これらのテンソルの積[b][size=150]C = A⊗B[/size][/b]は、 q + r 次元のテンソルになります。[br][b][br]・共通の方向（軸）の数を０とする[/b]([b]axes=0[/b])と、２つのテンソルの各要素に対する[b]直積[/b]になるので、[br]次元が和になります。[color=#0000ff]１次元と１次元のテンソル積は1＋1＝２次元。[br](例）[/color][br]Python[br][IN][br]#============[br]from numpy import array[br]from numpy import tensordot[br]A = array([1,2])[br]B = array([3,4])[br][color=#0000ff]C = [b]tensordot[/b](A, B, axes=0)[br][/color]print(C)[br]#=================[br][OUT][br][[3 4][br] [6 8]][br][br][b]・共通の方向（軸）の数を１とする[/b]([b]axes=1[/b])と、２つのテンソルの共通な軸で[b]内積[/b]がおきるので、[br][size=200][color=#0000ff]その軸が[b]つぶれます（縮約）[/b][/color][/size]。[br]軸数が１つかぶってそれぞれ、減ります。[br][color=#0000ff]２次元と２次元のテンソル積は行列の積と同じ2-1+2-1=2次元になります。[br][/color]どの軸をつぶして内積をとるかをそれぞれのテンソルの添え字の番号で指定します。[br][color=#0000ff](例）[/color][br]Python[br][IN][br]#============[br]import numpy as np[br][b]# (2, 3) と (3, 4) の行列。共通の「3」の部分（軸1と軸0）で積をとる[br][/b]A = np.array(range(1,7)) #１から６の数列[br]A.shape=(2,3)#2行3列の2軸[br]B = np.ones((3, 4))#１だけを3行4列の2軸[br][b][color=#0000ff]res = np.tensordot(A, B, axes=(1, 0)) [br]resre = np.tensordot(B,A,axes=(0,1))[br][/color][/b]# 結果は (2, 4) の行列。np.dot(A, B) と同じ。[br]print(A)[br]print(B)[br]print("np.tensordot(A,B,axes=(1,0))=")[br]print(res)[br]print(resre)[br]#=============[br][OUT][br][[1 2 3][br] [4 5 6]][br][[1. 1. 1. 1.][br] [1. 1. 1. 1.][br] [1. 1. 1. 1.]][br]np.tensordot(A,B,axes=(1,0))=[br][[ 6. 6. 6. 6.][br] [15. 15. 15. 15.]][br][[ 6. 15.][br] [ 6. 15.][br] [ 6. 15.][br] [ 6. 15.]][br][br]

[b]＜2つ以上の軸を同時につぶせるテンソル積＞[br][br]・AとBの２つの軸で、大きさが同じなら、2軸一緒に[b]（axes=2)[/b]つぶせます。Aの後ろ２つの軸と、Bの前２つの軸をつぶすのが、[b][color=#0000ff]np.tensordot(A, B, axes=2)です。[br][b][b][color=#0000ff]４次元と３次元で2軸をつぶす（縮約）するテンソル積は4-2+3-2=3次元になります。[br][/color][/b][/b]３次元と２次元で2軸をつぶす（縮約）するテンソル積は3-2+2-2=1次元になります。[br][/color][/b]つぶす軸を添え字の番号で指定することもできます。[br](例）[br][/b]# A: 4次元テンソル (2, 3, [b]4, 5[/b])[br]# B: 3次元テンソル ([b]4, 5,[/b] 2)[br]A = np.random.rand(2, 3, 4, 5)[br]B = np.random.rand(4, 5, 2)[br]# Aの「軸2, 3」(4, 5) と Bの「軸0, 1」(4, 5) を一気にぶつけて消す[br][b][color=#0000ff]res = np.tensordot(A, B, axes=([2, 3], [0, 1])) #axesのあとは添え字の位置の指定[br][/color]# 結果のテンソルの型は (2, 3, 2) と3次元になります。[br][/b]print(A)[br]print(B)[br]print("np.tensordot(A,B,axes=(1,0))=")[br]print(res)[br][OUT][br]略[br][br][b][color=#9900ff][u][size=150]課題：つぶれる前とつぶれた後のテンソルを視覚化するにはどうしたらよいでしょうか。[br][/size][/u][/color][/b]たとえば、[br]aが３軸 (2, 2, 2)のテンソルで成分が、[[[1 2], [3 4]], [[5 6], [7 8]]]だとします。[br]Aが２軸 (2, 2)のテンソルで成分が、[[1000 100], [10 1]]とするとき、[br][b][color=#0000ff]軸ごとにmatplotlib.pyplotの箱(voxel)を用意する[/color][/b]と、[br]テンソルaは3軸で、2×2×2＝８個の箱（ボクセル）でできた立方体になりますね。[br]テンソルAは2軸で、2×2＝4個の箱（ボクセル）が1段並びます。[br]これをツールとして使うことで、テンソル積の軸の選び方で、結果がどう変わるかを確かめましょう。[br][br]import numpy as np[br]import matplotlib.pyplot as plt[br]import numpy as np[br][br]#tensordotの視覚化。axes=1[br]a = np.array(range(1, 9))[br]a.shape = (2, 2, 2)[br]A = np.array((1000,100,10,1), dtype=object)[br]A.shape = (2, 2)[br]print("--- a ---")[br]print(f"形: {a.shape}") [br]print(a)[br]print(a[:,-1])[br]print("\n--- A ---")[br]print(f"形: {A.shape}")[br]print(A)[br]print(A[:,0])[br]print("\n--- res ---")[br][br]fig, ax = plt.subplots(1, 2, figsize=(11, 4), subplot_kw={'projection': '3d'})[br][br]def [b]drawTensor3D[/b](P,n):[br] z_size, y_size, x_size = P.shape[0], P.shape[1],P.shape[2][br] voxels = np.ones((x_size, y_size, z_size), dtype=bool)[br] ax[n].voxels(voxels, edgecolor='k', facecolors='green', alpha=0.1)[br] #テンソルの内容をその位置にかく[br] for pz in range(z_size):[br] for py in range(y_size):[br] for px in range(x_size):[br] contents = P[pz,py,px][br] ax[n].text(px, py, pz+0.5, f'{contents}', fontsize=15, verticalalignment="bottom")[br] ax[n].set_xlim([0, 2]);ax[n].set_ylim([0, 2]);ax[n].set_zlim([0, 2])[br] ax[n].set_xlabel('x');ax[n].set_ylabel('y');ax[n].set_zlabel('z')[br] ax[n].set_aspect('equal')[br]def [b]drawTensor2D[/b](P,n):[br] y_size, x_size = P.shape[0],P.shape[1][br] #テンソルの内容をその位置にかく[br] for py in range(y_size):[br] for px in range(x_size):[br] contents = P[py,px][br] ax[n].text(px, py, 0.5, f'{contents}', fontsize=15, verticalalignment="top")[br] ax[n].set_xlim([0, 2]);ax[n].set_ylim([0, 2]);ax[n].set_zlim([0, 2])[br] ax[n].set_xlabel('x');ax[n].set_ylabel('y');ax[n].set_zlabel('z')[br] ax[n].set_aspect('equal')[br]def [b]drawTensor1D[/b](P,n):[br] x_size = P.shape[0][br] #テンソルの内容をその位置にかく[br] for px in range(x_size):[br] contents = P[px][br] ax[n].text(px, 0, 0, f'{contents}', fontsize=15, verticalalignment="bottom")[br] ax[n].set_xlim([0, 2]);ax[n].set_ylim([0, 2]);ax[n].set_zlim([0, 2])[br] ax[n].set_xlabel('x');ax[n].set_ylabel('y');ax[n].set_zlabel('z')[br] ax[n].set_aspect('equal')[br][br][b][color=#0000ff]（例１）tensordot(a, A,axes=2)とすると、[br][/color][/b]aのｚ軸が残り、aのxy軸とAのxy軸がピッタリと重なってつぶれることで、2個の数値になります。[br]aの[[1 2], [3 4]]がAの[[1000 100], [10 1]]に要素ごとの内積になり[br]1×1000＋２×100＋3×10＋４×１＝1234ができます。段ｚを１にしてaの2段目から5678ができます。[br][b]#============ ここからを変える。================[br]res = np.tensordot(a, A,axes=2) #aの末２軸とAの前２軸をつぶす。[br]print(f"形: {res.shape}")[br]print(res)[br][br]drawTensor3D(a,0)[br]drawTensor2D(A,0)[br]drawTensor1D(res,1)[br]plt.show()[br][/b]#===========================================[br]--- a ---[br]形: (2, 2, 2)[br][[[1 2][br] [3 4]][br][br] [[5 6][br] [7 8]]][br][[3 4][br] [7 8]][br][br]--- A ---[br]形: (2, 2)[br][[1000 100][br] [10 1]][br][1000 10][br][br]--- res ---[br]形: (2,)[br][1234 5678]

[color=#0000ff][b]（例2）einsum('ijk,lk->ijl',a,A)[/b][br][/color]einsumはtensordotの添え字番号をijkのように文字にしたもです。[br]ijkはaの添え字で、lkはAの添え字で、共通の軸はkです。[br]テンソルaでのkはz軸,y軸,x軸の順のxなので、x軸方法のベクトルとします。[br]テンソルaは４本のベクトル(1,2),(3,4),(5,6),(7,8)がzｙ平面に立ってるイメージです。[br]テンソルAでのkはy軸,x軸の順のｘなので、ｘ軸方向の2本のベクトル(1000,100),(10,1)が寝てます。[br]これらの直積を内積によりつぶすと、１×1000＋２×100＝1200、1×10+2×1＝12、[br]3400、34、5600、56、7800、78ができますが。これらを2×2×2の箱に入れます。[br]np.tensordot(a, A, (2, 1))でも同じですね。[br][b]#============ ここからを変える。================[br][/b]res = np.einsum('ijk,lk->ijl',a,A)[br]print(f"形: {res.shape}")[br]print(res)[br]fig, ax = plt.subplots(1, 2, figsize=(11, 4), subplot_kw={'projection': '3d'})[br][br]drawTensor3D(a,0)[br]drawTensor2D(A,0)[br]drawTensor1D(res,1)[br]plt.show()[br]#=========================================[br][OUT][br]--- a ---[br]形: (2, 2, 2)[br][[[1 2][br] [3 4]][br][br] [[5 6][br] [7 8]]][br][[3 4][br] [7 8]][br][br]--- A ---[br]形: (2, 2)[br][[1000 100][br] [10 1]][br][1000 10][br][br]--- res ---[br]形: (2, 2, 2)[br][[[1200 12][br] [3400 34]][br][br] [[5600 56][br] [7800 78]]][br]

[color=#0000ff][b]（例3）einsum('ijk,li->jkl',a,A)[/b][br][/color]ijkはaの添え字で、liはAの添え字で、共通の軸はiです。[br]テンソルaでのiはz軸,y軸,x軸の順のzなので、z軸方法のベクトルとします。[br]テンソルaは４本のベクトル(1,5),(2,6),(3,7),(4,8)がxｙ平面に立っているイメージです。[br]テンソルAでのiはy軸,x軸の順のｘなので、ｘ軸方向の2本のベクトル(1000,100),(10,1)が寝てます。[br]これらの直積を内積によりつぶすと、１×1000＋5×100＝1500、1×10+5×1＝15、[br]3700、37、2600、26、4800、48ができますが。これらを2×2×2の箱に入れます。[br]np.tensordot(a, A, (0, 1))でも同じですね。[br][b]#============ ここからを変える。================[br][/b]res = np.einsum('ijk,li->jkl',a,A)[br]print("\n--- res ---")[br]print(f"形: {res.shape}")[br]print(res)[br][br]drawTensor3D(a,0)[br]drawTensor2D(A,0)[br]drawTensor3D(res,1)[br]plt.show()[br]#=========================================[br][OUT][br]--- res ---[br]形: (2, 2, 2)[br][[[1500 15][br] [2600 26]][br][br] [[3700 37][br] [4800 48]]][br]

[color=#0000ff][b]（例４）einsum('ijk,lk->ijl',a,A)[/b][br][/color]ijkはaの添え字で、liはAの添え字で、共通の軸はkです。[br]テンソルaでのiはz軸,y軸,x軸の順のzなので、z軸方法のベクトルとします。[br]テンソルaは４本のベクトル(1,5),(2,6),(3,7),(4,8)がxｙ平面に立っているイメージです。[br]テンソルAでのiはy軸,x軸の順のｘなので、ｘ軸方向の2本のベクトル(1000,100),(10,1)が寝てます。[br]これらの直積を内積によりつぶすと、１×1000＋5×100＝1500、1×10+5×1＝15、[br]3700、37、2600、26、4800、48ができますが。これらを2×2×2の箱に入れます。[br]np.tensordot(a, A, (0, 1))でも同じですね。[br][b]#============ ここからを変える。================[br][/b]res = np.einsum('ijk,lk->ijl',a,A)[br]print("\n--- res ---")[br]print(f"形: {res.shape}")[br]print(res)[br][br]drawTensor3D(a,0)[br]drawTensor2D(A,0)[br]drawTensor3D(res,1)[br]plt.show()[br]#=========================================[br][OUT][br]--- res ---[br]形: (2, 2, 2)[br][[[1200 12][br] [3400 34]][br][br] [[5600 56][br] [7800 78]]]

[b][size=150]＜テンソルでスカラー場、ベクトル場の視覚化をしよう＞[size=100]テンソルは軸の数が1以上であれば、数値データを自由に入れられた。[br]だから、座標データを入れたテンソルを作り、その点ごとのスカラー場、ベクトル場を[br]計算して、空間のようすを視覚化することができるはずだ。[br][/size][color=#9900ff][br]課題：ストークスの定理を視覚化するにはどうしますか。[br][/color][/size][/b][br]ガウスの定理、グリーンの定理、ストークスの定理、どれも次元下げの定理です。[br]だから、空間を扱うときに、空間全部ではなく、表面のようすを調べると、その中は打ち消されているので、全体のようすがわかるというものです。[br]部分部分をたし込むときに、次元を下げて積分してみましょう。[br][b][color=#0000ff]（例）[/color][/b][br]from sympy import symbols, cos, sin, Matrix,diff,integrate,pi,lambdify[br]import numpy as np [br]import matplotlib.pyplot as plt[br]#from sympy.plotting import plot3d_parametric_line,plot3d_parametric_surface[br][br][color=#38761d]# グリーンの定理を３Dにしたものがストークスの定理。曲面のふちだけで積分しても同じ量になる。[br]# 境界線Cの曲面Sでのベクトル場Fについて、曲面SでのFの回転の面積分G＝Fの各成分の境界線Cで線積分の和[br]# F=[F1,F2,F3] の回転の面積分。ｎは単位法線ベクトル,Rは曲面Sの領域,Cは領域Rの境界線[br]# G=∫∫R rotF・ｎ dA = ∫C(F1dx + F2dy + F3dz)=∫C F・dr/ds ds=∫C F・dr [br]# Sは 放物面 z=f(x,y)=1-(x^2+y^2)(z≧0)で、境界線Cは円周（x^2+y^2=1）F=[y,z,x][br][/color]def grad(w,vars): #スカラー場ｗの勾配はベクトル[br] return Matrix([diff(w, i) for i in vars])[br]def div(v,vars): #ベクトル場ｖの発散はスカラー[br] return sum([diff(v[i],vars[i]) for i in range(len(vars))])[br]def rot(F, vars):#回転はcurlとすることが多い。日本ではrot[br] F0,F1,F2 = F[br] return Matrix([diff(F2,vars[1]) - diff(F1,vars[2]),diff(F0,vars[2]) - diff(F2,vars[0]),diff(F1,vars[0]) - diff(F0,vars[1])])[br][br]f,xx,yy,zz,x,y,z,r,t,u,v,a,b= symbols('f xx yy zz x y z r t u v a b')[br]#領域Vで3重積分[br]F1= y;F2= z;F3= x[br][b]F = Matrix([F1,F2,F3])[br][/b]xx= cos(t)[br]yy= sin(t)[br]zz= 0[br][b]drdt = Matrix([diff(xx,t),diff(yy,t),diff(zz,t)])[br][/b]Ft=F.subs({x:xx,y:yy,z:zz})[br][color=#0000ff]Ftdotdr=Ft.[b]dot[/b](drdt)[br][/color]G=integrate(Ftdotdr,(t,0,2*pi))[br]print("線積分 ∫ F・(dr/dt) dt =",G)[br][br]Z=1-(x**2+y**2)[br][b]N=grad(z-(Z),[x,y,z])[br]xxr= r*cos(t)[br]yyr= r*sin(t)[br]zzr= 0[br][color=#0000ff]rotF = rot(F,[x,y,z])[br][/color][/b]FN=rotF.dot(N).subs({x:xxr,y:yyr,z:zzr})[br][b]FNr=FN*r[br][/b]G2=integrate(integrate(FNr,(r,0,1)),(t,0,2*pi))[br]print("面積分 ∫∫ rotF・N dxdy =∫∫ rotF・N(r,t) rdrdt =",G)[br][br][color=#45818e]#視覚化[br][/color]def [b]drawVecAt[/b](u0,v0):[br] #接点と法線の方向[br] x0, y0, z0 = np.array([b]func_pos[/b](u0, v0), dtype=float).flatten()[br] nx, ny, nz = np.array([b]func_normal[/b](x0, y0), dtype=float).flatten()[br] fx, fy, fz = np.array([b]func_F[/b](x0, y0, z0), dtype=float).flatten()[br] [br][color=#0000ff] ax.scatter(x0, y0, z0, color= "green", s=10)[br] #法線ベクトル[br] ax.quiver(x0, y0, z0, nx/10, ny/10, nz/10, length=1.0, color="green", linewidth=1) [br] #Fベクトル[br] ax.quiver(x0, y0, z0, fx, fy, fz, length=1.0, color="blue", linewidth=1) [br][/color][b]XXX= r*cos(t)[br]YYY= r*sin(t)[br]ZZZ= 1-r**2[br][/b]N=grad(z-(Z),[x,y,z])[br][br]fig = plt.figure(figsize=(10, 8))[br]ax = fig.add_subplot(111, projection='3d')[br]t_vals = np.linspace(0, 2*np.pi , 120);[br]r_vals = np.linspace(0, 1 , 120); [br][b]T, R = np.meshgrid(t_vals, r_vals)[color=#0000ff]#２パラメータのテンソル化[/color][br]func_pos [/b]= lambdify((t, r), [b][XXX,YYY,ZZZ][/b], "numpy") [br][b]func_normal [/b]= lambdify((x, y), [b]N,[/b] "numpy")[br][b]func_F [/b]= lambdify((x, y, z), [b]F[/b], "numpy")[br]X, Y, Z = func_pos(T, R)[br]ax.plot_surface(X, Y, Z, alpha=0.3, cmap='viridis', edgecolor='blue')[br][br]pT = np.random.uniform(0, 2*np.pi, size=(150, 1))[br]pR = np.random.uniform(0, 1.0 , size=(150, 1))[br]for i in range(len(pT)):[br] [b]drawVecAt[/b](pT[i].flatten(),pR[i].flatten()) [br]ax.set_xlim([-1.2, 1.2])[br]ax.set_ylim([-1.2, 1.2])[br]ax.set_zlim([0, 1])[br]ax.set_xlabel('x'); [br]ax.set_ylabel('y'); [br]ax.set_zlabel('z')[br]ax.set_aspect('equal')[br]ax.set_title("STOKES’S THEOREM")[br]plt.show()[br][br]#=====================================[br][OUT]線積分 ∫ F・(dr/dt) dt = -pi[br]面積分 ∫∫ rotF・N dxdy =∫∫ rotF・N(r,t) rdrdt = -pi

[u][color=#9900ff][b][size=150]課題：Geogebraで回転を計算して視覚化するにはどうしますか。[br][/size][/b][/color][/u][br]geogebraの数式に、1行ずつx、y、zの関数として入力します。[br]vx=x^2*y*z[br]vy=x*y^2*z[br]vz=x*y*z^2[br]rotV={Derivative(vz,y)-Derivative(vy,z),Derivative(vx,z)-Derivative(vz,x),Derivative(vy,x)-Derivative(vx,y)}[br]が回転です。[br][b][color=#0000ff]it=Sequence(t,t,-20,20,10)で、数列を作ります。これを3次元の直積にするために、[br]st=Flatten(Zip(Zip(Zip((p,q,r),p,it),q,it),r,it))とかきます。p,q,rが独立して動くので直積となり、[br]テンソルができます。３次元にしたので、stは3軸のテンソルになりますね。[br]a=Length(st)によって、テンソルを配列につぶしたときのサイズをaに入れときます。[br][br][/color][/b]このテンソルstは(x,y,z)の集まりともいえますが、[br]この[b][color=#0000ff][u]要素のi番目のｘ座標にアクセスするためにx(st(i))とかきます[/u][/color][/b]。[br]だから、stのi番目の点の座標は(x(st(i)),y(st(i)),z(st(i)))です。この点をPとして、[br]PからrotVの成分だけ加算した点Qまでつなぐと、ベクトルになりますね。[br]rotVのｘ成分、ｙ成分、ｚ成分は[br]Derivative(vz,y)-Derivative(vy,z)、[br]Derivative(vx,z)-Derivative(vz,x)、[br]Derivative(vy,x)-Derivative(vx,y)なので、[br]Qのｘ座標,y座標,z座標は[br]b(x,y,z)=x+(Derivative(vz,y)-Derivative(vy,z))/300[br]c(x,y,z)=y+(Derivative(vx,z)-Derivative(vz,x))/300[br]d(x,y,z)=z+(Derivative(vy,x)-Derivative(vx,y))/300[br]とおくと、x(st(i)),y(st(i)),z(st(i))をメッシュのテンソルと考えて、ベクトルの配列を作れますね。[br]それは、[br]Sequence(Vector((x(st(cnt)),y(st(cnt)),z(st(cnt))),(b(x(st(cnt)),y(st(cnt)),z(st(cnt))),c(x(st(cnt)),y(st(cnt)),z(st(cnt))),d(x(st(cnt)),y(st(cnt)),z(st(cnt))))),cnt,1,a)[br]です。[br]とても長ったらしいですが、構造は[br]Sequence(Vector(P,Q),cnt,1,a)です。[br]点Pをcntの関数として、(x(st(cnt)),y(st(cnt)),z(st(cnt)))とかいて動かし、[br]点Qもcntの関数(b(x(st(cnt)),y(st(cnt)),z(st(cnt))),c(x(st(cnt)),y(st(cnt)),z(st(cnt))),d(x(st(cnt)),y(st(cnt)),z(st(cnt))))[br]を動かしています。[br]できなくはないですが、大仕掛けですね。[br]300で発散の成分を割っているのは大きすぎてはみ出して、動きがわからなくなるからです。[br][color=#ff0000]※保存する前までは、表示されますが、構文の複雑さの限界を超えてしまうようです。[br]保存すると、座標メッシュとしての空間点テンソルstは表示されますが、[br]矢印のシークエンスは保存すると、Undefinedの？になってしまいます。[br]保存する前までの短い命です。[br]構文チェックの構造がgeogebraのコンパイラとインタプリタ-でちがっているのが原因かもしれません。[/color]