第239行,对x的方向导数 dgau2D=-x.*exp(-(x.*x+y.*y)/(2*ssq))/(pi*ssq); 是复合函数导数,对指数里面求导是不是少了个2*ssq,有matlab的帮忙看一下啊~~~~