欢迎光临散文网 会员登陆 & 注册

[快乐数学]一个"有趣"的公式

2023-08-03 07:00 作者:名浮半生  | 我要投稿

咳咳,少废话直接开始吧。

1.引子

这一次的讨论要从一个求和开始。

1%2B%5Cfrac%7B1%7D%7B1%C3%973%7D%20%2B%5Cfrac%7B1%7D%7B1%C3%973%C3%975%7D%20%2B%5Cfrac%7B1%7D%7B1%C3%973%C3%975%C3%977%7D%20%2B....

这个样貌一看就是一个级数的和。

抓到一个级数 ,应当先判断它收不收敛。

不过我都拿出来讨论了,当然是收敛的啦。这瓜保熟!

额,先说明一下符号。

n!!的“!”不是感叹号,两个“!”连用是双阶乘。

意思的话,你看前面给的式子应该也看懂了。

n!!指的是,不超过n且与n同奇偶性的正整数的乘积。

这类问题我想很多人都见过。常规的操作是构造幂级数然后求幂级数和函数的和。

即构造h(x)%3D%5Cfrac%7Bx%5E1%7D%7B1%7D%20%2B%5Cfrac%7Bx%5E3%7D%7B1%C3%973%7D%20%2B%5Cfrac%7Bx%5E5%7D%7B1%C3%973%C3%975%7D%20%2B%5Cfrac%7Bx%5E7%7D%7B1%C3%973%C3%975%C3%977%7D%20%2B...

求幂级数的和函数呢,往往利用逐项求积或者逐项求导的方法。

这里我们要做的应该是把分母消掉,所以应该逐项求导。

h'(x)%3D1%20%2B%5Cfrac%7Bx%5E2%7D%7B1%7D%20%2B%5Cfrac%7Bx%5E4%7D%7B1%C3%973%7D%20%2B%5Cfrac%7Bx%5E6%7D%7B1%C3%973%C3%975%7D%20%2B%5Cfrac%7Bx%5E8%7D%7B1%C3%973%C3%975%C3%977%7D%20%2B....

哦,这里我们可以看到,如果我们在h'(x)里提出一个x,那么h(x)就会再现。

h'(x)%3D1%20%2Bx(%5Cfrac%7Bx%5E1%7D%7B1%7D%20%2B%5Cfrac%7Bx%5E3%7D%7B1%C3%973%7D%20%2B%5Cfrac%7Bx%5E5%7D%7B1%C3%973%C3%975%7D%20%2B%5Cfrac%7Bx%5E7%7D%7B1%C3%973%C3%975%C3%977%7D%20%2B....)%3D1%2Bxh(x)

ok,这样的话,我们h(x)就是上面那个微分方程的解了。

那么问题来了,这个微分方程怎么解呢?

首先判断类型,这个微分方程是一个一阶非齐次线性常微分方程。

没错吧。这里的h(x)就是y,P(x)=-x,Q(x)=1

而这类微分方程的求解我们是有公式的。

就这玩意。

但是公式里有待定的数C(大多数微分方程的解都有很多很多个,所以才会有待定的数C)

但对我们而言,C是能够被确定的。

因为h(0)你看一下,是不是0啊?

所以,把上面的条件代入就可以了。

结果是

h(x)%3De%5E%5Cfrac%7Bx%5E2%7D%7B2%7D%5Cint_%7B0%7D%5E%7Bx%7D%20e%5E%5Cfrac%7B-t%5E2%7D%7B2%7D%20dt%20

你问我为什么不把积分号里面的玩意积出来?

拜托,那个结果无法用初等函数表示的。

继续咱们的思路,现在当x=1的时候h(1)就是答案了。这里你会发现h(1)%3De%5E%5Cfrac%7B1%7D%7B2%7D%20%5Cint_%7B0%7D%5E%7B1%7D%20e%5E%5Cfrac%7Bt%5E2%7D%7B2%7D%20dt

看过

的朋友或许会对这个结果表示熟悉。

因为这里又出现钟形曲线。不过对这期专栏来说,这个结果并不够有趣。

2.一个“有趣”连分式

第二部分咱们来考量一个和前面的微分方程有点像的玩意。

f'(x)%3Dxf(x)-1

为了方便,我们可以令y=f(x)

小标题提到了“连分式”,那我们就来试试用它构造一个连分式。

和前面差不多,我们对y“逐项求导”。

y'=xy-1

y''=xy'+1y

y'''=xy''+2y

y''''=xy'''+3y

......

这里你或许看不出什么,那如果我两边同时除以x边上的那个导数呢?

y'/y=x-1/y

y''/y=x+1y/y'

y'''/y''=x+2y'/y''

y''''/y'''=x+3y''/y'''

......

发现了吗?他们头尾相等了。

这不就可以迭代了。

第一个式子里的y格格不入,而且它可以承担前面类似和函数的作用。

所以我们变形一下全部的式子。把最右边那一列移到等号左边,第一列移到右边,然后分子分母颠倒以把y单独拎出来。

这样的话就是

y%3D%5Cfrac%7B1%7D%7Bx-%5Cfrac%7By'%7D%7By%7D%20%7D%20

%5Cfrac%7By'%7D%7By%7D%20%3D%5Cfrac%7B-1%7D%7Bx-%5Cfrac%7By''%7D%7By'%7D%20%7D%20

%5Cfrac%7By''%7D%7By'%7D%20%3D%5Cfrac%7B-1%7D%7Bx-%5Cfrac%7By'''%7D%7By''%7D%20%7D%20

%5Cfrac%7By'''%7D%7By''%7D%20%3D%5Cfrac%7B-1%7D%7Bx-%5Cfrac%7By''''%7D%7By'''%7D%20%7D%20

......

于是我们就可以迭代了。

y%3D%5Cfrac%7B1%7D%7Bx-%5Cfrac%7B-1%7D%7Bx-%5Cfrac%7B-2%7D%7Bx-%5Cfrac%7B-3%7D%7Bx-....%7D%20%7D%20%7D%7D%20

两个负号抵消就是

y%3D%5Cfrac%7B1%7D%7Bx%2B%5Cfrac%7B1%7D%7Bx%2B%5Cfrac%7B2%7D%7Bx%2B%5Cfrac%7B3%7D%7Bx%2B....%7D%20%7D%20%7D%7D%20

(你可以试着倒过来,由这个连分式得到前面的微分方程。如果你有什么自然的方法也欢迎发到评论区)

3.所以。。。"有趣"在哪里?

啊对啊,所以有趣在哪里呢?

前面那个微分方程我们依然能解的对吧。

解出来的结果有点恐怖哦。

f(x)%3De%5E%5Cfrac%7Bx%5E2%7D%7B2%7D%5Cint_%7Bx%7D%5E%7B%2B%E2%88%9E%7D%20e%5E%5Cfrac%7B-t%5E2%7D%7B2%7D%20dt%20

如果我们把h(x)和f(x)加起来,那么积分部分刚好就拼在一起了。

而这个积分在[0,+∞)的结果是可以利用转成二重积分算出来的。

结果是%5Cint_%7B0%7D%5E%7B%2B%E2%88%9E%7De%5E%5Cfrac%7B-t%5E2%7D%7B2%7D%20dt%3D%5Csqrt%7B%5Cfrac%7B%CF%80%7D%7B2%7D%20%7D%20%20

那么就有

h(x)%2Bf(x)%3De%5E%5Cfrac%7Bx%5E2%7D%7B2%7D%20%5Csqrt%7B%5Cfrac%7B%CF%80%7D%7B2%7D%20%7D%20

当x=1时,就有

h(1)%2Bf(1)%3De%5E%5Cfrac%7B1%7D%7B2%7D%20%5Csqrt%7B%5Cfrac%7B%CF%80%7D%7B2%7D%20%7D%20

而e的1/2次方就是根号e

h(1)和f(1)我们前面有结论,代入得

1%2B%5Cfrac%7B1%7D%7B1%C3%973%7D%20%2B%5Cfrac%7B1%7D%7B1%C3%973%C3%975%7D%20%2B%5Cfrac%7B1%7D%7B1%C3%973%C3%975%C3%977%7D%20%2B...%2B%5Cfrac%7B1%7D%7B1%2B%5Cfrac%7B1%7D%7B1%2B%5Cfrac%7B2%7D%7B1%2B%5Cfrac%7B3%7D%7B1%2B....%7D%20%7D%20%7D%7D%20%3D%5Csqrt%7B%5Cfrac%7Be%CF%80%7D%7B2%7D%20%7D%20

这个式子是由拉马努金(又是那个恐怖如斯的男人)提出的。

这个式子既有级数又有连分数,还有πe根号,属实是把各个要素凑齐了。


[快乐数学]一个"有趣"的公式的评论 (共 条)

分享到微博请遵守国家法律