User`s guide
Table Of Contents
- Preface
- Quick Start
- LTI Models
- Introduction
- Creating LTI Models
- LTI Properties
- Model Conversion
- Time Delays
- Simulink Block for LTI Systems
- References
- Operations on LTI Models
- Arrays of LTI Models
- Model Analysis Tools
- The LTI Viewer
- Introduction
- Getting Started Using the LTI Viewer: An Example
- The LTI Viewer Menus
- The Right-Click Menus
- The LTI Viewer Tools Menu
- Simulink LTI Viewer
- Control Design Tools
- The Root Locus Design GUI
- Introduction
- A Servomechanism Example
- Controller Design Using the Root Locus Design GUI
- Additional Root Locus Design GUI Features
- References
- Design Case Studies
- Reliable Computations
- Reference
- Category Tables
- acker
- append
- augstate
- balreal
- bode
- c2d
- canon
- care
- chgunits
- connect
- covar
- ctrb
- ctrbf
- d2c
- d2d
- damp
- dare
- dcgain
- delay2z
- dlqr
- dlyap
- drmodel, drss
- dsort
- dss
- dssdata
- esort
- estim
- evalfr
- feedback
- filt
- frd
- frdata
- freqresp
- gensig
- get
- gram
- hasdelay
- impulse
- initial
- inv
- isct, isdt
- isempty
- isproper
- issiso
- kalman
- kalmd
- lft
- lqgreg
- lqr
- lqrd
- lqry
- lsim
- ltiview
- lyap
- margin
- minreal
- modred
- ndims
- ngrid
- nichols
- norm
- nyquist
- obsv
- obsvf
- ord2
- pade
- parallel
- place
- pole
- pzmap
- reg
- reshape
- rlocfind
- rlocus
- rltool
- rmodel, rss
- series
- set
- sgrid
- sigma
- size
- sminreal
- ss
- ss2ss
- ssbal
- ssdata
- stack
- step
- tf
- tfdata
- totaldelay
- zero
- zgrid
- zpk
- zpkdata
- Index

9 Design Case Studies
9-42
LQG Design for the y-Axis
TheLQG design for the -axis (regulation of the thickness)followsthe exact
samestepsasforthe -axis.
% Specify model components
Hy = tf(7.8e8,[1 71 88^2],'inputn','u-y')
Fiy = tf(2e4,[1 0.05],'inputn','w-iy')
Fey = tf([1e5 0],[1 0.19 9.4^2],'inputn','w-ey')
gy = 0.5e–6 % force-to-gap gain
% Build open-loop model
Py = append([ss(Hy) Fey],Fiy)
Py = [–gy gy;1 1] * Py
set(Py,'outputn',{'y-gap' 'y-force'})
% State-feedback gain design
Pydes = append(lpf,1) * Py % Add low-freq. weigthing
set(Pydes,'outputn',{'y-gap*' 'y-force'})
ky = lqry(Pydes(1,1),1,1e–4)
% Kalman estimator design
esty = kalman(Pydes(2,:),eye(2),1e3)
% Form SISO LQG regulator for y-axis and close the loop
Regy = lqgreg(esty,ky)
cly = feedback(Py,Regy,1,2,+1)
y
y
x