用matlab写一段SG滤波的程序

SG（Savitzky-Golay）滤波是一种平滑信号的方法，可以用于去噪和信号平滑。下面是一个用MATLAB实现SG滤波的例子：

function y_filt = sg_filter(x, window_size, order, deriv)
% x: 输入信号
% window_size: 窗口大小
% order: 多项式阶数
% deriv: 求导数的阶数

% 构造系数矩阵
half_size = (window_size - 1)/2;
x_axis = -half_size:half_size;
A = zeros(window_size, order+1);
for i = 0:order
    A(:,i+1) = x_axis.^i;
end

% 计算伪逆矩阵
AtA = A' * A;
Aty = A' * x';
coef = AtA \ Aty;

% 计算导数系数
if deriv > 0
    coef = coef(end:-1:1);
    for i = 1:deriv
        coef = polyder(coef);
    end
end

% 应用滤波器
y_filt = filtfilt(coef, 1, x);
end

使用方法：将上述代码保存为.m文件，并在MATLAB命令行中输入以下命令：

y_filt = sg_filter(x, window_size, order, deriv)

其中，x为输入信号，window_size为滤波窗口大小，order为多项式阶数，deriv为求导数的阶数。返回值y_filt为滤波后的信号。