Code covered by the BSD License

# Feedback Control of Dynamic Systems, 4th Ed., G. F. Franklin, J. D. Powell, A. Emami-Naeini

### Abbas Emami-Naeini (view profile)

25 Apr 2003 (Updated )

.m and .mdl files for Feedback Control of Dynamic Systems

fig9_36.m
%  Figure 9.36      Feedback Control of Dynamic Systems, 4e
%                        Franklin, Powell, Emami
%
%  Fig. 9.36  Script for computation of the yaw damper Bode plot
%  for feedback through a washout circuit.

% the equations of the aircraft:

f =[-0.0558   -0.9968    0.0802    0.0415;
0.5980   -0.1150   -0.0318         0 ;
-3.0500    0.3880   -0.4650         0 ;
0    0.0805    1.0000         0] ;
g =[0.0073;
-0.4750;
0.1530;
0];
h = [0     1     0     0];
j =[0];
% the equations of the actuator:

na=[0 10];
da=[1 10];
[fa,ga,ha,ja]=tf2ss(na,da);

% the equations of the aircraft with actuator:

[fp,gp,hp,jp]=series(fa,ga,ha,ja,f,g,h,j);
% the washout circuit
nw=[1 0];dw=[1 1/3];
[fw,gw,hw,jw]=tf2ss(nw,dw);
%  Open loop equations with washout
sysw=ss(fw,gw,hw,jw);
sysp=ss(fp,gp,hp,jp);
syspw=series(sysp,sysw)
[fpw,gpw,hpw,jpw]=ssdata(syspw);
% the washout compensated rootlocus
hold off; clf
w=logspace(-1,1);
syspw1=ss(fpw,gpw,-hpw,-jpw);
[mag, ph]=bode(syspw1,w);
subplot(211); loglog(w,mag(:,:)); grid; hold on;