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


>>> 数学分析,奇异积分,几何,代数,微分方程,群与环,数论
Elinkage数学论坛基础数学 [返回] → 浏览:[求助]请求由点坐标求出函数表达式 标记论坛所有内容为已读 

 目前论坛总在线 4 人,本主题共有 1 人浏览。其中注册用户 0 人,访客 1 人。  [关闭详细列表]
发表一个新主题 回复贴子 开启一个新投票 ◆此帖被阅读 103 次◆  浏览上一篇主题  刷新本主题  树形显示贴子 浏览下一篇主题
 * 贴子主题: [求助]请求由点坐标求出函数表达式 不分页显示此帖  保存该页为文件  本贴有问题,发送短消息报告给版主  加入个人收藏&关注本贴  显示可打印的版本  把本贴打包邮递  把本贴加入收藏夹  发送本页面给朋友   
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [楼 主]
  各位老师和朋友们,请帮忙看看如何做?不必完全等于实际,近似值就行,略小于实际,不能大于实际。谢谢!如下是问题:[1s9A`
求如下这些点的函数9?4
(4,9),(5,28),(6,80),(7,221).(
我弄的公式:b = 1.74136140824861 ^ (1.44 * a) - 1.74136140824861 ^ (a + 0.92)。*
a依次输入4,5,6,7的值为:9.0903320129064,27.5776386262734,74.135055990122,187.135893810461[
朋友的公式r]Ox
y=0.1311exp(1.0653x)/K!v
您的公式的验证,依次输入4,5,6,7为9.294335344023,26.9694599213373,78.2575344579496,227.08062072062。最后一个值大于实际了。Gff
7往后都大于实际了。1[.
7往后的值分别是221,608,1694,4775,13606.6%


发贴时间2020/05/16 10:05am 此 IP 为代理服务器IP: 已设置保密[本文共624字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 2 楼]
  陆元鸿教授的精彩解答:(好象也是拉格朗日插值法吧)SGkTY'
y = ( 56 x^3 - 741 x^2 + 3367 x - 5142 ) / 6QA
©Elinkage数学论坛 -- Elinkage极酷超级论坛  ?P}
x = 4 ,y = 9 。MH
©Elinkage数学论坛 -- Elinkage极酷超级论坛  :
x = 5 ,y = 28 。8<rR[m
©Elinkage数学论坛 -- Elinkage极酷超级论坛  S
x = 6 ,y = 80 。W=>
©Elinkage数学论坛 -- Elinkage极酷超级论坛  i0gM#8
x = 7 ,y = 221 。E,
©Elinkage数学论坛 -- Elinkage极酷超级论坛  W
x = 8 ,y = 507 。_P1>mw
©Elinkage数学论坛 -- Elinkage极酷超级论坛  U^#
x = 9 ,y = 994 。".jB77
©Elinkage数学论坛 -- Elinkage极酷超级论坛  ==
x = 10 ,y = 1738 。[
©Elinkage数学论坛 -- Elinkage极酷超级论坛  G
x = 11 ,y = 2795 。y{.
©Elinkage数学论坛 -- Elinkage极酷超级论坛  ?Na
x = 12 ,y = 4221 。E_TC:


发贴时间2020/05/16 03:25pm 此 IP 为代理服务器IP: 已设置保密[本文共313字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 3 楼]
  其实前面那些整数点并非都是整数,而是取消了小数部分,依次输入4,5,6,7对应的带小数的值为9.1950831715873,28.9888206573201,80.4899477790795,221.074445101172。,i9AP
后面的依次是608.992317564808,1694.99578749547,4775.7128358109,13606.088463196。sIq"W*
©Elinkage数学论坛 -- Elinkage极酷超级论坛  t
应该能弄出更准确的公式,欢迎各位老师和朋友继续给予指导!BW1


发贴时间2020/05/16 03:28pm 此 IP 为代理服务器IP: 已设置保密[本文共328字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 153297 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2236
精华: 0
资料:  
在线: 46天5时14分15秒
注册: 2010/12/07 01:27pm
造访: 2020/05/24 11:12am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 4 楼]
  更精確的通項只能基於更精確,更多的實際數據. qP,m
因為除此之外都是缺乏根據的拼湊.Ok2M


发贴时间2020/05/17 05:56am IP: 已设置保密[本文共102字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 5 楼]
  谢谢!很对,谢谢老师!%dD'X


发贴时间2020/05/17 09:29am 此 IP 为代理服务器IP: 已设置保密[本文共46字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 6 楼]
  我已经把陆元鸿教授的公式用于计算素数个数的程序,结果与实际非常接近,如下为结果:ew{ST,
10内的素数个数下限为 4,实际为4个。100内的素数个数下限为 22,实际为25个。:cM
1000内的素数个数下限为 163,实际为168个。b*{
10000内的素数个数下限为 1227,实际为1229个。>wJf!
100000内的素数个数下限为 9564,实际为9592个。Z6/"~?
1000000内的素数个数下限为 78463,实际为78498个。R9
10000000内的素数个数下限为 664561,实际为664579个。J0E}K
(只用于10^7以内的数。)jW#~0!
大于10^7采用我们的近似公式:(下面的b1和b4,而b3为陆元鸿教授的)U>
b = Exp(1.033 * a)l/
b1 = 0.1578 * bYu
b3 = (56 * a ^ 3 - 741 * a ^ 2 + 3367 * a - 5142) / 6* X`@j
b4 = (1.74136140824861 ^ (1.8 * a)) / 4.9+B
b5 = (b4 + b1) / 25\A
©Elinkage数学论坛 -- Elinkage极酷超级论坛  p
a依次输入7,8,9,10,11时b1,b4的平均值b5分别为219.664329196356,606.570528142824,1675.45599266651,;
4629.29243022057,`h
12794.5765370479,(调整一下还能更接近实际。)Y(<3|H
如下是利用这个公式得出的素数个数:L?x
100000000内的素数个数下限为 5760140,实际为5761455.rh
1000000000内的素数个数下限为 50817717,实际为50847534._A/^Hm
1600000000内的素数个数下限为 79394377,实际为79451833个。10000000000内的素数个数下限为 454416615,实际为455052511.fzM<L
100000000000内的素数个数下限为 4107923014,实际为4118054813.u
160000000内的素数个数下限为 8972209,实际为8974458.%
1900000000内的素数个数下限为 93475776。&y1:z


发贴时间2020/05/17 09:38am 此 IP 为代理服务器IP: 已设置保密[本文共1298字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 153297 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2236
精华: 0
资料:  
在线: 46天5时14分15秒
注册: 2010/12/07 01:27pm
造访: 2020/05/24 11:12am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 7 楼]
  既然考虑的是素数,当然就要用到素数定理,^`
问题归结为如何得到更好的π(10x)π(10x)的渐近公式.],*Rk9
按此在新窗口浏览图片lD


发贴时间2020/05/19 03:18am IP: 已设置保密[本文共149字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 8 楼]
  由于当输入1时原公式的值是0,我们不从1开始是从4开始的,由于原公式是改进的,加了提前量4,若去掉这个提前量是否会好弄?请老师用下面的数据试试?*
可以试试这几个数,重新做出个函数公式吧。就是(4,13),(5,32),(6,84),(7,225),(8,612)。)Otg(
还有下面的数(9,1698),(10,4779),(11,13610).2_%$
电脑中有这样的自动给出函数结果的程序吗?o~g
©Elinkage数学论坛 -- Elinkage极酷超级论坛  _3 LiL
感谢老师的数据,我需要这个,再多算几项!3


发贴时间2020/05/19 11:19am 此 IP 为代理服务器IP: 已设置保密[本文共411字节]  
 ysr 




等级: 新手上路
信息: 该用户目前不在线
威望: 0 积分: 0
现金: 2898 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 50
精华: 0
资料:  
在线: 13 时 55 分 00 秒
注册: 2013/04/04 00:24am
造访: 2020/05/24 03:34am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 9 楼]
  对,所谓的素数定理就是这两个公式,好象是,我用的是第一个,就是下限公式,而且是改进了一下。$75g
不改进的话数据是上面发的值,规律是一样的。`9Gm
谢谢老师!*


发贴时间2020/05/19 11:31am 此 IP 为代理服务器IP: 已设置保密[本文共174字节]  
 elim 
 头衔: 论坛版主

 

等级: 新手上路
信息: 该用户目前不在线 此人为版主
威望: 0 积分: 0
现金: 153297 雷傲元
存款: 没开户
贷款: 没贷款
来自: 保密 blank
发帖: 2236
精华: 0
资料:  
在线: 46天5时14分15秒
注册: 2010/12/07 01:27pm
造访: 2020/05/24 11:12am
消息 查看 搜索 好友 引用 回复贴子回复 只看我 [第 10 楼]
  现在按照我的帖子(一般公式)写代码. 用计算机完成所需数据的近似计算.首先贴出代码:s#E

代码:
from math import *zVP
©Elinkage数学论坛 -- Elinkage极酷超级论坛  wZu80"
def fct(n):Col
   p = 1),5n
   for k in range(n):ftDA
       p *= (k+1)q2%>.
   return pe6U=R
©Elinkage数学论坛 -- Elinkage极酷超级论坛  +
©Elinkage数学论坛 -- Elinkage极酷超级论坛  DJ:1
def lmd(x,k,m,a,h):{zG
   v=(-1)**(m-k-1)*((0.0+a+k*h)/(x*h))**(m-1)=_
   v/= (fct(k)*fct(m-1-k))mNqA4Q
   for j in range(m):YJ1l
       if k != j:!
           v *= (x-a-h*j) /be/@
   return vwE&q
©Elinkage数学论坛 -- Elinkage极酷超级论坛  ^j
©Elinkage数学论坛 -- Elinkage极酷超级论坛  Q#"K
def px(x):8L
   w= 10**x:we
   return (1.0+45./w)*w/log(w)HF=I{
©Elinkage数学论坛 -- Elinkage极酷超级论坛  =$`/Cs
©Elinkage数学论坛 -- Elinkage极酷超级论坛  }
def g(x,L,funct,a=1,h=1):z0W
   m=len(L)~=s*^
   M=[L[k]/funct(a+h*k) for k in range(m)]9JtwG
   w = 02.ulzV
   for k in range(m):uF
       w += M[k]*lmd(x,k,m,a,h)$!O7
   return w*funct(x)

下面是对这段代码的简单注释:U3Jz1
$(1)\quad\text{ fct}(n)\,$是阶乘$\,n!\,$的python代码.xe^t
$(2)\quad\text{ lmd}(x,\,k,\,m,\,a,\,h)\,$是一般公式中的$\,\lambda_k(x,x_0,h,m).$XFeyi2
$(3)\quad\text{ px}(x)\,$是$\,(1+\frac{a}{w})\frac{w}{\log w},\;w=10^x\,$的python函数定义.[hL9
$(4)\quad\text{ g}(x,L,\text{funct},a=1,h=1)\,$是一般公式本身.其中$\text{L}$0{sP<(
$\quad$是插值数据(python List,数列$m$是数列长),funct是.?
$\quad$被逼近函数.可以是上面的$\,\text{px}(x)=(1+{\scriptsize\dfrac{45}{10^{\large x}}}){\scriptsize\dfrac{10^{\large x}}{\log(10^{\large x})}}$IuMt
$\quad\approx\pi(10^x).\;$其中$1+\frac{45}{10^{\,x}}\,$是矫正因子,大有改进余地.910
$\quad$最后$,\,x_k=a+kh\,$解释了参数$\,a,\,h:$MLpUSf
$\qquad L=[\varphi(a+0h),\ldots,\varphi(a+(m-1)h)].$q8|Rv
按此在新窗口浏览图片gMrE


发贴时间2020/05/19 02:17pm IP: 已设置保密[本文共1812字节]  

 2 9 7 [ 1 2 ] 8 :

快速回复主题: [求助]请求由点坐标求出函数表达式
您目前的身份是: 客人 ,要使用其他用户身份,请输入用户名和密码。未注册客人请输入网名,密码留空。
输入用户名和密码: 用户名: 没有注册? 密码: 忘记密码?
上传附件或图片 (最大容量 10000KB)
目前附件:(如不需要某个附件,只需删除内容中的相应 [UploadFile ...] 标签即可) [删除]
选项

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

使用字体转换?

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


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

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