No BSD License  

Highlights from
Signals, Systems, and Transforms

image thumbnail
impz(b,a);
function [h,n] = impz(b,a);
	n = 0:4/(1-max(abs(roots(a))));  
	delta = zeros(size(n)); delta(1) = 1;       % unit impulse
	h = filter(b,a,delta);

Contact us at files@mathworks.com