《自动控制原理课程设计(温度控制系统的分析与校正)最终版本》由会员分享,可在线阅读,更多相关《自动控制原理课程设计(温度控制系统的分析与校正)最终版本(20页珍藏版)》请在金锄头文库上搜索。
1、学号:帮程设护题目温度控制系统的分析与校正学院自动化学院专业自动化班级姓名指导教师2012 年 1 月 2 日课程设计任务书10自动化学院学生姓名: 专业班级:指导教师: 工作单位:题 目: 温度控制系统的分析与校正初始条件:某温箱的开环传递函数为Gp(s)=昭要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等 具体要求)1、试用 Matlab 绘制其波特图和奈奎斯特图,计算相角裕度和幅值裕度2、试设计超前校正装置,使系统的相角裕度增加 10度;3、用 Matlab 对校正后的系统进行仿真,画出阶跃响应曲线。时间安排:任务时间(天)审题、查阅相关资料1.5分析、计算2.5
2、编写程序2.5撰写报告1论文答辩0.5指导教师签名:系主任(或责任教师)签名目录引 言 11 系统开环传递函数分析 21.1 比例环节-1 21.2 积分环节-1/S 21.3 惯性环节-1/(5s+1) 31.4 延迟环节-e-2s 31.5 开环传递函数-Gp(s) 3p2 利用 Matlab 分析传递函数 42.1 绘制波特图 42.2 绘制奈奎斯特图 62.3 计算相角裕度 72.4 计算幅值裕度 73 设计超前校正装置 83.1 无源超前校正装置 83.2 确定校正函数 93.2.1 估算校正函数93.2.2 检验相角裕度93.2.3 增大补偿角后确定校正函数103.3 校正装置参数
3、设置 124 校正后系统的仿真以及其阶跃响应曲线 134.1 仿真校正后的系统 134.2 阶跃响应曲线的绘制 13结束语 15参考文献 16引言本次课程设计要求运用所学的理论知识去分析并设计校正温度控制系统的开环传递函 数,并通过软件 Matlab 辅助设计。自动控制原理如今已经运用到我们的各个领域了,如温度控制、气压控制、水位控制、 航天控制等等,通过自动控制原理的运用极大的改变着我们的生活,使我们的生活变得简 单而又丰富多样。而此次课程设计要求我们利用自动控制原理相关知识以及Matlab软件, 并通过查阅相关资料来分析并校正一个温箱的开环传递函数。这不仅要求我们学好书本上 的知识,还要求
4、我们能够灵活的运用所学的理论知识,将理论知识与实践相结合,通过使 用功能强大的软件Matlab分析并协助解决相关问题。这将很大程度的锻炼我们的自主学习 能力以及动手解决实际问题的能力。通过此次课程设计应该掌握一些简单的开环传递函数 的分析与校正,更重要的是学会一种解决问题的方法,为我们将来分析和校正更复杂的传 递函数打下坚实基础。学会综合分析自动控制原理类问题在当今社会更显得尤为重要,因此,这次课程设计 将是一个良好的开端,为进一步深入研究打下良好的基础。1 系统开环传递函数分析系统的开环传递为:G (s)= 亠,容易看出传递函数G(s)由比例、积分、惯性、 ps(5s +1)p延迟四个环节共
5、同组成。1.1 比例环节-1系统传递函数Gp(s)的比例环节为1,它的基本特性如下:比例环节的传递函数为:G(s)=1,频率特性为:G(jw)=1;幅值特性为:A(w)=IG(jw)l=1,相频特性为:叽w)=ZG(jw)=O;对数幅频特性为:L(w)=201gA(w)=201g1=0,对数相频特性为:叽w)=0;对数幅频特性L(w)是w轴线。1.2 积分环节-1/S系统传递函数Gp(s)的比例环节为1/s,它的基本特性如下:积分环节的传递函数为:G(s) = 1/s,频率特性为:G(jw) = 1/jw = e-j90/w ;幅频特性为:A(w) = 11/jw I = 1/w,相频特性为:
6、p(w) = Z(1/jw) = -90;对数幅频特性为:L(w) = 20lgA(w) = 20lg(1/w) = -20lgw,对数相频特性为:叽w)= -90。由于Bode图的横坐标按lgw刻度,故上式可视为自变量为lgw,因变量为L(w)的关 系式,因此该式在Bode图上是一个直线方程式。直线的斜率为-20dB/dec。当=1时,- 201gw=0,即L(1) = 0,所以积分环节的对数幅频特性是与w轴相交于w = 1,斜率为- 20dB/dec的直线。积分环节的相频特性是p(w) = -90,相应的对数相频特性是一条位于w 轴下方,且平行于 w 轴的水平直线。1.3 惯性环节-1/(
7、5s+1)系统传递函数Gp(s)的比例环节为l/(l+5s),它的基本特性如下:惯性环节的传递函数为:G(s) = 1/(1+5s),频率特性为:G(jw) = 1/(1+j5w); 幅值特性为:A(w)=l1/1+j5wl,相频特性为:p(w)=Z(1/1+j5w) = -arctan5w; 对数幅频特性为:L(w)=201gA(w),对数相频特性为p(w)= -arctan5w。绘制惯性环节的对数幅频特性曲线时,可以将不同的 w 值代入上式逐点计算求出L(w),但通常是用渐近线的方法先画出曲线的大致图形,然后再加以精确化修正。1.4 延迟环节-e-2s系统传递函数Gp(s)的比例环节为e-
8、2s,它的基本特性如下:延迟环节的的传递函数为:G(s)=e-2s,频率特性为:G(jw)=e-2jw=cos2w-jsin2w;幅频特性为: A(w)=IG(jw)l=lcos2w-jsin2wl=1,相频特性为:p(w)=ZG(jw)= Ze-2jw= -arctan(sin2w/cos2w) = -2w(rad)= -57.3*2w();对数幅频特性为:L(w)=201gA(w)=201g1=0,对数相频特性为:叽w)= -57.3*2w( )。 由此可知,延迟环节并不影响系统的幅频特性,而只是影响系统的相频特性。1.5 开环传递函数-G (s)p根据以上的分析可知,开环函数Gp(s)的
9、幅频特性为:A(w)=1/(w*l1+j5wl);相频特性 为:p(w) = -pi/2-arctan5w-2w。2 利用 Matlab 分析传递函数MATLAB软件提供的函数bode()和nyquist()不能直接绘制具有延迟环节系统的伯德图 和奈奎斯特图,由于延迟环节不影响系统的幅频特性而只影响系统的相频特性的缘故,因 此可以通过对相频的处理结合绘图函数的应用来绘制具有纯延迟环节系统的伯德图和奈 奎斯特图,从而完成相应的曲线绘制。2.1 绘制波特图绘制对数坐标图的程序如下:num=1;%开环传递函数的分子den=conv(1 0,5 1);%开环传递函数的分母w=logspace(-2,1
10、,100);%确定频率范围mag,phase,w=bode(num,den,w);%计算频率特性的幅值和相角%利用相频特性求加上延迟环节后的相频特性 phase1=phase-w*57.3*2;%绘制幅频特性 subplot(211),semilogx(w,20*log10(mag); v=0.01,10,-60,40;axis(v)grid%绘制相频特性 subplot(212),semilogx(w,phase1); v=0.01,10,-270,-90;axis(v) %设置坐标轴的标尺属性 set(gca, ytick ,-270 -240 -210 -180 -150 -120 -9
11、0) grid绘制的伯德图如图 2-1所示:图 2-1 开环传递函数的 bode 图手工绘制伯德图的步骤:绘制系统开环对数幅频渐近特性曲线方法如下: 第一步:系统开环传递函数可分解为比例、积分、惯性、延迟四个环节;第二步:确定各环节的交接频率。非最小相位惯性环节:w =0.2,斜率减小20dB/dec,1最小交接频率为w =w =0.2;min 1第三步:绘制低频段渐近特性曲线。因为v=l,则低频渐近线斜率k= -20dB/dec,过(0.2,13.98)点;第四步,绘制频段w w 渐近特性曲线。w w 时,k= -40 dB/dec;min min第五步,绘制整个曲线于图上。 绘制系统开环对
12、数相频曲线,取若干个频率点,列表计算各点的相角并标注在对数坐 标图中,最后将各点光滑连接。2.2 绘制奈奎斯特图在设绘制极坐标图的程序如下:num=1;den=conv(1 0,5 1);w=logspace(-1,2,100); mag,phase,w=bode(num,den,w);%开环传递函数的分子%开环传递函数的分母%确定频率范围%计算频率特性的幅值和相角%利用相频特性求加上延迟环节后的相频特性phase1=phase*pi/180-w*2;hold on %用极坐标曲线绘制函数画出奈奎斯特图 polar(phase1,mag) v=-2.5,1,-1,1;axis(v)grid图
13、2-2 开环传递函数的 nyquis 图原因分析:因系统传递函数具有延迟环节,所以当w从08变化时,幅角也会从0T 变化,所以奈氏曲线为螺旋线。2.3 计算相角裕度由 bode 图在-1 时的值可知截止频率 wc=0.45rad/scY=180 + (w) = 180 - 90 - arctan5w - 2w =90 - arctan2.25 巧1.57cc c= -27.6 0用 Matlab 函数 margin 验算:gm,pm,wcg,wcp=margin(mag,phase1,w)得出:wcp = 0.4254pm = -23.5870所以系统不稳定。2.4 计算幅值裕度由相频特性曲线
14、在-180的值可知相角穿越频率 wg=0.30rad/s。g则IG(jw )1 = 1/w *l1+j5w I = 1.849g g gh = 1/|G(jwg)| = 0.541g用 Matlab 函数 margin 验算:gm,pm,wcg,wcp=margin(mag,phase1,w)得出:wcg = 0.2965gm = 0.53053 设计超前校正装置3.1 无源超前校正装置如图 3-1 所示,图为由电阻和电容所组成的无源超前校正网络的电路图,其传递函数为: Gc(s) = Uc(s)/Ur(s) = (R2(l+RlCs)/(Rl+R2)/(Rl+R2+RlR2Cs)/(Rl+R2) = a-l(aTs+l)/(Ts+l) 其中T=R1R2C/(R1+R2),a=(R1+R2)/R2, a为超前网络的分度系数。RlcR2OO图 3-l 无源超前校正网络采用无源超前网络进行串联校正时,整个系统的开环增益要下降a倍,因此需要提高 放大器增益加以补偿,如图3-2所示,此时的传递函数为:Gc(s) = Uc(s)/Ur(s) =(aTs+l)/(Ts+l)Rlc图 3-2 带有附加放大器的无源超前校正网络由上式可得,无源超