24小时热门版块排行榜    

查看: 2032  |  回复: 7
本帖产生 1 个 数学EPI ,点击这里进行查看

笃学明志

至尊木虫 (知名作家)

[求助] 二阶变系数线性方程 求助递推数列通项求解

在验证本学科几率效应的时候遇到一个极限问题,该数列是这样的:

我想知道大侠们能否求出n—>infinite时S的极限值,我主要还是想了解通项的求法,我转换成二阶变系数方程后不知道怎么处理了?拼不出来函数?可能要变形,重新构造数列。谢谢!
回复此楼
一日红黑,终生红黑!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笃学明志

至尊木虫 (知名作家)

怎么没有一点动静啊?!
一日红黑,终生红黑!
2楼2012-09-23 17:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笃学明志

至尊木虫 (知名作家)

不是沉了吧?大神救救我啊
一日红黑,终生红黑!
3楼2012-09-25 13:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笃学明志

至尊木虫 (知名作家)


闲来无事  根据拟合得出
bingo  
一日红黑,终生红黑!
4楼2013-04-24 09:36:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 笃学明志 at 2013-04-24 09:36:33
闲来无事  根据拟合得出
bingo  

令f(x) = S_1 + S_2 x + S_3 x^2 + ... = Sum_{n >=1} S_n* x^{n-1}.

则 f'(x)= S_2 + 2*S_3 x + 3* S_4 x^2 +...= Sum_{n>=2} (n-1)* x^{n-2} * S_n

你的级数 当n>=3时,
(n-1) * S_n - (n-2) * S_{n-1} = 2 * S_{n-2}. 两边同时乘以x^{n-2}, 求和,
请直接验证
(f' - S_2) - x* f' = 2 * x * f.
因此 f' / f = 2x / (1-x).  积分得 f(x)= C* e^{-2x} / (1-x)^2.
带入 f(0)= S_1=1, 得到 C=1, 即 f(x)= e^{-2x} / (1-x)^2.

利用 e^{-2x}= Sum_{k>=0} (-2x)^k / k!
  (1-x)^{-2} = Sum_{n>=0} (-2 choose n) (-x)^n
得到
S_{n+1} = Sum_{k=0}^n (-2)^k / k! * (-1)^{n-k} * (-2 choose n-k).

关于进一步的分析, 我无能为力了。
We_must_know. We_will_know.
5楼2013-07-29 13:35:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 笃学明志 at 2013-04-24 09:36:33
闲来无事  根据拟合得出
bingo  

继续楼下的分析。

设f(x)= Sum_{n>=1} S_n *x^{n-1} = e^{-2x} / (1-x)^2.
楼主猜测 S_n = (n+2) e^{-2}, 所以我们避开通项公式,看看可否证明楼主的猜想。

比较 (1-x)^2 * f(x) = e^{-2x} 的 x^n 的系数, 得到
S_{n+1} - 2* S_n + S_{n-1} = (-2)^n / n!.
该式子对所有的n>=0 成立,如果我们取 S_0=S_{-1}=0的话.

故
S_{n+1} - S_n =Sum_{k=0}^n (-2)^k/ k! 趋于  e^{-2},
所以 楼主的拟合 是惊人的准确。

事实上,我个人认为 S_{n+1} = S_n + Sum_{k=0}^n (-2)^k / K!
会比通项公式 更好用,也更直观。
We_must_know. We_will_know.
6楼2013-07-29 15:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
笃学明志(sweety代发): 金币+10, 鼓励虫友认真的回复和解答 2014-10-14 21:23:38
sweety: 数学EPI+1 2014-10-14 21:23:44
引用回帖:
6楼: Originally posted by hank612 at 2013-07-29 15:18:57
继续楼下的分析。

设f(x)= Sum_{n>=1} S_n *x^{n-1} = e^{-2x} / (1-x)^2.
楼主猜测 S_n = (n+2) e^{-2}, 所以我们避开通项公式,看看可否证明楼主的猜想。

比较 (1-x)^2 * f(x) = e^{-2x} 的 x^n 的系 ...

我按楼主的程序画了图, 这根本就是一条直线, 我还以为程序出错了呢. 把它在0-100放大,就会发现
除了S1=1, S2=0, S3=1, S4=2/3, S5=1 以后都是单调递增的, 几乎是直线.

由于 S_{n+1} - S_n =Sum_{k=0}^n (-2)^k/ k!,
所以 S_{n+1}  = Sum_{m=0}^n Sum_{k=0}^m (-2)^k/ k!
               = Sum_{k=0}^n  (-2)^k/ k! * (n+1-k)
= (n+1)*  Sum_{k=0}^n  (-2)^k/ k! -  Sum_{k=0}^n  (-2)^k/ (k-1)!
---> (n+1)*e^{-2} - (-2)* e^{-2}
       = (n+3)* e^{-2}
这就是楼主断言的 S_n ~ (n+2) * e^{-2}.
We_must_know. We_will_know.
7楼2013-07-30 01:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

引用回帖:
6楼: Originally posted by hank612 at 2013-07-29 15:18:57
继续楼下的分析。

设f(x)= Sum_{n>=1} S_n *x^{n-1} = e^{-2x} / (1-x)^2.
楼主猜测 S_n = (n+2) e^{-2}, 所以我们避开通项公式,看看可否证明楼主的猜想。

比较 (1-x)^2 * f(x) = e^{-2x} 的 x^n 的系 ...

以前曾经得到母函数的表达式, 即.

由于的精确表达式颇为复杂,况且楼主只需要考察是否有, 我们可以直指问题,避开麻烦.

考虑在x=0处的Taylor展开式, 但我们来个明修栈道,暗度陈仓. 由于在x=1是非常容易展开的: .

注意到
以及显然的积分式对x 展开得到的


于是, 几乎没有经过任何运算, 我们就得到 等于
明显绝对值小于 趋于0. 这就是楼主看到的.

我翻老贴只是想说明, 想问题最好直指本心. 那些兜七兜八把大伙都快绕晕了的, 大体上是竺老师说的"小己玄思, 凿枘不容". 数学本该是简单而深刻的.
We_must_know. We_will_know.
8楼2016-07-01 01:59:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 笃学明志 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见