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


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

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

 

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



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

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 128509 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 1942
精华: 0
资料:  
在线: 904 时 50 分 59 秒
注册: 2010/12/07 06:27am
造访: 2018/12/15 02:39am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 2 楼]
  题:计算极限$\;\;\displaystyle{\underset{\,}{\;\lim_{n\to\infty}}\frac{1^n+2^n+\cdots+n^n}{n^n}}$=^NDm
解:令$\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)$fjkC$
$\qquad$由Taylor定理,;1)9
$\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)}$kUs
$\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)$}AgG
$\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)$B
$\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$J]nK



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

 

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



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

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 128509 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 1942
精华: 0
资料:  
在线: 904 时 50 分 59 秒
注册: 2010/12/07 06:27am
造访: 2018/12/15 02:39am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 4 楼]
  数值计算 $\displaystyle{R_n=\small\frac{1^n+2^n+\cdots+n^n}{n^n}}$ 颇具挑战性. ($n^n$ 增长极快!)Hi)
用Mathematica 直接计算,当$\,n\underset{\,}{>}10^6\,$时开始很辛苦,很快就变得几乎不可能. xBUy
为此作以下估算:对正整数$\underset{\,}{\,}d,\,m,\,n\,(n>m)\,$为使+,hf+
$\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.$Ur\t2
举例说来,对$\,d=73,\,$有$\,m> 168.547\ldots$ 于是有E1/
$\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\,$项和,3
就可得到误差小于$\,10^{-73}\,$的所求比值. 以下用 GP/Pari来践行我们的分析:

Bk%
$\qquad\qquad\qquad\qquad\qquad$GP/PARI CALCULATOR Version 2.9.2 (released)dL
............\CWC
$\qquad\qquad\qquad\qquad\qquad\quad$Copyright (C) 2000-2017 The PARI Group)y
©Elinkage数学论坛 -- Elinkage极酷超级论坛  5e
PARI/GP is free software, covered by the GNU General Public License, and comesjSz}]-
WITHOUT ANY WARRANTY WHATSOEVER.0j"
©Elinkage数学论坛 -- Elinkage极酷超级论坛  dE}yi
(10:05) gp > \p 720!@o:S
  realprecision = 77 significant digits (72 digits displayed)/z
(10:05) gp > sm(n,m) =5@H]
{ my(s=1);`{v^4
 for(j=1,m,:
   s += exp(n*log((n-j)/n))e7{.\
 );e(ye
 return(s);Jj$7
}©Elinkage数学论坛 -- Elinkage极酷超级论坛  ~;2)
(10:06) gp > v = 1/(1-exp(-1))]y&EI7
%2 = 1.58197670686932642438500200510901155854686930107539613626678705964804382&ud[D
(10:07) gp > sm(1000000,169)dgr[
%3 = 1.58197571072294110657980881326425291863135428213215632546731155979676445mx>
(10:07) gp > sm(10^80,169)&e*5
%4 = 1.58197670686932642438500200510901155854686930107539613626678705964804382u<~{9
(10:07) gp > sm(10^70,169)Y'<'ku
%5 = 1.58197670686932642438500200510901155854686930107539613626678705964804372g6
(10:07) gp > sm(10^60,169)b]
%6 = 1.581976706869326424385002005109011558546869301075396136266786063500660258yDKev
(10:07) gp > sm(10^50,169)I=kd
%7 = 1.58197670686932642438500200510901155854686930107538617479295143471107919e
(10:07) gp > sm(10^40,169)x+U_/
%8 = 1.58197670686932642438500200510901155854676968633703988689714075881739371TbGt
(10:07) gp > sm(10^30,169)(p6E>k
%9 = 1.58197670686932642438500200510801541116330680737893312796028699919981775;<:wq
(10:07) gp > sm(10^20,169)e
%10= 1.58197670686932642437504053127338662158233904257934893714275487447672291JWSW
(10:07) gp > sm(10^10,169)"
%11= 1.58197670676971168603873509235804837780261504742119104113010901864482951k
(10:07) gp > sm(10^72,169)+^Q
%12= 1.58197670686932642438500200510901155854686930107539613626678705964804382

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



发贴时间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数学论坛 立场无关