我假设，通过阅读你的代码，你得到了每个x1-y1对之间的距离。所以你应该用毕达哥拉斯定理，也就是a^2+b^2=c^2，而不是a^2-b^2=c^2。
将距离计算更改为 Dis1[i][j] = sqrt((x1[i] - x1[j])**2 + (y1[i] - y1[j])**2) 实际上，你试图取可能为负数的平方根，它返回nan。

如果你真的想计算矩阵 Dis1 ，则必须使用复数。

import numpy as np
x1=[0.0,0.0,1.0,0.5,2.0,10.0,11.0,12.0,8.0,8.9]
x1=np.array(x1, dtype=complex)
y1=[1.0,1.5,1.0,1.2,2.0,9.0,13.0,8.0,8.0,9.0]
y1=np.array(y1, dtype=complex)
Dis1=np.zeros((10,10), dtype=complex)
for i in range(10):
    for j in range(10):
        Dis1[i][j]=sqrt((x1[i]-x1[j])**2-(y1[i]-y1[j])**2)

在你的问题中，每次当你试图计算 sqrt 输出为负数 nan . 对于复数，虚部是非零的。

Dis1.imag
array([[ 0.        ,  0.5       , -0.        ,  0.        ,  0.        ,
         0.        ,  4.79583152,  0.        ,  0.        ,  0.        ],
       [ 0.5       ,  0.        , -0.        , -0.        ,  0.        ,
         0.        ,  3.35410197,  0.        ,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ,
         0.        ,  6.63324958,  0.        ,  0.        ,  1.26095202],
       [ 0.        ,  0.        , -0.        ,  0.        ,  0.        ,
         0.        ,  5.38423625,  0.        ,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ,
         0.        ,  6.32455532,  0.        ,  0.        ,  1.17898261],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ,
         0.        ,  3.87298335, -0.        ,  0.        ,  0.        ],
       [ 4.79583152,  3.35410197,  6.63324958,  5.38423625,  6.32455532,
         3.87298335,  0.        , -4.89897949,  4.        ,  3.40440891],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ,
         0.        ,  4.89897949,  0.        ,  0.        ,  0.        ],
       [ 0.        ,  0.        ,  0.        ,  0.        ,  0.        ,
         0.        ,  4.        , -0.        ,  0.        ,  0.43588989],
       [ 0.        ,  0.        ,  1.26095202,  0.        ,  1.17898261,
        -0.        ,  3.40440891, -0.        ,  0.43588989,  0.        ]])