《现代控制理论》实验指导书Experiment Instructor of Modern ControlTheory张存山吴兴华齐晓妹Zhang CunShan Wu Xinghua Qi Xiaomei山东理工大学电气与电子工程学院自动化系Dept. of Automation, School of Electrical and ElectronicsEngineering, SDUT2015年4月April, 2015
《现代控制理论》实验指导书 Experiment Instructor of Modern Control Theory 张存山 吴兴华 齐晓妹 Zhang CunShan Wu Xinghua Qi Xiaomei 山东理工大学电气与电子工程学院自动化系 Dept. of Automation, School of Electrical and Electronics Engineering, SDUT 2015年4月 April，2015
目录Contents一状态空间模型的建立与状态响应的分析实验一EXPERIMENT1 Modeling with State SpaceDescription andAnalysis ofTime Response by Using MATLAB实验二 控制系统的稳定性、能控性和能观测性分析EXPERIMENT2 Analysis of Stability、Controllability and Observabi-litybyUsingMATLAB10实验三控制系统状态反馈控制律的确定EXPERIMENT3 Designing the Feedback Control Law for a 2-OrderSystem15实验四控制系统综合设计EXPERIMENT4 System Synthesis for a 3-order System by Using MATLAB
目 录 Contents 实验一 状态空间模型的建立与状态响应的分析.1 EXPERIMENT1 Modeling with State Space Description and Analysis of Time Response by Using MATLAB 实验二 控制系统的稳定性、能控性和能观测性分析.5 EXPERIMENT2 Analysis of Stability、Controllability and Observability by Using MATLAB 实验三 控制系统状态反馈控制律的确定.10 EXPERIMENT3 Designing the Feedback Control Law for a 2-Order System 实验四 控制系统综合设计.15 EXPERIMENT4 System Synthesis for a 3-order System by Using MATLAB
实验一状态空间模型的建立与状态响应分析EXPERIMENT1 Modeling with State Space Description and Analysis ofTime Response by Using MATLABI.ExperimentalPurpose1. Master the methods of obtaining the state space description2. Master the methods of state transformation.3. Master the methods of obtaining the diagonal canonical form and Jordan canonicalform by state transformation4. Master the methods of calculation of the time response of the LTI systemIl.ExperimentalContent1.1 The state space description ofthe system is「2.25-5-1.25[4 0.56242.25-4.251.25-0.25X(t) =X(t)u(t)220.25-1-0.5-1.25021.251.75-0.250.750100u(t)y(t)10oL1.1.1 Modeling with the state space descriptionThis equation can be input to the MATLAB work space by the followingstatement.>>A=[2.25.-5,-1.25.-0.5:2.25,-4.25,-1.25,-0.25:-0.25,-0.5,-1.25-1:1.25.-1.75,0.25,-0.75];>> B=[4,6;2,4;2,2;0,2];>>C=[1,0,0,0;0,1,0,0];>> D=zeros(2,2);G=ss(A,B,C,D)1
1 实验一 状态空间模型的建立与状态响应分析 EXPERIMENT1 Modeling with State Space Description and Analysis of Time Response by Using MATLAB I. Experimental Purpose 1. Master the methods of obtaining the state space description. 2. Master the methods of state transformation. 3. Master the methods of obtaining the diagonal canonical form and Jordan canonical form by state transformation 4. Master the methods of calculation of the time response of the LTI system II. Experimental Content 1.1 The state space description of the system is 2.25 5 1.25 0.5 4 6 2.25 4.25 1.25 0.25 2 4 ( ) () () 0.25 0.5 1.25 1 2 2 1.25 1.75 0.25 0.75 0 2 X t X t ut −− − −−− = + − −− − −−− ( ) 0 1 0 0 1 0 0 0 y(t) u t = 1.1.1 Modeling with the state space description This equation can be input to the MATLAB work space by the following statement. >>A=[2.25,-5,-1.25,-0.5;2.25,-4.25,-1.25,-0.25;-0.25,-0.5,-1.25,-1;1.25,-1.75,- 0.25,-0.75]; >> B=[4,6;2,4;2,2;0,2]; >> C=[1,0,0,0;0,1,0,0]; >> D=zeros(2,2);G=ss(A,B,C,D)
1.1.2 Conversion of state space model established by 1.1.1 as the transferfunction modelThe model of transferfunction can be obtained by thefollowing statements:>>G1=tf(G)or>>[Numl,denl]=ss2tf(A,B,C,D,1)>>[Num2,den2]=ss2tf(A,B,C,D,2)1.2 The model of the system is$+5G(s) = ,$*+2s3 +3s? +4s+51.2.1 Set up the model of the transfer function of the systemTransfer function can be input to the MATLAB work space by the followingstatements.>> num=[1,5];den-[1,2,3,4,5],G=tf(num,den1.2.2 Conversion ofthe transfer function established by 1.2.1 as the state spacedescription>>GI=ss(G)or>>[A,B,C,D]=tf2ss(num,den)1.3Thestatespacedescriptionofthesystem is[1 0-1]1X=010x+0u[o02][2]y=[1 2 1]Transform the state space description into the diagonal canonical formThe transformation can be executed by the following statements.>> A=[1 0 -1;0 1 0;0 0 2];>> b=[1 0 2];2
2 1.1.2 Conversion of state space model established by 1.1.1 as the transfer function model The model of transfer function can be obtained by the following statements: >>G1=tf(G) or >> [Num1,den1]=ss2tf(A,B,C,D,1) >>[Num2,den2]=ss2tf(A,B,C,D,2) 1.2 The model of the system is 2 3 4 5 5 ( ) 4 3 2 + + + + + = s s s s s G s 1.2.1 Set up the model of the transfer function of the system Transfer function can be input to the MATLAB work space by the following statements. >> num=[1,5];den=[1,2,3,4,5],G=tf(num,den) 1.2.2 Conversion of the transfer function established by 1.2.1 as the state space description >> G1=ss(G) or >> [A,B,C,D]=tf2ss(num,den) 1.3 The state space description of the system is [1 2 1] 2 0 1 0 0 2 0 1 0 1 0 1 = + − = y X X u Transform the state space description into the diagonal canonical form. The transformation can be executed by the following statements. >> A=[1 0 -1;0 1 0;0 0 2]; >> b=[1 0 2]';
>> c=[1 2 1];>>d=0;>>[V,D]=eig(A);>>[Ap,Bp,Cp,Dp]=ss2ss(A,b,c,d,inv(V))1.4 The state space description of the system is[010][0]X=001x2301y=[1 0 1]XTransform the state space description into the Jordan canonical form and thetransfor mation can be executed by the following statements.>> A=[ 0 1 0;0 0 1;2 3 0];>> b=[0 0 1]';>> c=[1 0 1];>>d=0;>>[V,J]=Jordan (A);>>[Ap,Bp,Cp,Dp]=ss2ss(A,b,c,d,inv(V))1.5ConsideraLTI systemdescribedbythefollowingstatespacedescription[。]x(0+101[ 0 ()=X(t)=-2-3y(t)=[1 2]X(t)1.5.1Determine the matrixexponiential function e4tThe matrix exponiential function can be derived by the following statements:>> A=[0 1;-2 -3];>> syms t>>eAt=expm(A*t)1.5.2 Find the zero input response of the system in the period [0, 6s]It canbecalculated bythefollowingstatements%Zeroinputresponse3
3 >> c=[1 2 1]; >> d=0; >> [V,D]=eig(A); >> [Ap,Bp,Cp,Dp]=ss2ss(A,b,c,d,inv(V)) 1.4 The state space description of the system is y X u [1 0 1] 1 0 0 2 3 0 0 0 1 0 1 0 = + X = X Transform the state space description into the Jordan canonical form and the transformation can be executed by the following statements. >> A=[ 0 1 0;0 0 1;2 3 0]; >> b=[0 0 1]'; >> c=[1 0 1]; >> d=0; >> [V,J]=Jordan (A); >> [Ap,Bp,Cp,Dp]=ss2ss(A,b,c,d,inv(V)) 1.5 Consider a LTI system described by the following state space description 01 0 1 ( ) ( ) ( ), (0) 23 1 2 ( ) [1 2] ( ) X t X t ut X yt Xt =+ = − − = 1.5.1 Determine the matrix exponiential function t e A . The matrix exponiential function can be derived by the following statements: >> A=[0 1;-2 -3]; >> syms t >> eAt=expm(A*t) 1.5.2 Find the zero input response of the system in the period [0, 6s] It can be calculated by the following statements. %Zero input response