zir.m

很多matlab的源代码

function nzi = zir(d,ic,ty)
% ZIR	Zero Input Response 
% Support file for SYSRESP1 AND SYSRESP2


% ADSP Toolbox: Version 2.0 
% For use with "Analog and Digital Signal Processing", 2nd Ed.
% Published by PWS Publishing Co.
%
% Ashok Ambardar, EE Dept. MTU, Houghton, MI 49931, USA
% http://www.ee.mtu/faculty/akambard.html
% e-mail: [e-mail]
% Copyright (c) 1998


if nargin<3,ty='s';end,
ic=ic(:).';n=length(ic);
l=length(d);m=l-1;ic=[ic zeros(1,m-n)];
if ty=='s',d=d(m:-1:1);else,d(1)=[];d=-d;end
h=zeros(m);
for j=1:m,h(j,m-j+1:m)=ic(1:j)*d(j);end
nzi=sum(h);
if ty=='z',nzi=nzi(m:-1:1);end
if ty=='s',nzi=[0 nzi];else,nzi=[nzi 0];end