1

愚直に計算する定積分の自作問題

200
0

問題

次の定積分を求めてください。
(1) log|x+1|+log|x1|x2dx
(2) log|x+1|log|x1|x2dx

2問とも広義積分です。
愚直に原始関数を求めて極限を計算して定積分を求めます。
(長々と計算し続けるので、飛ばし飛ばし読んでも構いません。)

(1)の答え

不定積分を求める

不定積分log|x+1|+log|x1|x2dxを求めます。
部分積分や部分分数分解により、
log|x+1|+log|x1|x2dx=log|x+1|+log|x1|x(1x)(1x+1+1x1)dx()=log|x+1|+log|x1|x+(1x11x+1)dx=log|x+1|+log|x1|x+log|x1|log|x+1|=(x1)log|x1|(x+1)log|x+1|x
であることが分かります(どうせ定積分に使うので積分定数は省略)。

積分区間を分ける

被積分関数にはx=1,0,1の定義できないところがあるので、
求める定積分は(,1),(1,0),(0,1),(1,)4つの積分範囲による広義積分の和になりますが、
被積分関数は偶関数なので、
log|x+1|+log|x1|x2dx=2(01log(x+1)+log(1x)x2dx+1log(x+1)+log(x1)x2dx)
となり、
01log(x+1)+log(1x)x2dx1log(x+1)+log(x1)x2dx2つを求めればいいことが分かります。

極限を計算して定積分を求める

1つ目の01log(x+1)+log(1x)x2dxは、
01log(x+1)+log(1x)x2dx=lima0limb1ablog(x+1)+log(1x)x2dx=lima0limb1[(x1)log(1x)(x+1)log(x+1)x]ab=limb1(b1)log(1b)(b+1)log(b+1)blima0(a1)log(1a)(a+1)log(a+1)a

ここで、
limb1(b1)log(1b)(b+1)log(b+1)b=limb1(b1)log(1b)2log(2)=limb1log(1b)(b1)12log(2)=limb1(1b)1(b1)22log(2)()=limb1(1b)2log(2)=2log(2)

lima0(a1)log(1a)(a+1)log(a+1)a=lima0(log(1a)log(a+1))()=0

(1-4)(1-5)より、(1-2)は
limb1(b1)log(1b)(b+1)log(b+1)blima0(a1)log(1a)(a+1)log(a+1)a=2log(2)
となります。

2つ目の1log(x+1)+log(x1)x2dxも同様に計算します。
1log(x+1)+log(x1)x2dx=lima1limbablog(x+1)+log(x1)x2dx=lima1limb[(x1)log(x1)(x+1)log(x+1)x]ab=limb(b1)log(b1)(b+1)log(b+1)blima1(a1)log(a1)(a+1)log(a+1)a

ここで、
limb(b1)log(b1)(b+1)log(b+1)b=limb(log(b1)log(b+1)log(b1)blog(b+1)b)=limb(log(11/b1+1/b)log(b1)blog(b+1)b)=limb(log(b1)blog(b+1)b)=limb(1b11b+1)()=0

lima1(a1)log(a1)(a+1)log(a+1)a=lima1(a1)log(a1)2log(2)=lima1log(a1)(a1)12log(2)=lima1(a1)1(a1)22log(2)()=lima1(a1)2log(2)=2log(2)

(1-8)(1-9)より、(1-7)は
limb(b1)log(b1)(b+1)log(b+1)blima1(a1)log(a1)(a+1)log(a+1)a=2log(2)
となります。

従って、(1-6)(1-10)から(1-2)は
2(01log(x+1)+log(1x)x2dx+1log(x+1)+log(x1)x2dx)=2(2log(2)+2log(2))=0
となり、求めたい値は0であることが分かりました。

答え:log|x+1|+log|x1|x2dx=0

(2)の答え

積分区間を分ける

求める定積分は(,1),(1,0),(0,1),(1,)4つの積分範囲による広義積分の和になりますが、
被積分関数は偶関数なので、
log|x+1|log|x1|x2dx=2(01log(x+1)log(1x)x2dx+1log(x+1)log(x1)x2dx)=2(lima0limb1ablog(x+1)log(1x)x2dx+lima1limbablog(x+1)log(x1)x2dx)
となり、2つの広義積分を求めればいいことが分かります。

不定積分を求めようとするが・・・

部分積分と部分分数分解から、
log|x+1|log|x1|x2dx=log|x+1|log|x1|x(1x)(log|x1|x+1+log|x+1|x1)dx()=log|x+1|log|x1|x+(log|x1|x(x+1)+log|x+1|x(x1))dx=log|x+1|log|x1|x+(log|x1|xlog|x1|x+1+log|x+1|x1log|x+1|x)dx
なので、(2-1)は
2(lima0limb1ablog|x+1|log|x1|x2dx+lima1limbablog|x+1|log|x1|x2dx)=2(lima0limb1([log(x+1)log(1x)x]ab+ab(log(1x)xlog(1x)x+1+log(x+1)x1log(x+1)x)dx)(2-2)+lima1limb([log(x+1)log(x1)x]ab+ab(log(x1)xlog(x1)x+1+log(x+1)x1log(x+1)x)dx))
となります。

残念ながら(1)の問題とは違い、(2-2)の積分は、不定積分が初等関数で表せないらしいです。
被積分関数を見ると、以下の6つの積分をしないといけません。
log(1x)xdx (0<x<1)log(1x)x+1dx (0<x<1)log(x+1)x1dx (0<x<1,1<x<)log(x+1)xdx (0<x<1,1<x<)log(x1)xdx (1<x<)log(x1)x+1dx (1<x<)
これらを積分するために、「2重対数関数(Dilogarithm)」という特殊関数を使います。

積分結果を2重対数関数を使って表す

2重対数関数

2重対数関数Li2(x)とは、Li2(x)=n=1xnn2という冪級数で定義される関数です。
これは1x1で定義できます。

しかし、普通はこの冪級数を積分を使って0xlog(1t)tdtと表し、
これを改めてLi2(x)と定義します。
こうすることで、1x1のときの値はそのままに、
実数の範囲では定義域を<x1に拡張できます(解析接続というやつです)。

では(2-3)の①~⑥を求めてみましょう。(以下、積分定数は省略します。)

一番簡単なのは①で、定義からlog(1x)xdx=Li2(x)なので、
log(1x)xdx=Li2(x)です。
これ以降は、被積分関数をlog(1x)xの形にすることを目標にして置換などします。

次は④です。x=tと置換すればlog(1x)xの形になりそうなのが分かると思います。
dx=(1)dtであることに注意して、
log(x+1)xdx=log(1t)t(1)dt=log(1t)tdt=Li2(t)=Li2(x)
となります。

次は②をやってみましょう。分母のx+1が扱いにくいので、x+1=tで置換します(dx=dt)。
log(1x)x+1dx=log(2t)tdt
分子がlog(1x)の形ではなくなったのですが、log(2t)=log(2)+log(1t2)と分けて、t=2uと置換すれば大丈夫です。
log(1x)x+1dx (0<x<1)=log(2t)tdt(x+1=t,dx=dt)=log(2)+log(1t2)tdt=log(2)log|t|+log(1t2)tdt=log(2)log|t|+log(1u)2u2du(t=2u,dt=2du)=log(2)log|t|+log(1u)udu=log(2)log|t|Li2(u)=log(2)log|x+1|Li2(x+12)=log(2)log(x+1)Li2(x+12)(0<x<1)

同様に③を計算します。④と同様の置換をすると②と同じになります。
log(x+1)x1dx (0<x<1,1<x<)=log(1t)t+1dt(x=t,dx=dt)=log(2)log|t+1|Li2(t+12)=log(2)log|1x|Li2(1x2)

⑤が割と大変です。なかなかlog(1x)xの形になってくれません。
一旦、部分積分すると目的の形が見えてきます。
log(x1)xdx (1<x<)=log(x)log(x1)log(x)x1dx()=log(x)log(x1)log(1t)tdt(x=1t,dx=dt)=log(x)log(x1)+Li2(t)=log(x)log(x1)+Li2(1x)

最後に⑥です。⑤と同様に部分積分で③に帰着します。
log(x1)x+1dx (1<x<)=log(x+1)log(x1)log(x+1)x1dx()=log(x+1)log(x1)log(2)log(x1)+Li2(1x2)=log(x+12)log(x1)+Li2(1x2)

以上、(2-3)の①~⑥により、(2-2)は
2(lima0limb1([log(x+1)log(1x)x]ab+ab(log(1x)xlog(1x)x+1+log(x+1)x1log(x+1)x)dx)+lima1limb([log(x+1)log(x1)x]ab+ab(log(x1)xlog(x1)x+1+log(x+1)x1log(x+1)x)dx))=2(lima0limb1([log(x+1)log(1x)xLi2(x)log(2)log(x+1)+Li2(x+12)+log(2)log(1x)Li2(1x2)+Li2(x)]ab)+lima1limb([log(x+1)log(x1)x+log(x)log(x1)+Li2(1x)log(x+12)log(x1)+log(2)log(x1)2Li2(1x2)+Li2(x)]ab))=2(limb1(log(b+1)log(1b)bLi2(b)log(2)log(b+1)+Li2(b+12)+log(2)log(1b)Li2(1b2)+Li2(b))lima0(log(a+1)log(1a)aLi2(a)log(2)log(a+1)+Li2(a+12)+log(2)log(1a)Li2(1a2)+Li2(a))+limb(log(b+1)log(b1)b+log(b)log(b1)+Li2(1b)log(b+12)log(b1)+log(2)log(b1)2Li2(1b2)+Li2(b))(2-4)lima1(log(a+1)log(a1)a+log(a)log(a1)+Li2(1a)log(a+12)log(a1)+log(2)log(a1)2Li2(1a2)+Li2(a)))
となります。もうてんやわんやですね。

これから利用する2重対数関数の特殊値と恒等式

(2-4)は以下の4つの極限の計算でできています。
limb1(log(b+1)log(1b)bLi2(b)log(2)log(b+1)+Li2(b+12)+log(2)log(1b)Li2(1b2)+Li2(b))lima0(log(a+1)log(1a)aLi2(a)log(2)log(a+1)+Li2(a+12)+log(2)log(1a)Li2(1a2)+Li2(a))limb(log(b+1)log(b1)b+log(b)log(b1)+Li2(1b)log(b+12)log(b1)+log(2)log(b1)2Li2(1b2)+Li2(b))lima1(log(a+1)log(a1)a+log(a)log(a1)+Li2(1a)log(a+12)log(a1)+log(2)log(a1)2Li2(1a2)+Li2(a))
これらを計算するのにLi2(x)の特殊値と、xにおける挙動について確認しましょう。

まず、Li2(0)は明らかに0です。
次に、Li2(1)は定義からLi2(1)=n=11n2なので、「バーゼル問題」と同じで、値はπ26です。「バーゼル問題」の解き方はネットにいっぱいあるので調べてみてください。
続いてLi2(1)です。先程のLi2(1)を利用して、
Li2(1)=n=1(1)nn2=2n=11(2n)2n=11n2=12n=11n2n=11n2=12n=11n2=π212
と計算できます。
正の数tに対し、Li2(t1)を考えます。
2重対数関数の定義と微分積分学の基本定理、合成関数の微分公式から、
ddtLi2(t1)=log(1+t1)t1t2=log(t+1)tlog(t)t
になりますので、t1からx(>0)までで定積分することで、
Li2(x1)Li2(1)=1x(log(t+1)tlog(t)t)dt=[Li2(t)(log(t))22]1x=Li2(x)(log(x))22+Li2(1)
すなわち、
Li2(x)=Li2(x1)(log(x))22π26
この等式が任意の正の数xで成り立つのです。
この等式の何が嬉しいかというと、Li2(x)xで発散しますが、Li2(x1)0になります。代わりにlog(x)が発散しますが、2重対数関数よりも普通の対数関数の方が扱いやすいです。

纏めるとこうなります。

  • Li2(0)=0
  • Li2(1)=π26
  • Li2(1)=π212
  • Li2(x)=Li2(1x)(log(x))22π26(x>0)

極限を計算して定積分を求める

上で纏めた式とロピタルの定理などの道具を使って(2-5)の①~④を計算します。

limb1(log(b+1)log(1b)bLi2(b)log(2)log(b+1)+Li2(b+12)+log(2)log(1b)Li2(1b2)+Li2(b))=(log(2))2π212+limb1(log(b+1)log(1b)b+log(2)log(1b))=(log(2))2π212+limb1(log(2)b1log(b+1)(log(1b))1)=(log(2))2π212+limb1(b2log(b+1)b1(b+1)1(log(1b))2(1b)1)()=(log(2))2π212+(log(2)12)limb1((log(1b))2(1b)1)=(log(2))2π212+(log(2)12)limb1(2log(1b)(1b)1(1b)2)()=(log(2))2π212(2log(2)1)limb1(log(1b)(1b)1)=(log(2))2π212(2log(2)1)limb1((1b)1(1b)2)()=(log(2))2π212(2log(2)1)limb1((1b))=(log(2))2π212

lima0(log(a+1)log(1a)aLi2(a)log(2)log(a+1)+Li2(a+12)+log(2)log(1a)Li2(1a2)+Li2(a))=lima0(log(a+1)log(1a)a)=lima0(((a+1)1log(1a)log(a+1)(1a)1))()=0

③はちょっと大変です。
limb(log(b+1)log(b1)b+log(b)log(b1)+Li2(1b)log(b+12)log(b1)+log(2)log(b1)2Li2(1b2)+Li2(b))=limb(log(b+1)log(b1)b+log(b)log(b1)Li2(1b1)(log(b1))22π26log(b+12)log(b1)+log(2)log(b1)2(Li2(2b1)(log((b1)/2))22π26)Li2(1b)(log(b))22π26)((26))=(log(2))2+limb(log(b+1)log(b1)b+log(b)log(b1)log(b+1)log(b1)+(log(b1))22(log(b))22)=(log(2))2+limb(log(b+1)log(b1)b+log(b)log(b+1)(log(b1))1+log(b1)log(b)2(log(b1)+log(b))1)=(log(2))2+limb(log(b1)b+1log(b+1)b1+b1(b+1)1(log(b1))2(b1)1+(b1)1b14(log(b1)+log(b))2((b1)1+b1))()=(log(2))2+limb(log(b1)b+1log(b+1)b1(b1)(log(b1))2b(b+1)(log(b1)+log(b))24(2b1))=(log(2))2+limb(1b11b+1(log(b1))2+2log(b1)2b+1(log(b1)+log(b))((b1)1+b1)4)()=(log(2))2+limb((log(b1))22b+12log(b1)2b+1(log(b1)+log(b))(2b1)4b(b1))=(log(2))2+limb(log(b1)b11b1((b1)1+b1)(2b1)+2(log(b1)+log(b))4(2b1))()=(log(2))2+limb(log(b1)b11b1((b1)1+b1)4log(b1)2(2b1)log(b)2(2b1))=(log(2))2+limb(log(b1)b1log(b1)2(2b1)log(b)2(2b1))=(log(2))2+limb(1b114(b1)14b)()=(log(2))2

lima1(log(a+1)log(a1)a+log(a)log(a1)+Li2(1a)log(a+12)log(a1)+log(2)log(a1)2Li2(1a2)+Li2(a))=π212+lima1(log(a+1)log(a1)a+log(a)log(a1)log(a+12)log(a1)+log(2)log(a1))=π212+lima1(a1log(a+1)+log(a)log(a+1)+log(4)(log(a1))1)=π212+lima1(a2log(a+1)a1(a+1)1+a1(a+1)1(log(a1))2(a1)1)()=π212log(2)lima1((log(a1))2(a1)1)=π212log(2)lima1(2log(a1)(a1)1(a1)2)()=π2122log(2)lima1(log(a1)(a1)1)=π2122log(2)lima1((a1)1(a1)2)()=π2122log(2)lima1((a1))=π212

以上、(2-5)の①~④より、(2-4)は
2(limb1(log(b+1)log(1b)bLi2(b)log(2)log(b+1)+Li2(b+12)+log(2)log(1b)Li2(1b2)+Li2(b))lima0(log(a+1)log(1a)aLi2(a)log(2)log(a+1)+Li2(a+12)+log(2)log(1a)Li2(1a2)+Li2(a))+limb(log(b+1)log(b1)b+log(b)log(b1)+Li2(1b)log(b+12)log(b1)+log(2)log(b1)2Li2(1b2)+Li2(b))lima1(log(a+1)log(a1)a+log(a)log(a1)+Li2(1a)log(a+12)log(a1)+log(2)log(a1)2Li2(1a2)+Li2(a)))=2((log(2))2π212+(log(2))2+π212)=0
となり、求めたい値は0であることが分かりました。

答え:log|x+1|log|x1|x2dx=0

おわりに

log|x+1|+log|x1|x2dx=0
log|x+1|log|x1|x2dx=0
となり、足し算も掛け算も0になっちゃいました。自分で計算してて予想外な結果でした。

あと愚直に計算しまくって疲れました。。。
たぶんもっと効率のいい方法があると思ってるので、何か知ってる人は教えてください。

以上、読んでいただきありがとうございました。

投稿日:2022429
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

SunPillar
SunPillar
145
28827
自由気儘な社会人。 Twitterアカウントはフォローしてもそんなに数学ツイートはないので注意。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 問題
  2. (1)の答え
  3. (2)の答え
  4. おわりに