授業/制御工学特論2012 の履歴(No.8)

• 履歴一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• 授業/制御工学特論2012 へ行く。
  • 1 (2012-09-07 (金) 09:52:58)
  • 2 (2012-09-07 (金) 09:57:34)
  • 3 (2012-09-27 (木) 11:48:37)
  • 4 (2012-09-27 (木) 14:30:02)
  • 5 (2012-09-28 (金) 11:55:37)
  • 6 (2012-10-04 (木) 14:30:26)
  • 7 (2012-10-10 (水) 21:23:20)
  • 8 (2012-10-11 (木) 14:27:46)
  • 9 (2012-10-15 (月) 15:57:07)
  • 10 (2012-10-18 (木) 14:35:36)
  • 11 (2012-10-24 (水) 20:56:46)
  • 12 (2012-10-25 (木) 14:26:41)
  • 13 (2012-10-30 (火) 10:28:24)
  • 14 (2012-11-01 (木) 14:24:13)
  • 15 (2012-11-05 (月) 16:31:35)
  • 16 (2012-11-08 (木) 11:21:13)
  • 17 (2012-11-08 (木) 14:35:19)
  • 18 (2012-11-09 (金) 11:05:46)
  • 19 (2012-11-15 (木) 14:34:25)
  • 20 (2012-11-19 (月) 18:54:09)
  • 21 (2012-12-06 (木) 14:24:29)
  • 22 (2012-12-07 (金) 14:28:33)
  • 23 (2012-12-13 (木) 14:29:12)
  • 24 (2012-12-19 (水) 13:25:34)
  • 25 (2012-12-20 (木) 13:30:49)

---

Advanced Automation†

[lecture #1] 2012.9.6 review of classical control theory†

[lecture #2] 2012.9.20 given by Prof. Kimura†

[lecture #3] 2012.9.27 CACSD introduction with review of classical and modern control theory†

1. introduction of Matlab and Simulink
   • Basic usage of MATLAB and Simulink used for 情報処理演習及び考究II/Consideration and Practice of Information Processing II: Advanced Course of MATLAB
2. relationship between TF and SSR
   ex0927_1.m
3. open-loop stability can be checked by poles of TF and eigenvalues of A-matrix in SSR
   ex0927_2.m
   mod0927_2.mdl
4. closed-loop stability
   • graphical test by Nyquist stability criterion and Bode plot with GM(gain margin) and PM(phase margin)
     ex0927_3.m
     mod0927_3.mdl
   • numerical test by poles of closed-loop TF and eigenvalues of A-matrix in closed-loop SSR
     ex0927_4.m
   • example: stabilization of unstable system
     ex0927_5.m
     mod0927_5.mdl

%-- 9/27/2012 1:07 PM --%
a = 1
pwd
G1 = 1/(s+1)
ex0927_1
G1 = 1/(s+1)
A
ex0927_2
G2_tf
G2_tf.den
G2_tf.den{:}
G2_tf.num{:}
eig(G2_tf)
G2_ss.a
eig(G2_ss.a)
G2_ss
G2_ss.a
G2_ss.b
ex0927_3
nyquist(1.5*G3_tf)
nyquist(-1.5*G3_tf)
bode(-1.5*G3_tf)
bode(1.5*G3_tf)
K
K = 1
ex0927_4
ex0927_5
A
eig(A)
K = 0
K=1
ex0927_5 

[lecture #4] 2012.10.4 Intro. to robust control theory (H infinity control theory)†

1. Typical design problems
   1. robust stabilization
   2. performance optimization
   3. robust performance problem (robust stability and performance optimization are simultaneously considered)
2. H infinity norm
   • definition
   • example
3. H infinity control problem
   • definition
   • application example : reference tracking problem
     • relation to the sensitivity function S(s) (S(s) -> 0 is desired but impossible)
     • given control system
       ex1004_1.m
       please change the line `K = 1' to `K = ss(1)'

• controller design with H infinity control theory
  ex1004_2.m

%-- 10/4/2012 1:11 PM --%
s = tf('s')
G = 1/(s+1)
norm(G, 'inf')
help norm
ex1004_1
ex1004_2
K
K_hinf
size(K_hinf.a)

[lecture #5] 2012.10.11 Introduction to Robust Control (cont.)†

1. Typical design problems
   1. robust stabilization ... plant is given as class
   2. performance optimization; ... plant has no variation
   3. robust performance problem (robust stability and performance optimization are simultaneously considered)
2. connection between H infinity control problem and robust stabilization
   • normalized uncertainty \Delta
   • small gain theorem
   • sketch proof ... Nyquist stability criterion
3. How to design robust stabilizing controller with H infinity control problem ?
   • connection between \Delta and plant ?
   • generalized plant G ?
   • practical example of robust stabilization problem
     • unstable plant with perturbation
       ex1011_1.m
     • multiplicative uncertainty model
       ex1011_2.m
     • H infinity controller design
       ex1011_3.m
       mod1011.mdl

%-- 10/11/2012 12:56 PM --%
ex1011_1
ex1011_2
ex1011_3
mod1011
ex1011_1
P
ex1011_2
1i
j
j = 2
1i = 2
ex1011_3
mod1011
c = 0.8
c = 1.2
c = 1.3
c = 2

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! the followings are under construction !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

[lecture #6] 2011.10.27 norm (vector space, normed linear space), singular value, mixed sensitivity problem†

• H infinity norm ... {robust stabilization, performance improvement}
• norm

  #ref(): File not found: "norm.pdf" at page "授業/制御工学特論2012"

  • size of {number, signal, system}
  • defined on vector space (linear space)
    • ---> normed linear space
    • ... optimization
• H infinity norm : scalar (SISO) ---> matrix (MIMO)
  • absolute value ---> maximum singular value

    #ref(): File not found: "norm2.pdf" at page "授業/制御工学特論2012"

  • robust stabilization or performance optimization problem ---> mixed sensitivity problem
• reference
  • Feedback control theory / John C. Doyle, Bruce A. Francis and Allen R. Tannenbaum

    #ref(): File not found: "feedback_control_theory.jpg" at page "授業/制御工学特論2012"

    • 2. Norms for signals and systems
    • 3. Basic concepts
    • 4. Uncertainty and robustness ... p.62 Eq.(4.10) H infinity norm condition in mixed sensitivity problem
  • Robust and optimal control / Kemin Zhou, John C. Doyle and Keith Glover

    #ref(): File not found: "robust_and_optimal_control.jpg" at page "授業/制御工学特論2012"

    • 2. Linear algebra (2.7 Vector norms and matrix norms, 2.8 Singular value decomposition)
    • 4. Performance specifications (4.1 Normed spaces, 4.3 Hardy spaces H_2 and H_infinity)
    • 9. Model uncertainty and robustness
    • 16. H infinity control: simple case
    • 17. H infinity control: general case

#ref(): File not found: "2011.10.27-1.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "2011.10.27-2.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "2011.10.27-3.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "2011.10.27-4.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "2011.10.27-5.jpg" at page "授業/制御工学特論2012"

[lecture #8] 2011.11.10 state space representation of connected system, state space representation of generalized plant for various control problem, mixed sensitivity problem (by Prof. Kimura)†

see Prof. Kimura's page

#ref(): File not found: "aa-111110-1.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "aa-111110-2.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "aa-111110-3.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "aa-111110-4.jpg" at page "授業/制御工学特論2012"

[lecture #9] 2011.12.1 robust control design example: Robust Control System Synthesis for Pneumatic Systems (given by Prof. Kimura)†

[lecture #10] 2011.12.8 Speed control of two inertia system with servo motor (1/3)†

• Experimental setup

  #ref(): File not found: "apparatus.ppt" at page "授業/制御工学特論2012"

  #ref(): File not found: "setup.pdf" at page "授業/制御工学特論2012"

  • Objective for control system design
  1. speed tracking
  2. robust against inertia-load variation
• Modelling
  • frequency response experiment results:

    #ref(): File not found: "frdata_small_offset0.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "frdata_small_offset5.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "frdata_small_offset10.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "frdata_large_offset0.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "frdata_large_offset5.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "frdata_large_offset10.dat" at page "授業/制御工学特論2012"

  • example of m-file
    freqresp.m
    nominal.m
    weight.m

    >> freqresp
    >> nominal
    >> weight

• Controller design (mixed sensitivity problem)
  • example of m-file
    cont.m

    >> cont

  • design example
    cont.dat
    cont_order.dat
    cont.mat

    #ref(): File not found: "result_small.dat" at page "授業/制御工学特論2012"

    #ref(): File not found: "result_large.dat" at page "授業/制御工学特論2012"

• report

1. design your controller(s) so that the system performance is improved compared with the given example above
2. Draw the following figures and explain the difference between two control systems (your controller and the example above):
   1. bode diagram of controllers
   2. gain characteristic of closed-loop systems
   3. time response of control experiment
3. Why is the performance of your system improved(or unfortunately decreased)?
   • due date: 28th(Wed) Dec 17:00
   • submit your report(pdf or doc) by e-mail to kobayasi@nagaokaut.ac.jp
   • You can use Japanese
   • maximum controller order is 20
   • submit your cont.dat, cont_order.dat, and cont.mat to kobayasi@nagaokaut.ac.jp not later than 23th(Fri) Dec

• program sources for frequency response experiment
  freqresp.h
  freqresp_module.c
  freqresp_app.c
  • format of frdata.dat file
    • 1st column: frequency (Hz)
    • 2nd column: gain
    • 3rd column: phase (deg)
• program sources for control experiment
  hinf.h
  hinf_module.c
  hinf_app.c
  • format of result.dat file
    • 1st column: time (s)
    • 2nd column: motor speed (rad/s)
    • 3rd column: motor torque (Nm)
    • 4th column: reference speed (rad/s)
• configuration of control experiment
  • reference signal is generated as described in hinf_module.c:

    if((t > 3)&&(t < 7)){
      r = 10.0;
    }else{
      r = 5;
    }

  • two inertia-load discs (small and large) are used

• Difficulties of our plant: As motor speed is approximately calculated by using difference for one sampling period in hinf_module.c like

  thetaM_rad = (double)read_theta(0) / (double)Pn212 * 2 * M_PI;
  speedM_rad = (thetaM_rad - thetaM_rad_before) / msg->sampling_period;
  thetaM_rad_before = thetaM_rad

  sampling period should not become too small. On the other hand, sampling period should be chosen as small as possible so that desiged continuous-time controller could be closely implemented by its descretized version. Therefore, we have a dilemma to control our plant. The sampling period 0.25 msec was chosen by traial and error so that noise in measured speed is not too large. The gain in high frequency range of continuous-time controller should be small enough for discretization.

#ref(): File not found: "2011.12.08-1.jpg" at page "授業/制御工学特論2012"

%-- 12/8/2011 1:22 PM --%
freqresp
nominal
weight
load result_small.dat
result_small
plot(result_small(:,1),result_small(:,2))
plot(result_small(:,1),result_small(:,2), 'b', result_small(:,1),result_small (:,4),'r')

[lecture #11] 2011.12.15 Speed control of two inertia system with servo motor (2/3)†

• explanation of design example (cont. from the previous lecture)
• preparation of your own controller(s)

#ref(): File not found: "2011.12.15-1.jpg" at page "授業/制御工学特論2012"

#ref(): File not found: "2011.12.15-2.jpg" at page "授業/制御工学特論2012"

%-- 12/15/2011 1:06 PM --%
cont
help hinfsyn
K
Ktf = tf(K)
bode(K, 'r', Ktf, 'b')
legend('ss', 'tf')
help fitsys

[lecture #12] 2011.12.22 Speed control of two inertia system with servo motor (3/3)†

• preparation of your own controller(s)

participant list2011

supplementary lectures will be given by Prof. Kimura.†

• example for changing ws = 2 pi 0.1 -> 2 pi 0.01

  #ref(): The style ref(filename,pagename) is ambiguous and become obsolete. Please try ref(pagename/filename)

• example to compare time responces

  #ref(): File not found: "compare_result.m" at page "授業/制御工学特論2012"

%-- 12/22/2011 1:20 PM --%
nominal
ctrlpref
nominal
weight
cont
ks
cont
compare_result

related links†

• Prof. Kimura's page