6

我有以下用于查找四分位数的代码：

#include <stdio.h>
#include <stdlib.h>

typedef struct {
    double qrt[3];
    double *value;
    int count;
} t_data;

static void set_qrt(t_data *data, int qrt)
{
    int n, e;
    double d;

    d = qrt * 0.25 * data->count + 0.5;
    n = (int)d;
    e = n != d;
    data->qrt[qrt - 1] = data->value[n - 1];
    if (e) {
        data->qrt[qrt - 1] += data->value[n];
        data->qrt[qrt - 1] *= 0.5;
    }
}

static void set_qrts(t_data *data)
{
    set_qrt(data, 2);
    if (data->count > 1) {
        set_qrt(data, 1);
        set_qrt(data, 3);
    } else {
        data->qrt[0] = 0.0;
        data->qrt[2] = 0.0;
    }
}

static int comp(const void *pa, const void *pb)
{
    const double a = *(const double *)pa;
    const double b = *(const double *)pb;

    return (a > b) ? 1 : (a < b) ? -1 : 0;
}

int main(void)
{
    double values[] = {3.7, 8.9, 7.1, 5.4, 1.2, 6.8, 4.3, 2.7};
    t_data data;

    data.value = values;
    data.count = (int)(sizeof(values) / sizeof(double));
    qsort(data.value, data.count, sizeof(double), comp);
    set_qrts(&data);
    printf("Q1 = %.1f\nQ2 = %.1f\nQ3 = %.1f\n", data.qrt[0], data.qrt[1], data.qrt[2]);
}

是

d = qrt * 0.25 * data->count + 0.5;
n = (int)d;
e = n != d;

保证按预期工作？(e == isinteger(d))

4

4 回答 4

8

数字0.5、0.25等0.125表示 2 的负幂，因此可以在IEEE 754 类型中精确表示。使用这些数字不会导致表示错误。

于 2013-07-30T10:44:06.970 回答
3

值 0.5 和 0.25 本身将是准确的。您计算的中间值可能是也可能不是，取决于它们的范围。IEEE double 有一个 52 位尾数，因此它们将精确地表示为尾数中需要 50 位或更少的 0.25 个数字，即大约 15 个十进制数字。

因此，如果将 0.25 添加到 100000000000000 (10^14)，您将得到 100000000000000.25。但是如果你把 0.25 加到 10000000000000000 (10^16) 上，你就会失去分数。

于 2013-07-30T11:04:40.257 回答
1

dasblinkenlight is absolutely correct. Double/float and integer types are stored differently according to IEEE754. Watch this for an easy tutorial if you are curious about it.

于 2013-07-30T10:47:28.760 回答
1

双精度浮点格式的主要内容有 53 位，其中一位是隐含的。这意味着它可以表示 2^0 到 2^53-1 范围内的所有正整数和负整数。

0（零）是一种特殊情况，它有自己的格式。

当涉及 0.25 间距时，范围直接计算为 2^-2 到 2^51-0.25。这意味着相当多但绝不是所有 0.25 的倍数都可以在双精度格式中精确表示，就像相当多但并非所有整数都可以精确表示一样。

因此，如果您的可精确表示的间距为 2^x，则可表示的范围为 2^x 到 2^(53+x)-2^x。

于 2013-07-30T13:24:24.103 回答