双变量的精度仅为十点后的2个位置

Double variable has precision of only 2 places after decimal point

本文关键字：十点 2个 位置 精度 变量      更新时间：2023-10-16

我正在求解一个问题，该问题找到一个值k，该值k正确〜4个小数点。我将K初始化为双重。但是，在小数点之后只有2个位置。我想知道这种奇怪行为背后的原因。

我已经尝试了不同的数字，以尝试了解为什么会发生这种情况，但没有成功。

#include <bits/stdc++.h>
using namespace std;
double calc_err(double x, double y, double z, double a, double k)
{
    return ((y*y)/((z-k)*(z-k)) - (a*a*x*x)/(k*k) - (a*a) + 1.0);
}
int main()
{
    int t;
    cin >> t;
    while (t--)
    {
        double x, y, z, a;
        cin >> x >> y >> z >> a;
        double cerrp = 100.0;
        double hi = z;
        double lo = 0;
        while(1)
        {
            cout << "cerrp : " << cerrp << endl;
            double k  = (hi + lo)/2;
            if (cerrp  < 0.00005 && cerrp >= 0.0 )
            {
                cout << "k final : " << k << endl;
                break;
            }
            cout << "hi : " << hi << endl;
            cout << "lo : " << lo << endl;
            cout << "k : " << k << endl;
            cerrp = calc_err(x, y, z, a, k);
            if (cerrp >= 0) hi = k;
            else lo = k;

        }
    }
}

考虑以下输入：17390 8573 9287 285

在一定数量的迭代后，" hi"采用一个值9263.49和" lo"采用一个值9263.48。因此，我希望" k"会占据一个值（HI LO）/2 = 9263.485，但值为9263.49。另外，请注意，即使" k"没有更改，" CERRP"也会不断变化（所有其他参数都是恒定的）。为什么会发生？

编辑：这与浮点数学有所不同吗？因为我的问题与变量的精度有关，而链接中的问题与浮点操作的固有不准确有关。

只需使用std :: setprecision（）即可打印时更多数字

例如

int main()
{
    int t;
    cin >> t;
    while (t--)
    {
        double x, y, z, a;
        cin >> x >> y >> z >> a;
        double cerrp = 100.0;
        double hi = z;
        double lo = 0;
        while(1)
        {
            cout << "cerrp : " << std::setprecision(10) <<cerrp << endl;
            double k  = (hi + lo)/2;
            if (cerrp  < 0.00005 && cerrp >= 0.0 )
            {
                cout << "k final : " << std::setprecision(10) << k << endl;
                break;
            }
            cout << "hi : " << std::setprecision(10) << hi << endl;
            cout << "lo : " << std::setprecision(10) << lo << endl;
            cout << "k : " << std::setprecision(10) << k << endl;
            cerrp = calc_err(x, y, z, a, k);
            if (cerrp >= 0) hi = k;
            else lo = k;

        }
    }
}

输入1 7390 8573 9287 285

执行 
cerrp : 100
hi : 9287
lo : 0
k : 4643.5
cerrp : -286945.888
hi : 9287
lo : 4643.5
k : 6965.25
cerrp : -172643.8308
hi : 9287
lo : 6965.25
k : 8126.125
cerrp : -148345.0696
hi : 9287
lo : 8126.125
k : 8706.5625
...
hi : 9263.485111
lo : 9263.485111
k : 9263.485111
cerrp : -8.219334995e-05
hi : 9263.485111
lo : 9263.485111
k : 9263.485111
cerrp : 1.337719732e-05
k final : 9263.485111
相关文章：