滤除阶跃信号中的毛刺（matlab代码）

【滤除阶跃信号中的毛刺（matlab代码）】
function aph = smoothing(aphEW2,NS,NT)
%函数作用为滤除阶跃信号中的毛刺
%NS为矩阵的行数
%NT为矩阵的列数
%aphEW2为待滤波矩阵
%aph为滤波后矩阵
aphEW2=floor(aphEW2);
bb(1:NS,1:NT)=0;
bb_1(1:NS,1:NT)=1;
aa=1;
m_m=0;
for i=1:1:(NT)
aa=1;
m_m=1;
for mm=1:1:(NS-11)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(mm
else
if((aphEW2(mm,i)==aphEW2(mm+1,i))&&(aphEW2(mm,i)==aphEW2(mm+2,i))&&(aphEW2(mm,i)==aphEW2(mm+3,i))&&(aphEW2(mm,i)==aphEW2(mm+4,i))&&(aphEW2(mm,i)==aphEW2(mm+5,i))&&(aphEW2(mm,i)==aphEW2(mm+6,i))&&(aphEW2(mm,i)==aphEW2(mm+7,i))&&(aphEW2(mm,i)aphEW2(mm+8,i))&&(aphEW2(mm,i)aphEW2(mm+9,i))&&(aphEW2(mm,i)aphEW2(mm+10,i