11

二重数の解析学考えてみた!!

741
2
$$$$

はじめに

今回は, 二重数について考えていきます。二重数というのは, いわば複素数の兄弟のような数で, 自動微分などの考え方がプログラミングに応用されています。今日は, そんな数で微分積分やゼータ関数を考えてみましょう!

今回はあくまで研究成果を書くだけなので, 厳密な証明を省く場合が多々あります。後々記事にはするつもりですので悪しからず。
また, 研究段階なので「あれ?これで終わり」となることが多々あります
(´;ω;`) <研究することが多すぎる…

二重数とは

二重数とは$\begin{align} \varepsilon ^2=0(\varepsilon \neq 0) \end{align}$という数のことです。複素数と同じように計算できるのでちょっとやってみましょう。
$ (a+\varepsilon b)+(c+\varepsilon d)=(a+c)+\varepsilon(b+d)$
$(a+\varepsilon b)-(c+\varepsilon d)=(a-c)+\varepsilon(b-d)$
$(a+\varepsilon b)(c+\varepsilon d)=ac+\varepsilon(bc+ad)$
$$\frac{a+\varepsilon b}{c+\varepsilon d}=\frac ac +\varepsilon \frac{bc-ad}{c^2}$$

上の式に$c=0$を代入したらわかりますが, 純二重数は逆元を持ちません。
そのため, 除法には細心の注意を払う必要があります。
また, 重要な性質として, 以下のように微分できることが挙げられます。

二重数による自動微分

多項式$P(x)$に対して
$$P(a+\varepsilon b)=P(a)+\varepsilon bP'(a)$$

これは, 展開すればすぐにわかります。
$$P(x):=c_0+a_1x+c_2x^2+c_3x^3+...+c_nx^n$$
とすると,
$$(a+\varepsilon b)^k=a^k+ka^{k-1}\varepsilon b+{}_k\mathrm{C }_2a^{k-2}\varepsilon ^2b^2+{}_k\mathrm{C }_3a^{k-3}\varepsilon ^3b^3+...$$

3項目以降には$\varepsilon ^2$が入っているので0になり
$$(a+\varepsilon b)^k=a^k+\varepsilon ka^{k-1}b$$
$$P(a+\varepsilon b)=c_0+c_1(a+\varepsilon b)+c_2(a+\varepsilon b)^2+...+c_n(a+\varepsilon b)^n$$
$$     =c_0+c_1(a+\varepsilon b)+c_2(a^2+\varepsilon 2ab)+...+c_n(a^n+\varepsilon na^{n-1}b)$$
$$     =(c_0+c_1a+c_2a^2+...+c_na^n)+\varepsilon b (c_1+c_22a+...+c_nna^{n-1})$$
$$     =P(a)+\varepsilon b P'(a)$$
となります。

無限次の場合に適応していいのかについては議論すべきかもしれませんが, とりあえず一般の関数にも適用させてみます。すると
$$\sin(a+\varepsilon b)=\sin(a)+\varepsilon b\cos(a)$$
$$\cos(a+\varepsilon b)=\cos(a)-\varepsilon b\sin(a)$$
$$\tan(a+\varepsilon b)=\tan(a)-\frac{\varepsilon b}{\cos ^2(a)} $$

$$\exp(a+\varepsilon b)=\exp(a)(1+\varepsilon b)$$
$$\log(a+\varepsilon b)=\log(a)+\frac {\varepsilon b}{a}$$
であることがわかります。
では, 微分について考えてみましょう!

二重数関数の微分

二重数関数の微分

定義が不完全なので無視してください
$x,y \in \mathbb{R}, u(x,y),v(x,y)$を実関数とする
$f(z)=f(x+\varepsilon y:=u(x,y)+\varepsilon v(x,y)$と表せるとき, $f(z)$を二重数関数といい, 微分は
$$\frac{df(z)}{dz}:=\lim_{\sigma \rightarrow 0}\lim_{\tau\rightarrow 0}\frac{f(z+\sigma+\varepsilon\tau)-f(z)}{\sigma+\varepsilon\tau}=\lim_{\sigma \rightarrow 0}\lim_{\tau\rightarrow 0}\frac{(f(z+\sigma+\varepsilon\tau)-f(z))(\sigma+\varepsilon\tau)}{\sigma^2}$$
となる。

(※二重数の性質上, 極限をとるのが難しかったので, 実数の極限でお茶を濁しています)
「まあ, そうだよね」としか言いようがない定義ですね。
複素関数のときと全く同じです。しかし!次の定理で「二重数すげぇ!」ってなるはずです。

正則関数の定式化

正則関数(微分可能な二重数関数)はすべてこのように表せる。
$$f(x+\varepsilon y)=\varphi(x)+\varepsilon(y\varphi '(x)+\psi(x))$$
ただし, $\varphi,\psi$は実一変数関数

(概略)

$$\lim_{h\rightarrow 0}\frac{f(z+h)-f(z)}{h}$$
$$\lim_{h \rightarrow 0}\frac{f(z+h+\varepsilon h)-f(z)}{h+\varepsilon h}$$
を比較することによって, コーシー・リーマンの方程式の類似
$$\frac{\partial u}{\partial y}=0,\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y}$$
が得られます。そこから代入して色々すると元の式が得られます。

また, このことから$C^1$級関数$f$の, $\mathbb{R}$から$\mathbb{R}[\varepsilon]$への拡張は, $f(x+\varepsilon y)=f(x)+\varepsilon yf'(x)$ただ一つであるとわかります。
($\psi(x)=0$でないと定義域が実数なのに値域が二重数になるから)
そして, そこから複素数と同様にコーシーの積分定理が得られます。

二重数での積分

コーシーの積分定理(二重数ver)

$f(z)$を正則関数, $C$をジョルダン曲線とする。
$$\oint_C f(z)dz=0$$

$$\oint _Cf(z)dz=\int_a^b\left(\varphi+\varepsilon \left(y\varphi'+\psi\right)\right)\left(x'+\varepsilon y'\right)dt$$
$$    \ \ \ =\int_a^b\left(\varphi x'+\varepsilon(y\varphi'x'+\psi x'+\varphi y')\right)dt$$
$$    \ \ \ =\int_a^b(x'(\varphi(x)+\varepsilon \psi(x))dt+\varepsilon\int_a^b(y\varphi'x'+y'\varphi)dt$$
第一項は$C$は閉曲線なので, $x(a)=x(b),y(a)=y(b)$となり,
$$\int_a^b(x'(\varphi(x(t))+\varepsilon\psi(x(t)))dt=\int_{x(a)}^{x(a)}(\varphi(x)+\varepsilon\psi(x))dx=0$$
第二項は積の微分公式より,
$$\frac{d}{dt}(y(t)\varphi(x(t)))=y'(t)\varphi(x(t))+y(t)\frac{d}{dt}\varphi(x(t))$$
$$        =y'(t)\varphi(x(t))+y(t)x'(t)\varphi'(x(t))$$
よって
$$\varepsilon\int_a^b(y\varphi'x'+y'\varphi)dt=\varepsilon\left[y(t)\varphi(x(t))\right]_a^b$$
$$        \ \ \ =\varepsilon(y(b)\varphi(x(b))-y(a)\varphi(x(a)))$$
$$        \ \ \ =\varepsilon(y(a)\varphi(x(a))-y(a)\varphi(x(a)))$$
$$        \ \ \ =0$$
$$ \therefore \oint_Cf(z)dz=0$$

さて, これでコーシーの積分定理が示せたので, コーシーの積分公式(やグルサの定理)に...と言いたいところですが, ここで二重数の厄介な性質が浮き彫りになります。
それは有限個の点を除いて正則な関数が作れないという性質です。
例えば, $f(z)=\frac 1z$という関数を考えてみます。
するとこの関数は$z=0$だけでなく, $z=913\varepsilon$$z=\pi\varepsilon$,$z=- \sqrt[334]{114514} \varepsilon$などで発散してしまいます。これでは周回積分ができません。
ここでやめたらいいね貰えないので何か別の方法を考えましょう。

あ!複素数って便利じゃん!!

そう, 複素数なら積分がすごいやりやすいのです。
なのでここからは$\mathbb{R}[\varepsilon]$上の線積分ではなく$\mathbb{C}[\varepsilon]$上の面積分を考えます。
まず手始めに$\Gamma_k= \lbrace e^{i\theta}+\varepsilon k \vert 0 \leq \theta \lt 2\pi\rbrace $上で$f(z)=\frac 1z$を線積分します。
すると$$\oint_{\Gamma_k}\frac{dz}{z}=\int_{0}^{2\pi}\frac{ie^{i\theta}}{e^{i\theta}+\varepsilon k}d\theta=i\int_{0}^{2\pi}(1-\varepsilon ke^{-i\theta})d\theta=2i\pi$$
となり, 普通に複素積分したときと同じようになります。
つまり, 純二重数はある意味で0と同じように考えられそうだとわかります。
では, 今度は円柱の側面$D=\lbrace e^{i\theta}+\varepsilon t \vert 0 \leq \theta \lt 2\pi,a \leq t \lt b\rbrace$で面積分してみます。
$$\iint_{D} \frac {dz}{z}=\int_a^b\int_{0}^{2\pi}\frac{ie^{i\theta}・\varepsilon}{e^{i\theta}+\varepsilon t}d\theta dt=\varepsilon\int_a^b\left(\int_{0}^{2\pi}\frac{ie^{i\theta}}{e^{i\theta}+\varepsilon t}d\theta\right)dt=2i\pi\varepsilon(b-a)$$
となります。また, ここから応用してグルサの定理も導けます

二重数のグルサの定理(円柱)

$$\iint_D\frac{f(z)}{(z-\alpha)^{n+1}}dz=\frac{2i\pi \varepsilon(b-a)}{n!}f^{(n)}(\alpha)$$
ただし$D=\lbrace e^{i\theta}+\varepsilon t+\alpha \vert 0 \leq \theta \lt 2\pi,a \leq t \leq b\rbrace$

が...これって面積分の定義から考えたら当たり前なんですよね。
積分する領域が円柱形だから複素積分の結果に$\varepsilon(b-a)$をかけてあげればいいだけですので。
しかし, 残念ながらこれ以上書けるような一般論が見つかっていないので, 積分の話はここまでです。(トーラスでの積分も考えましたがめんどくさいのでまたの機会に)
次は, この二重数を拡張してみましょう!!!

多重数

ここまでは$2$乗して$0$になる数を考えていましたが, $n$乗して$0$になる数を考えてみましょう。

n重数

$$\varepsilon_n^n=0$$
ただし, $\varepsilon_n,\varepsilon_n^2,...,\varepsilon_n^{n-1}$は0でなく, 線形独立であるとする。

性質

  • $R[\varepsilon_n] \cong R[X]/(X^n)$
  • $A \subset R[\varepsilon_n] \Longrightarrow A\cong R[\varepsilon_d] (d|n)$
  • $$f(a+\varepsilon_nb)=\sum_{k=0}^{n}\frac{f^{(k)}(a)}{k!}\varepsilon_n^kb^k$$

$3$番目から伝わると思いますが, それなりにいい感じの性質を二重数から受け継いでいます。
$n$重数が与えられると, 今度は整数に二重数と三重数(みえじゃないよ)を添加した環など, $2$種類以上の$n$重数を添加したものを考えたくなりますよね。
しかし, ここで一つ問題が出てきます。
$$\varepsilon_{nm}^m=k\varepsilon_n$$としてよいのかという問題です。確かに両辺の数は$n$乗して初めて$0$になる数です。だからといって, 同一視しても良いものでしょうか...
じゃあもう同一視するものとしないもの, 二種類作っちゃいましょう!!

離散多重数

$R$の離散多重数環$\varepsilon_{\mathbb{N}}(R)$とは
$$\varepsilon_{\mathbb{N}}(R):=R[\varepsilon_2,\varepsilon_3,...,\varepsilon_n,...]$$
である。ただし, 任意の$n$に対して$\varepsilon_n\notin R[\varepsilon_2,\varepsilon_3,...,\varepsilon_{n-1},\varepsilon_{n+1},...]$
また,
$$\alpha=\prod_{i}\varepsilon_i$$
とするとき, $\alpha$の次数を
$$\mathrm{deg}(\alpha)=\sum_i\frac{1}{i}$$と定義する。

連続多重数

$n\in\mathbb{N},0< x<1,0< y<1$とする$$\varepsilon\left(\frac{1}{n}\right):=\varepsilon_n,\varepsilon\left(x\right)\varepsilon\left(y\right)=\varepsilon(x+y)$$
を添加した環
$$\varepsilon_{\mathbb{Q}}(R):=R[\varepsilon(q_1),\varepsilon(q_2),\varepsilon(q_3),...]   (q_k\in\mathbb{Q},0< q_k<1)$$$R$の連続多重数と定義する。

$\varepsilon\left(\frac{1}{n}\right)$のように逆数をとっているのは, $0$$\frac 1n$乗感を出させて, 指数法則もどきを違和感なく採用できるようにするためです。また, $\varepsilon(\mathrm{deg(\alpha)})=\alpha$となっていろいろ都合がいいです。連続多重数の基底は有理数なため稠密性があり, それをどう使っていくかが今後の課題です。
では最後にゼータ関数の話をしましょう!!

ゼータ関数

デデキント・ゼータ関数

$K$を環とする。$K$のイデアル$M$のノルム$N(M)=\vert K /M\vert $とするとき
$$\zeta_K(s):=\prod_{M\in \mathrm{Spec}( \mathcal{O} _K)}\frac{1}{1-{N(M)}^{-s}}$$

$K=\mathbb{Q}[\varepsilon]$とすると
$$\mathcal{O}_K=\mathbb{Z}+\varepsilon\mathbb{Q}$$
$\mathbb{Q}$が入っているから体だろ!と思われますが, 計算してみると除法に閉じていないことがわかります。
$\mathrm{Spec}(\mathcal{O}_K)= \lbrace (p) | p\mathrm{は素数}\rbrace $
となるため$\zeta_{\mathbb{Q}[\varepsilon]}(s)=\zeta(s)$となります。
これでは少々味気ないですね。少々定義を改良してみましょう!

どろぉどの$\zeta$関数

$\varepsilon_{x_1},\varepsilon_{x_2},...,\varepsilon_{x_n}$を離散多重数とする。

$$\zeta(\varepsilon_{x_1},\varepsilon_{x_2},...,\varepsilon_{x_n};s):=\prod_{\mathfrak{p}\in\mathrm{Spec}(\mathbb{Z}[\varepsilon_{x_1},\varepsilon_{x_2},...,\varepsilon_{x_n}])}\frac{1}{1-N(\mathfrak{p})^{-s}}$$

$$\zeta(\varepsilon_n;s)=\prod_{k=1}^{n}\zeta(ks)$$

($n=2$のみ)

$\mathbb{Z}[\varepsilon]$のイデアルを$I$とし, $a+\varepsilon b\in I, x+\varepsilon y\in \mathbb{Z}[\varepsilon]$とすると
$$(a+\varepsilon b)(x+\varepsilon y)=ax+\varepsilon(ay+bx)\in I$$
$$I=mn\mathbb{Z}+\varepsilon n\mathbb{Z}  (m,n\in\mathbb{N})$$
と表せる。よって
$$\zeta(\varepsilon_2;s)=\sum_I\frac{1}{N(I)^s}=\sum_{(m,n)\in\mathbb{N}^2}\frac{1}{(mn^2)^s}=\sum_{m=1}^{\infty}\frac{1}{m^s}\sum_{n=1}^{\infty}\frac{1}{n^{2s}}=\zeta(s)\zeta(2s)$$

もしかしたら, このゼータ関数の性質が分かればリーマンゼータ関数もわかるかも...?
では, 添加する数を増やしたらどうなるのでしょうか。
$2$個の場合は...

$$\zeta(\varepsilon_a,\varepsilon_b;s)=\prod_{\begin{align}1 \leq n \leq a \\ 1 \leq m \leq b\end{align}}\zeta(nms) $$

となります。複数個の場合でも同様のことが言えそうですね。

今後の研究内容

  • 二重複素積分における留数
  • トーラスの表面上での二重複素積分
  • 二重四元数
  • どろぉどの$\zeta$関数の一般論(連続多重数の場合など)
  • 多重数の行列表現
    とくに, 二重四元数に関しては, 元々あった成果が実はミスだと判明した結果, 完全に手つかずなので何とかしたいです。

あとがき

色々と粗雑でしたが何とか終わりました(怪しい部分は指摘してください)。にしても勝手に自分の名を冠した関数を作ってもよかったのだろうか...
それにしても前半に比べて後半が手抜き過ぎるな
この記事を読んで少しでも二重数に興味をもったなら, この未開の地に踏み込んでみてください!!!

いいねお願いします!!!

投稿日:20231224
更新日:20231230

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

かくことは あれども筆は 進まぬか コマンド打つのが 割とめんどい

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中