10

ラプラス変換を用いて解ける積分(積分botの問題を解いてみた4)

466
0

今回はラプラス変換を用いて解ける積分を解説します。
次の変形を事前に把握しておいてください.
(a,bの範囲はarg=arctanとできる範囲)

1(a+ib)s
=esln(a+ib)=es2ln(a2+b2)siarctanba=esiarctanba(a2+b2)s/2

問題11

0πarctanlnsinxxdx=πarctan2ln2π

0πarctanlnsinxxdx
=0ππ2arctanxlnsinxdx
0πarctanxlnsinxdx
=Im0πln(lnsinx+ix)dx=Imdds0πdx(lnsinx+ix)s=Imdds1Γ(s)0π0ts1e(lnsinx+ix)tdtdx=dds1Γ(s)0ts10πsintxsintxdxdt
ここで
0πsinu1xsinvxdx=πsinπv22u1uB(u+v+12,uv+12)
を用いる。
これは留数定理を使って証明できるが今回は省略する。
dds1Γ(s)0ts10πsintxsintxdxdt
=ddsπΓ(s)0ts1sinπt22tdt=ImddsπΓ(s)0ts1et(ln2iπ2)dt=Imddsπ(ln2iπ2)s=πarctanπ2ln2
0πarctanlnsinxxdx=πarctan2ln2π

これ実は
0πsin(sarctanxlnsinx)(ln2sinx+x2)s/2dx=πsin(sarctanπ2ln2)(ln22+π24)s/2
を求めて微分してるだけなんですよね。
解き方は私の一個目の解説

https://mathlog.info/articles/U1SticPVy5AVxLrB6cax
の問題2と同じなんですよ。
つまり

0πln(ln2sinx+x2)dx=π2ln(ln22+π24)
0π/2arctanlncosxxdx

も求められる訳です。

問題12

sin(rarctanx)(1+x2)r/2(eπx+1)dx=η(r)ζ(r)

sin(sarctanx)(1+x2)s/2(eπx+1)dx
=ImΓ(s)1eπx+10ts1et(1+xi)dtdx=ImΓ(s)0ts1etetixeπx+1dxdt=ImΓ(s)0ts1et0etixeπx1+eπx+etix1+eπxdxdt=ImΓ(s)n0(1)n0ts1et0ex(ti+π(n+1))+ex(ti+πn)dx=ImΓ(s)n0(1)n0ts1et(1π(n+1)+ti+1πnti)dt=1Γ(s)n0(1)n0tset(tπ2(n+1)2+t2tπ2n2+t2)dt=0tsetnZ(1)n1π2n2+t2dt=1Γ(s)0ts1etsinhtdt=2Γ(s)0ts1e2t1e2tdt=2Γ(s)n10ts1e2ntdt=12s1n11ns=η(s)ζ(s)

ちなみに、
多重対数関数にアベル·プラナの和公式を使うと(zを変数とする)
sin(sarctanxxlnz)(1+x2)s/2(e2πx1)dx
という積分が現れます。
z=1とすれば、今回といたような積分になりますしゼータ関数のアベル·プラナともなります。
分母を因数分解すれば
sin(sarctanx)(1+x2)s/2(eπx1)dx
も求められそうです。

問題13

0(1cosaxex)dxx1+r=Γ(1r)r(1+a2)r/2cos(rarctana)aR,0<r<1

0(1cosaxex)dxx1+r
=Re0(1e(ai+1)x)dxx1+r=Re(ai+1)001e(ai+1)xydydxxr=Re01Γ(1r)(ai+1)ry1rdy=Re(ai+1)rΓ(1r)r=Γ(1r)r(1+a2)r/2cos(rarctana)

これはラプラス変換を用いて解くと言うよりも計算の途中にラプラス変換があるだけですね。
特に話すことはないです。

問題14

cos(sarctan(ax))(1+x2)(1+(ax)2)s/2dx=π(1+a)sa,s>0

cos(sarctan(ax))(1+x2)(1+(ax)2)s/2dx
=ReΓ(s)11+x20ts1e(1+axi)tdtdx=1Γ(s)0ts1etcosaxt1+x2dxdt=πΓ(s)0ts1e(1+a)tdt=π(1+a)s

途中で有名な積分
cosax1+x2dx=πe|a|
が出てきましたね。
こういう解ける積分と
ラプラス変換をかけ合わせたものは
見た目複雑で中身は綺麗になったりして面白いんじゃないかなと思いました。

問題15

cos(r2ln(1+(ax)2))cosh(rarctan(ax))dx1+x2=πcos(rln(1+a))a>0,rR

最後に異形を持ってきました。
いままでの
三角関数の中に逆正接関数ではなく
双曲線関数の中に入っているものです。
ラプラス変換を知っていればこれも解けます。

0tu1eatcos(bt)dt
=Γ(u)(a2+b2)u/2cos(uarctanba)
u→iu
0tui1eatcos(bt)dt
=Γ(ui)eiu2ln(a2+b2)cosh(uarctanba)
cos(r2ln(1+a2x2))cosh(rarctan(ax))dx1+x2
=ReΓ(ir)11+x20tir1etcos(axt)dtdx=ReΓ(ir)0tir1etcos(axt)1+x2dxdt=ReπΓ(ir)0tir1e(1+a)tdt=Reπ(1+a)ir=πcos(rln(1+a))

こんな感じでした。
これらの多くは留数定理でも解けますがこのように綺麗に逐次積分に変換して解くこともできます。
良い経験になってくれれば嬉しいです。

投稿日:2023820
更新日:20231118
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

KBHM
KBHM
141
11134
怠惰

コメント

他の人のコメント

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