自动控制原理实验报告 (2)

合集下载

自动控制实验2实验报告

自动控制实验2实验报告

⾃动控制实验2实验报告:实验报告项⽬名称: MATLAB⽤于时域分析课程名称: ⾃动控制原理信息科学与⼯程学院通信⼯程系⼀、实验名称:MATLAB⽤于时域分析⼆、1)⼀阶系统响应sys1=tf([100],[1 0]);sys2=tf([0.1],[1]);sys=feedback(sys1,sys2);step(sys)1)⼆阶系统响应%Wn=1;t=0:0.1:12;num=[1];zetal=0;den1=[1 2*zetal 1]; zeta3=0.3; den3=[1 2*zeta3 1]; zeta5=0.5; den5=[1 2*zeta5 1]; zeta7=0.7; den7=[1 2*zeta7 1]; zeta9=1.0; den9=[1 2*zeta9 1]; [y1,x,t]=step(num,den1,t);[y3,x,t]=step(num,den3,t);[y5,x,t]=step(num,den5,t);[y7,x,t]=step(num,den7,t);[y9,x,t]=step(num,den9,t);plot(t,y1,t,y3,t,y5,t,y7,t,y9); grid on3)稳定性分析den=[1 1 2 24];roots(den)4)动态性能分析t=0:0.01:2;num=[1000];den=[1 34.5 1000];[y,x,t]=step(num,den,t);plot(t,y);%求超调量maxy=max(y);yss=y(length(t));pos=100*(maxy-yss)/yss%求峰值时间for i=1:1:201if y(i)==maxy,n=i;endendtp=(n-1)*0.01%求调节时间for i=n:1:201if(y(i)<1.05&y(i)>0.95),m=i;break;endendym=y(18)ts=(m-1)*0.015)稳态误差分析%-----------单位冲击-------t=0:0.1:15;[num1,den1]=cloop([1],[1,1]);[num2,den2]=cloop([1],[1,1,0]); [num3,den3]=cloop([4,1],[1,1,0,0]); y1=impulse(num1,den1,t); y2=impulse(num2,den2,t);y3=impulse(num3,den3,t);subplot(3,1,1);plot(t,y1);subplot(3,1,2);plot(t,y2);subplot(3,1,3);plot(t,y3);er1=0-y1(length(t))%0型系统稳态误差er2=0-y2(length(t))%1型系统稳态误差er3=0-y3(length(t))%2型系统稳态误差figure;%-----------单位阶跃-------t=0:0.1:20;[num1,den1]=cloop([1],[1,1]);[num2,den2]=cloop([1],[1,1,0]); [num3,den3]=cloop([4,1],[1,1,0,0]); y1=step(num1,den1,t);y2=step(num2,den2,t);y3=step(num3,den3,t);subplot(3,1,1);plot(t,y1);subplot(3,1,2);plot(t,y2);subplot(3,1,3);plot(t,y3);er4=0-y1(length(t))%0型系统稳态误差er5=0-y2(length(t))%1型系统稳态误差er6=0-y3(length(t))%2型系统稳态误差figure%-----------单位斜坡-------t=0:0.1:20;t1=0:0.1:20;[num1,den1]=cloop([1],[1,1]);[num2,den2]=cloop([1],[1,1,0]); [num3,den3]=cloop([4,1],[1,1,0,0]); y1=step(num1,[den1 0],t);y2=step(num2,[den2 0],t);y3=step(num3,[den3 0],t);subplot(3,1,1);plot(t1,y1,t1,t1); subplot(3,1,2);plot(t,y2,t,t); subplot(3,1,3);plot(t,y3,t,t);er7=t1(length(t1))-y1(length(t))%0型系统稳态误差er8=t(length(t))-y2(length(t))%1型系统稳态误差er9=t(length(t))-y3(length(t))%2型系统稳态误差6)实例分析:kp=[0.11 6];t=[0:0.01:1];num1=303.03*kp(1);den1=[0.00001 0.00633 0.20167 21.21*kp(1)+1]; y1=step(num1,den1,t);num2=303.03*kp(2);den2=[0.00001 0.00633 0.20167 21.21*kp(2)+1]; y2=step(num2,den2,t);subplot(211),plot(t,y1);subplot(212);plot(t,y2);gtext('kp=0.11');gtext('kp=6');。

自动控制原理2 实验报告

自动控制原理2 实验报告

中国石油大学(北京)实验报告实验课程:自动控制原理2实验名称:采样控制系统分析班级:学号: 姓名:实验台号:成绩:实验日期:年月日实验1采样控制系统一、实验目的考察连续时间系统的采样控制中,零阶保持器的作用与采样时间间隔Ts对系统稳定性的影响。

二、实验步骤1、典型单位负反馈连续时间系统的开环传递函数为G(s)=K/(s2+s),借助于Matlab 仿真,并分析并验证K对系统性能的影响。

步骤:Matlab相关命令:Gs=tf([1],[1 1 0]) ;pzmap(Gs);figure(1)rlocus(Gs);K值变化时的阶跃相应曲线for k=[0,0.01,0.05,0.10,0.15,0.20,0.25]num=[k];den=[1,1,0]Gs=tf(num,den);figure(1)margin(Gs);figure(2)t=0:0.001:500;step(Gs,t);grid;hold onend2、将上述连续系统离散化,成为带零阶保持器的采样系统。

借助于Matlab仿真,调整采样周期T 和增益K 的大小,观察T 和K 对系统稳定性和调节性能的影响。

调整系数,给出[1]p384-385习题7-24和7-26的答案。

实验步骤:(1) 确定有零阶保持器的开环系统脉冲传递函数G(z)。

))(1()1()(T T e z z z e K z G -----=Matlab 相关命令:for k=[0,0.01,0.05,0.10,0.15,0.20,0.25]num=[k*0.1,0];den=[1,-1.9,0.9];G1=tf(num,den);G=tf2zp(num,den);Gd=c2d(G,0.1,’zoh ’);G0=feedback(Gd,a);t=0:0.1:50;u=1;tsim(G0,u,t,0);gridfor k=[0,0.01,0.05,0.10,0.15,0.20,0.25]G=tf([5],[1 1 0]);Gd=c2d(G,0.1,'zoh');G0=feedback(Gd,1);t=0:0.1:50;step(G0,t); gridxlabel('t');ylable('c(t)');title(‘ramp response ’)hold onend当T=0.1,0.5,1,2时分别重复上面的命令习题7-247-24(1)求出脉冲传递函数:程序代码:rlocus(G)G0=tf([1],[1 10 0 ]);G=c2d(G0,0.1,'zoh')G =0.003679 z + 0.002642----------------------z^2 - 1.368 z + 0.3679Sample time: 0.1 secondsDiscrete-time transfer function.(2)求闭环系统的z特征方程feedback(G,1)ans =0.003679 z + 0.002642----------------------z^2 - 1.364 z + 0.3705Sample time: 0.1 secondsDiscrete-time transfer function.(3)计算使系统稳定的K的最大值rlocus(G)(4)K=78(5)求闭环脉冲传递函数并绘出单位阶跃响应曲线程序代码:G0=tf([78],[1 10 0 ]);G=c2d(G0,0.1,'zoh')Gd= feedback(G,1);t=0:0.1:6;step(Gd,t)Gd =0.2869 z + 0.2061---------------------z^2 - 1.081 z + 0.574Sample time: 0.1 seconds Discrete-time transfer function. 阶跃响应曲线:(6)系统闭环极点以及超调量程序代码:G0=tf([120],[1 10 0 ]);G=c2d(G0,0.1,'zoh');Gd=feedback(G,1);t=0:0.1:6;step(Gd,t)Transfer function:0.4415 z + 0.3171----------------------z^2 - 0.9264 z + 0.685 Sampling time: 0.1b = [0.4415 0.3171];a = [1 -0.9264 0.685]; [b,a] = eqtflength(b,a); [z,p,k] = tf2zp(b,a)z =-0.7182p =0.4632 + 0.6859i0.4632 - 0.6859i k =0.4415超调量为53.8%. (7) t=0:0.1:6;step(Gd,t)7-267-26.程序代码:G0=tf([1],[1 1 0]);G=c2d(G0,0.2,'zoh');Gd=feedback(G,1);t=0:0.2:20;step(Gd,t)hold onG0=tf([1],[1 1 0]);G=c2d(G0,0.4,'zoh');Gd=feedback(G,1);t=0:0.4:20;step(Gd,t)hold onG0=tf([1],[1 1 0]);G=c2d(G0,0.6,'zoh');Gd=feedback(G,1);t=0:0.6:25;step(Gd,t)hold onG0=tf([1],[1 1 0]);G=c2d(G0,0.8,'zoh');Gd=feedback(G,1);t=0:0.8:30;step(Gd,t)hold onG0=tf([1],[1 1 0]);G=c2d(G0,1.0,'zoh');Gd=feedback(G,1);t=0:1.0:30;step(Gd,t)hold onG0=tf([1],[1 1 0]);G=c2d(G0,1.2,'zoh');Gd=feedback(G,1);t=0:1.2:30;step(Gd,t)hold on实验图形记录:(1)T=0.2s%21%;8.38s T σ==(2)T=0.4s%26%;8.53s T σ==(3)T=0.6s%31%;11.4s T σ==(4)T=0.8ss(5)T=1.0s(6)%40%;15.3s T σ==(7)T=1.2ssT 从0.2s 到1.2s3、计算机控制系统如图5-7所示,采样周期T=0.1s ,试分析不同的PID 调节器及不同参数对系统性能的影响,并分析各种情况下PID 参数的选择方法。

自控原理实验报告

自控原理实验报告

一、实验目的1. 理解并掌握自动控制原理的基本概念和基本分析方法。

2. 掌握典型环节的数学模型及其在控制系统中的应用。

3. 熟悉控制系统的时间响应和频率响应分析方法。

4. 培养实验操作技能和数据处理能力。

二、实验原理自动控制原理是研究控制系统动态性能和稳定性的一门学科。

本实验主要涉及以下几个方面:1. 典型环节:比例环节、积分环节、微分环节、惯性环节等。

2. 控制系统:开环控制系统和闭环控制系统。

3. 时间响应:阶跃响应、斜坡响应、正弦响应等。

4. 频率响应:幅频特性、相频特性等。

三、实验内容1. 典型环节的阶跃响应- 比例环节- 积分环节- 比例积分环节- 比例微分环节- 比例积分微分环节2. 典型环节的频率响应- 幅频特性- 相频特性3. 二阶系统的阶跃响应- 上升时间- 调节时间- 超调量- 峰值时间4. 线性系统的稳态误差分析- 偶然误差- 稳态误差四、实验步骤1. 典型环节的阶跃响应- 搭建比例环节、积分环节、比例积分环节、比例微分环节、比例积分微分环节的实验电路。

- 使用示波器观察并记录各个环节的阶跃响应曲线。

- 分析并比较各个环节的阶跃响应曲线,得出结论。

2. 典型环节的频率响应- 搭建比例环节、积分环节、比例积分环节、比例微分环节、比例积分微分环节的实验电路。

- 使用频率响应分析仪测量各个环节的幅频特性和相频特性。

- 分析并比较各个环节的频率响应特性,得出结论。

3. 二阶系统的阶跃响应- 搭建二阶系统的实验电路。

- 使用示波器观察并记录二阶系统的阶跃响应曲线。

- 计算并分析二阶系统的上升时间、调节时间、超调量、峰值时间等性能指标。

4. 线性系统的稳态误差分析- 搭建线性系统的实验电路。

- 使用示波器观察并记录系统的稳态响应曲线。

- 计算并分析系统的稳态误差。

五、实验数据记录与分析1. 典型环节的阶跃响应- 比例环节:K=1,阶跃响应曲线如图1所示。

- 积分环节:K=1,阶跃响应曲线如图2所示。

自控原理课程实验报告

自控原理课程实验报告

一、实验目的1. 理解并掌握自动控制原理的基本概念和基本分析方法。

2. 熟悉自动控制系统的典型环节,包括比例环节、积分环节、比例积分环节、惯性环节、比例微分环节和比例积分微分环节。

3. 通过实验,验证自动控制理论在实践中的应用,提高分析问题和解决问题的能力。

二、实验原理自动控制原理是研究自动控制系统动态和稳态性能的学科。

本实验主要围绕以下几个方面展开:1. 典型环节:通过搭建模拟电路,研究典型环节的阶跃响应、频率响应等特性。

2. 系统校正:通过在系统中加入校正环节,改善系统的性能,使其满足设计要求。

3. 系统仿真:利用MATLAB等仿真软件,对自动控制系统进行建模和仿真,分析系统的动态和稳态性能。

三、实验内容1. 典型环节实验(1)比例环节:搭建比例环节模拟电路,观察其阶跃响应,分析比例系数对系统性能的影响。

(2)积分环节:搭建积分环节模拟电路,观察其阶跃响应,分析积分时间常数对系统性能的影响。

(3)比例积分环节:搭建比例积分环节模拟电路,观察其阶跃响应,分析比例系数和积分时间常数对系统性能的影响。

(4)惯性环节:搭建惯性环节模拟电路,观察其阶跃响应,分析时间常数对系统性能的影响。

(5)比例微分环节:搭建比例微分环节模拟电路,观察其阶跃响应,分析比例系数和微分时间常数对系统性能的影响。

(6)比例积分微分环节:搭建比例积分微分环节模拟电路,观察其阶跃响应,分析比例系数、积分时间常数和微分时间常数对系统性能的影响。

2. 系统校正实验(1)串联校正:在系统中加入串联校正环节,改善系统的性能,使其满足设计要求。

(2)反馈校正:在系统中加入反馈校正环节,改善系统的性能,使其满足设计要求。

3. 系统仿真实验(1)利用MATLAB等仿真软件,对自动控制系统进行建模和仿真,分析系统的动态和稳态性能。

(2)根据仿真结果,优化系统参数,提高系统性能。

四、实验步骤1. 搭建模拟电路:根据实验内容,搭建相应的模拟电路,并连接好测试设备。

自动控制原理实验报告

自动控制原理实验报告

学生实验报告PID 控制器是一种线性控制器,它根据给定值()t r 与实际输出值()t y 构成控制偏差()t e()()()t y t r t e -=(2.2.1)将偏差的比例()P 、积分()I 和微分()D 通过线性组合构成控制量,对被控对象进行控制,故称PID 控制器。

其控制规律为()()()()⎥⎦⎤⎢⎣⎡++=⎰dt t de T dt t e T t e K t u D tp 011(2.2.2)或写成传递函数的形式()()()⎪⎪⎭⎫ ⎝⎛++==s T s T K s E s U s G D p 111(2.2.3) 式中:p K ——比例系数;I T ——积分时间常数;D T ——微分时间常数。

在控制系统设计和仿真中,也将传递函数写成()()()sK s K s K s K s K K s E s U s G I p D D Ip ++=++==2(2.2.4) 式中:P K ——比例系数;I K ——积分系数;D K ——微分系数。

上式从根轨迹角度看,相当于给系统增加了一个位于原点的极点和两个位置可变的零点。

简单说来,PID 控制器各校正环节的作用如下:A 、比例环节:成比例地反映控制系统的偏差信号()t e ,偏差一旦产生,控制器立即产生控制作用,以减少偏差。

B 、积分环节:主要用于消除稳态误差,提高系统的型别。

积分作用的强弱取决于积分时间常数I T ,I T 越大,积分作用越弱,反之则越强。

C 、微分环节:反映偏差信号的变化趋势(变化速率),并能在偏差信号值变得太大之前,在系统中引入一个有效的早期修正信号,从而加快系统的动作速度,减小调节时间。

2、 PID 参数的确定方法 (1) 根轨迹法确定PID 参数 PID 的数学模型可化为:()s K s K s K s G IP D ++=2从仿真曲线看出未校正系统震荡不稳定。

设球杆系统PID 校正的结构图为如图2.2.5 示:要求采用凑试法设计PID校正环节,使系统性能指标达到调节时间小于令Kp=2.5,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:令Kp=2,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:令Kp=2.1,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:令Kp=2.4,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:令Kp=2.5,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:令Kp=2.6,Ki=1.2,Kd=1.5,SampleTime=-1,位移响应曲线如下:PID参数整定:Time Offset(s) Kp Ki Kd SampleTime sT(s) %5 2.5 0.9 1.5 -1 23 4%学生实验报告从仿真曲线看出未校正系统震荡不稳定。

自控实验报告实验二

自控实验报告实验二

自控实验报告实验二一、实验目的本次自控实验的目的在于深入理解和掌握控制系统的性能指标以及相关参数对系统性能的影响。

通过实验操作和数据分析,提高我们对自控原理的实际应用能力,培养解决实际问题的思维和方法。

二、实验设备本次实验所使用的设备主要包括:计算机一台、自控实验箱一套、示波器一台、信号发生器一台以及相关的连接导线若干。

三、实验原理在本次实验中,我们主要研究的是典型的控制系统,如一阶系统和二阶系统。

一阶系统的传递函数通常表示为 G(s) = K /(Ts + 1),其中 K 为增益,T 为时间常数。

二阶系统的传递函数则可以表示为 G(s) =ωn² /(s²+2ζωn s +ωn²),其中ωn 为无阻尼自然频率,ζ 为阻尼比。

通过改变系统的参数,如增益、时间常数、阻尼比等,观察系统的输出响应,从而分析系统的稳定性、快速性和准确性等性能指标。

四、实验内容与步骤1、一阶系统的阶跃响应实验按照实验电路图连接好实验设备。

设置不同的时间常数 T 和增益 K,通过信号发生器输入阶跃信号。

使用示波器观察并记录系统的输出响应。

2、二阶系统的阶跃响应实验同样按照电路图连接好设备。

改变阻尼比ζ 和无阻尼自然频率ωn,输入阶跃信号。

用示波器记录输出响应。

五、实验数据记录与分析1、一阶系统当时间常数 T = 1s,增益 K = 1 时,系统的输出响应呈现出一定的上升时间和稳态误差。

随着时间的推移,输出逐渐稳定在一个固定值。

当 T 增大为 2s,K 不变时,上升时间明显变长,系统的响应速度变慢,但稳态误差基本不变。

2、二阶系统当阻尼比ζ = 05,无阻尼自然频率ωn = 1rad/s 时,系统的输出响应呈现出较为平稳的过渡过程,没有明显的超调。

当ζ 减小为 02,ωn 不变时,系统出现了较大的超调,调整时间也相应变长。

通过对实验数据的分析,我们可以得出以下结论:对于一阶系统,时间常数 T 越大,系统的响应速度越慢;增益 K 主要影响系统的稳态误差。

自动控制原理实验报告

自动控制原理实验报告

自动控制原理实验报告本实验为基于微处理器的温度控制系统的设计与实现。

实验目的是通过实践掌握基于微处理器的控制系统设计和实现方法,了解数字信号处理的基本原理和应用。

本报告将分为实验原理,系统设计,实验步骤,实验结果和结论等几个部分进行详细阐述。

一、实验原理数字信号处理的基本原理是将模拟信号经过采样、量化和编码后转换为数字信号,并在数字领域中对其进行处理。

在本实验中,采用的是基于单片机控制的数字温度控制系统。

该系统的设计要求基于以往的温度控制系统,并具备更过的实用价值和工程性能。

系统的基本原理如下:1.数字信号采样该系统通过传感器来采集温度值,并将其转化为数字信号,实现了数字化控制。

系统在稳态时,通过采用PID控制方法来对温度进行控制。

2.温度控制方法对于本实验中开发的系统,采用的是基于PID控制算法的控制方法。

PID即比例积分微分控制算法,它是一种最常用的控制算法,具备响应速度快、稳态误差小等优点。

PID控制算法的主要原理是,通过比例、积分和微分三个控制系数对输出进行调节,使系统的响应速度更快,而且在稳态时误差非常小。

3.系统设计本实验系统的设计通过单片机的程序控制,主要包含三部分:硬件设计、软件设计和温控系统设计。

二、系统设计1.硬件设计本实验采用的是基于AT89S52单片机的数字温度控制系统,其硬件电路主要包括以下模块:(1)单片机控制器:采用AT89S52单片机;(2)温度传感器:采用DS18B20数字温度传感器;(3)电源模块:采用稳压电源,提供系统所需电压。

2.软件设计本实验采用的是基于C语言开发的程序控制系统,该软件具备以下功能模块:(1)数据采集:通过程序控制读取温度传感器数值;(2)控制算法:实现PID控制算法的程序设计;(3)控制输出:将PID算法结果通过程序输出到负载端。

3.温控系统设计本实验设计的数字温度控制系统,其温控系统设计主要包括以下几个方面:(1)温度检测:系统通过DS18B20数字温度传感器检测环境温度。

自动控制原理实验报告二

自动控制原理实验报告二

-500
-400
-1
-300
-200
-100
Real Axis (seconds )
Step Response 1 0.9 0.8 0.7
System: sys Settling time (seconds): 0.174 System: sys Rise time (seconds): 0.00405 System: sys Peak amplitude: 1 Overshoot (%): 0 At time (seconds): 5
2 1 0 -1 -2 -3 -4 -5 -2.5 -2 -1.5 -1 -0.5
-1
0
0.5
Real Axis (seconds )
Root Locus 5 0.707 4 3
System: Gb Gain: 4.32 Pole: -3.16 + 3.16i Damping: 0.707 Overshoot (%): 4.33 Frequency (rad/s): 4.47
5
Real Axis (seconds )
0.03
System: sys Peak amplitude: 0.0288 Overshoot (%): 15.4 Step Response At time (seconds): 0.44
System: sys 0.025 Rise time (seconds): 0.194 System: sys Settling time (seconds): 0.927
1.2
1
Amplitude
0.8
0.6
0.4
0.2
0
0
0.5
1
1.5
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验一 典型环节的模拟研究及阶跃响应分析1、比例环节可知比例环节的传递函数为一个常数:当Kp 分别为0.5,1,2时,输入幅值为1.84的正向阶跃信号,理论上依次输出幅值为0.92,1.84,3.68的反向阶跃信号。

实验中,输出信号依次为幅值为0.94,1.88,3.70的反向阶跃信号, 相对误差分别为1.8%,2.2%,0.2%. 在误差允许范围内可认为实际输出满足理论值。

2、 积分环节积分环节传递函数为:(1)T=0.1(0.033)时,C=1μf (0.33μf ),利用MATLAB ,模拟阶跃信号输入下的输出信号如图: T=0.1 T=0.033与实验测得波形比较可知,实际与理论值较为吻合,理论上T=0.033时的波形斜率近似为T=0.1时的三倍,实际上为8/2.6=3.08,在误差允许范围内可认为满足理论条件。

3、 惯性环节惯性环节传递函数为:if i o R RU U -=TS1CS R 1Z Z U U i i f i 0-=-=-=1TS K)s (R )s (C +-=K = R f /R 1,T = R f C,(1) 保持K = R f /R 1 = 1不变,观测T = 0.1秒,0.01秒(既R 1 = 100K,C = 1μf ,0.1μf )时的输出波形。

利用matlab 仿真得到理论波形如下: T=0.1时 t s (5%)理论值为300ms,实际测得t s =400ms 相对误差为:(400-300)/300=33.3%,读数误差较大。

K 理论值为1,实验值2.12/2.28,相对误差为(2.28-2.12)/2.28=7%与理论值较为接近。

T=0.01时t s (5%)理论值为30ms,实际测得t s =40ms 相对误差为:(40-30)/30=33.3%由于ts 较小,所以读数时误差较大。

K 理论值为1,实验值2.12/2.28,相对误差为(2.28-2.12)/2.28=7%与理论值较为接近(2) 保持T = R f C = 0.1s 不变,分别观测K = 1,2时的输出波形。

K=1时波形即为(1)中T0.1时波形K=2时,利用matlab 仿真得到如下结果: t s(5%)理论值为300ms,实际测得t s =400ms 相对误差为:(400-300)/300=33.3% 读数误差较大K 理论值为2,实验值4.30/2.28,相对误差为(2-4.30/2.28)/2=5.7%与理论值较为接近。

4、 二阶振荡环节令R 3 = R 1,C 2 = C 11KTSS T 1)s (R )s (C 22++=T = R 1C 1,K = R 2/R 1n ω= 1/T = 1/R 1C 1ξ= 1/2K = R 1/2R 2(1) 取R 1 = R 3 = 100K,C 1 = C 2 = 1μf 既令T = 0.1秒,调节R 2分别置阻尼 比ξ= 0.1,0.5,1○1R2=500k,ξ=0.1时, n ω=10;matlab 仿真结果如下:超调量M p 理论值为e^(-ξ*π/(1-ξ^2)^0.5)=73%,实验值为(3.8-2.28)/2.28=66.7%与理论值较为接近.过渡过程时间理论值(计算时的估计公式)t s =4/(ξ*n ω)=4s ,由matlab 仿真得t s =2.89s ,实验值为3.1s,与仿真得到的理论值相对误差为(3.1-2.89)/2.89=7.2%较为接近。

○2R2=100k, ξ=0.5,n ω=10 ;matlab 仿真结果如下: 超调量M p 理论值为e^(-ξ*π/(1-ξ^2)^0.5)=16%,实验值为(2.8-2.28)/2.28=22.8%与理论值较为接近过渡过程时间理论值(计算时的估计公式)t s =4/(ξ*n ω)=0.8s ,由matlab 仿真得t s =0.525s ,实验值为0.59,与仿真得到的理论值相对误差为(0.59-0.525)/0.525=12.4%较为接近。

○3 R2=50k, ξ=1,n ω=10;matlab 仿真结果如下:超调量M p 理论值为0,实验值为(2.28-2)/2.28=12.3%,与理论值吻合。

过渡过程时间理论值,由matlab 仿真得t s =0.48s ,实验值为0.40,与仿真得到的理论值相对误差为(0.48-0.40)/0.48=20%较为接近。

(2)取R 1 = R 3 = 100K,C 1 = C 2 =0.1μf 既令T = 0.01秒,重复进行上述测试。

○1R2=500k,ξ=0.1时, n ω=100;matlab 仿真结果如下: 超调量M p 理论值为e^(-ξ*π/(1-ξ^2)^0.5)=73%,实验值为(3.8-2.28)/2.28=66.7%与理论值较为接近.过渡过程时间理论值(计算时的估计公式)t s =4/(ξ*n ω)=0.4s ,由matlab 仿真得t s =0.29s ,实验值为0.30,与理论值相对误差为(0.30-0.29)/0.29=3.4%较为接近。

○2R2=100k,ξ=0.5时, n ω=100;matlab 仿真结果如下: 超调量M p 理论值为e^(-ξ*π/(1-ξ^2)^0.5)=16%,实验值为(2.8-2.28)/2.28=22.8%与理论值较为接近过渡过程时间理论值(计算时的估计公式)t s =4/(ξ*n ω)=0.08s ,由matlab 仿真得t s =0.0525s ,实验值为0.05,与仿真得到的理论值相对误差为(0.0525-0.05)/0.0525=4.8%较为接近。

○3 R2=50k, ξ=1,n ω=10;matlab 仿真结果如下:超调量M p 理论值为0,实验值为(2.28-2)/2.28=12.3%,与理论值吻合。

过渡过程时间理论值,由matlab 仿真得t s =0.048s ,实验值为0.04,与仿真得到的理论值相对误差为(0.048-0.04)/0.048=16.7%较为接近。

六、思考题1、根据实验结果,分析一阶系统t s 与T,K 之间的关系。

参数T 的物理意义? T 越大,ts 越大,ts 与K 无关。

T 反映了系统的瞬态响应速度。

2、根据实验结果,分析二阶系统t s ,M p ,与n ω,ξ之间的关系。

参数n ω,ξ的物理意义? 超调量只与ξ有关,ξ越小,超调量越大;调节时间与n ω*ξ有关,乘积越大,调节时间越小;n ω*ξ反映了系统阶跃响应的衰减程度,n ω反映了阶跃响应的振荡快慢程度。

3、对于图1-5所示系统,若将其反馈极性改为正反馈;或将其反馈回路断开,这时的阶跃响应应有什么特点?试从理论上进行分析(也可在实验中进行观察)变成正反馈或将其反馈回路断开,理论上阶跃响应的大小不断增加,实际中受制于运放的最大输出电压的影响,阶跃响应快速上升,最后达到一个很大的幅值。

4、根据所学习的电模拟方法,画出开环传递函数为)1S T 2S T )(1S T (K)s (G 22221+ξ++=的单位反馈系统的模拟线路图,并注明线路图中各元件参数(用R 、C 等字符表示)和传递函数中参数的关系。

易知将一个一阶惯性环节与图1-5所示电路串联起来后,再加一个单位反相比例环节即可实现,电路图如下其中应有R3=R1,C2=C1,于是K=Rf/R1,T1=Rf*C,T2=R1*C1,ζ=R1/(2*R2)。

实验二开环零点及闭环零点作用的研究实验电路图见附件(a)选择T=3.14s,K=3.14,T(S)=L(S)/1+L(S)=3.14/3.14S^2+S+3.14利用MATLAB仿真如下Mp:理论值1.6 实际值1.7 相对误差6.25%tp:理论值3.26 实际值 2.9 相对误差11.0%ts:理论值23 实际值 24.2 相对误差5.2%(b) Td=0.033T(S)=L(S)/1+L(S)=1.0362S+3.14/3.14S^2+4.1762S+3.14 利用MATLAB仿真Mp:理论值1.065 实际值1.15 相对误差8.0% tp:理论值3.68 实际值3.6 相对误差2.2%ts:理论值5.77 实际值6.0 相对误差4.0%(c) T(S)=L(S)/1+L(S)=3.14/3.14S^2+4.1762S+3.14利用MATLAB仿真Mp:理论值1.06 实际值1.08 相对误差2.0%tp:理论值4.12 实际值4.3 相对误差4.4%ts:理论值6.09 实际值6.2 相对误差1.8%比较实验二、三,知开环零点加快了瞬态响应;比较实验一、三,知闭环零点改善了整体的闭环性能,其主要原因是改变了阻尼比。

由实验结果可知,增加比例微分环节后系统的瞬态响应改善了,其根本在于增大了阻尼比。

而第二个实验中由于引进了开环零点,所以其性能与第三个不一样。

实验心得及体会提前预习,熟悉电路图,设计好参数对完成实验有很大的帮助,可以起到事半功倍的效果,要养成提前预习的习惯。

思考题为什么说系统的动态性能是由闭环零点,极点共同决定的?从时域和频域的关系来看,极点的位置决定了系统的响应模态,而零点的位置决定了每个模态函数的相对权重。

实验三控制系统稳定性研究一、实验数据本实验的线路图如下,其中R11=R12=R21=R31=100K,1.对于方案一,取R13=R22=1M,C1=1μ,C2=10μ,R3=100K,C3=1μ,由实验现象得知,对任意α∈(0,1),系统均稳定,且α越大,响应速度越快,幅值也越大。

对于方案二,C3=1μ,知对于任意α系统仍稳定,且α越大,响应速度越快,幅值也越大。

方案三中R32=1M,C3=1μ,当输出呈现等幅振荡时,α=0.0192. 对于第一组,由实验可知对任意α∈(0,1)系统均稳定,且α越大,响应速度越快,幅值也越大。

第二组中,当输出呈现等幅振荡时,α=0.5103. 仍选择以上电路,要使T=RC=0.5s,可选取R=500K,C=1μ。

而由以上传a=1时,R13=R22=R32=500K,C1=C2=C3=1μ。

实验测得当输出开始呈现缓慢衰减,K=809.1Hz。

a=2时,R13=1M,R22=500K,R32=250K,C1=C2=C3=1μ。

实验测得当输出开始呈现缓慢衰减,K=924.1Hz。

a=5时,R13=250K,C1=10μ,R22=500K,C2=1μ,R32=100K,C3=1μ。

此时发现对任意α∈(0,1)系统均稳定。

二、数据处理1.对于前三个方案,由Hurwitz判据易知α=1.22,11.1,0.0242时系统临界稳定。

而实验中α不可能大于1,故前两个实验中系统均稳定,而第三个实验中测得α=0.019,与理论值相对误差为(0.0242-0.019)/0.0242=21.4%。

相关文档
最新文档