【Python】高精度三角函数计算实现方法

瞎扯

这几天我又修仙了一下,查阅了大量资料和进行了不少于百次的编程和运行,终于有脸在这为大家提供新的东西——高精度三角函数的 Python 实现方案。

说来也是,本来只用取前几位值进行计算或直接调用 math 模块就可以实现的东西,我却硬生生直接手动做了出来。在旁人看起来用处可能不大,但是我却从中学到了不少东西。接下来看我慢慢介绍吧。

代码我开源到 Github 和 Gitee 上了。

I've Opened Source it onto Github & Gitee.

If you need English version, please use translation function in your browser or contact me.

正文

介绍

文件名虽然只包含三角函数,但是因为平时数学计算需要,我也夹了一些大杀器在里面:最主要的就是三种计算精确圆周率的方式。

使用

首先你得获取我开源的 HighPrecisionTrigonometricFunctions.py 文件,然后在你需要的地方通过以下命令导入:

import HighPrecisionTrigonometricFunctions as hptf

依赖

该项目依赖库有 decimalsys ,都是 Python 预装库,无需额外安装。

函数文档

想要查看各个函数的文档和常量的值,你可以直接以运行脚本的形式运行 HighPrecisionTrigonometricFunctions.py 文件,非英语用户懒得看就别看了。

Translation: If you want to view the document of functions, you could run the HighPrecisionTrigonometricFunctions.py file as a script directly.

函数

正弦,余弦,正切

这三个乃是三角函数之基础,它们分别代表单位圆中,任意弧度角与单位圆的交点的纵坐标,横坐标和纵坐标比横坐标。

调用

要调用模块中的这三个三角函数,你只需要在导入模块后分别调用 sin() cos()tan() 即可。

实现

这三个函数中,正弦函数和余弦函数的实现方法是使用泰勒展开式Taylor Expansion);正切函数则分别调用正弦和余弦再输出返回比值。

参数

对于每个三角函数(下同),你需要传入一个角度 x ,它是弧度制Radian Measure),若你需要计算角度制Degree Measure),请调用编写好的 to_radian() (计算公式 radian = (degree / 180) * pi)函数将其转换为弧度制。类型可以是 int|float|Decimal(即整数、浮点数、十进制小数)。

对于第二个参数 precision ,是精度,他决定了计算多少位和返回多少位数字(尽管程序用上了 decimal 库,但计算结果的后几位总是具有摆尾性的,你可以选择性忽略后几位,或直接在某一位上四舍五入)。传入类型必须是整数!

关于圆周率(π)的使用

弧度制使用过程中经常需要用到圆周率,对此我有两个靠谱的方式:第一个,自己手输圆周率(有多少输多少,我背了六十多位:3.141592653587293238462643383279502884197169399375105820974974592307...);另外一个请看下面:

模块内提供了三种比较准确的方式计算圆周率,同时也提供了常量 pi 以便计算(虽说是常量,但也是直接调用三种方式中最快的一种进行计算,默认取一万位有效数字)。详细的内容后面会说。

余割、正割、余切

这三个比较不常用,但是和上面的三个紧密相连(分别是正弦、余弦和正切的倒数),所以,,函数内部直接分别返回了 1 / (sin|cos|tan) 的值。

调用

分别调用 sec csccot 即可。

参数

同正弦余弦和正切。

sin2 & cos2

这两个方法分别适用于已知余弦求正弦、已知正弦求余弦。

参数

cosine 值(for sin2);sine 值(for cos2)

类型可为字符串、整数、浮点数、Decimal 对象。

返回

列表,组成元素分别为一正一负的相加等于零的数(谁家同一正弦/余弦值只有一个对应的余弦/正弦值??)

圆周率计算

终于讲到圆周率了,关于实现方法,我查了大量资料,也问了 ChatGPT,最终确认以下两种方式:手动实现参考代码实现梅钦公式;BBP 公式。

反正切函数~~

反三角函数之一,用于求已知正切值对应的角度,值域为 (-pi/2, pi/2)。

在模块内,我保留了 arctan 反正切函数的程序函数。

参数~~

还是那两样:正切值(x,整数/浮点数/ Decimal 对象)和精度(precision,整数)

返回~~

Decimal 对象,根据精度的弧度制角度数。

梅钦公式(Machin Formula)

此公式计算速度较快,较推荐使用。

程序内提供了 pi_Machin 函数进行计算圆周率。

仅需传入精度 precision (整数)。

返回的是指定位数(有效数字位数)的圆周率。

返回示例:

print(pi_Machin(1000))
# 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117067982148086513282306647093844609550582231725359408128481117450284102701938521105559644622948954930381964428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273724587006606315588174881520920962829254091715364367892590360011330530548820466521384146951941511609433057270365759591953092186117381932611793105118548074462379962749567351885752724891227938183011949129833673362440656643086021394946395224737190702179860943702770539217176293176752384674818467669405132000568127145263560827785771342757789609173637178721468440901224953430146549585371050792279689258923542019956112129021960864034418159813629774771309960518707211349999998372978049951059731732816096318595024459455346908302642522308253344685035261931188171010003137838752886587533208381420617177669147303598253490428755468731159562863882353787593751957781857780532171226806613001927876611195909216420198

注意:尽管您设置了精度,您也应该知道要舍弃结果的后几位,因为毕竟是个无理数,在不完全计算出来的情况下总是会有摆尾性的。当然,您也可以尝试设置精度为多几位,再选择性地屏蔽那多出地几位。

这个提示同样适用于后面的 BBP 公式。

贝利-波尔温-普劳夫公式(BBP Formula)

此公式较以上梅钦公式计算较慢,荐备用。

参数和返回同上,不再赘述,只是时间慢些。


END

本项目和文档为草草而作且未作多检查,如若有任何出错,请指出。

If you have troubles in using or find bus/mistakes, please point out!


阅读原文 | 查看更多

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 159,458评论 4 363
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,454评论 1 294
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 109,171评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 44,062评论 0 207
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,440评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,661评论 1 219
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,906评论 2 313
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,609评论 0 200
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,379评论 1 246
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,600评论 2 246
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,085评论 1 261
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,409评论 2 254
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,072评论 3 237
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,088评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,860评论 0 195
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,704评论 2 276
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,608评论 2 270

推荐阅读更多精彩内容