控制基础第二次试验——系统频率特性地测试

1、实用标准文案东南大学自动化学院自动控制原理实验实验报告实验四 系统频率特性的测试姓 名：学 号：专 业： 自动化 实验室：组 别：同组人员：设计时间：2014年10月28日评定成绩：审阅教师：精彩文档实用标准文案一、实验目的(1)明确测量幅频和相频特性曲线的意义。(2)掌握幅频曲线和相频特性曲线的测量方法。(3)利用幅频曲线求出系统的传递函数。二、预习与回答(1) 实验时，如何确定正弦信号的幅值？幅度太大会出现什么问题，幅度过小又会出现 什么问题？答：根据实验参数，计算正弦信号幅值大致的范围，然后进行调节，具体确定调节幅值时，首先要保证输入波形不失真，同时，要保证在频率较大时输出信号衰减后人能

2、够测量出来。如果幅度过大，波形超出线性变化区域，产生失真；如果波形过小，后续测量值过小， 无法精确的测量。(2) 当系统参数未知时，如何确定正弦信号源的频率？答：从理论推导的角度看，应该采取逐点法进行描述，即3从0变化到8,得到变化时幅度和相位的值。从实际操作来看，3值过小所取得的值无意义，因此我们选取1.0,100.0的范围进行测量。(3) 先对本系统进行机理建模，求出开环传递函数。答：G(s) =G(s)G2(s)G3(s)G4(s)令 R=100K 贝UG1(s) = -2R2RG2(s)=R/(1 RC1S)R11 0.1sG3(s)=R/(1 RC2S)R11 - 0.047sG4(

3、s)=2R/(1 2RC3S)2R11 0.02s所以:G(s)=G1(s)Gz(s)G3(s)G4(s) =1(1 0.1s)(1 0.02s)(1 0.047 s)三、实验原理在设计控制系统时，首先要建立系统的数学模型， 而建立系统的数学模型是控制系统设 计的前提和难点。建模一般有机理建模和辨识建模两种方法。机理建模就是根据系统的物理关系式，推导出系统的数学模型。 辨识建模主要是人工或计算机通过实验来建立系统数学模 型。两种方法在实际的控制系统设计中，常常是互补运用的。 辨识建模又有多种方法。 本实验采用开环频率特性测试方法，确定系统传递函数，俗称频域法。还有时域法等。准确的系统建模是很困

4、难的，要用反复多次，模型还不一定建准。模型只取主要部分，而不是全部参 数。Bode图设计控制系统就是另外，利用系统的频率特性可用来分析和设计控制系统，用精彩文档实用标准文案其中一种。幅频特性就是输出幅度随频率的变化与输入幅度之比，即A(0)=匕)，测幅频特性时，Ui改变正弦信号源的频率测出输入信号的幅值或峰峰值和输输出信号的幅值或峰峰值测相频有两种方法：(1)双踪信号比较法：将正弦信号接系统输入端，同时用双踪示波器的Y1和Y2测量系统的输入端和输出端两个正弦波， 示波器触发正确的话， 可看到两个不同相位的正弦波， 测出波形的周期T和相位差A t,则相位差=生又3600。这种方法直观，容易理解。

5、就模拟示波 ,T器而言，这种方法用于高频信号测量比较合适。(2)李沙育图形法：将系统输入端的正弦信号接示波器的X轴输入，将系统输出端的正弦信号接示波器的丫轴输入，两个正弦波将合成一个椭圆。通过椭圆的切、割比值； 椭圆所在的象限；椭圆轨迹的旋转方向三个要素来决定相位差。就模拟示波器而言，这种方法用于低频信号测量比较合适。若用数字示波器或虚拟示波器，建议用双踪信号比较法。利用幅频和相频的实验数据可以作出系统的波Bode图和Nyquist图四、实验设备THBDC-1实验平台、THBDC-1虚拟示波器五、实验线路图六、实验步骤(1)如图接线，用实验台上的U7、U9 U11、U13单元，信号源的输入接“

6、数据采集接口”AD1(兰色波形)，系统输出接“数据采集接口”AD2(红色波形)。(2)信号源选“正弦波”，幅度、频率根据实际线路图自定，要预习。(3)点击屏上THBDC-1示波器图标，直接点击“确定”，进入虚拟示波器界面，点“示波器(E) ”菜单，选中“幅值自动”和“时基自动”。在“通道选择”下拉菜单中选“通道(1-2) ”，“采样频率”调至“ 1”。点“开始采集”后，虚拟示波器可看到正弦波，再点“停止采集”，波形将被锁住，利用示波器“双十跟踪”可准确读出波形的幅度。改变信号源的频率，分别读出系统输入和输出的峰峰值，填入幅频数据表中。f=0.16时要耐心。(4)测出双踪不同频率下的A t和T填

7、相频数据表，利用公式 =生父360°算出相位差。T精彩文档实用标准文案七、实验数据（1）数据表格频率f(Hz)0.160.320.641.111.592.393.184.786.3711.115.9W1.02.04.07.010.015.020.030.040.070.0100.02Uim4.01514.02084.01164.01164.01164.01164.01164.01164.01164.01164.01162Uom3.98903.91173.62983.05452.46941.68601.17020.60000.33720.10910.044620Lg 2Uim-0.05

8、66-0.2389-0.8687-2.3675-4.2145-7.6582-10.7012-16.5033-21.5086-31.3099-39.0797At0.18430.16770.16200.15430.14660.12340.11320.08970.07640.05150.0382T6.31293.12721.56200.90250.62870.41650.31650.20930.15720.08970.06326/010.5119.3137.3461.5583.94106.66129.58154.29174.96206.69217.6042 1意统言行正itram 屯近讨赢 frh&

9、lt;M 充e. （器*M占上V.的国i*ii丽i/ii-nwi当f=1.11Hz （即3 =7.0）时，输入输出波形如上图，其中蓝色为输入信号，红色为输精彩文档实用标准文案出信号。A t= 0.1543。八、实验分析(1)画出系统的实际幅度频率特性曲线、相位频率特性曲线，并将实际幅度频率特性曲线 转换成折线式Bode图，并利用拐点在 Bode图上求出系统的传递函数。答：在MATLA腌令行下输入：>> plot(Lgw,LgU)> > grid on> > xlabel('Lgw/(rad/s)');ylabel('20Lg(Uom/

10、Uim) / dB');得到幅频特性图如下:-55 0 5 0 -1-22-3 HPFT5至口 nmlR-400020.40.60,811 21.41.61.3Lgw/(rad/s)在MATLA瑜令行下输入：> > plot(Lgw,sita)> > grid on> > xlabel('Lgw/(rad/s)');ylabel('0 / ° ');得到相频特性图如下：精彩文档实用标准文案Lg/(rad/s)近似求转折频率如下：用MATLAB以合工具cftool拟合幅频特性可以得出:Linear model:

11、f(x) = a*(sin(x-pi) + b*(x-10)A2) + cCoefficients (with 95% confidence bounds):a =-44.99(-46.89,-43.1)b =2.225(2.17, 2.28)c =-222.7(-228.4,-217)Goodness of fit:SSE: 1.049R-square: 0.9994Adjusted R-square: 0.9993RMSE: 0.3621精彩文档实用标准文案Lgw可得幅频特性曲线拟合方程为：f(x) =-44.99sin(x -二)+ 2.225(x-10)A2 -222.7求其导数，得：

12、.，f (x) = - 44.99cos(x-二)+ 4.45(x -10)分别令一，f (x) =-20dB/dec一一 .f (x) = -40dB/dec.，f (x) =-60dB/dec>> w = 1:0.01:10;>> dlgu = -44.99.*cos(w-pi)+4.45.*(w-10);Six1 = 1.12x = 1.63x3 =2.16过上述 x1、x2、x3三点，分别作-20dB/dec、-40dB/dec、-60dB/dec 的直线> > w = 0:0.01:2.50;> > Igu = -44.99.*sin(

13、w-pi)+2.225.*(w-10).A2)-222.7;> > plot(w,lgu)> > grid on> > xlabel('Lgw/(rad/s)');ylabel('20Lg(Uom/Uim) / dB');> > hold on> > w1 = 0.60:0.01:1.60;> > w2 = 1.25:0.01:2.15;精彩文档实用标准文案> > w3 = 1.75:0.01:2.50;> > LgU1 = Igu1 - 20.*(w1 - 1.12

14、);> > LgU2 = Igu2 - 40.*(w2 - 1.63);> > LgU3 = Igu3 - 60.*(w3 - 2.16);> > pIot(w1,LgU1,'b')> > hoId on> > pIot(w2,LgU2,'g')> > hoId on> > pIot(w3,LgU3,'r')> > hoId on>> pIot(0.7823,0:-0.01:-70,'r')> > hoId on

15、>> pIot(1.3821,0:-0.01:-70,'g')> > hoId on>> pIot(1.8889,0:-0.01:-70,'b')> > hoId on贝U该三条直线即为三条斜率分另为-20dB/dec、 -40dB/dec、 -60dB/dec的切线，最终的图如下：D11522.5Lgwradi's)代入Xi、X2、X3三点可得> > lgu1 = -44.99.*sin(1.12-pi)+2.225.*(1.12-10)A2)-222.7;> > lgu2 = -

16、44.99.*sin(1.63-pi)+2.225.*(1.63-10)A2)-222.7;> > lgu3 = -44.99.*sin(2.16-pi)+2.225.*(2.16-10)A2)-222.7;> > lgw1 = 1.12 - lgu1 / (-20);> > lgw2 = (-lgu1 - 1.12*20 + lgu2 + 1.63*40) / (20);> > lgw3 = (-lgu2 - 1.63*40 + lgu3 + 2.16*60) / (20);可得：lg 1 =0.7823 Ig 2 =1.3821 Ig 3 =

17、1.8889精彩文档实用标准文案根据上面的Bode图得到的转折频率分别是：1 =6.0576rad/s.2 =24.1046rad /s3 =77.4283rad /s将三条渐近线向上平移，第一条渐近线向上平移3.03dB,第二条平移0.97dB则lg 1 =0.7823 3.03 ： 20 = 0.9338lg 2 =1.3821-0.97-40 = 1.3578lg 3 =1.8889-3.03 ： 20-0.9740 =1.7131则修正后的转折频率分别是：1 =8.586rad /s2 -22.792Srad /s3 = 51.6535rad / s1可对应G(S)=求出：(T1s 1

18、)(T2s 1)(T3s 1)T1 =0.1164T2 -0.0439T3 =0.0194G(S)=1(0.1164s 1)(0.0439s 1)(0.0194s 1)（2）用文字简洁叙述利用频率特性曲线求取系统传递函数的步骤方法。 答：系统传递函数表示形式为：G(S)=K(T1s 1)(T2s 1)(T3s 1)(T4s 1)在对数频率特性曲线上分别画出斜率为40dB/dec、20dB/dec、0dB/dec、-20dB/dec、-40dB/dec、-60dB/dec等的渐近线，平移这些渐近线直至与对数频率特性曲线有切点，找1出斜率临近的两条渐近线的交点，即为一个转折频率点。求出相应的时间常

19、数t='，且通过斜率可以判断为惯性环节（在分母上）还是一阶微分环节（在分子上），在确定好各个精彩文档实用标准文案环节的时间常数后可以确定出常数Ko(3)利用上表作出 Nyquist图。答：在MATLA腌令行下输入:Num = 0,0,0,1;D1 = 0.1,1;D2 = 0.047,1;D3 = 0.02,1;D4 = conv(D1,D2);Den = conv(D3,D4);nyquist(Num,Den);得到奈奎斯特图如下:R Ml在MATLA瑜令行下输入：> > Gw = 1./(1+0.1.*w.*i).*(1+0.047.*w.*i).*(1+0.02.*w.*i);> &g