使用python 进行Sen‘s斜率计算

要使用Python计算Sen's斜率，可以使用以下代码：

import numpy as np

def sen_slope(y):
    n = len(y)
    x = np.linspace(1, n, n)
    x_bar = np.mean(x)
    y_bar = np.mean(y)
    sxy = np.sum((x - x_bar) * (y - y_bar))
    sxx = np.sum((x - x_bar) ** 2)
    return sxy / sxx

# 示例
y = np.array([8, 7, 6, 5, 4, 3, 2, 1])
slope = sen_slope(y)
print(slope)

在这个例子中，输入变量y是一个numpy数组，表示一个时间序列。sen_slope函数将计算Sen's斜率，并返回结果。

在函数中，我们首先计算了时间序列的长度n，并生成一个等差数列x，用于表示时间。然后，我们计算了x和y的均值x_bar和y_bar，并用它们来计算Sxy和Sxx。最后，我们将Sxy除以Sxx，得到Sen's斜率。

在示例中，我们使用了一个简单的时间序列，其中y的值从8到1递减。运行代码后，我们得到Sen's斜率为-1.0。