1万以内相邻质数之差分布图——史前巨兽“心电图”

质数很神秘,质数的分布是否有规律呢?

(今天在网上看到一篇关于微分几何可视化的文章,感觉可以用来绘制更为精准细致的质数分布图,但本人水平有限,只能粗略的画一下简单的差值分布)

类比斐波那契数列,利用Python计算了1000以内相邻两个质数之差分布图,有点像心跳图,在不断地振动,有人见过这种“波”吗?

做折线图的时候,1000个数字为一组,一共10张图,只选取2张展示,

其他列在文章末尾

质数,或者说素数,就像是一个史前巨兽;现在,让我们看一下他的“心电图”

Python代码如下:

import matplotlib as mpl
from matplotlib import pyplot as plt

### 1000以内相邻质数之差分布-折线图
### 看起来好像是某种”波“,或者”振动“

num1=[];
i=2
for i in range(2,1000): ##修改点1,这里修改范围,1000个数字一组
   j=2
   for j in range(2,i):
      if(i%j==0):
         break
   else:
      num1.append(i)
print(num1)

num2 = [num1[i] for i in range(1,len(num1)) ]
print(num2)
print(len(num2))

num3 = []
for i in range(len(num2)):
    num3.append(num2[i]-num1[i])
print(num3)

#x = range(2,26,2)
#y = [15,13,14.5,17,20,25,26,26,27,22,18,15]

x = range(len(num2))
y = num3

#指定字体为SimHei,用于显示中文,如果Ariel,中文会乱码
mpl.rcParams["font.sans-serif"]=["SimHei"]
mpl.rcParams["axes.unicode_minus"]=False

plt.title(u"2~1000以内相邻质数之差分布图:")  ##修改点2 ,设置图片标题

#绘图
plt.plot(x,y)
#展示图形
plt.show()

毕竟不是数学专业,就到这里为止吧,以免走火入魔。

参考了几篇文章,如下:

https://blog.csdn.net/hakesashou/article/details/134365984

https://zhuanlan.zhihu.com/p/140502200

https://juejin.cn/post/6893098337903575053

  • 11
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值