shock滤波器的m函数,可是直接使用的shock.m