Matlab smooth函数与Python savgol_filter的差异及等效性咨询

输出是否一致？ #

直接使用你给出的两段代码，两者输出大概率不一致，核心差异体现在三个方面：

• 窗口长度的奇偶要求：Matlab的smooth函数允许span为偶数或奇数；但Scipy的savgol_filter要求window_length必须是奇数，传入偶数会直接报错。且Matlab遇到偶数span时，会自动将其加1转为奇数使用。
• 端点处理逻辑：Matlab对信号两端的拟合采用动态缩短窗口的方式，而Scipysavgol_filter默认使用插值模式（mode='interp'），二者的端点计算细节存在差异。
• 参数边界处理：当degree大于等于span时，Matlab和Scipy的报错、处理逻辑不同，Matlab会限制degree必须小于span，Scipy也有同样限制，但报错提示和内部处理细节有区别。

Matlab smooth('sgolay')的Python等效函数 #

要实现和Matlab完全对齐的效果，需要手动处理窗口长度的奇偶转换，并调整端点模式，以下是等效实现：

from scipy.signal import savgol_filter

def matlab_smooth_sgolay(y, span, degree):
    # 匹配Matlab对偶数span的自动转换逻辑
    if span % 2 == 0:
        span += 1
    # 确保多项式次数小于窗口长度，和Matlab行为一致
    if degree >= span:
        raise ValueError("Degree must be less than span")
    # 使用'nearest'模式匹配Matlab的端点处理逻辑
    return savgol_filter(y, window_length=span, polyorder=degree, mode='nearest')

测试时，将原Matlab代码中的y、span、degree参数传入这个函数，输出结果会和Matlab的yy高度一致。

内容的提问来源于stack exchange，提问作者Rocky