文档章节

使用MicroPython计算任意位数圆周率

shaoziyang
 shaoziyang
发布于 2016/11/04 10:38
字数 446
阅读 1.7K
收藏 8

精选30+云产品,助力企业轻松上云!>>>

计算任意精度的圆周率是个有趣的主题,得益于python的强大计算能力,我们在MicroPython中也可以轻松的计算pi的数值。

先输入下面的代码:

"""
文件:pi.py
说明:用MicroPython计算任意精度圆周率计算
作者:未知
版本:
时间:
修改:邵子扬
      2016.5
      v1.1
http://bbs.micro-python.com/forum.php
"""
import time

def pi(places=10):
  # 3 + 3*(1/24) + 3*(1/24)*(9/80) + 3*(1/24)*(9/80)*(25/168)
  # The numerators 1, 9, 25, ... are given by (2x + 1) ^ 2
  # The denominators 24, 80, 168 are given by (16x^2 -24x + 8)
  extra = 8
  one = 10 ** (places+extra)
  t, c, n, na, d, da = 3*one, 3*one, 1, 0, 0, 24

  while t > 1: 
    n, na, d, da = n+na, na+8, d+da, da+32
    t = t * n // d
    c += t
  return c // (10 ** extra)

def pi_t(n=10):
    t1=time.ticks_us()
    t=pi(n)
    t2=time.ticks_us()
    print('elapsed: ', time.ticks_diff(t1,t2)/1000000, 's')
    return t

def pi2(n=10):
    r=6*(10**n)*1000
    p=0
    k=0
    c=r//2
    d=c//(2*k+1)
    while d>0:
        p=p+d
        k=k+1
        k2=2*k
        c=c*(k2-1)//(4*k2)
        d=c//(k2+1)
    return p//1000

def pi2_t(n=10):
    t1=time.ticks_us()
    t=pi2(n)
    t2=time.ticks_us()
    print('elapsed: ', time.ticks_diff(t1,t2)/1000000, 's')
    return t

 

我们将它复制到STM32F7DISC中(STM32F7DISC已经下载了MicroPython固件),然后测试一下计算速度。(在其它MicroPython开发板上也可以进行这个测试)

MicroPython v1.8 on 2016-05-15; F7DISC with STM32F746
Type "help()" for more information.
>>> import pi
>>> npi=pi.pi_t(1000)
elapsed:  0.221486 s
>>> npi=pi.pi_t(2000)
elapsed:  0.793141 s
>>> npi=pi.pi_t(5000)
elapsed:  4.981964 s
>>> npi=pi.pi_t(10000)
elapsed:  21.02012 s
>>>


从运行结果可以看出,虽然不是很快,但是考虑到STM32的资源和性能,结果已经出乎预料了,毕竟计算函数还不到10行代码,没有做深度优化。最后打印出1000位的圆周率,大家可以和标准结果比较看看。


>>> pi.pi(1000)
31415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337867831652712019091456485669234603486104543266482133936072602491412737245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094330572703657595919530921861173819326117931051185480744623799627495673518857527248912279381830119491298336733624406566430860213949463952247371907021798609437027705392171762931767523846748184676694051320005681271452635608277857713427577896091736371787214684409012249534301465495853710507922796892589235420199561121290219608640344181598136297747713099605187072113499999983729780499510597317328160963185950244594553469083026425223082533446850352619311881710100031378387528865875332083814206171776691473035982534904287554687311595628638823537875937519577818577805321712268066130019278766111959092164201989
>>>

 

shaoziyang
粉丝 16
博文 78
码字总数 28248
作品 0
武汉
程序员
私信 提问
加载中
此博客有 3 条评论,请先登录后再查看。
micropython esp8266 烧录

MicroPython 介绍 MicroPython 是一个 Python 3 语言的精简、高效实现,其包括 Python 标准库的一小部分,并经过优化,可以运行在微控制器和受限环境中运行。 MicroPython 包含了诸如交互式提...

osc_ox6w5nxo
2019/02/17
12
0
MicroPython与Python速度对比

首先说明,micropython跟python是没有任何可比性的,python作为一种通用的语言,在扩展性上不是micropython能比的,比如大量的库,可以方便的用C语言加模块提升速度,有pypy这样的带JIT的解释...

osc_s24g046g
2018/05/14
5
0
什么是 MicroPython ?它有什么优势?

MicroPython 和 Python 都是一类 Python 的动态语言解释器。 若是在系统里安装了 Python 解释器,例如在 Windows 或者 Linux 命令行输入 python,进入 python 解释器的命令行模式即可输入运行...

osc_qzrqhs9y
2019/09/13
8
0
TPYBoard开发板带你轻松玩转MicroPython

说到MicroPython,也许有人会感到陌生。而说到和它密切相关的Python,是否会恍然大悟呢?Python属于解释型语言,1989年才开发完成的Python显然非常年轻。然而经过数十年的磨砺,如今Python已...

bodasisiter
2019/03/28
512
0
玩转MicroPython开发板!你和Ta会有什么火花?

  说到MicroPython,也许有人会感到陌生。而说到和它密切相关的Python,是否会恍然大悟呢?Python属于解释型语言,1989年才开发完成的Python显然非常年轻。然而经过数十年的磨砺,如今Pyt...

bodasisiter
2017/06/03
1
0

没有更多内容

加载失败,请刷新页面

加载更多

Hacker News 简讯 2020-07-14

更新时间: 2020-07-14 00:00 I Know What You Download on BitTorrent - (iknowwhatyoudownload.com) 我知道你在BitTorrent上下载了什么 得分:69 | 评论:40 Show HN: Primo – all-in-one I......

FalconChen
47分钟前
99
0
绕过移动端系统限制的 dlopen 库 byOpen

byOpen是一个绕过移动端系统限制的增强版dlfunctions库。 支持特性 Android 支持App中加载和使用Android系统库接口(即使maps中还没有被加载也支持)。 Android 7以上dlopen, System.load都是...

shzwork
昨天
21
0
Golang学习系列第二天:变量、常量、数据类型和流程语句

继golang第一天后,今天学习下golang的变量、常量、数据类型和控制流语句。 做过其他编程语言(比如JavaScript,java,python)项目的话,其实很好理解变量、常量、数据类型和控制流。 变量也...

董广明
昨天
30
0
redis系列之——一致性hash算法

一致性hash算法你了解吗?什么时候使用?解决什么问题?redis集群模式使用了一致性hash算法了吗? 数据分片(sharding) 分布式数据存储时,经常要考虑数据分片,避免将大量的数据放在单表或...

诸葛小猿
昨天
15
0
IMDB是否提供API? [关闭] - Does IMDB provide an API? [closed]

问题: I recently found a movie organizer application which fetches its data from the IMDB database . 最近,我发现了一个电影管理器应用程序,该应用程序从IMDB数据库中获取其数据。 ...

fyin1314
昨天
14
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部