现代控制理论课程设计心得【模版】

宁波理工学院现代控制理论课程设计报告

题目打印机皮带驱动系统能控能观和稳定性分析项目成员史旭东童振梁沈晓楠

专业班级自动化112

指导教师何小其

分院信息分院

完成日期 2014-5-28

目录

1. 课程设计目的 (4)

2.课程设计题目描述和要求 (4)

3.课程设计报告内容 (4)

3.1 原理图 (4)

3.2 系统参数取值情况 (5)

3.3 打印机皮带驱动系统的状态空间方程 (5)

4. 系统分析 (8)

4.1 能控性分析 (8)

4.2 能观性分析 (8)

4.3 稳定性分析 (9)

5. 总结 (11)

项目组成员具体分工

打印机皮带驱动系统能控能观和稳定性

分析

课程设计的内容如下：

1.课程设计目的

综合运用自控现代理论分析皮带驱动系统的能控性、能观性以及稳定性，融会贯通并扩展有关方面的知识。加强大家对专业理论知识的理解和实际运用。培养学生熟练运用有关的仿真软件及分析，解决实际问题的能力，学会应用标准、手册、查阅有关技术资料。加强了大家的自学能力，为大家以后做毕业设计做很好的铺垫。

2.课程设计题目描述和要求

(1)环节项目名称：能控能观判据及稳定性判据

(2)环节目的：

①利用MATLAB分析线性定常系统的可控性和客观性。

②利用MATLAB进行线性定常系统的李雅普诺夫稳定性判据。

(3)环节形式：课后上机仿真

(4)环节考核方式：

根据提交的仿真结果及分析报告确定成绩。

(5)环节内容、方法：

①给定系统状态空间方程，对系统进行可控性、可观性分析。

②已知系统状态空间方程，判断其稳定性，并绘制出时间响应曲线验

证上述判断。

3.课程设计报告内容

3.1 原理图

在计算机外围设备中，常用的低价位喷墨式或针式打印机都配有皮带驱动器。它用于驱动打印头沿打印页面横向移动。图1给出了一个装有直流电机的皮

带驱动式打印机的例子。其光传感器用来测定打印头的位置，皮带张力的变化用于调节皮带的实际弹性状态。

图1 打印机皮带驱动系统

3.2 系统参数取值情况

表1打印装置的参数

3.3 打印机皮带驱动系统的状态空间方程

图2 打印机皮带驱动模型

状态空间建模及系统参数选择。图2为打印机皮带驱动器的基本模型。模型中记皮带弹性系数为k，滑轮半径为r，电机轴转角为θ，右滑轮的转角为θP，打印头质量为m，打印头的位移为y(t)。光传感器用来测量y(t)，光传感器的输出电压为v1，且v1=k1y。控制器输出电压为v2，对系统进行速度反馈，即有v2=k2dv1/dt。注意到y=rθp,可知皮带张力T1，T2分别为

于是作用在质量m上的皮带净张力为T1-T2=2k（rθ-y）=2kx1，其中x1=（rθ-y）为第一个状态变量，表示打印头实际位移y与预期位移rθ之间的位移差。则质量m的运动方程为

取第二个状态变量x2=dy/dt，于是有

定义第三个状态变量x3=dθ/dt，x1的导数

推导电机旋转的运动方程：当L=0时，电机电枢电流i=v2/r，而电机转矩为M m=K m i,于是有：

设作用在驱动皮带上的扰动转矩为Md，则电机驱动皮带的有效转矩为M=Mm-Md。显然，只有有效转矩驱动电机轴带动滑轮运动，因此有

由于

故得

在上式中代入（1-3）以及

得到

最后可得

式（1-1）、（1-2）、（1-3）构成了描述打印机皮带驱动系统的一阶运动微分方程组，其向量矩阵形式为

将表1打印装置的参数代入得

4. 系统分析

4.1 能控性分析

根据能控性的秩判据

经计算可控性判别阵为

由上可知，rank(Sc)=3=rank（A），故系统能控。

4.2 能观性分析

根据能观性的秩判据

经计算可观测性判别阵为

由上可知，rank(S0)=3=rank（A），故系统能观。

4.3 稳定性分析

由打印机皮带驱动器系统的状态方程，可利用MATLAB求出其特征值：

代码：

A=[0 -1 0.015;200 0 0;-60 -8 -25];

Q=eye(3);

p=lyap(A',Q) val=eig(A)

结果：p =

1.0e+004 *

-6.0523 -0.0008 -0.0028

-0.0008 -0.0303 0.0001

-0.0028 0.0001 0.0000

val =

0.0009 +14.1755i

0.0009 -14.1755i

-25.0018

解的特征值为0.0090+14.4724i,0.0090-14.4724i,-25.0180。三个特征值中存在两个正实部根，一个负根，这说明打印机皮带驱动器系统，即被控系统是不稳定的。

采用MATLAB对被控对象进行仿真，如下图所示为打印机皮带驱动器没有添加任何控制器下三个变量的单位阶跃响应。如图可知，系统不稳定，不能到达控制目的。

代码:

A=[0 -1 0.015;200 0 0;-60 -8 -25]; b=[0;0;-100];

c=[1,0,0];

d=0;

sys0=ss(A,b,c,d);

t=0:0.01:5;

[y,t,x]=step(sys0,t);

subplot(2,2,1);

plot(t,x(:,1));grid

xlabel('t(s)');ylabel('x(t)');

title('z');

subplot(2,2,2);

plot(t,x(:,2));grid;

xlabel('t(s)');ylabel('x(t)');

title('z的微分');

subplot(2,2,3);

plot(t,x(:,3));grid

xlabel('t(s)');ylabel('x(t)');

title('\theta')

仿真结果如下：