[c语言日寄]浮点数的排序

siy2333

发布于 2025-03-17 05:36:47

7700

代码可运行

文章被收录于专栏：来自csdn的博客来自csdn的博客

运行总次数：0

代码可运行

前言

在C语言的编程实践中，浮点数的排序是一个常见且重要的问题。无论是科学计算、数据分析还是图形处理，我们常常需要对浮点数数组进行排序。然而，浮点数的排序并不像整数排序那样简单直接，因为浮点数的精度问题和特殊值（如NaN和无穷大）需要特别处理。

今天，我们将通过一个具体的例子来探讨如何在C语言中对浮点数数组进行排序。我们将从基础的排序算法入手，逐步深入到浮点数排序的细节和注意事项，并提供一些拓展应用的思路。

题目引入

对以下浮点数进行排序

double arr_double[12] = {
		0.0,               // 零
		-0.0,              // 负零
		1.0,               // 正数
		-1.0,              // 负数
		3.14159,           // 正数（π）
		-2.71828,          // 负数（e）
		INFINITY,          // 正无穷
		-INFINITY,         // 负无穷
		NAN,               // 非数字值（NaN）
	};

知识点分析

1. 浮点数的表示和精度问题

浮点数在计算机中是以二进制形式存储的，其表示精度是有限的。根据IEEE 754标准，单精度浮点数（float）有32位，其中1位用于符号位，8位用于指数位，23位用于尾数位；双精度浮点数（double）有64位，其中1位用于符号位，11位用于指数位，52位用于尾数位。

这种表示方式导致了浮点数的精度问题。例如，0.1 在计算机中无法精确表示，其二进制表示为一个无限循环小数，存储时会被截断或舍入。因此，直接比较两个浮点数是否相等可能会因为微小的误差而失败。

2. 浮点数的特殊值

浮点数有几种特殊值，包括：

NaN（Not a Number）：表示未定义或不可表示的数值。例如，0.0 / 0.0 和 sqrt(-1.0) 的结果都是NaN。
正无穷（INFINITY）和负无穷（-INFINITY）：表示数值过大或过小，超出了浮点数的表示范围。

这些特殊值在排序时需要特别处理，因为NaN与任何数（包括自身）比较的结果都是false，而无穷大值需要放在数组的开头或结尾。

3. 浮点数排序算法

在C语言中，可以使用标准库函数qsort来对浮点数数组进行排序。qsort需要一个比较函数来确定排序规则。对于浮点数排序，比较函数需要处理浮点数的精度问题和特殊值。

以下是一个简单的浮点数比较函数：

#define EPS 1e-6

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

    if (isnan(*pa)) {
        return isnan(*pb) ? 0 : 1; // 如果a是NaN，b也是NaN则相等，否则a大于b
    }
    if (isnan(*pb)) {
        return -1; // 如果b是NaN，a不是NaN，则a小于b
    }

    if (fabs(*pa - *pb) < EPS) {
        return 0; // 相等
    } else {
        if (*pa < *pb) {
            return -1; // a小于b
        } else {
            return 1; // a大于b
        }
    }
}

这个比较函数首先检查NaN值，然后使用一个非常小的阈值EPS来比较两个浮点数是否相等。

注意事项

1. 阈值的选择

在比较浮点数时，阈值EPS的选择非常重要。如果EPS过大，可能会导致两个明显不同的浮点数被误判为相等；如果EPS过小，可能会因为浮点数的精度问题导致比较失败。

通常，EPS的选择需要根据具体的应用场景来决定。对于一般的科学计算，1e-6是一个常用的阈值。

2. 特殊值的处理

在排序时，需要特别处理NaN和无穷大值。NaN值通常被放在数组的最后，而无穷大值可以放在数组的开头或结尾。

以下是一个处理特殊值的比较函数：

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

    if (isnan(*pa)) {
        return isnan(*pb) ? 0 : 1; // 如果a是NaN，b也是NaN则相等，否则a大于b
    }
    if (isnan(*pb)) {
        return -1; // 如果b是NaN，a不是NaN，则a小于b
    }

    if (isinf(*pa) && isinf(*pb)) {
        return (*pa > *pb) ? -1 : 1; // 如果a和b都是无穷大，比较它们的符号
    }
    if (isinf(*pa)) {
        return -1; // 如果a是无穷大，a小于b
    }
    if (isinf(*pb)) {
        return 1; // 如果b是无穷大，a大于b
    }

    if (fabs(*pa - *pb) < EPS) {
        return 0; // 相等
    } else {
        if (*pa < *pb) {
            return -1; // a小于b
        } else {
            return 1; // a大于b
        }
    }
}