>> 欢迎您,客人登录 按这里注册 忘记密码 在线 搜索 论坛风格  帮助  插件   


>>> 数学分析,奇异积分,几何,代数,微分方程,群与环,数论
Elinkage数学论坛基础数学 [返回] → 浏览:计算极限$\displaystyle{\;\lim_{n\to\infty}\frac{1^n+2^n+\cdots+n^n}{n^n}}$ 标记论坛所有内容为已读 

 目前论坛总在线 9 人,本主题共有 1 人浏览。其中注册用户 0 人,访客 1 人。  [关闭详细列表]
发表一个新主题 回复贴子 开启一个新投票 ◆此帖被阅读 56 次◆  浏览上一篇主题  刷新本主题  树形显示贴子 浏览下一篇主题
 * 贴子主题: 计算极限$\displaystyle{\;\lim_{n\to\infty}\frac{1^n+2^n+\cdots+n^n}{n^n}}$ 不分页显示此帖  保存该页为文件  本贴有问题,发送短消息报告给版主  加入个人收藏&关注本贴  显示可打印的版本  把本贴打包邮递  把本贴加入收藏夹  发送本页面给朋友   
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 133823 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2000
精华: 0
资料:  
在线: 944 时 31 分 14 秒
注册: 2010/12/07 06:27am
造访: 2019/03/17 08:50am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [楼 主]
  题:计算极限$\;\;\displaystyle{\;\lim_{n\to\infty}\frac{1^n+2^n+\cdots+n^n}{n^n}}$aU712


发贴时间2017/09/11 10:57am IP: 已设置保密[本文共106字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 133823 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2000
精华: 0
资料:  
在线: 944 时 31 分 14 秒
注册: 2010/12/07 06:27am
造访: 2019/03/17 08:50am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 2 楼]
  题:计算极限$\;\;\displaystyle{\underset{\,}{\;\lim_{n\to\infty}}\frac{1^n+2^n+\cdots+n^n}{n^n}}$#'C`BV
解:令$\displaystyle{\;R_n := \small\frac{1^n+2^n+\cdots+n^n}{n^n}=\sum_{k=0}^{n-1}\frac{(n-k)^n}{n^n} = \sum_{k=0}^{n}\big(1-\frac{k}{n}\big)^n}\quad(1)$1QP
$\qquad$由Taylor定理,G
$\underset{\,}{\qquad} f(x):=e^x-(1+x) = {\large\frac{x^2}{2}}f\,''(x\theta)={\large\frac{x^2}{2}}e^{\theta x}>0\;{\small(x\ne 0<\theta<1)}$*}`
$\underset{\,}{\therefore}\quad 1-\frac{k}{n}< e^{-\frac{k}{n}},\quad \big(1-\frac{k}{n}\big)^n < e^{-k}\;\;\small(k\le n)$v
$\qquad$故对$\,m< n,\;$有$\,\small\displaystyle{\;\sum_{k=0}^m\big(1-\frac{k}{n}\big)^n < R_n < \sum_{k=0}^{\infty}e^{-k}=\underset{\,}{\frac{e}{e-1}}}\qquad\quad(2)$0:dw
$\qquad$令$\,n\to\infty\,$得$\quad\small\displaystyle{\sum_{k=0}^m e^{-k}\le\lim_{\overline{n\to\infty}}R_n\le \overline{\lim_{n\to\infty}}R_n\le\underset{\,}{\frac{e}{e-1}}}.\;$再令$\;m\to\infty$?
$\qquad$即得$\displaystyle{\;\;\small\lim_{n\to\infty}\frac{1^n+2^n+\cdots+n^n}{n^n}=\frac{e}{e-1}}=1.5819767068693\ldots\;\;\square$2)%o9y


发贴时间2017/09/11 00:27pm IP: 已设置保密[本文共1098字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 133823 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2000
精华: 0
资料:  
在线: 944 时 31 分 14 秒
注册: 2010/12/07 06:27am
造访: 2019/03/17 08:50am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 3 楼]
  按此在新窗口浏览图片Y-#C


发贴时间2017/09/13 06:30pm IP: 已设置保密[本文共60字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 133823 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2000
精华: 0
资料:  
在线: 944 时 31 分 14 秒
注册: 2010/12/07 06:27am
造访: 2019/03/17 08:50am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 4 楼]
  数值计算 $\displaystyle{R_n=\small\frac{1^n+2^n+\cdots+n^n}{n^n}}$ 颇具挑战性. ($n^n$ 增长极快!)=
用Mathematica 直接计算,当$\,n\underset{\,}{>}10^6\,$时开始很辛苦,很快就变得几乎不可能. ?"bA
为此作以下估算:对正整数$\underset{\,}{\,}d,\,m,\,n\,(n>m)\,$为使6
$\displaystyle{\underset{\,}{\small\sum_{k=1}^{n-m}}{\small\frac{k^n}{n^n}=\sum_{k=m}^{n-1}\big(1-\frac{k}{n}\big)^n< \sum_{k=m}^{\infty}e^{-k}=\frac{e}{e-1}e^{-m}}< 10^{-d}},\;$只需$\;\ln\big({\large\frac{e}{e-1}}10^d\big)< m.$Ob.u+A
举例说来,对$\,d=73,\,$有$\,m> 168.547\ldots$ 于是有w91U
$\displaystyle{0< {\small\frac{1^n+\cdots+n^n}{n^n}-\sum_{k=0}^{169}\big(\frac{n-k}{n}\big)^n }< 10^{-73}\;(n > 169)}.\;$即只需计算$\,169\,$项和,2jnr
就可得到误差小于$\,10^{-73}\,$的所求比值. 以下用 GP/Pari来践行我们的分析:

#dL]o
$\qquad\qquad\qquad\qquad\qquad$GP/PARI CALCULATOR Version 2.9.2 (released)T&\v
............%+sx
$\qquad\qquad\qquad\qquad\qquad\quad$Copyright (C) 2000-2017 The PARI GroupzeY<#Z
©Elinkage数学论坛 -- Elinkage极酷超级论坛  g
PARI/GP is free software, covered by the GNU General Public License, and comes=
WITHOUT ANY WARRANTY WHATSOEVER.yH
©Elinkage数学论坛 -- Elinkage极酷超级论坛  7iE"
(10:05) gp > \p 72CA^(}
  realprecision = 77 significant digits (72 digits displayed)|\L%
(10:05) gp > sm(n,m) =}AL{)v
{ my(s=1);7OH`11
 for(j=1,m,|
   s += exp(n*log((n-j)/n)){
 );#RIh
 return(s);Y0rp
}©Elinkage数学论坛 -- Elinkage极酷超级论坛  >.G
(10:06) gp > v = 1/(1-exp(-1))8XGSQ
%2 = 1.58197670686932642438500200510901155854686930107539613626678705964804382ngUD`
(10:07) gp > sm(1000000,169)T%Mf
%3 = 1.58197571072294110657980881326425291863135428213215632546731155979676445wv^?!.
(10:07) gp > sm(10^80,169)>F[
%4 = 1.58197670686932642438500200510901155854686930107539613626678705964804382@g>2LQ
(10:07) gp > sm(10^70,169)iT9?
%5 = 1.58197670686932642438500200510901155854686930107539613626678705964804372tZL s}
(10:07) gp > sm(10^60,169)|gO
%6 = 1.58197670686932642438500200510901155854686930107539613626678606350066025<{
(10:07) gp > sm(10^50,169){@G^]
%7 = 1.58197670686932642438500200510901155854686930107538617479295143471107919x]zU|
(10:07) gp > sm(10^40,169)Z
%8 = 1.58197670686932642438500200510901155854676968633703988689714075881739371: &U
(10:07) gp > sm(10^30,169)O{%
%9 = 1.58197670686932642438500200510801541116330680737893312796028699919981775Y;
(10:07) gp > sm(10^20,169)-kDx+
%10= 1.58197670686932642437504053127338662158233904257934893714275487447672291V>
(10:07) gp > sm(10^10,169)4oYKd'
%11= 1.581976706769711686038735092358048377802615047421191041130109018644829518aU
(10:07) gp > sm(10^72,169)pj
%12= 1.58197670686932642438500200510901155854686930107539613626678705964804382

_Bqg~_
©Elinkage数学论坛 -- Elinkage极酷超级论坛   u@
用蓝色标示计算结果与极值的吻合部分. 我们看到$\,n\,$为一百万只能得到小数点后pV^
五位的准确度. 另一方面, 虽然需要$\,n=10^{72}\,$才能得到所欲的精度,却只要计算LZ
$\,169\,$项. 并且完全避免了$(10^{72})^{10^{72}}$之类疯狂的取幂运算. 以至全部计算时间几乎(
可以不计.,+Y


发贴时间2017/09/16 03:22am IP: 已设置保密[本文共3298字节]  

 该主题只有一页

快速回复主题: 计算极限$\displaystyle{\;\lim_{n\to\infty}\frac{1^n+2^n+\cdots+n^n}{n^n}}$
您目前的身份是: 客人 ,要使用其他用户身份,请输入用户名和密码。未注册客人请输入网名,密码留空。
输入用户名和密码: 用户名: 没有注册? 密码: 忘记密码?
上传附件或图片 (最大容量 10000KB)
目前附件:(如不需要某个附件,只需删除内容中的相应 [UploadFile ...] 标签即可) [删除]
选项

使用 LeoBBS 标签?
显示您的签名?
有回复时使用邮件通知您?

使用字体转换?

    快速引用第 楼层的回复
 顶端 加到"个人收藏夹" 主题管理总固顶 取消总固顶 区固顶 取消区固顶 固顶 取消固顶 提升 沉底
加重 取消加重 精华 取消精华 锁定 解锁 删除 删除回复 移动


© 中文版权所有: 雷傲科技
程序版权所有: 雷傲超级论坛  版本: LeoBBS X Build051231
 

本论坛言论纯属发表者个人意见,与 Elinkage数学论坛 立场无关