cos给出了错误的价值

时间:2015-03-17 15:35:12

标签: c++ cos

#include <iostream>
#include <math.h>
#include "boost/math/constants/constants.hpp"
const double pi = boost::math::constants::pi<double>();

int main() {
    double x = 2;
    double y = 1;
    double angle = 90; //in degrees
    double rad { angle * pi / 180 }; //converting to rads
    std::cout << "rad: " << rad << std::endl; //1.5708 - OK
    double c = cos( rad );
    std::cout << "cos(rad): " << c << std::endl; //6.12303e-017 - Huh? Should be ~ -3.2051033e-9
    double s = sin( rad );
    std::cout << "sin(rad): " << s << std::endl; //1 - OK
    x = x * c - y * s;
    y = x * s + y * c;
    std::cout << "(" << x << "," << y << ")" << std::endl;

    return 0;
}

所以这很奇怪,我正在尝试轮换,但cos给了我错误的价值。它应该接近-3.2051033e-9,但我最终得到6.12303e-017。我已将角度转换为弧度,然后检查出来。 sin( rad )也会返回正确的值。 是什么给了什么?

编辑: 代码有错误但与问题无关。在y = x * s + y * c;上,x应该引用顶部声明的值,但是前面的等式会覆盖该值,最后的y将是错误的。

感谢对指数的澄清。我有点困惑。

3 个答案:

答案 0 :(得分:4)

cos(90)(度)大约为零,如0.0中所示。担心获得e-17 vs e-9的数字是没有意义的,它们基本上都是零。

结果实际的接近程度取决于它们代表pi90.0180.0的方式以及cos的实施情况math.h中的函数。

答案 1 :(得分:2)

  

它应该接近-3.2051033e-9但我最终得到的是6.12303e-017。

我不知道为什么你认为它应该是-3.2051033e-9&#34;因为cos(90º)为零,-0.0000000032051033远远超出了计算误差范围,

您计算的错误范围内,您获得的值0.0000000000000000612303 ,因此正确

答案 2 :(得分:0)

您的代码中没有错误..

在cos函数(cos)中,返回值是实数。以下是一些例子;

(cos -1) returns 0.540302
(cos 0.0) returns 1.0
(cos (* pi 0.5)) returns 6.12303e-017 which is practically equal to **0**

你也可以用 cos <6.1> cos(6.12303e-017)用 cos检查它,你会发现 90 度< / p>