我收到一个错误 <string>:149: RuntimeWarning: 在生成列表时在 sqrt 中遇到无效值

Question

def ellipse(numPoints, genX=np.linspace, HALF_WIDTH=10, HALF_HEIGHT=6.5):
    xs = 10.*genX(-1,1,numPoints)
    ys = 6.5*np.sqrt(1-(xs**2))
    return(xs, ys, "-")

我收到一条错误消息，指出在平方根中遇到了无效值。我看不到它是什么。

sqrt(0) = 0
6.5*sqrt(1- (-1**2)) = 0

它们应该可以工作，但是 y 值有问题，它们返回“nan”

Answer 1

可能xs**2 返回一个数字> 1 带有负数的 sqrt 将返回 nan（不是数字）

>>> import numpy as np
>>> np.sqrt(-1)
nan

如果我是对的，numpy 提供了复数功能，我认为这是表示 x 的 sqrt(x) 的唯一方法