1

我必须编写一个函数，使用以下系列的前 10 项计算并返回角度的余弦： cosx = 1 - (x**2)/2! + (x**4)/4! - (x**6)/6!....

我不能使用阶乘函数，但我可以使用这样一个事实，即如果前一个分母是n!，那么当前分母就是n!(n+1)(n+2)。我正在尝试使用累加器循环，但我很难接受它从正到负交替并且分母也有问题的事实。

这就是我迄今为止所拥有的。对分母和累加器循环有任何帮助吗？

def factorial(x):

    if(x == 0): 
        return 1
    return x * factorial(x-1)


def cosine(angle):

    cosx = 1
    sign = -1
    for i in range(2, 20, 2):
        cosx = cosx + (sign*(angle**i))/factorial(i)
        sign = -sign
    return cosx
4

4 回答 4

2

注意：如果你使用的是Python2.x，你应该使用

from __future__ import division

作为文件的第一行

一种方法是像这样交替标志

def cos(angle):
    sign = -1
    cosx = 1
    for i in range(2,20,2):
         cosx += sign*(x**i)/(i)
         sign = -sign

您仍然需要正确处理阶乘部分

这是一个更简单的版本，它根据前一个术语计算每个术语。

def cos(x):
    res = 0
    term = 1
    for i in range(1, 20, 2):
        res += term
        term *= -x * x/ i /(i + 1)
    return res
于 2013-11-25T20:59:59.540 回答
2

也许是这样的：

#! /usr/bin/python3.2

def cos (a):
    d = 1
    c = 1
    for i in range (2, 20, 2):
        d *= i * (i - 1)
        sign = -1 if i % 4 else 1
        print ('adding {} * a ** {} / {}'.format (sign, i, d) )
        c += sign * a ** i / d
        print ('cosine is now {}'.format (c) )
    return c

cos (1.0)

基本上d（如分母）是你的累加器。

于 2013-11-25T20:55:02.710 回答
1

为什么不为阶乘创建一个函数？（我为重复添加了一个可选参数）：

编辑：

根据评论，这是一个“学习”阶乘函数，它将确保不会重新计算任何值（当然，这最多只能工作 199！，我假设 cos 函数的重复次数永远不会超过 99 次，因为它可能已经在 90 左右抛出溢出错误）：

facts = [0]*200
facts[0] = 1
def factorial(x): 
    if(facts[x] != 0): return facts[x]
    return x * factorial(x-1)

def cos(x, reps=10):
    final_val = 1
    neg = -1
    for n in range(2, reps*2, 2):
        final_val += neg*(x**n)/factorial(n)
        neg *= -1
    return final_val

print cos(3.14/2)
于 2013-11-25T21:29:34.007 回答
0

对于计算系列，请使用列表推导的力量。

In [1]: a = range(1,11)

In [2]: a
Out[2]: [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]

该mul函数将给定列表的元素相乘：

In [3]: def mul(r):
   ...:     rv = 1
   ...:     for j in r:
   ...:         rv *= j
   ...:     return rv
   ...:

现在您可以生成阶乘：

In [7]: b = [mul(a[:i]) for i in a]

In [8]: b
Out[8]: [1, 2, 6, 24, 120, 720, 5040, 40320, 362880, 3628800]

从列表中选择正确的元素；

In [10]: a[1::2]
Out[10]: [2, 4, 6, 8, 10]

In [11]: b[1::2]
Out[11]: [2, 24, 720, 40320, 3628800]

要产生优点和缺点：

In [12]: [(-1)**k for k in a[:5]]
Out[12]: [-1, 1, -1, 1, -1]

结合所有元素；

In [14]: data = zip(a[1::2], b[1::2], [(-1)**k for k in a[:5]])

然后对某一个进行计算x；

In [21]: x = 12

In [22]: sum([s*x**t/float(b) for t, b, s in data])
Out[22]: -9753.737142857142
于 2013-11-25T22:28:46.127 回答