找回密码
 欢迎注册
查看: 25891|回复: 13

[讨论] n!+2^n

[复制链接]
发表于 2014-6-22 16:32:10 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
$T_0=2,T_1=3,T_2=6$

$T_n=(n+4)T_{n-1}-4nT_{n-2}+(4n-8)T_{n-3}$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-6-22 20:48:21 | 显示全部楼层
这个可以用母函数来做,或者说,Z变换。
设 $T_n$的母函数为 $f(x)$,则化简得:
\[f(x) -a_2x^2-a_1x-a_0 = (f'(x)-a_1)x^2 +5x(f(x)-a_1x-a_0)  - x^2(xf'(x)+2f(x)-2a_0) +4x^3(xf'(x)+f(x))\]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-6-22 21:13:36 | 显示全部楼层
设`T_n`的母函数为`G(x)=\sum_{n=0}^{\oo}T_nx^n`
于是有
$$\begin{align}  G(x)    &= T_0+T_1x+T_2x^2+ T_3x^3+ \cdots +T_{n-1}x^{n-1}+T_nx^n+ \cdots \tag{1} \\
\frac{(x^5G(x))'}{x^3} &= \mathrel{{\phantom{T_0}}}5T_0x+6T_1x^2+ 7T_2x^3 + \cdots +\quad(n+4)T_{n-1}x^n+ \cdots \tag{2} \\
-4x(x^2G(x))'              &=  \mathrel{{\phantom{T_0+T_1}}}-8T_0x^2- 12T_1x^3- \cdots - \qquad \enspace 4nT_{n-2}x^n - \cdots \tag{3} \\
4x^3(xG(x))'               &=  \mathrel{{\phantom{T_0+T_1x+T_2x^2+}}}4T_0x^3+ \cdots +(4n-8)T_{n-3}x^n+ \cdots \tag{4}
\end{align}$$
将`(2)+(3)+(4)-(1)`,并利用1楼的递推关系以及`T_0=2,T_1=3,T_2=6`,得到线性常微分方程$$(1-2x)^2x^2G(x)'+(1-2x)^2(x-1)G(x)=-4x^2+7x-2$$解得$$G(x)=\frac{c_1\text{e}^{-1/x}}{x}+\frac{\text{e}^{-1/x}\,\text{Ei}(1/x)}{x}+\frac{1}{1-2x}$$
其中,`\text{Ei}(x)`是指数积分函数。
因为初始条件为`\D\lim_{x->0^+}G(x)=T_0=2`,可知通解中的`c_1`为任意数均满足。取`c_1=0`有$$G(x)=\frac{\text{e}^{-1/x}\,\text{Ei}(1/x)}{x}+\frac{1}{1-2x} \quad(x>0)$$
根据指数积分在Re(z)>0上的级数定义可知$$\frac{\text{e}^{-1/x}\,\text{Ei}(1/x)}{x}=\sum_{n=0}^{\oo}n!x^n$$而$$\frac{1}{1-2x}=\sum_{n=0}^{\oo}2^nx$$故原递推关系的通项公式就是$$T_n=n!+2^n\quad(n \geqslant 0)$$
验证:
  1. Table[SeriesCoefficient[(E^(-1/x) ExpIntegralEi[1/x])/x + 1/(
  2.    1 - 2 x), {x, 0, n}, Assumptions -> x > 0], {n, 0, 10}]

  3. RecurrenceTable[{a[n] == (n + 4) a[n - 1] -
  4.     4 n a[n - 2] + (4 n - 8) a[n - 3], a[0] == 2, a[1] == 3,
  5.   a[2] == 6}, a, {n, 1, 10}]
复制代码

结果是都是:
{2,3,6,14,40,152,784,5168,40576,363392,3629824}

点评

@倪举鹏 非线性的递归方程大多没有显式解。即使是线性的,大部分也不一定有有限的显式解。  发表于 2014-6-23 11:22
厉害呀,这样的也可以算。a(1)=1,a(n+1)=a(n)+1/a(n)这个有没有解?  发表于 2014-6-22 22:32
非齐次线性常微分方程有通解公式,其中有关于因子exp(x)的积分,所以会出现非初等函数。这从杜哈梅尔(Duhamel)原理也可看出。  发表于 2014-6-22 21:38
汗,对对,一时犯2了  发表于 2014-6-22 21:25
是线性的,只不过系数不是常系数。你可以再看看线性微分方程的定义。  发表于 2014-6-22 21:19

评分

参与人数 2威望 +5 金币 +5 贡献 +5 经验 +3 鲜花 +5 收起 理由
fungarwai + 2 + 2 + 2 + 2
wayne + 3 + 3 + 3 + 3 + 3 赞一个!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2014-6-23 13:02:00 | 显示全部楼层
$a_n=(n+4)a_{n-1}-4na_{n-2}+4(n-2)a_{n-3}$
$a_n-na_{n-1}=4a_{n-1}-4(n-1)a_{n-2}-4a_{n-2}+4(n-2)a_{n-3}$
$a_n-na_{n-1}=4(a_{n-1}-(n-1)a_{n-2})-4(a_{n-2}-(n-2)a_{n-3})$
$b_n=a_n-na_{n-1},b_1=1,b_2=0$
$b_n=4b_{n-1}-4b_{n-2}$
$b_n-2b_{n-1}=2b_{n-1}-4b_{n-2}$
$c_n=b_n-2b_{n-1},c_2=-2$
$c_n=2c_{n-1}=-2^{n-1}$
$b_n=2b_{n-1}-2^{n-1}=-(n-2)2^{n-1}$
$a_n-na_{n-1}=-(n-2)2^{n-1}$
$a_n-2^n=n(a_{n-1}-2^{n-1})$
$a_n-2^n=n!$

点评

嗯,这个方法有取巧的成分,或者说很可能是逆向出来的。  发表于 2014-6-23 19:29
关键是关于n的变动系数给的非常恰当,以至于能凑成这样的数列,如果1楼中右边的n+4随便改成n+3,或者4n-8改成4n-7就没这样的结果。  发表于 2014-6-23 14:27
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2015-12-28 14:49:54 | 显示全部楼层
\(T_n-(n+4)T_{n-1}+4nT_{n-2}-(4n-8)T_{n-3}=0\)
\([E^3-(n+4)E^2+4nE-(4n-8)I]T_{n-3}=0\)
\([E^3-E^2(n+2)+E(4n-4)-(4n-8)I]T_{n-3}=0\)
\(8-4(n+2)+2(4n-4)-(4n-8)=0\)
\((E-I2)[E^2-En+I(2n-4)]T_{n-3}=0\)
\((n-2)^2-(n-2)n+2(n-2)=(n-2)(n-2-n+2)=0\)
\((E-I2)(E-I2)[E-I(n-2)]T_{n-3}=0\)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-22 16:36 , Processed in 0.026562 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表