% PDSi TPclass18MAR2021
%

b=[1 2 1]; a=[1 -0.5 -0.5];
h=impz(b,a, 30);
stem([0:29], h)
pause

myh=8/3+1/3*(-0.5).^([0:29]);
myh(1)=myh(1)-2;
hold on
stem([0:29], myh, 'g')
hold off
pause

stem(h-myh.')