Question

我有这个用C语言编写的学术项目。练习之一是不使用math.h来计算cos（x）。为此，我们得到了cos（x）的一个数列（我认为这是泰勒/麦克劳林）。对于标准输入，我们有x个角度和k次迭代。总和为（（-1）^ n * x ^（2n））/（2n）的n = 0至k-1！

我尝试过改变周期和不进行周期，并且摆弄了无济于事的变量。

以下是代码（已更新）：

#include <stdio.h>

int fat(int num) {
    int fat_num=1;
    for (int cnt=1; cnt<=num; cnt++) fat_num*=cnt;
    return fat_num;
}

float expn(int x, int y) {
    int x_y=x;
    for (int cnt=1; cnt<y; cnt++) x_y*=x;
    return x_y;
}

int main() {
    const float pi = 3.1415926;
    float x; 
    scanf("%f",&x);
    int k; 
    scanf("%d", &k);
    float cosx = 0;
    int n = 0;
    x *= pi/180;
    while (n <= k-1) {
        if (n%2 == 0) 
            cosx +=    expn(x,2*n)/fat(2*n);
        else 
            cosx += -1*expn(x,2*n)/fat(2*n);
        n++;
    }
    printf("%f", cosx);
    return 0; 
}

输入96可以得到0.54的余弦值。 哪个不对。

已解决：该错误发生在expn中，必须对其进行更新才能浮动！

Answer 1

2个问题

expn(x,2*n)/fat(2*n)使用整数除法。

float expn(int x, int y)使用int x。需要浮点x。

// int expn(int x, int y) {
float expn(float x, int y) {
  float x_y = 1;
  for (int cnt = 0; cnt < y; cnt++)
    x_y *= x;
  return x_y;
}

代码具有各种效率。

不需要k的替代方法-需要一些思考。

static double my_cos_helper(double xx, double term, unsigned n) {
  if (term + 1.0 == 1.0) {
    return term;
  }
  return term - xx * my_cos_helper(xx, term / ((n + 1) * (n + 2)), n + 2);
}

double my_cos(double x) {
  return my_cos_helper(x * x, 1.0, 0);
}

条款汇总未返回正确结果

1 个答案: