使用泰勒级数近似罪恶

Question

我正在尝试使用泰勒级数而不使用阶乘来计算 sin(x)。

import math, time
import matplotlib.pyplot as plot

def sin3(x, i=30):
    x %= 2 * math.pi
    n = 0
    dn = x**2 / 2
    for c in range(4, 2 * i + 4, 2):
        n += dn
        dn *= -x**2 / ((c + 1) * (c + 2))
    return x - n

def draw_graph(start = -800, end = 800):
    y = [sin3(i/100) for i in range(start, end)]
    x = [i/100 for i in range(start, end)]

    y2 = [math.sin(i/100) for i in range(start, end)]
    x2 = [i/100 for i in range(start, end)]

    plot.fill_between(x, y, facecolor="none", edgecolor="red", lw=0.7)
    plot.fill_between(x2, y2, facecolor="none", edgecolor="blue", lw=0.7)
    plot.show()

当你运行 draw_graph 函数时，它使用 matplotlib 绘制图形，红线是我的 sin3 函数的输出，蓝线是我的函数的正确输出math.sin 方法。

如您所见，曲线不太正确，不够高或不够低（似乎在 0.5 处达到峰值），并且还有奇怪的行为，它在 0.25 左右产生一个小峰值然后再次下降。如何调整我的函数以匹配 math.sin 的正确输出？

Answer 1

你的 sin(x) 方程有误，而且循环不变量也搞砸了。

sin(x)的公式是x/1! - x^3/3! + x^5/5! - x^7/7!...，所以我真的不知道你为什么要将dn初始化为涉及x^2的东西。

你也想问问自己：我的循环不变量是什么？当我到达循环开始时 dn 的值是多少。从您更新 dn 的方式可以清楚地看出您希望它涉及 x^i / i!。然而在循环的第一次迭代中，i=4，但 dn 涉及 x^2.

这是您要写的内容：

def sin3(x, i=30):
    x %= 2 * math.pi
    n = 0
    dn = x
    for c in range(1, 2 * i + 4, 2):
        n += dn
        dn *= -x**2 / ((c + 1) * (c + 2))
    return n

使用泰勒级数近似罪恶

Approximating sin using the Taylor series

python

trigonometry

taylor-series