%
% Example 4.26 and Figure 4.10
%
num=[1.151 0.1774]
den=[1  0.739 0.921 0]
[cnum,cden]=feedback(num,den,1,1,-1)
nf=3
df=[1 0]
cnumf=conv(cnum,nf)
cdenf=conv(cden,df)
t=0:0.01:30;
y=impulse(cnumf,cdenf,t);
plot(t,y)
xlabel('time in [sec]')
ylabel('pitch angle zero-state feedback response')
grid
print -deps fig4_11.eps