Advanced Automation

[lecture #1] 2010.9.2 review of classical control theory (given by Prof. Kimura)

[lecture #2] 2010.9.9 CACSD introduction

  • introduction of Matlab and Simulink

[lecture #3] 2010.9.16 CACSD introduction

  • review of clasical control theory
    1. transfer function
    2. Bode diagram
    3. characteristics of 2nd-order system
    4. Nyquist stability criterion
    5. gain margin and phase margin
  • review of modern control theory
    1. state-space representation

[lecture #4] 2010.9.30 Intro. to Robust Control

  • H infinity norm
    1. robust stabilization
    2. performance improvement
s = tf('s')
G = 1/(s+1)
norm(G, 'inf')
P = 1/(s^2 + 0.2*s + 1)
K = 1
nyquist(P*K)
help impulse
impulse(P)
T = 1/(1 + P*K)
impulse(T)
sqrt(2)/(2*pi)
1/ans
step(T)
ex21
ex22

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

  • H infinity norm
    1. robust stabilization <--
    2. performance improvement

[lecture #6] 2010.10.14 norm, vector space, normed linear space

  • H infinity norm ... {robust stabilization, performance improvement}
  • norm
    • size of {number, signal, system}
    • defined on vector space (linear space)
      • ---> normed linear space
      • ... optimization

[lecture #7] 2010.10.21 eigenvalue, eigenvector, singular value decomposition

  • H infinity norm : scalar (SISO) ---> matrix (MIMO)
    • absolute value ---> maximum singular value
    • background: mixed sensitivity problem {robust stabilization and performance improvement are simultaneously considered}
A = [2, 1; 0, 1]
[X, L] = eig(A)
X*L/X - A
[U, S, V] = svd(A)
U*U'
V*V'
U*S*V'
U*S*V' - A
sqrt(eig(A'*A))
S

[lecture #8] 2010.10.28 state space representation of connected system, state space representation of generalized plant for various control problem, mixed sensitivity problem

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

see RubustControlOfPneumatic-e.pdf in Prof. Kimura's homepage for detail

[lecture #10] 2010.11.11 robust control design example (cont.)

  1. Design H infinity controller with Eqs.(41)-(47) and generalized plant depicted in Fig.3 in the pdf file.
    1. Confirm the following m-file for design:
    2. Derive the generalized plant by hand and correct the m-file.
    3. Run the m-file to find controller.
  2. Simulation
simu_pneum
plot(t, y, 'r', t, r, 'b');

simu_pneum_noise
plot(t, y, 'r', t, n, 'b')

[lecture #11] 2010.11.18 robust control design example (cont.) (given by Prof. Kimura)

Exercise

  1. step response
  2. noise response (noise type: step, sinusoidal wave 1Hz, 50Hz, 500Hz)
  3. initial response (x(0) = [1;0;0])

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

  • simulation (cont. from Nov.18)
  • Problem setup
  • Modelling (frequency response experiment)

[lecture #13] 2010.12.2 Speed control of two inertia system with servo motor (2/3)

  • report
  1. design your controller 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:
    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(Tue) 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 and cont_order.dat to kobayasi@nagaokaut.ac.jp not later than 24th Dec
  • program sources for frequency response experiment
    • format of frdata.dat file
      • 1st column: frequency (Hz)
      • 2nd column: gain
      • 3rd column: phase (deg)
  • program sources for control experiment
    • 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(1) / (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.

[lecture #14] 2010.12.9 Speed control of two inertia system with servo motor (3/3)

Due to a dimension change of the driven shaft, frequency response experiment has been re-conducted. Please use the following fixed files instead of the ones introduced on the previous lecture.

[lecture #15] 2010.12.16 Speed control of two inertia system with servo motor (cont.)

  • preparation of your own controller(s)

participant list2010

related links


添付ファイル: file2010.12.9-1.jpg 389件 [詳細] fileresult_small2_fixed.dat 489件 [詳細] fileresult_small1_fixed.dat 573件 [詳細] fileresult_large2_fixed.dat 539件 [詳細] fileresult_large1_fixed.dat 604件 [詳細] filefrdata_offset10_small_fixed.dat 451件 [詳細] filefrdata_offset10_large_fixed.dat 440件 [詳細] filefrdata_offset5_small_fixed.dat 438件 [詳細] filefrdata_offset5_large_fixed.dat 423件 [詳細] filecont_order2_fixed.dat 434件 [詳細] filecont_order1_fixed.dat 430件 [詳細] filecont2_fixed.mat 365件 [詳細] filecont2_fixed.dat 482件 [詳細] filecont1_fixed.mat 363件 [詳細] filecont1_fixed.dat 424件 [詳細] filecont_fixed.m 458件 [詳細] fileweight_fixed.m 436件 [詳細] filenominal_fixed.m 452件 [詳細] filefreqresp_fixed.m 433件 [詳細] file2010.12.2-2.jpg 362件 [詳細] file2010.12.2-1.jpg 390件 [詳細] filefreqresp_module.c 470件 [詳細] filefreqresp_app.c 535件 [詳細] filefreqresp.h 452件 [詳細] filehinf_module.c 438件 [詳細] filehinf_app.c 459件 [詳細] filehinf.h 414件 [詳細] filecompare.m 450件 [詳細] filecont2.mat 441件 [詳細] filecont1.mat 390件 [詳細] fileresult_small2.dat 480件 [詳細] fileresult_small1.dat 448件 [詳細] fileresult_large2.dat 514件 [詳細] fileresult_large1.dat 548件 [詳細] filecont_order2.dat 384件 [詳細] filecont_order1.dat 446件 [詳細] filecont2.dat 398件 [詳細] filecont1.dat 446件 [詳細] filecont.m 534件 [詳細] fileweight.m 497件 [詳細] filenominal.m 442件 [詳細] filefreqresp.m 428件 [詳細] filefrdata_offset10_small.dat 452件 [詳細] filefrdata_offset10_large.dat 448件 [詳細] filefrdata_offset5_small.dat 453件 [詳細] filefrdata_offset5_large.dat 437件 [詳細] filefrdata_offset0_small.dat 457件 [詳細] filefrdata_offset0_large.dat 450件 [詳細] filesetup.pdf 444件 [詳細] file2010.11.18-1.jpg 381件 [詳細] file2010.11.25-1.jpg 385件 [詳細] fileapparatus.pptx 449件 [詳細] file2010.11.11-2.jpg 418件 [詳細] file2010.11.11-1.jpg 426件 [詳細] filesimu_pneum_noise.mdl 505件 [詳細] filesimu_pneum.mdl 500件 [詳細] filepneum_ans.m 441件 [詳細] filepneum.m 464件 [詳細] file2010.11.4-2.jpg 364件 [詳細] file2010.11.4-1.jpg 393件 [詳細] file2010.10.28-2.jpg 407件 [詳細] file2010.10.28-1.jpg 444件 [詳細] fileex26.m 463件 [詳細] file2010.10.21-5.jpg 422件 [詳細] file2010.10.21-4.jpg 427件 [詳細] file2010.10.21-3.jpg 412件 [詳細] file2010.10.21-2.jpg 418件 [詳細] file2010.10.21-1.jpg 431件 [詳細] file2010.10.14-4.jpg 441件 [詳細] file2010.10.14-3.jpg 416件 [詳細] file2010.10.14-2.jpg 427件 [詳細] file2010.10.14-1.jpg 435件 [詳細] filenorm2.pdf 512件 [詳細] filenorm.pdf 530件 [詳細] file2010.10.7-4.jpg 436件 [詳細] file2010.10.7-3.jpg 418件 [詳細] file2010.10.7-2.jpg 431件 [詳細] file2010.10.7-1.jpg 437件 [詳細] filemod2.mdl 566件 [詳細] fileex25.m 640件 [詳細] fileex24.m 593件 [詳細] fileex23.m 541件 [詳細] file2010.9.30-3.jpg 435件 [詳細] file2010.9.30-2.jpg 435件 [詳細] file2010.9.30-1.jpg 417件 [詳細] fileex22.m 564件 [詳細] fileex21.m 552件 [詳細] file2010.9.16-2.jpg 421件 [詳細] file2010.9.16-1.jpg 424件 [詳細] filemod1.mdl 585件 [詳細] fileex4.m 537件 [詳細] fileex3.m 551件 [詳細] fileex2.m 556件 [詳細] fileex1.m 576件 [詳細] file2010.9.9-2.jpg 484件 [詳細] file2010.9.9-1.jpg 463件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2010-12-16 (木) 13:21:01 (3487d)