こんにちは. 今回も
超微分
に関する内容です.
ラグさんの
超・McLaughlin展開を構成的に導出してみる
にて導入された, 超・McLaughlin Diamond(SMD)の特殊な場合の収束性と, その系になる
超指数関数・超三角関数
(の一部)の超・McLaughlin展開(SM展開)を求めてみようと思います.
2026/10/26 19:00 主定理の主張及び証明を微修正・補足しました
超微分の操作を記号でとらえるために, 少し定義をしようと思います.
また, E関数, ダイヤモンド演算子の表記を少し変えます.
$U\subset\mathbb{R}^\times,\ n\in\mathbb{N}$とする.
$f:U\to\mathbb{R}$で$n$階超微分可能で, $n$階導関数が連続な関数全体を$\mathscr{P}^n(U)$とかく.
$f:U\to\mathbb{R}$で何回でも超微分可能な関数全体を$\mathscr{P}^\infty(U)$とかく.
($n+1$階超微分をするには$n$階超導関数$f^{[n]}$が零点を持たない必要があることに注意.)
自分は一般の議論では返り値が正と仮定しないので零点を持たないという形になります.
$U=(0,\infty)$とする.
$f(x)=e^x$とすると
$f^{`}(x)=x,\ f^{``}(x)=1,\ f^{```}(x)=0$
であるから,
$f\in \mathscr{P}^3(U)\setminus \mathscr{P}^2(U)$
$f^{`}\in \mathscr{P}^2(U)\setminus \mathscr{P}^1(U)$
$f^{``}\in \mathscr{P}^1(U)\setminus \mathscr{P}^0(U)$
となる.
(あとであらためて定義するが, )$\sxp\in\mathscr{P}^\infty((0,e))$である.
$U\subset\mathbb{R}^\times,\ n\in\mathbb{N}\cup\set{\infty}$とする.
超微分作用素$Q:\mathscr{P}^{n+1}(U)\to\mathscr{P}^{n}(U)$を$Q(f):=f^{`}$で定義する.
また, 超積分作用素$T:\mathscr{P}^{n}(U)\to\mathscr{P}^{n+1}(U)$を$\displaystyle T(f)(x):=\sint_1^x \r x^{f(x)}\left(=\exp{\left(\int_1^x \frac{f(t)}{t}dt\right)} \right)$で定義する.
($n,U$は暗黙的な表示となる.)
また, $Q,T$については誤解のおそれがない限り括弧を省略し右結合とする.
そして繰り返す場合, 指数を用いて表記する.
(例えば, $Q^2 Tf(x)$は$(Q(Q(T(f))))(x)$を意味する.)
$\r x$の意味については 自分の前の記事 を見てください. (見なくても別に定義してるのでわかると思いますし, 以降では使いませんが…)
$TQf=f,\ QTf=f(1)^{-1}f$が成り立つのが簡単にわかると思います.
$a\in\mathbb{R}^\times,\ n\in\mathbb{N}$について, $E_n(a)\in\mathscr{P}^{n+1}((0,\infty))$を次で定義する.
$E_0(a)=a$(定数関数)
$E_{m+1}(a)=TE_{m}(a)$
つまり, $\displaystyle E_{m+1}(a)(x)=\exp(\int_1^x\frac{E_m(a)(t)}{t})dt$
$U\subset\mathbb{R}^\times,\ n\in\mathbb{N}$について$\diamond_n:\mathscr{P}^n(U)\times\mathscr{P}^n(U)\to\mathscr{P}^{\min\set{1,n}}(U)$を
$f\diamond_0 g:=fg$ ただし$(fg)(x)=f(x)g(x)$
$f\diamond_{n+1} g:=T((Qf)\diamond_{n} (Qg))$
積を用いているので, あまり高階の超微分可能性はわかりません.
次の命題は定義から直ちにわかると思います.
$a\in\mathbb{R}^\times,\ n\in\mathbb{N}$について, $E_n(a)=T^n a$
$f\diamond_n g=T^n(Q^nf)(Q^ng)$
ただし, $T$より右側の関数の積は先に計算するものとする.
この表記を使ってSMDを有限個で打ち切ったものの簡単な表記を得ます.
$n\in\mathbb{N}_{>0}$とする. 各項が零でない実数列$(a_i)_{i=0}^n$について,
$M_n:=E_0(a_0)\diamond_0 (E_1(a_1)\diamond_1 \cdots (E_{n-1}(a_{n-1})\diamond_{n-1}E_n (a_n)))$
は$a_0Ta_1T\ldots a_{n-1}Ta_n$と表せる.
慣用的に数列を明示するために$M_n(a_i)$と書く. 第$i$項のみによるわけではない. (この記事では, $i$はダミー変数以外に用いない. )
$E_{n-1}(a_{n-1})\diamond_{n-1}E_n (a_n)$
$=T^{n-1}(Q^{n-1}E_{n-1}(a_{n-1}))(Q^{n-1}E_n (a_n))$
$=T^{n-1}(Q^{n-1}T^{n-1}a_{n-1})(Q^{n-1}T^n a_n)$
$=T^{n-1}a_{n-1}T a_n$
$E_{n-2}(a_{n-2})\diamond_{n-2}(T^{n-1}a_{n-1}T a_n)$
$=T^{n-2}(Q^{n-2}T^{n-2}a_{n-2})(Q^{n-2}T^{n-1}a_{n-1}T a_n)$
$=T^{n-2}a_{n-2}Ta_{n-1}Ta_n$
これを続けると目標の式を得る. $\square$
超微分方程式$y^{`}=y$の解を超指数関数,
超微分方程式$y^{``}=y$の解を超双曲線関数という.
初期値$y(1)=a$を満たす超指数関数を$\sxp_a$, 特に$\sxp:=\sxp_1$とする.
初期値$y(1)=a, y^{`}(1)=b$を満たす超双曲線関数を$\hb_{a,b}$とする.
$1$を中心に考えるので$\mathrm{scosh}$,$\mathrm{ssinh}$について直接は議論しません.
具体的な解の表示は次に示します.
超微分方程式$y^{`}=y,\ y(1)=a$の解は$\displaystyle y=\frac{1}{\frac{1}{a}-\ln (x)}$
$y^{`}(x)=y(x)$
$\iff xy'(x)=y(x)^2$
$$\iff \frac{y'(x)}{y(x)^2}=\frac{1}{x}$$
$$\iff \int_1^x\frac{y'(t)}{y^2(t)}dt=\int_1^x\frac{1}{t}dt$$
$$\iff \frac{1}{a}-\frac{1}{y(x)}=\ln (x)$$
$$\iff y(x)=\frac{1}{\frac{1}{a}-\ln (x)}$$
$\square$
(初期値が$1$で, $x=0$での超微分を考え(たく)ないため, $x>0$のみを考えた. 以下も同様.)
超微分方程式$y^{``}=y,\ y(1)=a,\ y^{`}(1)=b$の解は
$a=b$なら$\displaystyle y=\frac{1}{\frac{1}{a}-\ln (x)}$
$a\neq b$なら$\displaystyle y=\frac{(b-a)x^{b-a}}{\frac{b}{a}-x^{b-a}}$
$$y^{``}(x)=y(x)$$
$$\iff x(y^{`})'(x)=y^{`}(x)y(x)$$
$$\iff x(y^{`})'(x)=xy'(x)$$
$$\iff (y^{`})'(x)=y'(x)$$
$$\iff \int_1^x (y^{`})'(t) dx=\int_1^x y'(t)dx$$
$$\iff y^{`}(x)-b=y(x)-a$$
$a=b$なら超指数関数に等しい.
以下, $a\neq b$とする. $y$が定数でないことに注意する. また, 符号の変わる点以上は考えない.
$$y^{`}(x)-b=y(x)-a$$
$$\iff y^{`}(x)=y(x)+b-a$$
$$\iff xy'(x)=y(x)(y(x)+b-a)$$
$$\iff \frac{y'(x)}{y(x)(y(x)+b-a)}=\frac{1}{x}$$
$$\iff \int_1^x\frac{y'(t)}{y(t)(y(t)+b-a)}dt=\int_1^x\frac{1}{t}dt$$
$$\iff \frac{1}{b-a}\left(\ln\left({\frac{y(x)}{a}}\right)-
\ln\left({\frac{y(x)+b-a}{a+b-a}}\right)\right)=\ln(x)$$
$$\iff \ln\left({\frac{y(x)}{a}}\right)-
\ln\left({\frac{y(x)+b-a}{b}}\right)=\ln(x^{b-a})$$
$$\iff \left({\frac{y(x)}{a}}\right)
\left({\frac{b}{y(x)+b-a}}\right)=x^{b-a}$$
$$\iff y(x)\left({\frac{b}{a}}\right)=(y(x)+b-a)x^{b-a}$$
$$\iff y(x)\left({\frac{b}{a}}-x^{b-a}\right)=x^{b-a}(b-a)$$
$$\iff y(x)=\frac{(b-a)x^{b-a}}{{\frac{b}{a}}-x^{b-a}}$$
$\square$
$[1,a]$上で定義される連続関数$f,g$について,
$\forall x\in [1,a], 0\leq f(x)\leq g(x) $とすると
$Tf(a)\leq Tg(a)$
また, $[a,1]$($a>0$)上で定義される連続関数$f,g$について,
$\forall x\in [a,1], 0\leq f(x)\leq g(x) $とすると
$Tf(a)\geq Tg(a)$
$f,g$を$[1,a]$上で定義される連続関数とし,
$\forall x\in [1,a], 0\leq f(x)\leq g(x) $とする.
$$Tf(a)=\exp(\int_1^a \frac{f(t)}{t}dt)\leq\exp(\int_1^a \frac{g(t)}{t}dt)=Tg(a)$$
$ $
$f,g$を$[a,1]$上で定義される連続関数とし,
$\forall x\in [a,1], 0\leq f(x)\leq g(x) $とする.
$$Tf(a)=\exp(\int_1^a \frac{f(t)}{t}dt)=\exp(-\int_a^1 \frac{f(t)}{t}dt)$$
$$\geq\exp(-\int_a^1 \frac{g(t)}{t}dt)=\exp(\int_1^a \frac{f(t)}{t}dt)=Tg(a)$$
$$\square$$
$$Tf(e^{-x})=\exp\left(-\int_0^x f(e^{-s})ds\right)$$
$t=e^{-s}$と置換すると, $s=-\ln t,\ -ds=\frac{1}{t}dt$だから,
$$Tf(e^{-x})=\exp\left(\int_1^{e^{-x}} \frac{f(t)}{t}dt\right)$$
$$=\exp\left(-\int_0^x f(e^{-s})ds\right)$$
次の命題は簡単に示せますが, 一応紹介します.
$x,y\in\mathbb{R}$について, $|e^{-x}-e^{-y}|\leq |x-y| $
また$x>0$について, $0\leq 1-e^{-x}\leq x$
WLOGで$x< y$と仮定できる.
平均値の定理より, $x< c< y$で$e^{-y}-e^{-x}=-e^{-c}(y-x)$なる$c$がとれる.
$0< c$より, $e^{-c}\leq 1$
$|e^{-y}-e^{-x}|=e^{-c}|y-x|\leq |y-x|$
$x>0$のとき, $0<1-e^{-x}$より下の主張も従う. $\square$
数列$(a_n)_{n=0}^\infty$が
$$\forall n\in\mathbb{N},\ 0< a_n$$
$$\exists N\in \mathbb{N},\ \exists c>0,\ \forall n\in\mathbb{N},\ n\geq N\Longrightarrow a_n\leq c$$
を満たすならば, $M_n(a_n)$は$(0,e^{-c})$で(広義一様)収束する.
関数列$(f_n)_n$が$(0,e^{-c})$で広義一様収束するなら, 積分の一般論および$\exp$の連続性から,$p>0$を定数として, $(pTf_n)_n$も広義一様収束する.
よって, すべての$n$で$0< a_n\leq c$の場合について示せばいい.
★$[1,e^{-c})$において考える.
まず, $M_1(a_i)(x)=a_0 x^{a_1}\geq a_0=M_0(a_i)(x)$である.
また, $n$が「任意の正項数列$(a_i)_i$について$M_{n+1}(a_i)\geq M_{n}(a_i)$」
を満たすとすると$T$の単調性から,
$M_{n+2}(a_i)=a_0TM_{n+1}(a_{i+1})\geq a_0TM_{n}(a_{i+1})=M_{n+1}(a_i)$
より$n+1$でも成り立つ. よって帰納法より, $M_n(a_i)$は$n$について単調増加する. (ちなみに, ここまで有界性は用いていない.)
次に, $M_0(a_i)=a_0\leq c\leq \sxp_c$であり,
$n$が「任意の$n$で$0< a_n\leq c$となるような任意の数列$(a_i)_i$について$M_n(a_i)\leq \sxp_c$」とすると, $T$の単調性から, $M_{n+1}(a_i)=a_0TM_{n}(a_{i+1})\leq cT\sxp_c = \sxp_c$
よって, $M_n(a_i)$は$n$について単調な有界数列となるため, 広義一様収束する.
★$(0,1)$において考える.
数列$(a_i)_i$は固定し, $M_n(a_i)$を$M_n$と書く.
$\tilde{M}_{k,n}(u)=a_k T\cdots Ta_n(e^{-u})$とする.
$\tilde{M}_{n,n}(u)=a_n$であり,
$\tilde{M}_{k,n}(u)=a_k T (a_{k+1} T\cdots Ta_n)(e^{-u})$
$\displaystyle=a_k \exp\left(-\int_0^u(a_{k+1} T\cdots Ta_n)(e^{-s})ds\right) $ (補題7)
$\displaystyle=a_k \exp\left(-\int_0^u \tilde{M}_{k+1,n}(s) ds\right) $
が成り立ち,
$\displaystyle M_n(x)=\tilde{M}_{0,n}(-\ln x)=a_0 \exp\left(-\int_0^{-\ln x} \tilde{M}_{1,n}(s) ds\right)$が成り立つ.
($-\ln x>0$であることに注意.)
$n\geq 1,\ u=-\ln x\ (>0)$とする.
$0\leq k\leq n-1,\ s\in[0,u]$について
$\displaystyle |\tilde{M}_{k,n}(s)-\tilde{M}_{k,n-1}(s)|\leq \frac{c^{n-k+1} s^{n-k}}{(n-k)!}$を示す.
$\tilde{M}_{n-1,n}(s)=a_{n-1}e^{-a_n s},\ \tilde{M}_{n-1,n-a}(s)=a_{n-1}$だから,
$|\tilde{M}_{n-1,n}(s)-\tilde{M}_{n-1,n-1}(s)|=a_{n-1}(1-e^{-a_n s})\leq a_{n-1}a_n s\leq c^2 s$
$k\leq n-2$, $\displaystyle \forall s\in[0,u],\ |\tilde{M}_{k+1,n}(s)-\tilde{M}_{k+1,n-1}(s)|\leq \frac{c^{n-k} s^{n-k-1}}{(n-k-1)!}$とする.
$$\displaystyle |\tilde{M}_{k,n}(s)-\tilde{M}_{k,n-1}(s)|$$
$$ =a_k\left| \exp\left(-\int_0^s\tilde{M}_{k+1,n}(\sigma)d\sigma\right) - \exp\left(-\int_0^s\tilde{M}_{k+1,n-1}(\sigma)d\sigma\right)\right|$$
$$\leq a_k \left|\int_0^s\tilde{M}_{k+1,n}(\sigma)d\sigma - \int_0^s\tilde{M}_{k+1,n-1}(\sigma)d\sigma\right|$$
$$\leq a_k \int_0^s |\tilde{M}_{k+1,n}(\sigma)-\tilde{M}_{k+1,n-1}(\sigma)|d\sigma\leq c \int_0^s \frac{c^{n-k} \sigma^{n-k-1}}{(n-k-1)!}d\sigma$$
$$=\frac{c^{n-k+1}s^{n-k}}{(n-k)!}$$
よって任意の$k$で成立. 特に$k=1$とすると,
$\displaystyle |\tilde{M}_{1,n}(s)-\tilde{M}_{1,n-1}(s)|\leq \frac{c^{n} s^{n-1}}{(n-1)!}$
よって,
$$|M_n(x) - M_{n-1}(x)|$$
$$=a_0\left|\exp\left(-\int_0^u \tilde{M}_{1,n}(s)ds\right)-\exp\left(-\int_0^u \tilde{M}_{1,n-1}(s)ds\right)\right|$$
$$\leq a_0\left|\int_0^u \tilde{M}_{1,n}(s)ds-\int_0^u \tilde{M}_{1,n-1}(s)ds\right|$$
$$\leq a_0\int_0^u\left| \tilde{M}_{1,n}(s)ds- \tilde{M}_{1,n-1}(s)\right|ds$$
$$\leq c\int_0^u\frac{c^{n} s^{n-1}}{(n-1)!}ds=\frac{c^{n+1} u^{n}}{n!}$$
$$\sum_{k=1}^{N}|M_{k}(x)-M_{k-1}(x)|\leq \sum_{k=1}^{N}\frac{c^{n+1} u^{n}}{n!}\leq c(e^{cu}-1)$$
より, $\displaystyle M_N-M_0=\sum_{k=1}^{N}(M_{k}-M_{k-1})$の右辺の級数は絶対収束する.
よって$M_N$は$(0,1]$で広義一様収束する.
$\square$
よって, O(1)な正数列について, SMDは収束することがわかりました.
最後にSMDを計算するための命題と例を紹介しようとおもいます.
数列$(a_i)_{i=0}^\infty$について, ある$U\subset \mathbb{R}$上で$(M_n(a_i))_n,\ (M_n(a_{i+1}))_n$が広義一様収束するとする.
このとき, $M(a_i)^{`}=M(a_{i+1})$である.
仮定から,
$M_{n+1}(a_i)^{`}=M_n(a_{i+1})$は$M(a_{i+1})$に広義一様収束する
$$M_{n+1}(a_i)^{`}(x)=\frac{xM_{n+1}(a_i)'(x)}{M_{n+1}(a_i)(x)}$$
であるから, $\displaystyle M_{n+1}(a_i)'(x)=\frac{M_{n+1}(a_i)^{`}(x)\cdot M_{n+1}(a_i)(x)}{x}$は$\displaystyle\frac{M(a_{i+1})(x)\cdot M(a_i)(x)}{x}$に広義一様収束する.
よって, 微分と極限が交換でき,
$$M(a_i)'(x)=\frac{M(a_{i+1})(x)\cdot M(a_i)(x)}{x}$$
$$\therefore M(a_{i+1})(x)=\frac{xM(a_i)'(x)}{M(a_i)(x)}=M(a_i)^{`}(x)\ \square$$
特に, 定義域の各点$a$の近傍である数列$(a_i)_{i=0}^\infty$が存在して$f(x)=M(a_i)(\frac{x}{a})$の形で表せる関数を$\mathscr{P}^\omega$級と表しましょう.
$\sxp_c(x)$は$(0,e^{\frac{1}{c}})$で$\mathscr{P}^\omega$級です.
ここまでお読みいただきありがとうございました.
最後に, 予想などについて残したいと思います.
以上になります.
極限ダイヤモンドについてはおもったより難しかったのでおそらく記事にしないか, するにしてもだいぶ先になりそうです. (だれかがやりたければもちろん全然譲ります.)
改めて, なにか意見等があればコメント,
$\mathbb{X}$
,
マシュマロ
等お待ちしています!
特にマシュマロはアカウントがなくとも匿名で長文が送れるのでおススメです!
(画像や数式を含む場合は他二つのほうがやりやすいかもしれないです. )
ちなみにSM展開の表記が危ういのはわざとです.