QQ登录

只需要一步,快速开始

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

MATLAB 求和与对数之间的关系

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

1192

主题

4

听众

2946

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2024-9-26 17:19 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
这段MATLAB代码用于计算一个极限,具体是涉及到求和与对数之间的关系。以下是对这段代码的详细解释:
5 c6 Z/ w9 q' p8 n" c# C! f- }5 s& z& ]% w2 n. F
### 1. 定义符号变量* e8 H. e3 o1 P- V
```matlab+ G9 N" {( H; o# y8 P
syms m n;
: i3 h* S( m2 v  D4 @1 T0 D```2 r! V% O; c: W; z% E! D
- 使用 `syms m n` 定义了两个符号变量 `m` 和 `n`,这两个变量将用于后续的符号运算。% M* N1 Y* x1 x: o6 x
7 t. ?% N4 f! V/ x7 p) d' C& H, e
### 2. 计算求和和对数的差3 r- z2 b% W. ?3 X8 n# M9 j
```matlab
  1. limit(symsum(1/m, m, 1, n) - log(n), n, inf)
复制代码
```
" C4 s8 _! N/ h( [! E% [5 V- `symsum(1/m, m, 1, n)`:& }; Q8 A  b6 j/ g  o
  - `symsum` 函数计算从 `m=1` 到 `m=n` 的级数和,这里具体是求 `1/m` 的和。
; m, g7 _% m( }  a+ F1 C5 F* Z% `  - 结果是哈默尼克级数,表示为 \( H_n = \sum_{m=1}^{n} \frac{1}{m} \)。
6 q  f9 l: l4 w4 V! l7 ?& m
/ ~$ I4 g  M3 D" s. v- `log(n)`:
. Q9 H2 p" J/ a$ O* k8 g/ w  - 这是以自然对数为底的对数函数,表达 `n` 的对数。
( S7 T. y4 b) H5 |) @4 ]4 w2 J
/ s: `+ m' J. i: U7 a- `limit(..., n, inf)`:' h  c9 I+ a: q2 X& k3 M: I
  - `limit` 函数用于计算当 `n` 趋近于无穷大时,`(H_n - \log(n))` 的极限。
! I. `* j% Z9 W, ?( t: T: @2 X0 T0 S$ C  - 根据调和级数的性质,我们知道 \( H_n \) 的增长速率与 \( \log(n) \) 相关,且 \( H_n \) 与 \( \log(n) \) 的差收敛于一个常数。
' }4 d+ ?+ f( B& M" x" A* k& t2 G# O3 `( F3 }; W! J
### 3. 显示结果
  1. vpa(ans, 70)  % 显示 70 位有效数字
复制代码
- `vpa(ans, 70)`:1 [7 b4 h5 P# F. T7 F+ K' f( q
  - `vpa` 表示“可变精度算术”,用于以高精度显示计算结果。5 B7 `3 G$ L1 ?& @) m9 `* m2 i& M
  - `ans` 是 MATLAB 中的默认变量,它保存上一个计算的结果。9 z9 L! Y. \1 |  d
  - 该函数将结果显示为70位有效数字。' M% m" I% C6 B, H) ]: X( q

3 Z. t+ x, Q4 t6 i' L$ {3 }* E### 总结. R' E: `+ ~6 D  J- Z
这段代码首先计算出哈默尼克级数的和与自然对数之间的差,当 `n` 趋于无穷时的极限。然后,结果将以70位有效数字的形式输出。这个极限的值实际上是著名的常数——欧拉–马歇罗尼常数(Euler–Mascheroni constant),通常记作 \( \gamma \),即:
4 _4 g6 N6 W* P4 i9 h\[  R/ l. u" i; f1 v
\gamma = \lim_{n \to \infty} \left( H_n - \log(n) \right)/ m% A, F& Q% e: i! F& F6 _' Q; T
\]
7 \2 ]5 W" C, b2 P' C6 V" B此常数的值大约为 0.577215664901532。但是,通过 `vpa` 能够提供更多的有效位数,使结果更为精确。
  L9 g! ^/ p4 g7 y1 |% Y" D: m! T9 ?: F% e

  K4 M  t. N- a: t. x* O$ o" D- @/ M$ S
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-8-25 16:39 , Processed in 0.319047 second(s), 51 queries .

回顶部