stepres1.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

% step response of an ideal lowpass filter

p=0.001;
t0=-2;
t1=6;
t=t0:p:t1;
n=1:((t1-t0)/p)+1;
f=5*sin(5*pi*(t-2)+eps)./(5*pi*(t-2)+eps);
r(1)=f(1);
for i=2:(8/p)+1
    r(i)=r(i-1)+p*f(i);
end
subplot(2,1,1)
plot(t,f);grid on;title('The impulse response h(t)')
subplot(2,1,2)
plot(t,r);grid on;title('The step response r(t)')