使用python帮我写实现下述功能的代码：对于给定的相位连通域，使用加权最小二乘拟合直线，得到直线的两端坐标，相位联通域是一个Box2D结构的矩形

假设你已经安装了 NumPy、SciPy 和 Matplotlib 库，可以使用以下 Python 代码实现你所需的功能：
import numpy as np
from scipy.stats import linregress
import matplotlib.pyplot as plt

# 定义相位连通域的左上角和右下角坐标
x1, y1 = 0, 0
x2, y2 = 10, 10

# 构造相位连通域内的坐标点
x = np.arange(x1, x2+1)
y = np.arange(y1, y2+1)
xx, yy = np.meshgrid(x, y)
points = np.column_stack([xx.ravel(), yy.ravel()])

# 构造点的权重，这里使用均匀分布的权重
weights = np.ones(len(points)) / len(points)

# 对点进行加权最小二乘拟合直线
slope, intercept, _, _, _ = linregress(points, np.arange(len(points)), weights)

# 计算直线的两个端点
x_start = x1
y_start = slope * (x_start - intercept)
if y_start < y1:
y_start = y1
x_start = (y_start + intercept) / slope

x_end = x2
y_end = slope * (x_end - intercept)
if y_end > y2:
y_end = y2
x_end = (y_end + intercept) / slope

# 绘制相位连通域和拟合的直线
fig, ax = plt.subplots()
ax.plot([x_start, x_end], [y_start, y_end], 'r')
ax.set_xlim(x1, x2)
ax.set_ylim(y1, y2)
plt.show()

这段代码首先根据相位连通域的左上角和右下角坐标生成一个坐标点网格，然后使用这些点进行加权最小二乘拟合直线。最后，计算出直线的两个端点并绘制相位连通域和拟合的直线。