MATLAB中阶跃信号，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)))
bb(aa,i)=mm;
aa=aa+1;
for kk=mm:1:NS
if(aphEW2(mm,i)aphEW2(kk,i))
if(NS-7kk)
m_m=kk;
end
else
m_m=kk;
break;
end
end
end
end
%%%%%%%%%%%%%%%%%
end
end
a_a=0;
for i=1:1:(NT)
for mm=1:1:NS
for kk=1:1:NS
if(bb(kk,i)~=0)
if(mm
if(mmbb(kk-1,i))
break;
else
aphEW2(mm,i)=aphEW2(bb(kk-1,i),i);
break;
end
end
end
end
end
end
aph=aphEW2;
—————————————————————————————————————————————————
【MATLAB弹出对话框并读取对话框内容】
prompt={‘输入纵向切片方位角’};%设置提示字符串
name=‘纵向切片’;%设置标题
numlines=1;%指定输入数据的行
defAns={‘90’};%设定默认
Resize=‘on’;%设定对话框尺寸可调节
answer=inputdlg(prompt,name,numlines,defAns,‘on’);%创建输入对
M=cell2mat(answer);
[m,n]=size(M);
m1=0;
poi=n;
for k=1:1:n
if M(k)’.’
poi=k;
end
end
if poin
for k=1:1:n
m1=(M(k)-48.0)*10.0^(n-k)+m1;
end
else
for k=1:1:n
if M(k)~=’.’
m1=(M(k)-48.0)10.0^(n-k)+m1;
else
end
end
end
m1=m110^(-(n-poi));
————————————————————————————————————————————————