自适应平滑滤波matlab

function Z = adaptsmooth_filter（Xmn）
% 函数对输入的二维图像矩阵进行自适应平滑滤波
% input：
%        X：输入的二维图像矩阵
%        m：m行的滤波模板
%        n：n列的滤波末班
% output: 
%        Z：输出对m x n的二维图像矩阵的运算结果
if nargin < 2
    m = 3;  % 滤波模板尺寸
    n = 3;
end
if size（X3）~=1
    error（‘图像应该为2维矩阵‘）
end
if ~isa（X‘double‘）
    X = double（X）/255;  % 数据类型 
end
m = mn（1）;  n = mn（2）;  % 模板大小
[n1n2] =size（X）;
X（n1:n1+m n2:n2+n） = 0;   % 扩充，使得图像能够被m x n模板全部覆盖
% 计算掩膜系数
for i =1:size（X1）-m
    for j=1