资源描述：
《自动控制原理实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验报告课程自动控制原理实验日期年月日专业班级测控1201班姓名XXXX学号XXXXX同组人实验名称系统的能控性与能观测性分析及状态反馈极点配置批阅教师签字一、实验目的1、学习掌握MATLAB控制工具箱中的基本命令的操作方法；2、加深理解能观测性、能控性等观念；3、掌握状态反馈系统的极点配置方法，研究不同配置对系统动态特性的影响。二、实验内容1、学习掌握MATLAB控制工具箱中的基本命令的操作方法，特别是数学模型的描述。（自选控制对象模型，应用以下命令，并写出结果）1）step,damp,pzmap,rlocus,rlocfind,bode,margin,
2、nyquist；2）tf2ss,ss2tf,tf2zp,zp2ss；3）ss2ss,jordan,canon,eig。2、能控性、能观测性分析（a）了解以下命令的功能；自选对象模型，进行运算，并写出结果。gram,ctrb,obsv,ctrbf,obsvf（b）已知连续系统的传递函数模型，，当a分别取-1，0，1时，判别系统的能控性与能观测性；（c）已知系统矩阵为，，，判别系统的能控性与能观测性；3、状态反馈极点配置原系统如图3-2所示。图中，X1和X2是可以测量的状态变量。图3-2系统结构图试设计状态反馈矩阵,使系统加入状态反馈后其动态性能指标满足给定的
3、要求:(1)已知：K=10，T=1秒，要求加入状态反馈后系统的动态性能指标为：σ%≤20%，ts≤1秒。(12)已知：K=1，T=0.05秒，要求加入状态反馈后系统的动态性能指标为：σ%≤5%，ts≤0.5秒。状态反馈后的系统，如图3-3所示：图3-3状态反馈后系统结构图分别观测状态反馈前后两个系统的阶跃响应曲线，并检验系统的动态性能指标是否满足设计要求三、实验环境MATLAB2014四、实验原理（或程序框图）及步骤1、系统能控性、能观性分析设系统的状态空间表达式如下：能控性直接判据：系统能控的充分必要条件是间接判别法：通过线性非奇异变换之后，根据中是否有
4、全零行判断状态能观测性判别式为：2、状态反馈极点配置一个受控系统只要其状态是完全能控的，则闭环系统的极点可以任意配置。极点配置有两种方法：①采用变换矩阵T，将状态方程转换成可控标准型，然后将期望的特征方程和加入状态反馈增益矩阵K后的特征方程比较，令对应项的系数相等，从而决定状态反馈增益矩阵K；②基于Carlay-Hamilton理论，它指出矩阵状态矩阵A满足自身的特征方程，改变矩阵特征多项式的值，可以推出增益矩阵K，这种方法推出增益矩阵K的方程式叫Ackermann公式。四、程序源代码1、熟悉指令文件：Stepnum=[0,0,0,2000];den=[1
5、,205,1000,20000];t=0:0.05:2.5;sys=tf(num,den);step(sys,t);grid;文件：Dampnum=[1];den=[1,10,10];sys=tf(num,den);[wzeta]=damp(sys)%输出结果w=1.12708.8730zeta=11文件：PzmapH=tf([251],[123]);pzmap(H)文件：Rlocusclearnum=[0001];den=conv([10],[13+sqrt(-1)]);den=conv(den,[13-sqrt(-1)]);rlocus(num,den
6、);v=[-82-44];axis(v);文件:Rlocfindrlocfind(sys)Selectapointinthegraphicswindowselected_point=-3.7239+0.0062ians=5.4641文件：Bodenum=[1000,1000];den=[10,10,50];bode(num,den)文件：Nyquistclearnum=[0202010];den=conv([110],[110]);nyquist(num,den)tf2ss,ss2tf,tf2zp,zp2ss；文件：Zp2ssclearZ=[-1;-3];
7、P=[0;-2;-4;-6];K=4;[A,B,C,D]=zp2ss(Z,P,K)A=-10.0000-4.8990004.8990000-6.0000-4.2866-2.00000001.00000B=[1010]TC=0004D=0文件：Ss2tfclearA=[01;1-2];B=[0;1];C=[13];D=[1];[num,den]=ss2tf(A,B,C,D)printsys(num,den,'s')[z,p,k]=ss2zp(A,B,C,D)num=150den=1.00002.0000-1.0000num/den=s^2+5s-------
8、------s^2+2s-1z=0-5p=-2.41420.414
