找回密码
 欢迎注册
楼主: liangbch

[提问] 黎曼函数平凡零点

[复制链接]
发表于 2013-8-20 16:49:04 | 显示全部楼层
本帖最后由 mathematica 于 2013-8-20 16:50 编辑

通过这个来定义负数的zeta函数的.
sin(pi*s/2)在s=-2,-4,-6,-8处明显等于零!
QQ截图20130820164751.jpg
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-20 16:53:05 | 显示全部楼层
这个论坛上的人的水平还没高到可以搞这种问题的时候,当然民科除外!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-20 17:06:43 | 显示全部楼层
黎曼函数是常用的函数之一,自然需要了解,谈不上研究。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-21 17:04:59 | 显示全部楼层
本帖最后由 云梦 于 2013-8-21 17:10 编辑

大于800时的非整数及整数的黎曼函数(正偶数、正奇数及负奇数-∞<0<∞)计算程序已经完成。

点评

我觉得你可以考虑把你的代码贴出来,mathematica有现成的计算函数吧  发表于 2013-8-21 19:56
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-21 17:26:11 | 显示全部楼层
本帖最后由 云梦 于 2013-8-21 20:12 编辑

小于800的非整数如何计算还没找到。精度1024位。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-21 20:17:45 | 显示全部楼层
本帖最后由 云梦 于 2013-8-21 20:39 编辑

计算黎曼函数常用到k^n,如果用对数计算可以这样:
x=n*Log[k]
x=Exp[x]
但Log[x]和Exp[x]在高精度时很浪费时间,尤其是当n较小时。
因此采用乘法代替乘方。
比如计算2^n可以这样:乘法次数<2log[n]/log[2]


Dim n, i, j As Long
Dim m(10) As Integer

n = Val(Left(Ax.st, Ax.zs + 1))
For i = 0 To 100
    If n Mod 2 = 1 Then
         n = (n - 1) / 2
         m(i) = 1
    Else
         n = n / 2
         m(i) = 0
    End If
    If n = 1 Then Exit For
Next
j = i
Ax.Bz = ""
Ax.st = "2"
Ax.zs = 0
Ex = Ax
Px = Ex

For i = 0 To j
   Ax = Px
   Cx = Px
   Call Mult_
   If m(j - i) = 1 Then
      Ax = Ex
      Call Mults_
      Px = Cx
   Else
      Px = Cx
   End If
Next

乘法:CX*AX ->CX   Mult_() 高精度乘法(指定精度,最大到1024),Mults_() 双精度乘法(12位)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-21 20:25:31 | 显示全部楼层
云梦 发表于 2013-8-21 17:04
大于800时的非整数及整数的黎曼函数(正偶数、正奇数及负奇数-∞


代码即使贴出来也不好看懂,只有我才懂得代码的意义。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-22 11:20:47 | 显示全部楼层
n=999/1000
z=960
N[(z/(n-1)+1/2)/z^n+Sum[BernoulliB[2 k] Gamma[n+2 k-1]/Gamma[n]/(2 k)! z^(1-n-2 k),{k,1,400}]+Sum[1/k^n,{k,1,z-1}],1024]
当n不是整数时Mathmatica的计算速度也很慢,可为什么计算Zeta[n]时速度会很快?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-22 15:18:40 | 显示全部楼层
除了速度问题,黎曼函数的计算已经全部完成。对所有实数均计算通过(-∞<n<+∞)。
下一步将进行全面优化。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-22 17:31:00 | 显示全部楼层
j计算难度更大的是不完全黎曼函数。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-3-29 14:08 , Processed in 0.046673 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表