直接用距离公式等于半径判定圆上点几乎总会失败,因浮点舍入误差使理论在圆上的点计算距离不等于r;应改用容差比较,如abs(dist - r)
浮点数比较导致的“看似在圆上却判定失败”
直接用
sqrt((x - cx)*(x - cx) + (y - cy)*(y - cy)) == r判定几乎总会失败——浮点运算存在舍入误差,哪怕点理论上在圆上,计算出的距离也大概率不等于r。哪怕用abs(dist - r) ,仍可能因尺度失配误判(比如半径是 <code>1e6时,1e-9远小于机器精度相对误差)。正确做法是避免开方,改用平方距离与半径平方比较,并引入相对容差:
- 计算平方距离:
(x - cx)*(x - cx) + (y - cy)*(y - cy)- 计算
r * r(注意:若r是浮点数,确保不重复计算)- 判断:
abs(sq_dist - r_sq) ,其中 <code>epsilon取1e-10或更小,取决于坐标量级C++中推荐的健壮判定函数写法
封装成函数时,需处理边界情况(如
r 、无穷坐标),并统一使用 <code>double避免单精度累积误差:bool onCircle(double x, double y, double cx, double cy, double r) { if (r double dx = x - cx; double dy = y - cy; double sq_dist = dx*dx + dy*dy; double r_sq = r * r; double diff = std::abs(sq_dist - r_sq); double scale = std::max({sq_dist, r_sq, 1.0}); return diff <p>}</p><p>关键点:</p>
- 不用
sqrt,避免额外误差和性能损耗- 容差随尺度自适应(
1e-12 * scale),比固定阈值更可靠- 显式检查
std::isfinite,防止 NaN 或 inf 传播当需要更高精度或整数坐标时的特殊处理
如果所有输入都是整数(如像素坐标、网格点),且半径平方也是整数,可完全避免浮点误差:
- 用
long long存储平方距离和r*r,防止溢出(注意:坐标范围大时需检查溢出)- 直接判断
sq_dist == r_sq,无需容差- 示例:
onCircleInt(3, 4, 0, 0, 5)→(3-0)^2 + (4-0)^2 == 25→true但一旦涉及除法、旋转、缩放等操作,坐标大概率转为浮点,此时必须回归带容差的平方比较。
常见误判场景与调试建议
实际调试时,常发现“明明画出来就在圆上,代码却返回 false”。典型原因包括:
- 传入的
r是直径而非半径(检查变量命名和单位一致性)- 坐标系 y 轴方向反了(如 OpenGL 与图像坐标系差异),导致
y - cy符号错- 使用了
float类型,中间计算溢出或精度丢失(尤其大坐标时)- 调试打印用了
std::cout 默认只显示 6 位小数,掩盖了真实误差值建议打印完整精度:
std::cout真正难的不是公式本身,而是意识到“精确位于周界线”在计算机里从来不是数学意义上的相等,而是一个尺度敏感的近似区间。
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!












