很棒的改写，对比Python版本发现了点儿新问题，可以在md更新一下

[fuerpy]

我用的是Python版本filtfilt是scipy.signal.filtfilt函数，然后我发现输出和你的对应不上，然后我就从网上搜索看到了Matlab和Python在滤波过程中的差异，会在函数计算中差0.3个相位差，所以在Python函数中 output = signal.filtfilt(b,a,x,padlen=3*(max(len(b),len(a))-1)
就可以得到一样的输出了。Matlab中不支持更改padlen（参考http://dong.sh/posts/differencefiltfiltfunction/）