9
現代数学議論
文献あり

マクローリン展開を離散から連続へ1(アプローチ編)

561
1

厳密性を嫌ほど欠いている可能性があります。アレルギーの方はこの記事の閲覧に関してご遠慮願います。覚悟ができていれば見ていってください。
参考文献ありと書いてありますが、wikiによる分数階微分や反復積分の説明のみです。連続化についての具体的な文献はまだ見つかっておりません。

つぶやき

4/9に突然思いついたので書き記します。

マクローリン展開

f(x)=n=0xnn!f(n)(0)

まずマクローリン展開とはこういうものです。f(n)(0)fxn階微分しx=0としたものです。相変わらず美しい。
ここで思いついたやつを書きます。

0xtΓ(t+1)f(t)(0)dt

無限級数を広義積分に、階乗をガンマ関数に、n次導関数を実数階微分へ拡張しました。解ける気がしません。
ガンマ関数、実数階微分(実数次導関数)を積分として広げると四重積分になりますが...書いてみましょう。

Γ(t+1)=0uteudu
f(t)(0)=1Γ(t)c0f(v)(v)t+1dv=(kt1)kΓ(kt)c0f(v)(v)t+1dv=(kt1)k0skt1esdsc0f(v)(v)t+1dv  (c>0,k>t)

まずい予感がしますが、清書。厳密性欠きまくります。

0xtΓ(t+1)f(t)(0)dt=0xt(0uteudu)(kt1)k(0skt1esdsc0f(v)(v)t+1dv)dt=000c01v(uxv)te(u+s)skt1(kt1)kf(v)dvdudsdt

吐き気のする形になりました。何か別のアプローチが必要そうですね。

相反公式による簡易化

なぜか途中まで思いつかなかった。

0xtΓ(t+1)f(t)(0)dt=0xtΓ(t+1)1Γ(t)c0f(u)(u)t+1du
ここで、1Γ(t+1)Γ(t)=1Γ(1(t))Γ(t)=sin(πt)π=sinπtπ
0xtΓ(t+1)f(t)(0)dt=1π0xtsinπtc0f(u)(u)t+1du
=1π0c01u(xu)tf(u)sinπt dudt  (c>0)

発散しない場合はかなり限られますが、こっちのほうがまだ意味を持ちそうな形です。あとc>0定数をどうにかしたい。
ということで、一旦別のことを考えてみましょう。
形を保つのではなく、関数の視点から。

係数c(s)の定義

例えば、
f(x)=n=0cnxn
を満たすcn
cn=fn(0)n!
となるんでしたね。(特にf(x)cnの母関数なんて呼んだり)
つまるところ、
f(x)=0c(s)xsds
なる係数(の関数)c(s)を定義したくなってきたわけです。
知識がないので一概に言いたくはないのですが、chatGPT君に聞いたところメリン変換やラプラス変換を例に挙げてくれました。
...察したでしょうか。
これを満たすc(s)は、基本的に存在しません。(存在しないそうです)
分布(超関数的なの)として現れることはあるそうなので、一旦ディラックのデルタ関数みたいに定義しましょう。

f(x)=0c(s)xsds
を満たす仮想上(かもしれない)の関数c(s)積分展開係数と定める。
また、積分展開係数から定まる関数f(x)父関数と定める。

安直な名前であったほうが覚えやすい。
少し考えるとわかるのですが、収束させたくはあるもののsに関して積分するのでxsが指数関数のオーダーになってしまうため、c(s)にガンマ関数が含まれる可能性が非常に高いのです。つまり嬉しい関数でもあるし、病的でもある。
ということで、c(s)はなかなか見つけにくいので簡単な父関数側から見ていきましょう。

指数関数の積分展開係数

たい焼き氏の記事 から一部いただきます。

a>0 , Re(z)>0なるaR, zCに対し、
az=1Γ(z)0sz1easds

この式、変形できそうですよね。

az=01Γ(z)sz1zas1lnzds=0sz1Γ(z)z...ds...

...別アプローチを探ります。

積分展開係数から定める

こちらは私のフォロワー様からの案を一つ提示させていただきます。aを正の定数としてc(s)=saとしてみましょう。

0<x<1とする。
I=0saxsds=0saesln1xds
sln1x=t とおけば、
I=0(tln1x)aet1ln1xdt=Γ(a+1)(ln1x)a+1
よって、
Γ(a+1)(ln1x)a+1=0saxsds
特にa=1の場合、0<x<1で、
1(ln1x)2=0sxsds

最後の積分がめっちゃ重要そうに見える。発散が怖いので範囲は絞りました。あとマイナスがついてるとあれなのでlnを(実際正にはなるんですが)常に正にするために逆数にしました。
一種のラプラス変換にも見えるところがいい。そこでここのパートのおまけですが、x=Ωが面白そうなのでやってみましょう。

0sΩsds=1(ln1Ω)2=1Ω2

すっげぇ(語彙力皆無)
二個目紹介します。

x<aなるaに対し、c(s)=asとすると、
J=0asxsds=0(ax)sds=0eslnaxds
slnax=t とおけば、
J=0et1lnaxdt=1lnax
よって、
1lnax=0asxsds

公式的なのが量産できそうなやつですね。ノートに量産したやつを適当にまとめて、今度奮発してみよう。

父関数から定める

最高の関数(?)を導入すると、展開が容易になります。

0<x<1に対し、
11x=0(n=0δ(sn))xsds

f(x)=11x (0<x<1)とする。
f(x)=n=0xn
xn=0δ(sn)xsds , f(x)=n=00δ(sn)xsds
ここで、積分と無限和は交換できるので、
f(x)=0n=0δ(sn)xsds=0(n=0δ(sn))xsds
すなわち、父関数11xの積分展開係数はn=0δ(sn)である。

少し思考実験しただけなので積分無限和交換は一旦許してほしうございます。

Conclusion

これからのアプローチを個人的に確定させたので、シリーズ化します。シリーズ化好きなので。
2以降これから更新入れていきます。協力してくださった方々に感謝。誤りなどありましたら書き込んでいただきたい。ではまた。
ところでこの積分展開係数...結構...有用に見えますよね?

追記

'25/04/30: 誤りを訂正しました。
'25/05/03: 「追記」枠を導入し、参考文献に関する注意を加えました。
'25/05/03: 別アプローチとして、「積分展開係数」なるものを追記しました。
'25/05/04: 積分展開係数から定義することによって関数を求め、そのおまけを追記しました。
'25/05/05: 公式として一部証明を追加、指数の説明を追加しました。

参考文献

投稿日:427
更新日:55
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

関数をつくろう(掛詞)

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. つぶやき
  2. 相反公式による簡易化
  3. 係数c(s)の定義
  4. 指数関数の積分展開係数
  5. 積分展開係数から定める
  6. 父関数から定める
  7. Conclusion
  8. 追記
  9. 参考文献