《DSP using MATLAB》Problem 4.18

代码：

%% ------------------------------------------------------------------------
%%            Output Info about this m-file
fprintf('\n***********************************************************\n');
fprintf('        <DSP using MATLAB> Problem 4.18 \n\n');
 
banner();
%% ------------------------------------------------------------------------
 
% ---------------------------------------
%       1   H1(z)
% ---------------------------------------
 
b = 1/4*[1, 2, 0, 1]; a = [1];               %  
 
[R, p, C] = residuez(b,a)
 
Mp = (abs(p))'
Ap = (angle(p))'/pi
 
%% ----------------------------------------------
%%   START a    determine H(z) and sketch
%% ----------------------------------------------
figure('NumberTitle', 'off', 'Name', 'P4.18.1 H(z) its pole-zero plot')
set(gcf,'Color','white');
zplane(b,a);
title('pole-zero plot'); grid on;
 
%% ----------------------------------------------
%%    END
%% ----------------------------------------------
 
% ------------------------------------
%               y1(n)  Output
% ------------------------------------
yb1 = conv(b, [2]);
ya1 = conv(a, [1,-0.9]);
 
[R1, p1, C1] = residuez(yb1, ya1)
 
Mp = (abs(p1))'
Ap = (angle(p1))'/pi
 
[delta, n] = impseq(0, 0, 7);
y1_check = filter(yb1, ya1, delta)                                             % check sequence
 
y1_answer = -1.7970*impseq(0,0,7) - 0.6173*impseq(1,0,7) - 0.5556*impseq(2,0,7) + 2.2970*(0.9).^n .* stepseq(0,0,7)         % answer sequence

　　运行结果：

只做第1小题了，后面的多少有些出入之处，以后有机会补上。