cutoff.m

来自「滤波器的设计以及零极点对应频率响应分析」· M 代码 · 共 20 行

function is=cutoff(s,td)
% s is the signal
% td is the treshold when the signal smaller than it, the signal must be interpolat;
sp=find(abs(s)<td);
l=length(sp);
l1=length(s);
j=1;
l1-l
for i=1:l1
    
  if abs(s(i))>=td
   is(j)=s(i);
   j=j+1;
  end
  

end
%wn=rectwin(3);
%isp=conv(is,wn);