如何排查NumPy中arccos的RuntimeWarning报错？为何正态分布会触发？

为什么正态分布会触发警告？ #

np.arccos()的输入必须落在**[-1, 1]**范围内，超出这个区间的数值会导致计算无效，返回nan并抛出警告。
你用的np.random.normal(-1, 1, N)是生成均值为-1、标准差为1的正态分布随机数，这种分布的取值范围是无限延伸的——虽然99.7%的数值集中在[-3, 1]区间内，但仍会有部分数值小于-1或者大于1，这些超出范围的数传入arccos()就会触发警告。
而你使用均匀分布时，应该是限定在了[-1, 1]区间内生成数值，自然不会超出定义域，也就不会出现警告。

排查与解决方法 #

• 先验证随机数的范围：添加代码打印生成数值的极值，确认确实存在超出[-1,1]的情况：

  vals = np.random.normal(-1, 1, N)
  print(f"最小值: {vals.min()}, 最大值: {vals.max()}")
  

• 处理超出范围的数值：
  1. 截断法：把超出[-1,1]的数强制拉到边界上，不会丢失数据点：

     vals = np.random.normal(-1, 1, N)
     vals = np.clip(vals, -1, 1)  # 将数值限制在[-1,1]区间内
     tau = np.arccos(vals)
     

  2. 过滤法：直接删除超出范围的数，适合对数据完整性要求不高的场景：

     vals = np.random.normal(-1, 1, N)
     valid_vals = vals[(vals >= -1) & (vals <= 1)]
     tau = np.arccos(valid_vals)
     

内容的提问来源于stack exchange，提问作者Lavender_square_2