python求过圆心直线与圆的两个交点

假设圆的圆心坐标为(x0, y0)，半径为r，直线的一般式方程为Ax+By+C=0，我们可以通过以下步骤求解交点：

1. 将直线的一般式方程转化为斜截式方程，即y = -A/B * x - C/B。

2. 将斜截式方程代入圆的方程，得到一个关于x的二次方程，解出x的两个解。

3. 将x的两个解代入斜截式方程，得到对应的y坐标。

代码如下：

import math

# 圆的圆心坐标和半径
x0 = 0
y0 = 0
r = 1

# 直线的一般式方程Ax+By+C=0
A = 1
B = -1
C = 0

# 将直线的一般式方程转化为斜截式方程
k = -A / B
b = -C / B

# 将斜截式方程代入圆的方程，得到一个关于x的二次方程
a = k ** 2 + 1
b = 2 * x0 * k - 2 * y0 * k - 2 * b * k
c = x0 ** 2 + y0 ** 2 + b ** 2 / (4 * a) - r ** 2

# 解出x的两个解
x1 = (-b + math.sqrt(b ** 2 - 4 * a * c)) / (2 * a)
x2 = (-b - math.sqrt(b ** 2 - 4 * a * c)) / (2 * a)

# 将x的两个解代入斜截式方程，得到对应的y坐标
y1 = k * x1 + b
y2 = k * x2 + b

print("交点1坐标：", (x1, y1))
print("交点2坐标：", (x2, y2))

注意，以上代码只适用于直线和圆有两个交点的情况。如果直线和圆没有交点或者有一个交点，代码会出现异常。