控制系统数字仿真.doc

现代工程控制理论实验报告实验名称:控制系统数字仿真技术实验时间:2015/5/3目录一、 实验目的 3二、 实验内容 3三、 实验原理 3四、 实验方案 61、 分别离散法; 62、 整体离散法; 73、 欧拉法 94、 梯形法 105、 龙格——库塔法 11五、 实验结论 12小结: 14实验目的探究多阶系统状态空间方程的求解;探究多种控制系统数字仿真方法并对之进行精度比较;实验内容对上面的系统进行仿真,运用分别离散法进行分析;对上面的系统进行仿真,运用整体离散法进行分析;对上面的系统进行仿真,运用欧拉法进行分析;对上面的系统进行仿真,运用梯形法进行分析;对上面的系统进行仿真,运用龙泽——库塔法进行分析;对上面的几种方法进行总计比较,对他们的控制精度分别进行分析比较;实验原理控制系统状态空间方程整体离散法的求解;控制系统的传递函数一般为有两种控制框图简化形式如下:KI控制器可以用框图表示如下:惯性环节表示如下:高阶系统的框图如下对于上面的框图可以简写传递函数根据各环节间的关系可以列写出式子中出现的系数A、B、C和D,下面进行整体离散法求传递函数的推导这样,如果知道系数,就可以知道高阶系统的传递函数和状态空间方程。在控制系统的每一个环节都加一个采样开关,构成分别离散法求解系统的状态空间方程;采样开关其实是一个零阶保持器比例环节:积分环节:惯性环节:实验方案分别离散法;系统框图根据上面提到的分别离散法得到仿真的公式已知系数:K1=;K2=;T1=;T2=;n1=2;n2=4;kp1=;ki1=;kp2=2;ki2=;惯性环节的系数:fai1=exp(-dt/T1);faiM1=1-fai1;fai2=exp(-dt/T2);faiM2=1-fai2;PID控制环节:up1=e*kp1;x(1)=x(1)+ki1*dt*e;up2=e1*kp2;x(2)=x(2)+ki2*dt*e1;惯性环节:x(3)=fai1*x(3)+K1*faiM1*u1;x(4)=fai1*x(4)+faiM1*x(3);x(5)=fai2*x(5)+K2*faiM2*x(4);x(6)=fai2*x(6)+faiM2*x(5);x(7)=fai2*x(7)+faiM2*x(6);x(8)=fai2*x(8)+faiM2*x(7);整体离散法;将系统框图拆开系统的状态空间方程为:可以得到此时状态方程的系数由上面的推导可知求出就可以得到系统的状态空间方程在Matlab中仿真时为fori=1:n1*n2faiM=faiM+(dt^i)*(a^(i-1))/factorial(i);endfai=faiM*a+eye(n1*n2);faiM=faiM*b;forj=1:lpx=fai*x+faiM*r;y=c*x+d*r;y1=[y1y];t=[tj*dt];end欧拉法由上面已经求出系统的状态空间方程,所以这里直接引用,欧拉法的求解过程如下:在Matlab中的仿真程序如下:fori=1:lpxk=a*x+b*r;x=x+xk*dt;y=c*x+d*r;y1=[y1y];t=[tdt*i];end梯形法类似于欧拉法,梯形法的推导如下在Matlab中仿真的程序如下:fori=1:lpxk=a*x+b*r;xk1=x+dt*xk;xk2=a*xk1+b*r;E=(xk+xk2)/2;x=x+dt*E;y=c*x+d*r;y1=[y1y];t=[tdt*i];end龙格——库塔法推导如下: