コンテンツにスキップ

Topic:微分積分学/研究

提供:ウィキバーシティ

変分法について

[編集]
本論文の主な著者:--すじにくシチュー (トーク) 2014年9月14日 (日) 08:15 (UTC) (日付は署名時の日付)

以下の理論は、変分法を意識した仮説(かも?)である。

曲線の集団は曲面と見なすべきかも?

[編集]

 たとえば、曲線族 (曲線の集まりのことを曲線族という)

y(x) = Cx (Cは任意定数) ・・・①

のような一変数の一次関数のあつまりをかんがえる。①式を説明すると、たとえば任意定数とは、定数関数なら、C=1でもC=2でも、C=125.35でも何でもいいので、①式は

{ y = 0x 、 y=1x 、 y=1.3x 、 y=-42x 、・・・ }

のような、一変数の一次関数の集まりと考えることができる。一変数の一次関数は、常微分が可能なので、微分すると

dy/dx = C = y´(x) 

である。Cが任意の値を取れるので、y´もCと同じ値であれば、任意の定数となる。

C = y´ より、
y = Cx = y´x 

である。 一方、y = Cx はyの値が、Cの値とxの値によって決まるので、yはCとxを変数とする二変数関数 y(x、C) とも考えられる。二変数関数はグラフにあらわすと曲面になるので、つまり、 「曲線の集合は曲面を表せる」 と考えるべきである。 あるいは 「一変数関数の集合は、二変数関数として表現できる」 とかんがえるべきである。y(x、C)は

 y = Cx = y´x  より、 y(x、y´)という二変数関数である。Cの値が、xとは独立に任意に取れて、Cを独立変数とみなしたのだから、C = y´ より、
y´(x)もxとは独立な変数とみなさなければならない。つまり、「関数の集合では、xと、xを変数とする導関数は、お互いに独立である。」といえる。または
y = y´x ⇔  y´= x/y 

より、y´(x、y)という二変数関数を考えることもできる。 このように、曲線の集合を扱う計算では、1本の曲線を扱う計算とは違った現象が起こる。曲線の集合は曲面とみなせるので、偏微分も定義できる。 たとえば、y = y´x をy´で偏微分すると、xとなる。 以上の議論は、一時間数の集合に限らず、2次関数の集合 y=c・x^2 、 三次間数の集合でも同類の議論ができて、一変数の関数の集合を二変数関数とみなすこともできる。 また、任意定数が2つ以上の場合も、同様の議論によって、一変数関数の集合を、他変数関数とみなせる。

変分法を意識して、作用素を仮定

[編集]

グラフ上で関数 y(x)の形を動かして、別の形の関数を作ることを考える。 形を変えて新しく作った関数を Ny とかくことにする。このとき、関数の形を変える作用素を N とおく。  たとえば  y=2x を変形させて y=x-5 という関数を作った場合は、

Ny=x-5 と書くする。

この場合の変動量は Ny-y である。これを

(N-1)y

書くとしよう。

y=2x から y=x-5 への変形の場合、変動量の計算には、代入をすればいいだけである。

Ny-y=(x-5)-(2x)=-3x+5

なので、つまりyを -3x+5 だけ動かせば、Nyをつくれる。この関数を動かした量、Ny-y を δyとおくことにする。

つまり、定義より 

δ=N-1

である。

いっぽう、 y=2xを微分すれば、y´=2   Ny を微分すれば、この問題でのNy=x-5 の条件の場合、

このとき、導関数は

(Ny)´ -y´=-1-2=-3

よって導関数はy´から-3だけ動いたことになる。  このとき、D(Ny)とは、

y=2x に N を作用させ、変化させた新関数 Ny =-x+5 を(x、y)平面上で作ってから、

Dで微分して、 微分されたD(Ny)=(-x+5)´= -1を、(x、y´)上に持ってきたものである。

ここで、仮に、 

DNy=NDy 
にするためには、どう設定すればよいか。

NDyの意味を考えると、これは、 微分と新関数作成の順序を変えてみて、 yを微分したDy=y´=(2x)´=2 を(x、y´)グラフ上に持ってきて、 この(x、y´)上の関数 Dy に Nを作用させて、

N(Dy)=-1 を得たものとも見ることができる。

つまり、関数yの導関数y´を変化させるとき、yを変化させずに y´ だけ動かすことは、普通はできないから、(x、y´)上の変化は、(x、y)平面上の変化と

DN=ND の関係で結ばれるとすればよいのである。


ともかく、作用素Nと微分作用素Dが、可換だとする。つまり、DN=NDだとする。 このとき δ=N-1 と D も可換である。実際に

D(δy) = D(N-1)y = DNy - Dy = NDy - Dy 
=(N-1)Dy = δDy

よって

 Dδ = δD
(説明、おわり)

(・・・という事らしい。昔、自分が書いた自主レポート(未発表)を元にしてるので、詳しいアイデアは忘れた。)

テイラー展開を元に、組み合わせ公式

[編集]
本論文の主な著者:--すじにくシチュー (トーク) 2014年9月14日 (日) 08:15 (UTC) (日付は署名時の日付)

たとえば4個のものの中から、2個を選ぶ組み合わせの場合を

C(4、2)

とかくことにする。

この記法で組み合わせ公式を用いるとき、

C(j-1、i-1)・(a+b)^(j-i)
=Σ(i≦k≦j){C(k-1、i-1)a^(k-i)}{C(j-1、k-1)b^(j-k)}

が成り立ちそうである。

  • 説明

 まず、テイラー展開の普通の関数を考える。 たとえば、

Y=C0+C1・x+C2・x^2

をベクトル (C0.C1,C2)であらわす。 つまり、べき級数の係数をベクトルの成分にする。 このときYにたいする微分Dや差分⊿、代入、シフトEが、行列で表せる。 シフトEとは、

E・Y(x) =Y(x+1)

のことである。 差分⊿は

⊿Y(x)=Y(x+1)-Y(x)=EY(x)-Y(x)=(E-1)・Y(x)

である。 さて、シフト E の行列 M(i,j) は、

M(i,j)=C(j-1、i-1) ・・・(i≦j のとき)

その他

M(i,j)=0 ・・・(その他のとき)

となる。 これは、パスカルの3角形と同じ配置である。

 G(A)[Y(x)]=Y(x+A)

という演算子G(A)を考えると、このG(A)をあらわす行列 M(i、j) は、

M(i,j)=C(j-1、i-1)・A^(j-i)・・(i≦j のとき)

その他

M(i,j)=0 ・・・(その他のとき)

となる。

A=1のときは、G(1)はシフトEと、同一になる。

さて、G(A+B)・Y(x) = G(A)・G(B)・Y(x) なので、この関係を、行列で表したのが、冒頭に掲げた組み合わせの公式である。

テイラー展開からテンソルに

[編集]
本論文の主な著者:--すじにくシチュー (トーク) 2014年9月14日 (日) 08:15 (UTC) (日付は署名時の日付)

(既存研究が、あるかも?)

  • テンソル積

 関数の係数をベクトルと考えると、

f(x) = A0 + A1・x + A2・x^2 + A3・x^3 + ・・・
g(y) = B0 + B1・y + B2・y^2 + B3・y^3 + ・・・

の積は、

f(x)・g(y) = A0・B0 + A0・B1・x + A0・B2・y^2 + ・・・
 +A1・B0・x + A1・B1・x・y +A1・B2・x・y^2 + ・・・
         +A2・B0・x^2 + A2・B1・x^2・y + ・・・

・・・・・ のようになる。 ここで、たとえば一変数関数f(x)の係数をベクトルの成分と見たように、2変数関数f(x)・g(y)の係数を行列の成分と見る。この場合、行列の一番左上は0行0列目と見る。 たとえば

A1・B2・x・y^2  は1行2列目に対応し、x・y^2 を基底とする成分A1・B2とする。

これはベクトル

a = (A0、A1、A2、A3、・・・)
b = (B0,B1、B2、B3、・・・)

を使って、行列の転地tを用いて、

(ta)・b の形にかける。

これは、テンソル代数でいうところのベクトルa,bのテンソル積a〇bの計算と同じである。

f(x)・g(y)・h(z)のような3変数以上の場合も、係数の積によって、定義する。

テンソルT(I,J,K)の第(I,j,K)成分には、x^I・y^J・z^Kの係数Ai・Bj・Ck を使う。

つまり、3つのベクトル a,b,c のテンソル積 T = a〇b〇c は
T(I,J,K)=Ai・Bj・Ck である。
  • まとめ

テンソルのもともとの由来は材料力学の計算からだが、とりあえず、それは無視して以上のようにテンソル積の計算を考えたほうが手っ取り早い。テーラー級数の積によるテンソル積の定義では、無限次元のベクトルや行列、テンソルを自然に導入できる。