24小时热门版块排行榜    

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

笃学明志

至尊木虫 (知名作家)

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

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

我想知道大侠们能否求出n—>infinite时S的极限值,我主要还是想了解通项的求法,我转换成二阶变系数方程后不知道怎么处理了?拼不出来函数?可能要变形,重新构造数列。谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一日红黑,终生红黑!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

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的回帖
普通回帖

笃学明志

至尊木虫 (知名作家)

怎么没有一点动静啊?!
一日红黑,终生红黑!
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

至尊木虫 (著名写手)

引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 笃学明志 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 17:22 by skyjudy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 基金不中,共勉 +3 eulota 2026-08-26 3/150 2026-08-26 16:17 by 路yyyy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 我不理解! +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见