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閉運算，先膨脹，再腐蝕%%