自动控制原理课程设计说明书

  • 格式:doc
  • 大小:450.50 KB
  • 文档页数:19

《自动控制原理》课程设计说明书日期:2012.5.21-2012.6.1自动控制原理课程设计评阅书摘要MATLAB 是一种直观、高效的计算机语言,同时也是一个科学计算平台。

MATLAB 软件技术进行控制系统建模,时域、负域和频域分析,前馈校正设计,离散系统与非线性系统分析,以及现代控制理论中的可控性、可观性和稳定性的判别。

在MATLAB中,常用的系统建模方法有传递函数模型,零极点模型以及状态空间模型等。

在这次课程设计中,就使用了MATLAB软件对我们过去学过的一些简单的函数进行仿真。

此次学会了使用简单的MA TLAB语言编写简单的传递函数,还有使用Simulink进行仿真。

我们可将综合性和设计性实验项目通过MATLAB 在计算机上仿真,使系统的观察实验的动态过程。

目前,MATLAB 已经成为我们当代大学生必须掌握的基本技能,在设计研究单位和工业部门,MATLAB 已经成为研究和解决各种具体工程问题的一种标准软件。

在完成了验证性、综合性和设计性实验后,课程设计必不可少。

课程设计是工科实践教学的一个重要的环节,目的是培养我们综合运用理论知识分析和解决实际问题的方法和能力,实现由知识向技能的初步化。

所以课程设计是培养我们思维创造能力最有效的图径。

关键词:控制系统;MATLAB;模拟仿真;思维创造目录1 课题描述 (1)2 仿真过程 (2)2.1 控制系统建模 (2)2.2 线性系统时域分析 (2)2.3 线性系统根轨迹分析 (4)2.4 线性系统频域分析 (5)2.5 线性系统校正 (6)3 Simulink仿真 (7)总结 (11)参考文献 (12)1 课题描述MATLAB 是一种直观、高效的计算机语言,同时也是一个科学计算平台。

它的伴随工具Simulink 是用来对真实世界的动力学系统建模、模拟仿真和分析的软件。

我们可将综合性和设计性实验项目通过MATLAB 在计算机上仿真,使系统的观察实验的动态过程。

目前,MATLAB 已经成为我们当代大学生必须掌握的基本技能,在设计研究单位和工业部门,MATLAB 已经成为研究和解决各种具体工程问题的一种标准软件。

在完成了验证性、综合性和设计性实验后,课程设计必不可少。

课程设计是工科实践教学的一个重要的环节,目的是培养我们综合运用理论知识分析和解决实际问题的方法和能力,实现由知识向技能的初步化。

所以课程设计是培养我们思维创造能力最有效的途径。

在学习自动控制原理这个学科之后,这次课程设计又针对这个学科做了MATLAB 软件的仿真学习。

这次课程设计的基本要求是学会使用MATLAB软件并进行仿真,学会建立简单的数学模型,例如,传递函数模型,零极点模型以及状态空间模型等。

然后对这个软件进行一定的了解与认识。

而且还要学会使用简单的MARLAB语言和Simulink构建简单的数学模型。

通过学习培养理论联系实际的科学态度,训练综合运用经典控制理论和相关课程知识的能力。

掌握自动控制原理的时域分析法,根轨迹法,频域分析法,以及各种(矫正)装置的作用及用法,能够利用不同的分析方法对给定系统进行性能分析,能根据不同的系统性能指标要求进行合理的系统设计,并调试满足系统的指标。

学会使用MATLAB语言及Simulink动态仿真工具进行系统的仿真与调试。

锻炼独立思考和动手解决控制系统实际问题的能力。

2 设计过程 2.1 控制系统建模已知多回路反馈系统的结构如图,求闭环系统的传递函数C(S)/R (S)32)(1+=s s G ,11)(1+=s s G ,441)(222+++=s s s G S ,61)(1++=s S s G ,21)(1++=s S s H ,2)(2=s H ,1)(3=s H 。

解 MATLAB 文本如下G1=tf([1],[1 10]);G2=tf([1],[1 1]);G3=tf([1 0 1],[1 4 4]); numg4=[1 1];deng4=[1 6];G4=tf(numg4,deng4); H1=zpk([-1],[-2],1);numh2=[2];denh2=[1];H3=1;nh2=conv(numh2,deng4);dh2=conv(denh2,numg4); H2=tf(nh2,dh2); sys1=series(G3,G4);sys2=feedback(sys1,H1,+1); sys3=series(G2,sys2); sys4=feedback(sys3,H2); sys5=series(G1,sys4); sys=feedback(sys5,H3) 程序执行结果 Zero/pole/gain:0.083333 (s+1)^2 (s+2) (s^2 + 1)---------------------------------------------------------- (s+10.12) (s+2.44) (s+2.349) (s+1) (s^2 + 1.176s + 1.023)2.2 线性系统时域分析已知系统的闭环传递函数为1682^16++=s s ζφ,其中&=0.707,求二阶系统的单位脉冲响应,单位脉冲响应和单位斜坡响应。

解 MATLAB 文本如下zeta=0.707;num=[16];den=[1 8 * zeta 16]; sys=tf(num,den);p=roots(den)t=0:0.01:3;figure(1)impulse(sys,t);gridxlabel('t');ylabel('c(t)');title('impulse response'); figure(2)step(sys,t);gridxlabel('t');ylabel('c(t)');title('step response'); figure(3)u=t;lsim(sys,u,t,0);gridxlabel('t');ylabel('c(t)');title('ramp response');程序执行结果单位脉冲响应单位阶跃响应单位斜坡响应分析系统对输入信号导数的响应,就等于系统对该信号响应的导数;或者,系统对输入信号的积分的响应,就等于系统对该输入信号响应的积分,而积分常数由零输入初始条件确定。

2.3 线性系统根轨迹分析已知单位负反馈系统的开环传递函数为))(4(20)(K s s s G ++=试画出K 从零点变化到无穷时的根轨迹图,并求出系统临界阻尼时对应的K 值及其闭环极点。

解 系统闭环特征多项式为D (s)=2s +4s+Ks+20=2s +4s+20+K(s+4)=0等效开环传递函数204)4()(2*+++=s s s K s G MATLAB 文本如下G=TF([1 4],[1 4 20]); figure(1) pzmap(G); figure(2) rLocus(G); rlocfind(G) 结果得到如下运行结果 零极点分布图根轨迹图分析根轨迹的起点和终点,根轨迹起于开环极点,终于开环零点;根轨迹的分支数,对称性和连续性,根轨迹的分支数与开环有限零点数m 和有限极点数n 中大者相等,它们是连续并且对称于实轴;当开环有限极点数n 大于有限零点数m 时,有n-m 条根轨迹分支沿着与实轴交角为mn k a -+=∂π)12(; k=0,1,2,...n-m-12.4 线性系统频域分析已知单位负反馈系统的开环传递函数1624.32081.16042.246401280)(234+++++=s S S s s s G 试绘制其伯德图和奈奎斯图,并判断闭环系统的稳定性解 MATLAB 文本如下G=tf([1280 640],[1 2402 1604.81 320.24 16]); figure(1) margin(G);figure(2) nyquist(G); axis equal 运行结果 伯德图奈奎斯特图分析由于系统无右半边开环极点,从图可以看出,开环幅相曲线不包括(-1,j0)点,闭环系统稳定,另外在伯德图可以得出系统的幅值裕度h=29.5dB,相角裕度09.72=γ,相应的截止频率904.0=c ω,穿越频率9.39=x ω。

有奈式判知,系统闭环稳定2.5 线性系统校正设单位负反馈系统的开环传递函数为)1()(0+=s s Ks G 若要求系统在单位斜坡输入信号作用时,位置输出稳态误差rad e ss 1.0)(≤∞,开环系统截止频率s rad c/4.4≤''ω,相角裕度045≥''r ,幅值裕度dB h 10≥'',式设计串联无缘超前网络。

解1)根据稳态误差要求,确定开环增益K 。

2)利用以确定的开环增益,计算待校正系统的幅值裕度、相角裕度及其对应的截止频率、穿越频率。

3)根据截止频率cω''的要求,计算超前网络参数a 和T 。

为保证正系统的响应速度,并充分利用网络的橡胶超前特性,可选择最大超前角频率等于截止频率,即cm ωω''=,其中a 由a L L m C C lg 10)()(=='''-ωω确定,然后再由aT m ω1=,确定T 值。

4)确定无源超前网络和最大超前角m ∂;Ts aTs s aG c ++=11)(;11arcsin +-=∂a a m5)验算已校正系统的幅值裕度,相角裕度以及对应的截止频率,穿越频率。

若验算结果不满意指标要求,需要新选择)(c m ωω''=,然后重复以上的设计步骤。

MTALAB 文本如下 k=1/0.1;G0=zpk([],[0 -1],k); [h0,r,wx,ec]=margin(G0) wm=4.4;L=bode(G0,wm); Lwc=20*log10(L) a=10^(-0.1*Lwc) T=1/(wm*sqrt(a)); phi=asin((a-1)/(a+1))Gc=(1/a)*tf([a*T 1],[T 1]); Gc=a*Gc;G=Gc*G0;bode(G,'r',G0,'b--');grid; [h,f,wx,wc]=margin(G) 运行结果对数幅频特性曲线分析图中虚线部分为系统校正前的对数幅频特性曲线,实线部分为系统校正后的对数幅频特性曲线,如果存在误差时,重新选取截止频率"c ω,从而得到准确的实验数据3 Simulink 仿真已知系统的闭环传递函数为16&816)(2++=s s s ϑ,其中&=0.707,求二阶系统的单位脉冲响应,单位脉冲响应和单位斜坡响应。