宁波理工学院
现代控制理论课程设计报告
题
目 打印机皮带驱动系统能控能观和稳定性分析
项目成员 史旭东 童振梁 沈晓楠
专业班级 自动化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預頌圣鉉儐歲龈讶骅籴。
项目组成员具体分工
项目组成员 史旭东 沈晓楠 童振梁
具体完成的工作 设计课题，制作Word 内容仿真 PPT制作，答辩 备注
打印机皮带驱动系统能控能观和稳定性
分析
课程设计的内容如下： 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，而电机转矩为Mm=Kmi,于是有：
设作用在驱动皮带上的扰动转矩为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('\heta')
仿真结果如下：
5.总结
本次实验以打印机皮带和打印头为研究对象，讨论了将极点配置在期望的区域内的状态反馈控制方法。从仿真的结果可以看出，该方法可以保证系统具有一定的动态特性要求，夜间鼓励抑制外部扰动对系统的影响。由此可知，极点配置控制方法可以实现打印机皮带驱动系统平衡控制。从实验的研究结果还可看出，打印机皮带驱动系统是研究各种控制理论的一个不错的实验装置。驅踬髏彦浃绥譎饴憂锦。
参考书目：
［1］ ［2］ ［3］ ［4］ ［5］ 薛定宇，控制系统计算机辅助设计，清华大学出版社，2012年。 胡寿松，自动控制原理简明教程，科学出版社，2005年。 方水良，现代控制理论及其MATAB实践，2005年。
王晓凯，基于简化模型的倒立摆控制实验研究[J]，计算机技术与自动化，1997年。 王海英,袁丽英,吴勃.控制系统的MATLAB仿真与设计(第1版) [M]，高等教育出版社.2009年。猫虿驢绘燈鮒诛髅貺庑。 ［6］ 刘卫国，MATLAB程序设计与应用（第二版），2006
附录
源代码：
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('\heta')
因篇幅问题不能全部显示，请点此查看更多更全内容