温度控制系统滞后校正

-z.题目:温度控制系统的滞后校正初始条件：*温箱的开环传递函数为要求完成的主要任务:〔包括课程设计工作量及其技术要求，以及说明书撰写等具体要求〕试用Matlab绘制其波特图和奈奎斯特图，计算相角裕度和幅值裕度；试设计滞后校正装置，使系统的相角裕度增加15度。用Matlab对校正后的系统进展仿真，画出阶跃相应曲线时间安排：任务时间〔天〕审题、查阅相关资料1分析、计算1.5编写程序1撰写报告1论文辩论0.5指导教师签名：年月日系主任〔或责任教师〕签名：年月日-z.温度控制系统的滞后校正1系统传递函数分析该传递函数由比例环节，延迟环节，积分环节，惯性环节组成。1.1比例环节比例环节的传递函数和频率特性:幅值特性和相频特性:对数幅频特性和对数相频特性:所以对数幅频特性L(ω)是ω轴线。1.2延迟环节延迟环节的传递函数和频率特性:幅频特性和相频特性:对数幅频特性和对数相频特性:由以上可知延迟环节不影响系统的幅频特性，只影响系统的相频特性。1.3积分环节积分环节的传递函数和频率特性:幅频特性和相频特性:积分环节的对数幅频特性和对数相频特性:由于Bode图的横坐标按lgω刻度，故上式可视为自变量为lgω，因变量为L(ω)的关系式，因此该式在Bode图上是一个直线方程式。直线的斜率为−20dB/dec。当ω=1时，−20lgω=0，即L(1)=0，所以积分环节的对数幅频特性是与ω轴相交于ω=1，斜率为−20dB/dec的直线。积分环节的相频特性是ω=−90°，相应的对数相频特性是一条位于ω轴下方，且平行于ω轴的水平直线。1.4惯性环节惯性环节的传递函数和频率特性:幅值特性和相频特性：对数幅频特性和对数相频特性：1.5开环传递函数将jω代替s带入传递函数求得：2用Matlab绘制波特图和奈奎斯特图，计算相角裕度和幅值裕度MATLAB提供的函数bode()和nyquist()不能直接绘制有延迟环节系统的波特图和奈奎斯特图。由于延迟环节不影响系统的幅频特性而只影响系统的相频特性，因此可通过对相频的处理结合绘图函数的应用来绘制具有纯延迟环节系统的波特图和奈奎斯特图。2.1绘制波特图绘制对数坐标图的程序如下：num=[1];%开环传递函数的分子den=conv([10],[41]);%开环传递函数的分母w=logspace(-2,1,100);%确定频率围[mag,phase,w]=bode(num,den,w);%计算频率特性的幅值和相角%利用相频特性求加上延迟环节后的相频特性phase1=phase-w*57.3*3;subplot(211),semilog*(w,20*log10(mag));%绘制幅频特性v=[0.01,10,-60,40];a*is(v)grid%绘制相频特性subplot(212),semilog*(w,phase1);v=[0.01,10,-270,-90];a*is(v)%设置坐标轴的标尺属性set(gca,'ytick',[-270-240-210-180-150-120-90])grid绘制的波特图如图2-1图2-1开环传递函数的bode图2.2绘制奈奎斯特图绘制极坐标图的程序如下：num=[1];%开环传递函数的分子den=conv([10],[41]);%开环传递函数的分母w=logspace(-1,2,100);%确定频率围[mag,phase,w]=bode(num,den,w);%计算频率特性的幅值和相角%利用相频特性求加上延迟环节后的相频特性phase1=phase*pi/180-w*3;holdon%用极坐标曲线绘制函数画出奈奎斯特图polar(phase1,mag)v=[-2.5,1,-1,1];a*is(v)grid绘制的奈奎斯特图如图2-2图2-2开环传递函数的nyquis图分析：因为有延迟环节，所以当ω从0→∞变化时，幅角也从0→∞变化，所以奈奎斯特曲线为螺旋线。2.3计算幅值裕度由相频特性曲线可知相角穿越频率2.4计算相角裕度由bode图可知截止频率所以系统不稳定。3设计滞后校正装置3.1无源滞后校正装置图3-1是由电阻和电容组成的无源滞后校正网络的电路图，其传递函数为其中,,b为滞后网络的分度系数。图3-1无源滞后校正网络3.2计算校正函数使系统的相角裕度增加15度。未校正系统，校正后又由得得到校正函数3.3检验相角裕度校正后系统的开环传递函数为绘制对数坐标图的程序如下：num=[25.81];%开环传递函数的分子den=conv([conv([10],[41])],[36.451]);%开环传递函数的分母w=logspace(-2,1,100);%确定频率围[mag,phase,w]=bode(num,den,w);%计算频率特性的幅值和相角%利用相频特性求加上延迟环节后的相频特性phase1=phase-w*57.3*3;subplot(211),semilog*(w,20*log10(mag));%绘制幅频特性v=[0.01,10,-40,60];a*is(v)grid%绘制相频特性subplot(212),semilog*(w,phase1);v=[0.01,10,-270,-90];a*is(v)%设置坐标轴的标尺属性set(gca,'ytick',[-270-240-210-180-150-120-90])grid绘制的波特图如图3-2所示由波特图可知，校正后函数的截止频率其相角裕度为满足给定指标要求。图3-2校正后开环传递函数bode图3.4校正装置参数设置经分析得T=(R1+R2)C,b=R2/(R1+R2)因此可设R2=2577ΩR1=3645Ω4校正后的系统4.1校正后的仿真电路在Matlab中输入simulink命令，进入仿真界面。仿真电路如图4-1图4-1校正后开环传递函数的Matlab仿真电路4.2校正后系统阶跃响应曲线阶跃响应曲线如图4-2图4-2校正后系统的阶跃响应曲线5心得体会通过本次?自动控制原理?课程设计，我稳固了对课本知识的理解，更熟悉了对Matlab的使用，并开场将理论运用于实践。为期5天的课程设计，时间不长，也算是对期末考试提前进展了笼统的复习，加深了我对自动控制原理的理解，扩展了自己的认识。另外，此次课程设计让我了解到，我们以后在学习生活中要培养好的习惯，做事情要按照规则办