多项式平滑：Savitzky-Golay滤波

多项式平滑

%五点三次平滑法
n = length(x);
a = x;
for k = 1:50
    b(1) = (69*a(1)+4*(a(2)+a(4))-6*a(3)-a(5))/70;
    b(2) = (2*(a(1)+a(5))+27*a(2)+12*a(3)-8*a(4))/35;
    for j = 3:n-2
        b(j) = (-3*(a(j-2)+a(j+2))+12*(a(j-1)+a(j+1))+17*a(j))/35;
    end
    b(n-1) = (2*(a(n)+a(n-4))+27*a(n-1)+12*a(n-2)-8*a(n-3))/35;
    b(n) = (69*a(n)+4*(a(n-1)+a(n-3))-6*a(n-2)-a(n-4))/70;
    a = b;
end

%七点三次
 n = length(x);
c = x;
for k = 1:50
    d(1) = (39*c(1) + 8*c(2) - 4*(c(3) + c(4)) + c(5) + 4*c(6) - 2*c(7))/42;
    d(2) = (8*c(1) + 19*c(2) + 16*c(3) + 6*c(4) - 4*c(5) - 7*c(6) + 4*c(7))/42