31点泰勒加权幅度 代码

31点泰勒加权幅度的代码如下：

import numpy as np

def taylor_weights(n, order):
    # 构造泰勒加权系数矩阵
    weights = np.zeros((n, order+1))
    for i in range(n):
        mid = n//2
        for j in range(order+1):
            if i < mid:
                weights[i][j] = (-1)**(j) * np.math.factorial(order) / \
                                (np.math.factorial(j) * np.math.factorial(order-j)) * \
                                (i-mid)**(order-j)
            else:
                weights[i][j] = (-1)**(j) * np.math.factorial(order) / \
                                (np.math.factorial(j) * np.math.factorial(order-j)) * \
                                (i-mid+1)**(order-j)
    # 归一化
    weights /= np.sum(weights, axis=0)
    return weights

# 使用31点泰勒加权幅度来平滑一个信号
def smooth_signal(signal):
    n = len(signal)
    order = 4
    weights = taylor_weights(n, order)
    smoothed_signal = np.zeros(n)
    for i in range(order//2, n-order//2):
        smoothed_signal[i] = np.dot(weights[i,:], signal[i-order//2:i+order//2+1])
    return smoothed_signal

这里的`taylor_weights`函数用于计算31点泰勒加权系数矩阵，`smooth_signal`函数则使用该系数矩阵来平滑一个输入信号。需要注意的是，这里的泰勒加权系数矩阵是针对31个点的情况进行计算的，如果需要使用其他点数的泰勒加权，需要相应地修改代码。