c9ex5.m

来自「数字信号处理Matlab实现?中藕糯鞰atlab实现数字信号处理Matlab」· M 代码 · 共 9 行

num = 0.5*[1 4];                % Define numerator
den = conv([1 1],[1 2]);        % Define denominator
Fs = 20/(2*pi);                 % Define samp. freq. in Hz
[bz,az] = impinvar(num,den,Fs); % Determine imp. invar. filter
[r,p,k] = residuez(bz,az);      % Place in parallel form
r                               % Display residues
p                               % Display poles
k                               % Display k vector