一、二阶系统频率特性测试与分析

1、【实验目的】掌握测量典型一阶系统和二阶系统的频率特性曲线的方法；掌握软件仿真求取一、二阶系统的开环频率特性的方法；学会用Nyquist判据判定系统的稳定性。【实验设备与软件】labACT实验台与虚拟示波器MATLAB软件【实验原理】系统的频率特性测试方法对于现行定常系统，当输入端加入一个正弦信号X(t)=Xsin(t)时，其稳m态输出是一个与输入信号频率相同，但幅值和相位都不同的正弦信号Y(s)=Ysin(t+,)=XG(j)sin(t+,)。mm幅频特性：G(j)=Y/X，即输入与输出信号的幅度比值，通常转换成mm20lgG(j)形式。相频特性：()=argG(j)，可以直接基于虚拟示波器读

2、取，也可以用“李沙育图行”法得到。可以将用Bode图或Nyquist图表示幅频特性和相频特。在labACT试验台采用的测试结构图如下:定系统对于实验就是有源放大电路模拟的一、二阶稳定系统。系统的频率测试硬件原理1)正弦信号源的产生方法频率特性测试时，一系列不同频率输入正弦信号可以通过下图示的原理产生。按照某种频率不断变化的数字信号输入到DAC0832,转换成模拟信号，经一级运放将其转换为模拟电压信号，再经过一个运放就可以实现双极性电压输出。根据数模转换原理，知0128R再根据反相加法器运算方法，得(1)V=0R1V+R2V(R01RR,222xN28N128V128R(2)由表达式可以看出输出

3、时双极性的：当N大于128时，输出为正；反之则为负；当输入为128时，输出为0.在labACT实验箱上使用的参考电压时5V的，内部程序可以产生频率范围是对一阶系统是0.5H64H、对二阶系统是0.5H16H的信号，并由B2单元ZZZZ的OUT2输出。2）被测对象输出信号的采样方法对被测对象的输出信号夏阳，首先将其通过LM324与基准电压进行比较嵌位，再通过CD14538进行脉冲整形，一保证有足够的IRQ采样时间，最后将信号送到处理器的IRQ6脚，向处理器申请中断，在中断中对模拟量V进行采样并模y数转换，进而进行处理与计算幅值与相位。途中采用ADC089采集模拟量，以单极性方式使用，所以在出现振

4、荡的情况下需要加入一个二极管，将V出现负值y时将其直接拉倒0。3)实验对象的描述与计算(1)水箱液位对象模型水箱液位对象模型可以抽象成一个一阶惯性环节G(s)=K,K0,T0，这Ts+1里可以选取不同的模拟不同水箱的情况。它的频率特性图可以用Nyquist图和Bode图只管的表示。(2)直流电机空载对象模型直流电机空载对象在忽略粘性摩擦时的模型开环出传递函数为G()，Ts(Ts+1)01。设电机时间常数T，1，电磁时间常数T，0.1，放大倍数K，25,得开环传递函数为125G(s)，叭2nTs(Ts+1)s(0.1s+1)s(s+2匚)其中01n，自然频率，K/TT，25/0.1，15.81r

5、ad/s阻尼比n101匚，0.5T/KT，0.316。011系统的开环频率特性为ZL(i)，一90一arctan(/10)由此可依据Nyquist曲线画法与Bode图画法得到相应图形。令ZL(i)，1，则得幅值交越频率，1+4匚2-2匚2,14.186cn将其代入相角表达式，得到相角裕度180+()arctan=34.93mc-2匚2+1+4匚4令ZL（i）=-180。，则得到相位交越频率=gg表明裕度为g。直流机对象模型传递函数W(S)2S2+2匚+2nn250s2+10s+250由闭环传递函数可计算谐振频率和谐振峰值L()=20lgr=4.44dB=1，2g2=14.14rad/srn测试

6、数据观察用labACT软件自选择自动控制菜单下的线性控制系统的频率响应分析-实验项目，分别选择一阶或二阶系统，再选择开始实验就会弹出虚拟示波器的频率特性界面，电机开始，试验机将自动产生幅值为1，频率为0.5H64HZZ（一阶）、0.5H16H（二阶）的多个频率信号，测试被测系统的频率特性，ZZ等待奖金十分钟，测试结束。测试结束后，可点击界面下方的“频率特性”选择框的任一项进行切换，将显示被测系统的对数幅频、相频特性曲线，同时在界面上方将显示该系统用户取频率点的L、Im、Re等相关参数。【实验内容】搭建单容水箱液位对象模拟模型（一阶惯性环节），这里的K取1,T任选，但由于设备本身的限定，惯性环节

7、开环增益不要大于1，并且转折点频率应在0.5之后，所以要根据试验箱上的资源选择合适的电阻和电容值。要求画出原理图（需要经过Multtisim软件测试通过）后在labACT实验箱上做实验,记录其开环频率特性曲线（Bode图）；选择不同的频率测试点，填写表1编制程序将实测的数据与MATLAB计算图画在一幅图中进行比较。表1惯性环节测试数据与理论计算数据表(rad/s)3.146.2812.5625.1250.24100.48200.96401.92f(Hz)测0.51.02.04.08.016.032.064.0-0.35-1.4-4.24-8.72-14.2-20.09-26.11-30.542

8、0lgU0量(dB)理-0.41-1.44-4.11-8.68-14.9-20.08-26.07-32.09论测-19-32-52-69-80-83-84-87量p()理-17.4-32.1-51.5-68.3-80-84.32-87.15-88.6论阶系统：搭建直流电机空载对象模拟模型（二阶环节的闭环形式），绘制其闭环频率特性曲线（Bode图）和开环频率特性曲线（Bode图）；选择不通过的频率测试点，填写下表2和表3.需要注意，测试二阶系统的开环频率特性曲线也要在闭环的状态下测试，然后再反求，这一工作实验平台已做好，可以在示波器界面左上角的红色开环或闭环字上双击完成切换显示。编制程序将实测数

9、据线与MATLAB计算图画在一幅图中。表2二阶闭环系统测试数据与理论计算数据表(rad/s)3.146.2812.5618.825.150.362.875.4f(Hz)0.51.02.03.04.08.010.012测CC1TT0.110.973.961.33-5.36-20.77-26.19-31.920lgU0量(dB)理0.271.124.071.33-5.20-19.41-23.51-26.83论（）测量-8-16.6-54-119.3-148-173.8-180-180理论7.4516.6-54.09-118.82-146.5-167.56-170.4-172.1表3二阶开环系统测试

10、数据与理论计算数据表(rad/s)3.146.2812.5618.825.150.362.875.4f(Hz)0.51.02.03.04.08.010.01220lgU0测量17.619.872.13-4.00-9.01-20.31-24.91-29.32(dB)理论17.6110.561.88-3.83-8.66-20.08-24.10-27.2（）测量-99-106-141-151.2-159-178.1-179-180理论-107.4-122.1-141.5-152.0-158.3-168.8-171.0-172.4二阶系统：用MATLAB软仿真求取一、二阶系统开环幅相Nyquist频率

11、特性曲图。1、num=1;den=0.11;kp=tf(num,den);margin(kp);grid;_QP-opnluBGsja_1iGau1u21cFrequency(radZsBodeDaaram-rajuencyi.rao/s):3.izPhase(d$q:-17.32、num=25;den=0.110;kp=tf(num,den);margin(kp1);gridmp)mpnuuBmplBadeDiagramGm=hfdBfatInf,Pm=34.Sdeg(at14.3rad/sr&ep)eeEMJFreeliency(rad.3、num1=25;den1=0.1125;kp1=tf(num1,den1);margin(kp1);gridBodsDiaoram-fflp】mpnHLCCsmap】wwELId【实