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


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

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

 

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



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

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 124046 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 1882
精华: 0
资料:  
在线: 874 时 40 分 59 秒
注册: 2010/12/07 06:27am
造访: 2018/09/22 09:26am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 2 楼]
  题:计算极限$\;\;\displaystyle{\underset{\,}{\;\lim_{n\to\infty}}\frac{1^n+2^n+\cdots+n^n}{n^n}}$iET
解:令$\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)$3O
$\qquad$由Taylor定理,("7t/v
$\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)}$P
$\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)$b|uF+P
$\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)$dF
$\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$j
$\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$He/-y



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

 

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



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

 

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

(m|Z
$\qquad\qquad\qquad\qquad\qquad$GP/PARI CALCULATOR Version 2.9.2 (released)!-u
............#
$\qquad\qquad\qquad\qquad\qquad\quad$Copyright (C) 2000-2017 The PARI Group<~2
©Elinkage数学论坛 -- Elinkage极酷超级论坛  yz<
PARI/GP is free software, covered by the GNU General Public License, and comesEWg}
WITHOUT ANY WARRANTY WHATSOEVER.&p7$Q
©Elinkage数学论坛 -- Elinkage极酷超级论坛  p]9
(10:05) gp > \p 724
  realprecision = 77 significant digits (72 digits displayed)u!wi
(10:05) gp > sm(n,m) =4
{ my(s=1);1Y#R
 for(j=1,m,`$G_
   s += exp(n*log((n-j)/n))|>6$5N
 );'
 return(s);j~K[[
}©Elinkage数学论坛 -- Elinkage极酷超级论坛  \
(10:06) gp > v = 1/(1-exp(-1))OL
%2 = 1.581976706869326424385002005109011558546869301075396136266787059648043825M1Tnw
(10:07) gp > sm(1000000,169)d)
%3 = 1.58197571072294110657980881326425291863135428213215632546731155979676445B&*l
(10:07) gp > sm(10^80,169)2K3IS
%4 = 1.58197670686932642438500200510901155854686930107539613626678705964804382@KOt
(10:07) gp > sm(10^70,169)ftxg
%5 = 1.58197670686932642438500200510901155854686930107539613626678705964804372g_k,W
(10:07) gp > sm(10^60,169)asu6
%6 = 1.58197670686932642438500200510901155854686930107539613626678606350066025&@ ,8/
(10:07) gp > sm(10^50,169)zwCc{
%7 = 1.581976706869326424385002005109011558546869301075386174792951434711079191r@8?U
(10:07) gp > sm(10^40,169)'c
%8 = 1.58197670686932642438500200510901155854676968633703988689714075881739371)Dp
(10:07) gp > sm(10^30,169)]
%9 = 1.58197670686932642438500200510801541116330680737893312796028699919981775yT
(10:07) gp > sm(10^20,169)x^<\
%10= 1.58197670686932642437504053127338662158233904257934893714275487447672291S
(10:07) gp > sm(10^10,169)IZElm
%11= 1.581976706769711686038735092358048377802615047421191041130109018644829512
(10:07) gp > sm(10^72,169)e&#&z
%12= 1.58197670686932642438500200510901155854686930107539613626678705964804382

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



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