华北水利水电学院计算机控制技术课程设计

计算机控制技术课程设计学院：电力学院专业：自动化班级：2009班姓名：学号：2009指导老师:徐俊红华北水利水电学院目录―、设计内容1.1对象模型1.2设计目的1.3设计要求二、设计方法及步骤2.1设计方法2.2设计具体步骤2.3matlab仿真结果展示三、仿真设计程序心得体会控制系统的状态空间设计1.1【设计对象】系统的对象模型为：G&)=s(s+4)(s+8)1.2【设计目的】A:试确定一个状态负反馈阵K,使相对于单位阵阶跃参考输入的输出过渡过程，满足如下的期望指标：超调量<=20%,峰值时间<=0.4soB:如果系统的状态变量在实际上无法测量，试确定一个状态观测器(全维状态观测器)，使得通过基于状态观测器的状态反馈，满足上述期望的性能指标。1.3【设计要求】要求学生掌握当Gc(s)设计好后如何将其变换为离散算法Gc(Z)以及如何将Gc(Z)转换在计算机上可完成计算的迭代方程。要求学生能掌握工业中常用的基本PID算法。掌握一阶向前，向后差分及双线性变换离散化的具体做法及应用场合。熟悉PID两种基本算法的计算公式：位置算法和增量算法。熟练使用MATLAB软件，掌握其仿真的方法、步骤及参数设置。

了解计算机控制系统的组成及相应设备的选用等问题。二.【设计方法及步骤】1、 利用Simulink进行仿真，判断是否满足期望的性能指标。系统仿真方框图如下：系统仿真结果如下：

有图可知，系统不满足期望的性能指标，需要进行配置。判断系统的能控能观性，确定系统是否能够通过状态反馈实现极点任意配置。能控性判别，能控性判别矩阵为[0 0 1]Q=0 1 -12口-12-112/系统的可控矩阵阶数为3，为满秩，则系统是能控的。由期望的性能指标求出闭环系统的期望极点。第一步：由典型二阶系统性能指标与系统参数之间的关系，确定系统参数，然后确定系统的主导极点和非主导极点。要求确定的主导极点使得闭环系统至少满足前述性能指标。由系统的性能指标：超调量＜=20%，峰值时间＜=0.4s。可以求出Z=0.4559Wn=8.82rad/s。因此可以选取Z=0.46Wn=8.9rad/s。可以求得：g*s)= 1 S2+8.2s+79.21第二步：由系统的特征方程可以求出系统的特征根：S=—4.1+j7.9S=—4.1-j7.92可以令非主导极点为：S3=-90确定状态反馈矩阵：受控系统的特征方程为F(s)=s3+12s2+32s；a=12a=32a=0；由期望的极点构成的特征多项式：f*(s)=(s+90)(s2+8.2s+79.21)=s3+98.2s2+817.21s+7128.9a*=98.2a*=817.21a*=712891 2 3 .所以状态反馈矩阵为：K』a；—a3,气—a2,a1*—a」=G128.9,785.21，86.2]确定输入放大系数对应的闭环传递函数为：、(s)= L s3+98.2s2+817.21s+7128.9要求跟踪阶跃信号的误差%=0，有(1_W(s)) 7128.9-LeP=0=limss\=limit-、(s)]=7128,9s—0 s—0所以L=7128.9对上面的初步结果，再用对跟踪速度信号的误差要求来验证，即-1W(s)一s一— —e=limLs2 s2」=0.115所以L满足要求。画出对应的能控规范型的闭环系统方块图L已知Wk(s)s3+98.2s2+817.21s+7128.97128.9=s3+98.2s2+817.21s+7128.97128.9s3+(12+86.2)s2+(32+785.21)s+7128.9其中可设W'(s)= 1k s3+(12+86.2)s2+(32+785.21)s+7128.9对应的规范型状态方程为：f〕•尤1-172£t3J>1=:0"7128.910-(32+785.21)0 :1-(12+86.2)/V尤1|x2>尤13J能控规范型的闭环系统方框图：

StepR&sponse1.20.4oO0*0.与3、配置Q.2由于期StepR&sponse1.20.4oO0*0.与3、配置Q.2由于期P=[-6.3237+10.4720j,-6.3237-1Q.4720j,-94.856];又由于希望观测器的响应要快于原系统的响应，配置状态观测器的极点应尽量离原极点距离远一些，取的点为大于十倍距离。判断系统的能观性由前可知系统能观测，可以进行配置。求状态观测器增益矩阵：L=1.0e+005*(3.2441-0.0172)求状态观测器设G=g2根据状态观测后的表达式写出特征式，与配置极点的特征式进行比较，求出G=1.0e+002*(-0.2213+0.3665j-0.2213-0.3665j-1.4228)t 绘制出带状态观测器的状态空间表达式如下图：由上述图运行结果可知，系统的阶跃响应比原系统有较快的响应，且满足所设定的条件，即所构建的状态观测器是完全可行的，满足题目要求。三.仿真设计程序Num=1;Den=[112321];Tf(num,den)[A,B,C,D]=tf2ss(num,den);N=length(A);Q=ctrb(A,B)N=rank(Q);Ifm==nDisp(‘systemstatevariablecanbetotallycontrolled(系统可控)’)P=[-4.1+j7.9,-4.1-j7.9,-90];K=acker(A,B,P)ElseEndS=size(A,1);Z=[zeros([1,s])1];N=inv([A,B,C,D])*Z;Nx=N(1:s);Nu=N(1+s);Nbar=Nu+K*NxAt=A-B*KBt=B*NbarCt=CDt=Dstep(At,Bt,Ct,Dt);holdon;num=1;den=[112321];tf(num,den)[A,B,C,D]=tf2ss(num,den)n=length(A);Qo=obsv(A,C)m=rank(Qo);ifm==ndisp('systemstatevariablecanbetotallycontrolled(系统可观)')P=[-10.5+7.85j-10.5-7.85j-100];K=acker(A,B,P);L=place(A',C',P)';elseends=size(A,1);Z=[zeros([1,s])1];N=inv([A,B;C,D])*Z';Nx=N(1:s);Nu=N(1+s);Nbar=Nu+K*Nx;At=[A-B*KB*Kzeros(size(A))A-L*C];Bt=[B*Nbarzeros(size(B))];Ct=[Czeros(size(C))];Dt=D;ss=sys(At,Bt,Ct,Dt);step(ss);holdon;grid;经过程序的仿真运行，与所求的系统的运行结果如下：

1■DTDno-d与StepResponse1■DTDno-d与StepResponse课程设计是每一个大学生在大学生涯中都不可或缺的，它使我们在实践中了巩固了所学的知识、在实践中锻炼自己的动手能力；实习又是对每一位大学生所学专业知识的一种拓展手段，它让我们学到了很多在课堂上根本就学不到的知识，不仅开阔了自己的视野，增战了自己的见识，也为我们以后进一步走向社会打下了坚实的基础，是我们走向以后走向工作岗位的奠基石。计算机控制技术是一门以电子技术、自动控制技术、计算机应用技术为基础，以计算机控制技术为核心，综合可编程控制技术、单片机技术、计算机网络技术，从而实现生产技术的精密化、生产设备的信息化、生产过程的自动化及机电控制系统的最佳化的专门学科。计算机控制系统技术的应用领域非常广泛，不但国防、航空航天等高精尖学科必不可少的组成部分尖学科必不可少的组成部分，而且在现代化的工、农、医等领域也发挥着重要作用。计算机控制作为一门理论与实践相结合并且具有极强的实用性的课