QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 1573|回复: 0
打印 上一主题 下一主题

MATLAB 求和与对数之间的关系

[复制链接]
字体大小: 正常 放大

1186

主题

4

听众

2922

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2024-9-26 17:19 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
这段MATLAB代码用于计算一个极限,具体是涉及到求和与对数之间的关系。以下是对这段代码的详细解释:
, F& ?, N( G! R; E4 V  l9 r7 |; _. S7 f1 d- C( G, O/ e5 M5 i
### 1. 定义符号变量
0 d: B7 e: G. L```matlab( b5 S5 E8 N7 j) ^
syms m n;
5 v* A1 s: Q) w; t, C7 ````
9 C! V4 V, s2 F1 w- 使用 `syms m n` 定义了两个符号变量 `m` 和 `n`,这两个变量将用于后续的符号运算。! P( G/ k0 r, `4 G  X. K
, a8 n' \( x9 c2 p& a2 [6 J
### 2. 计算求和和对数的差
& w3 c2 x( Z5 P```matlab
  1. limit(symsum(1/m, m, 1, n) - log(n), n, inf)
复制代码
```
( L* G# }- X( l- x2 x. R+ L; C% A- `symsum(1/m, m, 1, n)`:  W' b) Z* ]! g, u. |( X6 u+ m5 S, b4 V
  - `symsum` 函数计算从 `m=1` 到 `m=n` 的级数和,这里具体是求 `1/m` 的和。9 B1 j1 [3 Q% m
  - 结果是哈默尼克级数,表示为 \( H_n = \sum_{m=1}^{n} \frac{1}{m} \)。" k# Y1 a+ l8 t, V, R& @

: U& f7 J( V8 g& e7 N7 S: f$ {- `log(n)`:
2 A2 T: j7 s: y  ^% c3 S  - 这是以自然对数为底的对数函数,表达 `n` 的对数。
$ A( d7 E; V3 j0 m2 q* b% i; |2 P9 u) h% r
- `limit(..., n, inf)`:/ b% Z" c: R) R% ?7 m  d8 y$ q8 G, x
  - `limit` 函数用于计算当 `n` 趋近于无穷大时,`(H_n - \log(n))` 的极限。. J7 H9 p6 H0 f+ E( I
  - 根据调和级数的性质,我们知道 \( H_n \) 的增长速率与 \( \log(n) \) 相关,且 \( H_n \) 与 \( \log(n) \) 的差收敛于一个常数。
: q: M) b! @$ \& p- P6 ^; e
2 p, L. O3 G* E) X### 3. 显示结果
  1. vpa(ans, 70)  % 显示 70 位有效数字
复制代码
- `vpa(ans, 70)`:2 `/ l3 r8 C# O5 B8 i# z6 e) M
  - `vpa` 表示“可变精度算术”,用于以高精度显示计算结果。
2 A# R# e' t' {6 |9 R+ T2 n  - `ans` 是 MATLAB 中的默认变量,它保存上一个计算的结果。
* n1 |. A# C5 o1 W' g  - 该函数将结果显示为70位有效数字。5 `" E- B8 u7 W3 K3 x

2 k# y& P3 o3 g  p4 g### 总结! t. z% Q4 h# Z, V4 q0 P$ g
这段代码首先计算出哈默尼克级数的和与自然对数之间的差,当 `n` 趋于无穷时的极限。然后,结果将以70位有效数字的形式输出。这个极限的值实际上是著名的常数——欧拉–马歇罗尼常数(Euler–Mascheroni constant),通常记作 \( \gamma \),即:
5 G6 m) z7 C2 L" j\[9 }0 P( ^0 s% u9 A. P6 K# T- A
\gamma = \lim_{n \to \infty} \left( H_n - \log(n) \right)
7 J0 e' Y/ A6 V) D8 f# z/ D\]9 R$ y9 C+ ?# \* @. U# x, J
此常数的值大约为 0.577215664901532。但是,通过 `vpa` 能够提供更多的有效位数,使结果更为精确。( n& f; ^$ s" c

% \" w/ c8 l  l5 T" q$ x9 V. m3 B) M$ i" T5 x8 |

) P. d9 a4 {; X$ r! j) A& l
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-4-10 16:44 , Processed in 0.352084 second(s), 51 queries .

回顶部