8
現代数学解説
文献あり

sinxの分解にチャレンジ

425
0
$$$$

今回の問題

$\mathrm{rin}(\mathrm{rin}~x)=\sin x$を満たす$\mathrm{rin}x$は存在するか。
また、存在するならどのような関数か。

$\mathrm{rin}(0)$の値を求める

まずは一つ値を求めてみましょう。
前回導出した、合成不動点定理を利用して$\mathrm{rin}(0)$を求めます。
これからは$g(g(x))=f(x)$を満たす関数$g$$f_{\tiny[\frac{1}{2}]}(x)$と書きます。
前回の記事→ 合成関数の逆操作について

合成不動点定理

方程式$x=f(x)$の唯一の解が$\alpha$であるとき、$f_{\tiny[\frac{1}{2}]}(\alpha)=\alpha$

方程式$x=\sin x$を満たすのは$x=0$のときだけなので、$\mathrm{rin}(0)=0$

$\mathrm{rin}(0)=0$

一つ値が分かったので、微分してテイラー展開での表現を目指します。

微分しまくる

$\mathrm{rin}'(0)$の値

$\mathrm{rin}(\mathrm{rin}~x)=\sin x$の両辺を微分すると、左辺は合成関数の微分なので、
$$\mathrm{rin}'(x)\cdot\mathrm{rin}'(\mathrm{rin}~x)=\cos x$$ $x=0$を代入して、\begin{eqnarray}\mathrm{rin}'(0)\cdot\mathrm{rin}'(\mathrm{rin}~0)&=&\cos 0 \\[6pt] \mathrm{rin}'(0)\cdot\mathrm{rin}'(0)&=&1 \\[6pt] \left(\mathrm{rin}'(0)\right)^2&=&1\\[6pt]\mathrm{rin}'(0)&=&1 \end{eqnarray}

$\mathrm{rin}'(0)=1$

$\mathrm{rin}''(0)$の値

$\mathrm{rin}'(x)\cdot\mathrm{rin}'(\mathrm{rin}~x)=\cos x$をもう一度微分すると、
$$(\operatorname{rin}'(x))^2 \cdot \operatorname{rin}''(\operatorname{rin}(x)) + \operatorname{rin}''(x) \cdot \operatorname{rin}'(\operatorname{rin}(x)) =-\sin x$$


 導出
\begin{eqnarray}(\mathrm{rin}'(x)\cdot\mathrm{rin}'(\mathrm{rin}~x))'&=&\mathrm{rin}'(x)\cdot(\mathrm{rin}'(\mathrm{rin}~x))'+(\mathrm{rin}'(x))'\cdot\mathrm{rin}'(\mathrm{rin}~x) \\[7pt] &=&\operatorname{rin}'(x)\cdot \operatorname{rin}'(x)\cdot\operatorname{rin}''(\operatorname{rin}(x))+\operatorname{rin}''(x)\cdot\mathrm{rin}'(\mathrm{rin}~x) \\[7pt] &=&(\operatorname{rin}'(x))^2 \cdot \operatorname{rin}''(\operatorname{rin}(x)) + \operatorname{rin}''(x) \cdot \operatorname{rin}'(\operatorname{rin}(x)) \end{eqnarray}

$x=0$を代入して、\begin{eqnarray} (\operatorname{rin}'(0))^2 \cdot \operatorname{rin}''(\operatorname{rin}(0)) + \operatorname{rin}''(0) \cdot \operatorname{rin}'(\operatorname{rin}(0))&=&0 \\[6pt] 2\operatorname{rin}''(0)&=&0 \\[7pt] \operatorname{rin}''(0)&=&0 \end{eqnarray}

$\operatorname{rin}''(0)=0$

$\operatorname{rin}'''(0)$の値

結構めんどくさくなってきたのでwolfram alphaに計算してもらいました。

$(\operatorname{rin}'(x))^2 \cdot \operatorname{rin}''(\operatorname{rin}(x)) + \operatorname{rin}''(x) \cdot \operatorname{rin}'(\operatorname{rin}(x)) =-\sin x$をもう一度微分すると、
$$ \operatorname{rin}'''(\operatorname{rin}(x)) (\operatorname{rin}'(x))^3 + 3\operatorname{rin}''(\operatorname{rin}(x)) \operatorname{rin}'(x) \operatorname{rin}''(x) + \operatorname{rin}'(\operatorname{rin}(x)) \operatorname{rin}'''(x)=-\cos x $$
$x=0$を代入すると、
\begin{eqnarray} \operatorname{rin}'''(\operatorname{rin}(0)) (\operatorname{rin}'(0))^3 + 3\operatorname{rin}''(\operatorname{rin}(0)) \operatorname{rin}'(0) \operatorname{rin}''(0) + \operatorname{rin}'(\operatorname{rin}(0)) \operatorname{rin}'''(0)&=&-\cos 0 \\[7pt] 2\operatorname{rin}'''(0) &=& -1 \\[7pt] \operatorname{rin}'''(0)&=&-\frac{1}{2} \end{eqnarray}

$\operatorname{rin}'''(0)=-\frac{1}{2}$

同様にして、$\operatorname{rin}^{(4)}(0), \operatorname{rin}^{(5)}(0)...$の値も求める

あとはもう同様に計算していくだけです。
$\operatorname{rin}^{(10)}(0)$までの結果を記します。

$$ \operatorname{rin}^{(4)}(0)=0$$
$$ \operatorname{rin}^{(5)}(0)=-\frac{3}{4}$$
$$\operatorname{rin}^{(6)}(0)=0 $$
$$\operatorname{rin}^{(7)}(0)=-\frac{53}{8}$$
$$\operatorname{rin}^{(8)}(0)=0$$
$$\operatorname{rin}^{(9)}(0)=-\frac{1863}{16}$$
$$\operatorname{rin}^{(10)}(0)=0$$
微分した回数が偶数回の時は0になっていて、奇数回の時は分母が2の累乗のかたちになっています。

テイラー展開

$\operatorname{rin}(x)$をテイラー展開で近似した結果は次のようになります。

$$\operatorname{rin}(x)\sim x-\frac{1}{2\cdot3!}x^{3}-\frac{3}{2^2\cdot5!}x^{5}-\frac{53}{2^3\cdot7!}x^{7}-\frac{1863}{2^4\cdot9!}x^{9}$$

これを見ればわかるように、分母は規則的な形をしていますが、分子の規則性はあまりよくわかりません。
!FORMULA[44][-1276447451][0]の近似 (!FORMULA[45][-249001301][0]位まで近似できている) $\operatorname{rin}(x)$の近似 ($-\frac{\pi}{2}\le x\le \frac{\pi}{2}$位まで近似できている)
青:上式の2回合成 赤:!FORMULA[46][1755256236][0] 青:上式の2回合成 赤:$\sin x$
グラフを見る限り、正しく$ \operatorname{rin}(x)$を近似できてそうです。

テイラー展開の分子の形

検索してみたところ オンライン整数列大辞典のA098932 $\operatorname{rin}(x)$のテイラー展開の分子が載っていました。
最初の20項を載せておきます。
$$1, -1, -3, -53, -1863, -92713, -3710155, 594673187, 329366540401, 104491760828591, 19610322215706989, -5244397496803513989, -7592640928150019948759, -2156328049189410651012985, 3923796638128806973444887205,4835995123161649256573006676435, -4268483751449486260917713386079007,-18581347500475204057583673499997006497,-2144424268151426688705328023020488421667,108688458572837735664509183201007928496768811$$
符号が変化しながら、急速に数が巨大になっています。
一般項についてはよくわかりませんでした。わかったことがあればまた書きます。

結論

$ \operatorname{rin}(x)$は次のように近似される。
$$\operatorname{rin}(x)\sim x-\frac{1}{2\cdot3!}x^{3}-\frac{3}{2^2\cdot5!}x^{5}-\frac{53}{2^3\cdot7!}x^{7}-\frac{1863}{2^4\cdot9!}x^{9}-\frac{92813}{2^5\cdot11!}x^{11}-\frac{3710155}{2^6\cdot13!}x^{13}+\frac{594673187}{2^7\cdot15!}x^{15}$$

参考文献

投稿日:215
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。

投稿者

cbrtx
cbrtx
49
2355
何の変哲もないただの学生。

コメント

他の人のコメント

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