robust control - scl.ipc.i.u-tokyo.ac.jp€¦ · robust control spring, 2019 instructor: prof....
TRANSCRIPT
![Page 1: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/1.jpg)
Robust Control
Spring, 2019Instructor: Prof. Masayuki Fujita (S5-303B)
8th classTue., 4th June, 2019, 10:45~11:30,
S423 Lecture Room1
![Page 2: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/2.jpg)
8.1 HiMAT: Control(Highly Maneuverable Aircraft Technology)
8. Design Example
Reference:[SP05] S. Skogestad and I. Postlethwaite,
Multivariable Feedback Control; Analysis and Design,Second Edition, Wiley, 2005.
[SLH81] M. G. Safonov, A. J. Laub and G. L. Hartmann,Feedback Properties of Multivariable Systems: The Role and Use of the Return Difference Matrix,IEEE Transactions on Automatic Control, Vol. 26, No. 1, pp. 47-65, 1981.
[ES84] M. B. Evans and L. J. Schilling,The Role of Simulation in the Development and Flight Test of the HiMAT Vehicle,NASA Technical Memorandum 84912, 1984.
[CS98] R. Y. Chiang and M. G. Safonov,Robust Control Toolbox User’s Guide,The MathWorks, 1998. 2
![Page 3: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/3.jpg)
HiMAT: Attitude and flight path control6.86m
4.75m
Output
ObjectiveControl attitude ( ) and flight path ( ) separately
MIMO System
Input〜1 2
[Ex.]
Control attitude without changing flight pathSpecification
Performance Spec.: Minimize the sens. function as much as possible
Robustness Spec.: dB/dec roll-off and dB at 100rad/s(Because of the unmodeled low-damped structual mode)
3
Max takeoff weight: 1,542kg
![Page 4: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/4.jpg)
Modeling HiMAT
Frequency ResponseInteraction(Coupling)
Unstable
Poles
Zeros4
〜3 5
![Page 5: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/5.jpg)
(Roll-off)
MIMO Loop Shaping [SP05, p. 343]
Time Delay: 0.025s40 rad/s
Performance Weight1.43 rad/s
Uncertainty Weight30 rad/s
Actuator dynamics
Robustness
Performance
〜6 9
5
Unstable Pole:
Disturbance (Gust):
![Page 6: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/6.jpg)
6
HiMAT: Uncertainty WeightwM = tf([1 0 0], [1000]);WM = eye(2)*wM;
MATLAB Command
Uncertain Factors:Actuator Dynamics
30 rad/s0.033,
Roll-off: dB at 100rad/sdB/dec
Frequency:
Target Loop
2 Oscillation mode (2nd order system)
+Unmodeled dynamics (High freq.)
Natural angular frequencyDamping factor
![Page 7: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/7.jpg)
HiMAT: Performance Weight
,1.5
Ms = 2; A = 0.01; wb = 1.5;wP = tf([1/Ms wb], [1 wb*A]);WP = eye(2)*wP;
MATLAB Command
Update Larger
(Trade-off)
+
−
7
〜10 14
![Page 8: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/8.jpg)
−
%Generalized Plant%systemnames = 'Pnom WP WM';inputvar = '[w(2);u(2)]';outputvar = '[WP;WM;-w-Pnom]';input_to_Pnom= '[u]';input_to_WP = '[w+Pnom]';input_to_WM = ' [Pnom]';G = sysic;
MATLAB Command
HiMAT: Control Problem Formulation
−
Nominal Performance
−
Robust Stability
sysic
−
Generalized Plant
Mixed Sensitivity 8
![Page 9: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/9.jpg)
HiMAT: -iteration to obtain Controller
[SV,w]=sigma(CLhi);figure; semilogx(w,SV)
MATLAB Command
Find such that
[SP05, p. 358]
Appropriately sub-optimalCheck 1)
(Default settings)
Closed-loop Transfer Function
9
There is No Controller
No SolutionCheck 2) Data Structure
Resetting value of Gamma min based on D_11, D_12, D_21 terms
Test bounds: 0.1667 < gamma <= 0.9538
gamma hamx_eig xinf_eig hamy_eig yinf_eig nrho_xy p/f0.954 2.1e-02 1.8e-05 1.6e-01 1.0e-22 0.9915 p 0.560 2.1e-02 -1.1e+06# 1.6e-01 -6.4e-20 0.2088 f0.757 2.1e-02 -8.2e+06# 1.6e-01 -1.4e-20 0.6112 f0.855 2.1e-02 -5.0e+07# 1.6e-01 -4.1e-20 2.9142# f 0.905 2.1e-02 1.9e-05 1.6e-01 -4.4e-23 3.1083# f0.929 2.1e-02 1.8e-05 1.6e-01 -4.7e-20 1.5096# f0.941 2.1e-02 1.8e-05 1.6e-01 -2.5e-20 1.1979# f0.948 2.1e-02 1.8e-05 1.6e-01 -1.1e-19 1.0852# f
Gamma value achieved: 0.9538
[Khi,CLhi,ghi,hiinfo] = …hinfsyn(G,nmeas,ncon,’Display’, ‘on’);
MATLAB Command
![Page 10: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/10.jpg)
HiMAT: Controllerfiguresigma(Khi)
MATLAB Command
Zeros
8OrderModel Reduction
Poles
〜15 20
[Ex.]
10
![Page 11: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/11.jpg)
HiMAT: Open-loop Frequency Response
Loop Transfer Function
figuresigma(Fhi.Lo,WP, inv(WM),WP/ghi,ghi*inv(WM))
MATLAB Command
[SV,w]=sigma(WP*Fhi.So);figure; semilogx(w,SV)[SV,w]=sigma(WM*Fhi.To);figure; semilogx(w,SV)
MATLAB Command
Robust Stability (RS)
Nominal Performance (NP)NP/RS Test
corresponding maximum stability margin
12
![Page 12: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/12.jpg)
HiMAT: Closed-loop PerformanceNominal Stability (NS)
Poles of
Nominal Performance (NP) Robust Stability (RS)
Zeros of
pole(CLhi)zero(CLhi)figure; pzmap(CLhi)
MATLAB Command
11
![Page 13: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/13.jpg)
HiMAT: Step Response (closed loop)
−
ReferenceNominal ModelPerturbed Model
13
![Page 14: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/14.jpg)
8. Design Example
Reference:[SP05] S. Skogestad and I. Postlethwaite,
Multivariable Feedback Control; Analysis and Design,Second Edition, Wiley, 2005.
[SLH81] M. G. Safonov, A. J. Laub and G. L. Hartmann,Feedback Properties of Multivariable Systems: The Role and Use of the Return Difference Matrix,IEEE Transactions on Automatic Control, Vol. 26, No. 1, pp. 47-65, 1981.
[ES84] M. B. Evans and L. J. Schilling,The Role of Simulation in the Development and Flight Test of the HiMAT Vehicle,NASA Technical Memorandum 84912, 1984.
[CS98] R. Y. Chiang and M. G. Safonov,Robust Control Toolbox User’s Guide,The MathWorks, 1998.
8.1 HiMAT: Control(Highly Maneuverable Aircraft Technology)
14
![Page 15: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/15.jpg)
Class evaluation by students
Class Name : Robust ControlCode of Class : SCE.C402
Instructor : Prof. Masayuki Fujita
How to fill out questionnaire
・Use Pencil only (Do not use Pen)
・Mark Carefully
15
![Page 16: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/16.jpg)
Crew: NoneLength: 6.86mWingspan: 4.75mHeight: 1.31mMax takeoff weight: 1,542kg
General characteristics
PerformanceMaximum speed: Mach 1.6Maneuvering: 8G performance
1
16
NASA HiMAT
Real Product
![Page 17: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/17.jpg)
NASA HiMAT
17Classic Aircraft
Real ProductTesting Product
Testing high technology
UAV(Unmanned Aerial Vehicle)
HiMAT
Teleoperation
2
20
![Page 18: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/18.jpg)
18
Modeling HiMAT
: elevon actuator: canard actuator
: attitude angle: angle of attack
states to represent actuator dynamics
6.86m
4.75m
Control InputOutput 3
![Page 19: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/19.jpg)
19
HIMAT: Nominal Plant Model
% NASA HiMAT model G(s) ag =[ -2.2567e-02 -3.6617e+01 -1.8897e+01 -3.2090e+01 3.2509e+00 -7.6257e-01;
9.2572e-05 -1.8997e+00 9.8312e-01 -7.2562e-04 -1.7080e-01 -4.9652e-03; 1.2338e-02 1.1720e+01 -2.6316e+00 8.7582e-04 -3.1604e+01 2.2396e+01;
0 0 1.0000e+00 0 0 0; 0 0 0 0 -3.0000e+01 0; 0 0 0 0 0 -3.0000e+01];
bg = [ 0 0; 0 0; 0 0; 0 0; 30 0; 0 30]; cg = [ 0 1 0 0 0 0; 0 0 0 1 0 0]; dg = [ 0 0; 0 0]; G=ss(ag,bg,cg,dg);
MATLAB Command
State Space Form (Matrix Representation) 4
![Page 20: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/20.jpg)
20
-5.6757, -0.2578, -30, -30,0.6898 0.2488i
-0.0210Unstable
HiMAT: Nominal Plant ModelControllabilityObservabilityPoles (Stability)
Zeros
Frequency ResponseStep Response for Nominal Plant Model
0.5
5
![Page 21: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/21.jpg)
Actuators: Hydraulic serve system
Block diagrams of elevon
Block diagrams of canard
20°
-20°
Ratelimit [deg]
Bode diagram of actuators
30rad/s30°
-20°Saturation
Ratelimit [deg]
Saturation21
6
![Page 22: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/22.jpg)
Actuators: Hydraulic serve systemBlock diagrams of elevonSymmetric input
Asymmetric input
Saturation
0.5
Scaling
+
―
+
+
30°
-20°
30°
-20°
Saturation
Saturation
Ratelimit
Ratelimit[deg]
[deg] [deg]
[deg]
Block diagrams of canard
20°
-20°
+0
Symmetric canard pulse
Ratelimit
+ [deg]
Bode diagram of actuators
30rad/s
7
24
![Page 23: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/23.jpg)
23
Loop Time: Control Input time
Processing Time 25 ms (40Hz)
Onboard computer
Uplink53.3times/s
Downlink110kbits/s
8
![Page 24: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/24.jpg)
Models: Understand Performance Specifications
Process[Nm] [rad]
STEP 5STEP 6
STEP 7
Process
Controller
A/DD/A
Actuators Sensors
ControllerDecision Information
: Process Noise
24
9
![Page 25: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/25.jpg)
HiMAT: Closed-loop Performance (without update )Nominal Stability (NS)
Poles of
Nominal Performance (NP) Robust Stability (RS)
〇
Zeros of 〇⇒ Pole/Zero Cancellations Nothing
25
1.5rad/s
10
![Page 26: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/26.jpg)
HiMAT: Controller (without update )figuresigma(Khi)
MATLAB Command
Poles
Zeros
8Order
Difficult to implementNumerical problems or inaccuracies may be caused too high order
26
11
![Page 27: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/27.jpg)
27
HiMAT: Open-loop Frequency Response (without update )
Loop Transfer Function
figuresigma(Fhi.Lo,WP, inv(WM),WP/ghi,ghi*inv(WM))
MATLAB Command
[SV,w]=sigma(WP*Fhi.So);figure; semilogx(w,SV)[SV,w]=sigma(WM*Fhi.To);figure; semilogx(w,SV)
MATLAB Command
Robust Stability (RS)
Nominal Performance (NP)NP/RS Test
corresponding maximum stability margin
12
![Page 28: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/28.jpg)
HiMAT: -iteration to obtain Controller (without update )
[SV,w]=sigma(CLhi);figure; semilogx(w,SV)
MATLAB Command
Find such that [SP05, p. 358]
No Solution
Appropriately sub-optimal
[Khi,CLhi,ghi,hiinfo] = …hinfsyn(G,nmeas,ncon,’Gmax’,100,’Gmin’,100);
MATLAB Command
Check 1)
Check 2)
(Default settings)
Test) Tune and sub-optimality
Closed-loop Transfer Function
“interested” frequency range
28
There is No Controller
Data Structure
13
![Page 29: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/29.jpg)
29
HiMAT: Step Response of closed loop (without update )
0.514
![Page 30: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/30.jpg)
30
Approximately a bandwidth of 10 rad/sDesired Loop Shape
Target loop
HIMAT: Specifications and Open-loop (Loop Shaping Synthesis)
Gd= tf( 10, [1 .001] ) ;[K,CL,GAM,INFO]=loopsyn(G,Gd); sigma(Gd,‘b',G*K,'r',Gd/GAM,‘g:',Gd*GAM,‘g:',{.1,100})
MATLAB Command
15
![Page 31: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/31.jpg)
31
HIMAT: Controller Model Reduction (Loop Shaping Synthesis)
size(K)
hsv = hankelsv(K); semilogy(hsv,'*--'), grid title('Hankel singular values of K'), xlabel('Order')
Kr = reduce(K,9); order(Kr)
sigma(K,'b',K-Kr,'r-.') legend('K','error K-Kr')
MATLAB Command
Order: 16 → 9
Step responses
16
![Page 32: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/32.jpg)
32
HIMAT: Sensitivity (Loop Shaping Synthesis)
T = feedback(G*K,eye(2)); S = eye(2)-T;sigma(inv(S),'m',T,'g',L,'r--',Gd,'b',Gd/GAM,'b:',Gd*GAM,'b:',{.1,100})
MATLAB Command
-3
17
![Page 33: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/33.jpg)
33
HIMAT: Designed Weight (Loop Shaping Synthesis)
sigma(INFO.Gs,‘b',G,'r',INFO.W,‘g',{.01,100});MATLAB Command
18
![Page 34: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/34.jpg)
34
HIMAT: Controller (Loop Shaping Synthesis)
sigma(K);MATLAB Command
Numerical problems or inaccuracies may be caused too high orderDifficult to implement
Order: 16
19
![Page 35: Robust Control - scl.ipc.i.u-tokyo.ac.jp€¦ · Robust Control Spring, 2019 Instructor: Prof. Masayuki Fujita (S5-303B) 8th class. Tue., 4th June, 2019, 10:45 ~11:30, S423 Lecture](https://reader034.vdocument.in/reader034/viewer/2022051606/6018ee61925ac432c16c0304/html5/thumbnails/35.jpg)
35
0.5
HIMAT: Step Response of closed loop (Loop Shaping Synthesis)
20