形态学滤波matlab程序

function [last1last2]=morphological2（f1ML）
% %%%%%%%%%%原始数据f开运算，先腐蚀，再膨胀%%%%%%%%
% M=20;
g=zeros（1M）;
ero_z=[];N=length（f1）;
for i=0:N-M
    for j=0:M-1
        ero_z（j+1）=f1（i+j+1）-g（j+1）;
    end
    ero（i+1）=min（ero_z）;
end
x=ero;exp_z=[];N=length（x）;
for i=M-1:N-1
    for j=0:M-1
        exp_z（j+1）=x（i-j+1）+g（j+1）;
    end
    exp（i-M+2）=max（exp_z）;
end
%%%%%%%%%%在开运算的基础上进行闭运算，先膨胀，再腐蚀%%%%%%%%

x=0;
K=[];
% M=10;
K=zeros（1L）;
x=exp;exp_z=[];N=length（x）;
for i=L-1:N-1
    for j=0:L-1
        exp_z（j+1）=x（i-j+1）+K（j+1）;
    end
    exp1（i-L+2）=max（exp_z）;
end
x=0;
x=exp1;ero_z=[];N=length（x）;
for i=0:N-L
    for j=0:L-1
        ero_z（j+1）=x（i+j+1）-K（j+1）;
    end
    ero1（i+1）=min（ero_z）;
end
%%%%%%%%%%原始数据f闭运算，先膨胀，再腐蚀%%