自动控制原理》实验2(线性系统时域响应分析
- 格式:doc
- 大小:309.00 KB
- 文档页数:12
系别:机电工程学院专业:课程名称:自动控制原理实验班级:姓名:学号:组别:实验名称:典型系统的时域响应和稳定性分析实验时间:学生成绩:教师签名:批改时间:一、目的要求1.研究二阶系统的特征参量 (ξ、ωn) 对过渡过程的影响。
2.研究二阶对象的三种阻尼比下的响应曲线及系统的稳定性。
3.熟悉 Routh 判据,用 Routh 判据对三阶系统进行稳定性分析。
二、实验设备PC机一台,TD—ACC教学实验系统一套三、实验原理及内容1.典型的二阶系统稳定性分析(1) 结构框图:如图 1.2-1 所示。
图1.2-2(2) 对应的模拟电路图:如图 1.2-2 所示。
图1.2-2系别:机电工程学院专业:课程名称:自动控制原理实验班级:姓名:学号:组别:实验名称:实验时间:学生成绩:教师签名:批改时间:(3) 理论分析系统开环传递函数为:;开环增益:(4) 实验内容先算出临界阻尼、欠阻尼、过阻尼时电阻 R 的理论值,再将理论值应用于模拟电路中,观察二阶系统的动态性能及稳定性,应与理论分析基本吻合。
在此实验中(图 1.2-2),系统闭环传递函数为:其中自然振荡角频率:2.典型的三阶系统稳定性分析(1) 结构框图:如图 1.2-3 所示。
系别:机电工程学院专业:课程名称:自动控制原理实验班级:姓名:学号:组别:实验名称:实验时间:学生成绩:教师签名:批改时间:图 1.2-3(2)模拟电路图:如图1.2-4 所示。
图 1.2-4(3)理论分析:系统的特征方程为:(4)实验内容:实验前由Routh 判断得Routh 行列式为:系别:机电工程学院专业:课程名称:自动控制原理实验班级:姓名:学号:组别:实验名称:实验时间:学生成绩:教师签名:批改时间:为了保证系统稳定,第一列各值应为正数,所以有五、实验步骤1.将信号源单元的“ST”端插针与“S”端插针用“短路块”短接。
由于每个运放单元均设臵了锁零场效应管,所以运放具有锁零功能。
竭诚为您提供优质文档/双击可除线性系统时域分析实验报告篇一:自动控制原理实验报告《线性控制系统时域分析》实验一线性控制系统时域分析1、设控制系统如图1所示,已知K=100,试绘制当h 分别取h=0.1,0.20.5,1,2,5,10时,系统的阶跃响应曲线。
讨论反馈强度对一阶系统性能有何影响?图1答:A、绘制系统曲线程序如下:s=tf(s);p1=(1/(0.1*s+1));p2=(1/(0.05*s+1));p3=(1/(0.02*s+1) );p4=(1/(0.01*s+1));p5=(1/(0.005*s+1));p6=(1/(0.002 *s+1));p7=(1/(0.001*s+1));step(p1);holdon;step(p2); holdon;step(p3);holdon;step(p5);holdon;step(p6);hol don;step(p7);holdon;b、绘制改变h系统阶跃响应图如下:stepResponse1.41.21Amplitude0.80.60.40.200.050.10.150.20.250.30.350.40.450.5Time(seconds)结论:h的值依次为0.1、0.2、0.5、1、2、5、10做响应曲线。
matlab曲线默认从第一条到第七条颜色依次为蓝、黄、紫、绿、红、青、黑,图中可知随着h值得增大系统上升时间减小,调整时间减小,有更高的快速性。
2?n?(s)?22,设已知s?2??ns??n2、二阶系统闭环传函的标准形式为?n=4,试绘制当阻尼比?分别取0.2,0.4,0.6,0.8,1,1.5,2,5等值时,系统的单位阶跃响应曲线。
求出?取值0.2,0.5,0.8时的超调量,并求出?取值0.2,0.5,0.8,1.5,5时的调节时间。
讨论阻尼比变化对系统性能的影响。
答:A、绘制系统曲线程序如下:s=tf(s);p1=16/(s^2+1.6*s+16);p2=16/(s^2+3.2*s+16);p3=16/(s^ 2+4.8*s+16);p4=16/(s^2+6.4*s+16);p5=16/(s^2+8*s+16) ;p6=16/(s^2+12*s+16);p7=16/(s^2+16*s+16);p8=16/(s^2 +40*s+16);step(p1);holdon;step(p2);holdon;step(p3); holdon;step(p4);holdon;step(p5);holdon;step(p6);hol don;step(p7);holdon;step(p8);holdon;b、绘制系统阶跃响应图如下:c、?取值为0.2、0.5、0.8、1.5、5时的参数值。
实验2——时域系统分析和线性系统的稳定性研究1. 研究性教学目的① 学习二阶系统阶跃响应曲线的实验测试方法;② 研究二阶系统的两个重要参数n ωζ,对阶跃瞬态响应指标的影响;③ 研究线性系统的开环比例系数K 对稳定性的影响; ④ 研究线性系统的时间常数T 对稳定性的影响。
2. 知识点训练① 自行设计二阶及三阶系统电路。
② 选择好必要的参数值, 计算出相应的阶跃响应数值, 进行仿真分析。
3. 研究性要求① 自行设计二阶及三阶仿真电路, 可以使用Proteus\Multisim\EWB\Matlab 等仿真软件;② 针对各环节选择好必要的参数值, 理论计算环节的输出响应; ③ 仿真分析环节的在典型输入信号作用下的输出波形。
4. 研究性内容实验一:典型二阶系统方块图和实现电路如图1-1所示。
图1-1 二阶系统闭环传递函数如下:, (T 是时间常数)。
各运算放大器运算功能: OP1, 积分, ; OP2, 积分, ; OP9, 反相, (-1);OP6, 反相比例, 。
可以得到:31010021211⨯====f n R k RCT ζω实验一步骤, 使, , 取, , 使T=0.47s, , 加入单位阶跃扰动, 记录响应曲线, 记作[1]。
仿真结果实验结果通过游标读图可以得到以下数据:保持不变, 单位阶跃扰动不变, 取, , 使T=1.47s, , 记录响应曲线, 记作[2]。
仿真结果实验结果保持不变, 单位阶跃扰动不变, 取, ,使T=1.0s, , 记录响应曲线, 记作[3]。
仿真结果实验结果保持不变, 单位阶跃扰动不变, 取, 使k=0.8, , 记录响应曲线, 记作[4]。
仿真结果实验结果保持不变, 单位阶跃扰动不变, 取, 使k=2.0, , 记录响应曲线, 记作[5]。
仿真结果实验结果要求: 将曲线[1]、[2]、[3]进行对比, [3]、[4]、[5] 进行对比, 将[3]中的和理论值进行比较。
线性系统的时域分析实验报告线性系统的时域分析实验报告引言:线性系统是控制理论中的重要概念,它在工程领域中有广泛的应用。
时域分析是研究线性系统的一种方法,通过对系统输入和输出的时域信号进行观察和分析,可以得到系统的动态特性。
本实验旨在通过对线性系统进行时域分析,探究系统的稳定性、阶数和频率响应等特性。
实验一:稳定性分析稳定性是线性系统的基本性质之一,它描述了系统对于不同输入的响应是否趋于有界。
在本实验中,我们选取了一个简单的一阶系统进行稳定性分析。
首先,我们搭建了一个一阶系统,其传递函数为H(s) = 1/(s+1),其中s为复变量。
然后,我们输入了一个单位阶跃信号,观察系统的输出。
实验结果显示,系统的输出在输入信号发生变化后,经过一段时间后稳定在一个有限的值上,没有出现发散的情况。
因此,我们可以判断该系统是稳定的。
实验二:阶数分析阶数是线性系统的另一个重要特性,它描述了系统的动态响应所需的最小延迟时间。
在本实验中,我们选取了一个二阶系统进行阶数分析。
我们搭建了一个二阶系统,其传递函数为H(s) = 1/(s^2+2s+1)。
然后,我们输入了一个正弦信号,观察系统的输出。
实验结果显示,系统的输出在输入信号发生变化后,经过一段时间后才稳定下来。
通过进一步分析,我们发现系统的输出波形具有两个振荡周期,这表明系统是一个二阶系统。
实验三:频率响应分析频率响应是线性系统的另一个重要特性,它描述了系统对于不同频率输入信号的响应情况。
在本实验中,我们选取了一个低通滤波器进行频率响应分析。
我们搭建了一个低通滤波器,其传递函数为H(s) = 1/(s+1),其中s为复变量。
然后,我们输入了一系列不同频率的正弦信号,观察系统的输出。
实验结果显示,随着输入信号频率的增加,系统的输出幅值逐渐减小,表明系统对高频信号有较强的抑制作用。
这一结果与低通滤波器的特性相吻合。
结论:通过以上实验,我们对线性系统的时域分析方法有了更深入的了解。
实验二 控制系统稳定性分析和时域响应分析一、实验目的与要求1、熟悉系统稳定性的Matlab 直接判定方法和图形化判定方法;2、掌握如何使用Matlab 进行控制系统的动态性能指标分析;3、掌握如何使用Matlab 进行控制系统的稳态性能指标分析。
二、实验类型设计三、实验原理及说明1. 稳定性分析 1)系统稳定的概念经典控制分析中,关于线性定常系统稳定性的概念是:若控制系统在初始条件和扰动共同作用下,其瞬态响应随时间的推移而逐渐衰减并趋于原点(原平衡工作点),则称该系统是稳定的,反之,如果控制系统受到扰动作用后,其瞬态响应随时间的推移而发散,输出呈持续震荡过程,或者输出无限偏离平衡状态,则称该系统是不稳定的。
2)系统特征多项式以线性连续系统为例,设其闭环传递函数为nn n n mm m m a s a s a s a b s b s b s b s D s M s ++++++++==----11101110......)()()(φ 式中,n n n n a s a s a s a s D ++++=--1110...)(称为系统特征多项式;0...)(1110=++++=--n n n n a s a s a s a s D 为系统特征方程。
3)系统稳定的判定对于线性连续系统,其稳定的充分必要条件是:描述该系统的微分方程的特征方程具有负实部,即全部根在左半复平面内,或者说系统的闭环传递函数的极点均位于左半s 平面内。
对于线性离散系统,其稳定的充分必要条件是:如果闭环系统的特征方程根或者闭环传递函数的极点为n λλλ,...,21,则当所有特征根的模都小于1时,即),...2,1(1n i i =<λ,该线性离散系统是稳定的,如果模的值大于1时,则该线性离散系统是不稳定的。
4)常用判定语句2.动态性能指标分析系统的单位阶跃响应不仅完整反映了系统的动态特性,而且反映了系统在单位阶跃信号输入下的稳定状态。
实验二 线性系统时域响应分析一、实验目的1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。
2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。
3.熟练掌握系统的稳定性的判断方法。
二、基础知识及MATLAB 函数(一)基础知识时域分析法直接在时间域中对系统进行分析,可以提供系统时间响应的全部信息,具有直观、准确的特点。
为了研究控制系统的时域特性,经常采用瞬态响应(如阶跃响应、脉冲响应和斜坡响应)。
本次实验从分析系统的性能指标出发,给出了在MATLAB 环境下获取系统时域响应和分析系统的动态性能和稳态性能的方法。
用MATLAB 求系统的瞬态响应时,将传递函数的分子、分母多项式的系数分别以s 的降幂排列写为两个数组num 、den 。
由于控制系统分子的阶次m 一般小于其分母的阶次n ,所以num 中的数组元素与分子多项式系数之间自右向左逐次对齐,不足部分用零补齐,缺项系数也用零补上。
1.用MATLAB 求控制系统的瞬态响应1)阶跃响应求系统阶跃响应的指令有:step(num,den) 时间向量t 的范围由软件自动设定,阶跃响应曲线随即绘出step(num,den,t) 时间向量t 的范围可以由人工给定(例如t=0:0.1:10)[y ,x]=step(num,den) 返回变量y 为输出向量,x 为状态向量在MATLAB 程序中,先定义num,den 数组,并调用上述指令,即可生成单位阶跃输入信号下的阶跃响应曲线图。
考虑下列系统:25425)()(2++=s s s R s C 该系统可以表示为两个数组,每一个数组由相应的多项式系数组成,并且以s 的降幂排列。
则MATLAB 的调用语句:num=[0 0 25]; %定义分子多项式den=[1 4 25]; %定义分母多项式step(num,den) %调用阶跃响应函数求取单位阶跃响应曲线grid %画网格标度线xlabel('t/s'),ylabel('c(t)') %给坐标轴加上说明title('Unit-step Respinse of G(s)=25/(s^2+4s+25)') %给图形加上标题名则该单位阶跃响应曲线如图2-1所示:注意:在figure 中点鼠标右键,在右键菜单中选择“Characteristics”,其中包括四个系统性能指标:“Peak Response 峰值”、“Settling Time 调节时间”、“Rise Time”和“Steady State 稳态值”,选中其中的任何一个指标后,都会用大点点在图上标出指标对应的位置。
实验二 线性系统时域响应分析一、实验目的1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。
2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。
二、基础知识及MATLAB 函数(一)基础知识时域分析法直接在时间域中对系统进行分析,可以提供系统时间响应的全部信息,具有直观、准确的特点。
为了研究控制系统的时域特性,经常采用瞬态响应(如阶跃响应、脉冲响应和斜坡响应)。
本次实验从分析系统的性能指标出发,给出了在MATLAB 环境下获取系统时域响应和分析系统的动态性能和稳态性能的方法。
用MATLAB 求系统的瞬态响应时,将传递函数的分子、分母多项式的系数分别以s 的降幂排列写为两个数组num 、den 。
由于控制系统分子的阶次m 一般小于其分母的阶次n ,所以num 中的数组元素与分子多项式系数之间自右向左逐次对齐,不足部分用零补齐,缺项系数也用零补上。
1.用MATLAB 求控制系统的瞬态响应1)阶跃响应 求系统阶跃响应的指令有:step(num,den) 时间向量t 的范围由软件自动设定,阶跃响应曲线随即绘出step(num,den,t) 时间向量t 的范围可以由人工给定(例如t=0:0.1:10)[y ,x]=step(num,den) 返回变量y 为输出向量,x 为状态向量在MATLAB 程序中,先定义num,den 数组,并调用上述指令,即可生成单位阶跃输入信号下的阶跃响应曲线图。
考虑下列系统:25425)()(2++=s s s R s C 该系统可以表示为两个数组,每一个数组由相应的多项式系数组成,并且以s的降幂排列。
则MATLAB的调用语句:num=[0 0 25]; %定义分子多项式den=[1 4 25]; %定义分母多项式step(num,den) %调用阶跃响应函数求取单位阶跃响应曲线grid %画网格标度线xlabel(‘t/s’),ylabel(‘c(t)’) %给坐标轴加上说明title(‘Unit-step Respinse of G(s)=25/(s^2+4s+25)’) %给图形加上标题名则该单位阶跃响应曲线如图2-1所示:为了在图形屏幕上书写文本,可以用text命令在图上的任何位置加标注。
例如:text(3.4,-0.06,’Y1’) 和text(3.4,1.4,’Y2’)第一个语句告诉计算机,在坐标点x=3.4,y=-0.06上书写出’Y1’。
类似地,第二个语句告诉计算机,在坐标点x=3.4,y=1.4上书写出’Y2’。
若要绘制系统t在指定时间(0-10s)内的响应曲线,则用以下语句:num=[0 0 25];den=[1 4 25];t=0:0.1:10;step(num,den,t)即可得到系统的单位阶跃响应曲线在0-10s间的部分,如图2-2所示。
2)脉冲响应图2-1 二阶系统的单位阶跃响应图2-2 定义时间范围的单位阶跃响应①求系统脉冲响应的指令有:impulse (num,den) 时间向量t的范围由软件自动设定,阶跃响应曲线随即绘出impulse (num,den,t) 时间向量t的范围可以由人工给定(例如t=0:0.1:10)[y,x]=impulse(num,den) 返回变量y为输出向量,x为状态向量[y,x,t]=impulse(num,den,t)向量t 表示脉冲响应进行计算的时间例:试求下列系统的单位脉冲响应:12.01)()()(2++==sssGsRsC在MATLAB中可表示为num=[0 0 1];den=[1 0.2 1];impulse(num,den)gridt itle(‘Unit-impulse Response of G(s)=1/(s^2+0.2s+1)’)由此得到的单位脉冲响应曲线如图2-3所示:②求脉冲响应的另一种方法应当指出,当初始条件为零时,G (s)的单位脉冲响应与sG(s)的单位阶跃响应相同。
考虑在上例题中求系统的单位脉冲响应,因为对于单位脉冲输入量,R(s)=1所以sssssssGsCsRsC112.012.01)()()()(22⨯++=++===图2-3 二阶系统的单位脉冲响应因此,可以将G(s)的单位脉冲响应变换成sG(s)的单位阶跃响应。
向MATLAB 输入下列num 和den ,给出阶跃响应命令,可以得到系统的单位脉冲响应曲线如图2-4所示。
num=[0 1 0];den=[1 0.2 1];step(num,den) gridtitle(‘Unit-step Response ofsG(s)=s/(s^2+0.2s+1)’)3)斜坡响应MATLAB 没有直接调用求系统斜坡响应的功能指令。
在求取斜坡响应时,通常利用阶跃响应的指令。
基于单位阶跃信号的拉氏变换为1/s ,而单位斜坡信号的拉氏变换为1/s 2。
因此,当求系统G(s)的单位斜坡响应时,可以先用s 除G(s),再利用阶跃响应命令,就能求出系统的斜坡响应。
例如,试求下列闭环系统的单位斜坡响应。
11)()(2++=s s s R s C 对于单位斜坡输入量,R(s)=1/s 2 ,因此 ss s s s s s s C 1)1(1111)(222⨯++=⨯++=在MATLAB 中输入以下命令,得到如图2-5所示的响应曲线: num=[0 0 0 1]; den=[1 1 1 0];step(num,den)title(‘Unit -Ramp Response Cuve for System G(s)=1/(s^2+s+1)’)图2-4 单位脉冲响应的另一种表示法2. 特征参量ζ和n ω对二阶系统性能的影响 标准二阶系统的闭环传递函数为:2222)()(nn n s s s R s C ωζωω++= 二阶系统的单位阶跃响应在不同的特征参量下有不同的响应曲线。
1)ζ对二阶系统性能的影响设定无阻尼自然振荡频率)/(1s rad n =ω,考虑5种不同的ζ值:ζ=0,0.25,0.5,1.0和2.0,利用MATLAB 对每一种ζ求取单位阶跃响应曲线,分析参数ζ对系统的影响。
为便于观测和比较,在一幅图上绘出5条响应曲线(采用“hold ”命令实现)。
num=[0 0 1]; den1=[1 0 1]; den2=[1 0.5 1]; den3=[1 1 1]; den4=[1 2 1]; den5=[1 4 1];t=0:0.1:10; step(num,den1,t)gridtext(4,1.7,’Zeta=0’); holdstep(num,den2,t) text (3.3,1.5,’0.25’) step(num,den3,t) text (3.5,1.2,’0.5’) step(num,den4,t) text (3.3,0.9,’1.0’) step(num,den5,t) text (3.3,0.6,’2.0’)title(‘Step -Response Curves for G(s)=1/[s^2+2(zeta)s+1]’)由此得到的响应曲线如图2-6所示:2),n ω对二阶系统性能的影响 同理,设定阻尼比25.0=ζ时,当n ω分别取1,2,3时,利用MATLAB 求取单位阶跃响应曲线,分析参数n ω对系统的影响。
num1=[00 1]; den1=[10.5 1];t=0:0.1:10; step(num1,den1,t); grid; hold on text(3.1,1.4,’wn =1’)num2=[0 0 4]; den2=[1 1 4]; step(num2,den2,t); hold on text(1.7,1.4,’wn=2’)num3=[0 0 9]; den3=[1 1.5 9]; step(num3,den3,t); hold on text(0.5,1.4,’wn=3’) 由此得到的响应曲线如图2-7所示:三、实验内容图2-6 ζ不同时系统的响应曲线图2-7 n ω不同时系统的响应曲线1.观察函数step( )调用格式,假设系统的传递函数模型为S 146473)(2342++++++=s s s s s s s G 绘制出系统的阶跃响应曲线。
2,典型二阶系统:2222)(nn n s s s G ωζωω++= 分别绘出)/(2s rad n =ω,ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线,分析参数ζ对系统的影响绘制出当ζ=0.25, n ω分别取1,2,4,6时单位阶跃响应曲线,分析参数n ω对系统的影响。
四、实验报1.根据内容要求,写出调试好的MATLAB语言程序,及对应的MATLAB运算结果。
2. 记录各种输出波形,根据实验结果分析参数变化对系统的影响。
3.写出实验的心得与体会。
五、预习要求1. 预习实验中基础知识,运行编制好的MATLAB语句,熟悉MATLAB指令及step( )和impulse( )函数。
ω对二阶系统性能的影响。
2. 结合实验内容,提前编制相应的程序。
3.思考特征参量ζ和n六,实验心的与体会通过这次实验,自己觉得对matlab这款软件有了点新的认识。
初步学会了对该软件的基本操作,对相关数据以及方法的处理,使得自己的知识得到了一点拓展,能过独立的进行相关实验的基本操作,以及数据的处理。
不过,同时也发现了自己的相关知识的不足,比如说软件的功能的具体操作不够,以及相关基础知识的清晰认识不足,特别是英语基础的不足,不能过灵活的选择相关操作及处理。
所以在今后的学习深造中要学习和掌握的还有许多,有待更加努力。
(注:本资料素材和资料部分来自网络,仅供参考。
请预览后才下载,期待您的好评与关注!)。