数学研发论坛

 找回密码
 欢迎注册
查看: 301|回复: 18

[求助] 用mma怎么编程

[复制链接]
发表于 2023-1-13 20:57:41 | 显示全部楼层 |阅读模式

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

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

x
20230113141431(1).png
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-13 22:25:48 | 显示全部楼层
  1. Limit[n (n*RSolve[{a[1] == 1, a[n + 1] == Log[1 + a[n]]}, a[n], n] - 2)/Log[n], n -> Infinity]
复制代码

没解出来
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-14 00:46:03 | 显示全部楼层
本帖最后由 lsr314 于 2023-1-14 01:06 编辑

假设极限是r,将a(n)用n和r表示出来,由于x=1/n→0,我们可以将a(n+1)-a(n)在x=0处展开为级数:
  1. Series[r Log[1+n]/(1+n)^2+2/(1+n)- Log[1+r Log[n]/n^2+2/n]/.n->1/x,{x,0,3}]
复制代码

结果是(-2/3+r) x^3+O[x]^4,分析x^3的系数,极限应该是r=2/3.

评分

参与人数 2威望 +18 金币 +18 贡献 +18 经验 +18 鲜花 +18 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 赞一个!
uk702 + 10 + 10 + 10 + 10 + 10 学习了!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-1-14 11:54:55 | 显示全部楼层

你好,请问先生这个程序能否改进一下?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-14 13:49:03 | 显示全部楼层
本帖最后由 uk702 于 2023-1-14 15:18 编辑

确实是 2/3 ,我用的是类似的方法 。


  1. Series[2/(n + 1) + r Log[n + 1] / (n + 1)^2, {n, Infinity, 4}]
复制代码

得: \( \frac{2}{n}+\frac{r \log (n)-2}{n^2}+\frac{-2 r \log (n)+r+2}{n^3}+\frac{6 r \log (n)-5 r-4}{2 n^4}+O\left(\left(\frac{1}{n}\right)^5\right) \)

  1. Series[Log[1 + 2/n + r Log[n] /n^2], {n, Infinity, 4}]
复制代码

得: \( \frac{2}{n}+\frac{r \log (n)-2}{n^2} + \frac{-6 r \log (n) + 8}{3 n^3}+\frac{-r^2 \log ^2(n)+8 r \log (n)-8}{2 n^4}+O\left(\left(\frac{1}{n}\right)^5\right) \)

因此, \( \frac{1}{n^3} \) 系数应匹配,就必须有 2 + r - 2 r Log[n] = (8 - 6 r Log[n])/3, ∴ 2 + r = 8/3, 得 r = 2/3

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-1-14 16:08:21 | 显示全部楼层
uk702 发表于 2023-1-14 13:49
确实是 2/3 ,我用的是类似的方法 。

以2楼的编程为基础,请问2楼的程序怎么改进

点评

我是没那个能耐,MMA 有没有这个能耐,不知。  发表于 2023-1-14 16:42
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-14 17:02:37 | 显示全部楼层
本帖最后由 uk702 于 2023-1-14 17:04 编辑

3# 的方法就是大招。更进一步,似有
\[
a_n = \frac{2}{n} + \frac{2 Log[n]}{3 n^2} + \frac{2 Log[n]^2 - 2 Log[n]}{9 n^3} + o(n^{-3})
\]


  1. Series[(2/9  Log[1 + n]^2 - 2/9  Log[1 + n])/(1 + n)^3 +  2/3 Log[1 + n]/(1 + n)^2 + 2/(1 + n)  -  Log[1 + (2/9  Log[n]^2 - 2/ 9 Log[n])/n^3 + 2/3 Log[n]/n^2 + 2/n], {n, Infinity, 5}]
复制代码


结果:
\[
\frac{1}{9 n^4}+\frac{20 \log ^3(n)-80 \log ^2(n)+90 \log (n)-69}{135 n^5}+O\left(\left(\frac{1}{n}\right)^6\right)
\]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-14 22:35:31 | 显示全部楼层
  1. ClearAll["Global`*"]; Limit[
  2. RSolve[{a[1] == 1, a[n + 1] == a[n] + 1/(n + 1)^4}, a[n], n],
  3. n -> Infinity]
复制代码


$\underset{n\to \infty }{\text{lim}}\{a(n)\to \frac{1}{90} \left(\pi ^4-15 \psi ^{(3)}(n+1)\right)}$

点评

这个程序可以再改进一下吗  发表于 2023-1-14 23:12
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-1-14 22:51:14 | 显示全部楼层
本帖最后由 笨笨 于 2023-1-14 23:15 编辑
northwolves 发表于 2023-1-14 22:35
$%underset{n\to \infty }{\text{lim}}\{a(n)\to \frac{1}{90} \left(\pi ^4-15 \psi ^{(3)}(n+1)\righ ...


请问2楼的程序可以再改进吗,还是没有直接得到答案啊2/3
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-1-14 23:15:48 | 显示全部楼层
笨笨 发表于 2023-1-14 22:51
请问2楼的程序可以再改进吗,还是没有直接得到答案啊2/3

似乎有的能解出来,有的不能

点评

谢谢,很遗憾  发表于 2023-1-14 23:36
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2023-1-30 06:12 , Processed in 0.082658 second(s), 25 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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