找回密码
 欢迎注册
查看: 23768|回复: 12

[讨论] 另一类差分方程求解

[复制链接]
发表于 2016-4-2 11:15:56 | 显示全部楼层 |阅读模式

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

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

×
在下面的问题中,mathe提供了非常绝妙的求解方案得到相应的结果:

http://bbs.emath.ac.cn/forum.php ... 71&fromuid=1455

我们现在讨论另一类差分方程:

1.\(x_{n+1}=\sin(x_n)\)

2.\(\sin(\tan(x_{n+1}))=\tan(\sin(x_n))\)

尽可能给出精确\(x_n\)关于\(n\)的渐近的表达式?

进一步给出更一般差分方程:(\(x_n\)有界)

\(x_{n+1}=ax_n+bx_n^2+cx_n^3+dx_n^4+ex_n^5+fx_n^6+gx_n^7+\dots\)

求\(x_n\)的渐近表达式?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2016-4-2 11:28:08 | 显示全部楼层
在<估计基础>中已求出第一问的渐近表示:

\(x_n=\sqrt{\frac{3}{n}}(1-\frac{3\ln(n)}{10n}-\frac{c}{2n}+\frac{\alpha\ln(n)^2+\beta\ln(n)+\gamma}{n^2}+O(\frac{\ln(n)^3}{n^3}))\)

\(\alpha=\frac{27}{200},\beta=\frac{9c}{20}-\frac{9}{50},\gamma=\frac{3c^2}{8}-\frac{3c}{10}+\frac{79}{100}\),且c为依赖于\(x_0\)的常数.

我们希望给出更精确的渐近结果包含c的具体渐近表达式?

点评

记$y_n=1/{x_n}$,$y_{n+1}=y_n+1/{6y_n}+7/{360y_n^3}+...$转化为那边的问题  发表于 2016-4-2 11:38
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-4-2 11:50:02 | 显示全部楼层
第二问同样设$y_n=1/{x_n}$,有$y_{n+1}=y_n-1/{30y_n^5}-41/{1890y_n^7}-88/{4725y_n^9}+...$
对应链接中k=6的问题

点评

不知\(c\)的求解公式如何推导?  发表于 2016-4-3 17:58

评分

参与人数 1金币 +32 贡献 +24 经验 +24 鲜花 +24 收起 理由
数学星空 + 32 + 24 + 24 + 24 神马都是浮云

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2016-4-2 12:34:44 | 显示全部楼层
当然,我们可以根据2#表达式的特点利用待定系数法求解得到第1问题更精确的渐近展开:

\(x_n=\sqrt{\frac{3}{n}}(1-\frac{3(\ln(n)+c)}{10n}+\frac{\frac{27}{200}(\ln(n)+c)^2-\frac{9}{50}(\ln(n)+c)+\frac{79}{700}}{n^2}+\frac{-\frac{27}{400}(\ln(n)+c)^3+\frac{27}{125}(\ln(n)+c)^2-\frac{1941}{7000}(\ln(n)+c)+\frac{411}{3500}}{n^3}+\frac{\frac{567}{16000}(\ln(n)+c)^4-\frac{1917}{1000}(\ln(n)+c)^3+\frac{1683}{4000}(\ln(n)+c)^2-\frac{7227}{17500}(\ln(n)+c)+\frac{1168171}{7480000}}{n^4}+\cdots)\)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2016-4-3 18:24:04 | 显示全部楼层
发现下面链接中方法不太适用此超越函数形式的差分方程:

http://bbs.emath.ac.cn/forum.php ... 71&fromuid=1455

例如:对于第二题

\(\sin(\tan(x_{n+1}))=\tan(\sin(x_n))\)

首先做级数展开,然后令\(x_n=\frac{1}{y_n}\)

1.\(y_{n+1}=y_n-\frac{1}{30y_n^5}-\frac{41}{1890y_n^7}-\frac{88}{4725y_n^9}-\frac{2477}{138600y_n^{11}}+\dots \)

然后左右各6次方后展开

\(y_{n+1}^6-y_n^6+\frac{1}{5}+\frac{41}{315y_n^2}+\frac{176}{1575y_n^4}+\frac{523}{5775y_n^6}+\dots=0\)

下一步无法做代换\(y_n^6=z_n\)???

并且\(c\)计算也同样存在问题?


毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-4-3 18:46:14 | 显示全部楼层
这个石子可以假设$y_n^6=n(1/5+\sum_{k=1}^{+\infty} {c_k}/{n^{k/3}})$
然后$y_n^2=n^{1/3}(1/5+\sum_{k=1}^{+\infty} {c_k}/{n^{k/3}})^{1/3}$
代入比较系数即可,结果应该比第一题更加简单一些

点评

对于分数指数形式的级数展开及运算还是麻烦哈.  发表于 2016-4-3 20:17
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-4-4 12:40:07 | 显示全部楼层
上面问题本质上就是展开$(1+1/n)^{-k/3}$,利用泰勒级数可以非常容易就展开了
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-4-4 15:31:31 | 显示全部楼层
可以得出$y_n^6=n(1/5+41/{210n^{1/3}}+176/{525n^{2/3}}+...)$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2016-4-6 09:58:56 | 显示全部楼层
数学星空 发表于 2016-4-2 11:28
在中已求出第一问的渐近表示:

\(x_n=\sqrt{\frac{3}{n}}(1-\frac{3\ln(n)}{10n}-\frac{c}{2n}+\frac{\a ...


你说的《阶的估计基础》是指哪本书

点评

潘承洞的<阶的估计基础>  发表于 2016-4-6 13:17
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-6-16 20:28 , Processed in 0.049333 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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