饱和函数需要自己构造，程序如下： function M=sat(x) d=abs(x/k)

; if d