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

东南大学自动化学院《自动控制原理实验》实验报告实验四系统频率特性的测试姓名：学号：专业：自动化实验室：组别：同组人员：设计时间：2014年10月28日评定成绩：审阅教师：一、实验目的（1）明确测量幅频和相频特性曲线的意义。（2）掌握幅频曲线和相频特性曲线的测量方法。（3）利用幅频曲线求出系统的传递函数。二、预习与回答（1）实验时，如何确定正弦信号的幅值？幅度太大会出现什么问题，幅度过小又会出现什么问题？答：根据实验参数，计算正弦信号幅值大致的范围，然后进行调节，具体确定调节幅值时，首先要保证输入波形不失真，同时，要保证在频率较大时输出信号衰减后人能够测量出来。如果幅度过大，波形超出线性变化区域，产生失真；如果波形过小，后续测量值过小，无法精确的测量。（2）当系统参数未知时，如何确定正弦信号源的频率？答：从理论推导的角度看，应该采取逐点法进行描述，即ω从0变化到∞，得到变化时幅度和相位的值。从实际操作来看，ω值过小所取得的值无意义，因此我们选取[1.0,100.0]的范围进行测量。（3）先对本系统进行机理建模，求出开环传递函数。答：)()()()()(4321sGsGsGsGsG令R=100K，则122-)(1RRsGsRSRCRsG1.011)1/()(12sRSRCRsG047.011)1/()(23sRSRCRsG02.0112)21/(2)(34所以：)047.01)(02.01)(1.01(1)()()()()(4321ssssGsGsGsGsG三、实验原理在设计控制系统时，首先要建立系统的数学模型，而建立系统的数学模型是控制系统设计的前提和难点。建模一般有机理建模和辨识建模两种方法。机理建模就是根据系统的物理关系式，推导出系统的数学模型。辨识建模主要是人工或计算机通过实验来建立系统数学模型。两种方法在实际的控制系统设计中，常常是互补运用的。辨识建模又有多种方法。本实验采用开环频率特性测试方法，确定系统传递函数，俗称频域法。还有时域法等。准确的系统建模是很困难的，要用反复多次，模型还不一定建准。模型只取主要部分，而不是全部参数。另外，利用系统的频率特性可用来分析和设计控制系统，用Bode图设计控制系统就是Y1XorY2100K100K200K200K100K100K200K200K0.47μF0.1μF1μF正弦信号源虚拟示波器AD1AD2-+-+-+-+其中一种。幅频特性就是输出幅度随频率的变化与输入幅度之比，即)()(ioUUA，测幅频特性时，改变正弦信号源的频率测出输入信号的幅值或峰峰值和输输出信号的幅值或峰峰值测相频有两种方法：（1）双踪信号比较法：将正弦信号接系统输入端，同时用双踪示波器的Y1和Y2测量系统的输入端和输出端两个正弦波，示波器触发正确的话，可看到两个不同相位的正弦波，测出波形的周期T和相位差Δt，则相位差0360Tt。这种方法直观，容易理解。就模拟示波器而言，这种方法用于高频信号测量比较合适。（2）李沙育图形法：将系统输入端的正弦信号接示波器的X轴输入，将系统输出端的正弦信号接示波器的Y轴输入，两个正弦波将合成一个椭圆。通过椭圆的切、割比值；椭圆所在的象限；椭圆轨迹的旋转方向三个要素来决定相位差。就模拟示波器而言，这种方法用于低频信号测量比较合适。若用数字示波器或虚拟示波器，建议用双踪信号比较法。利用幅频和相频的实验数据可以作出系统的波Bode图和Nyquist图四、实验设备THBDC-1实验平台、THBDC-1虚拟示波器五、实验线路图六、实验步骤（1）如图接线，用实验台上的U7、U9、U11、U13单元，信号源的输入接“数据采集接口”AD1(兰色波形)，系统输出接“数据采集接口”AD2(红色波形)。（2）信号源选“正弦波”，幅度、频率根据实际线路图自定，要预习。（3）点击屏上THBDC-1示波器图标，直接点击“确定”，进入虚拟示波器界面，点“示波器（E）”菜单，选中“幅值自动”和“时基自动”。在“通道选择”下拉菜单中选“通道（1-2）”，“采样频率”调至“1”。点“开始采集”后，虚拟示波器可看到正弦波，再点“停止采集”，波形将被锁住，利用示波器“双十跟踪”可准确读出波形的幅度。改变信号源的频率，分别读出系统输入和输出的峰峰值，填入幅频数据表中。f=0.16时要耐心。（4）测出双踪不同频率下的Δt和T填相频数据表，利用公式0360Tt算出相位差。七、实验数据（1）数据表格频率f(Hz)0.160.320.641.111.592.393.184.786.3711.115.9ω1.02.04.07.010.015.020.030.040.070.0100.02imU4.01514.02084.01164.01164.01164.01164.01164.01164.01164.01164.01162omU3.98903.91173.62983.05452.46941.68601.17020.60000.33720.10910.044620Lg-0.0566-0.2389-0.8687-2.3675-4.2145-7.6582-10.7012-16.5033-21.5086-31.3099-39.0797Δt0.18430.16770.16200.15430.14660.12340.11320.08970.07640.05150.0382T6.31293.12721.56200.90250.62870.41650.31650.20930.15720.08970.06320/10.5119.3137.3461.5583.94106.66129.58154.29174.96206.69217.60（2）图当f=1.11Hz（即ω=7.0）时，输入输出波形如上图，其中蓝色为输入信号，红色为输出信号。Δt=0.1543。八、实验分析（1）画出系统的实际幅度频率特性曲线、相位频率特性曲线，并将实际幅度频率特性曲线转换成折线式Bode图，并利用拐点在Bode图上求出系统的传递函数。答：在MATLAB命令行下输入：plot(Lgw,LgU)gridonxlabel('Lgw/(rad/s)');ylabel('20Lg(Uom/Uim)/dB');得到幅频特性图如下：在MATLAB命令行下输入：plot(Lgw,sita)gridonxlabel('Lgw/(rad/s)');ylabel('θ/°');得到相频特性图如下：近似求转折频率如下：用MATLAB拟合工具cftool拟合幅频特性可以得出：Linearmodel:f(x)=a*(sin(x-pi))+b*((x-10)^2)+cCoefficients(with95%confidencebounds):a=-44.99(-46.89,-43.1)b=2.225(2.17,2.28)c=-222.7(-228.4,-217)Goodnessoffit:SSE:1.049R-square:0.9994AdjustedR-square:0.9993RMSE:0.3621可得幅频特性曲线拟合方程为：222.7-10)^2-2.225(x+)-44.99sin(x-=f(x)求其导数，得：10)-4.45(x+)-44.99cos(x-=(x)f’分别令0dB/dec2-=(x)f’0dB/dec4-=(x)f’0dB/dec6-=(x)f’w=1:0.01:10;dlgu=-44.99.*cos(w-pi)+4.45.*(w-10);解得.121=1x.631=2x.162=x3过上述1x、2x、3x三点，分别作0dB/dec2-、0dB/dec4-、60dB/dec-的直线w=0:0.01:2.50;lgu=-44.99.*sin(w-pi)+2.225.*((w-10).^2)-222.7;plot(w,lgu)gridonxlabel('Lgw/(rad/s)');ylabel('20Lg(Uom/Uim)/dB');holdonw1=0.60:0.01:1.60;w2=1.25:0.01:2.15;w3=1.75:0.01:2.50;LgU1=lgu1-20.*(w1-1.12);LgU2=lgu2-40.*(w2-1.63);LgU3=lgu3-60.*(w3-2.16);plot(w1,LgU1,'b')holdonplot(w2,LgU2,'g')holdonplot(w3,LgU3,'r')holdonplot(0.7823,[0:-0.01:-70],'r')holdonplot(1.3821,[0:-0.01:-70],'g')holdonplot(1.8889,[0:-0.01:-70],'b')holdon则该三条直线即为三条斜率分别为0dB/dec2-、0dB/dec4-、60dB/dec-的切线，最终的图如下：代入1x、2x、3x三点可得lgu1=-44.99.*sin(1.12-pi)+2.225.*((1.12-10)^2)-222.7;lgu2=-44.99.*sin(1.63-pi)+2.225.*((1.63-10)^2)-222.7;lgu3=-44.99.*sin(2.16-pi)+2.225.*((2.16-10)^2)-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);可得：7823.0lg11.3821lg21.8889lg3根据上面的Bode图得到的转折频率分别是：srad/6.05761srad/24.10462srad/77.42833将三条渐近线向上平移，第一条渐近线向上平移3.03dB，第二条平移0.97dB则9338.02003.37823.0lg1‘1.3578400.97-1.3821lg21.7131400.97-203.03-1.8889lg3则修正后的转折频率分别是：srad/586.81srad/22.79292srad/6535.153可对应)1)(1)(1(1)(321sTsTsTSG求出：0.11641T0.04392T0.01943T)10194.0)(10439.0)(11164.0(1)(sssSG（2）用文字简洁叙述利用频率特性曲线求取系统传递函数的步骤方法。答：系统传递函数表示形式为：)1)(1)(1()1()(4321sTsTsTsTKSG在对数频率特性曲线上分别画出斜率为40dB/dec、20dB/dec、0dB/dec、-20dB/dec、-40dB/dec、-60dB/dec等的渐近线，平移这些渐近线直至与对数频率特性曲线有切点，找出斜率临近的两条渐近线的交点，即为一个转折频率点。求出相应的时间常数1T，且通过斜率可以判断为惯性环节（在分母上）还是一阶微分环节（在分子上），在确定好各个环节的时间常数后可以确定出常数K。（3）利用上表作出Nyquist图。答：在MATLAB命令行下输入：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);得到奈奎斯特图如下：在MATLAB命令行下输入：Gw=1./((1+0.1.*w.*i).*(1+0.047.*w.*i).*(1+0.02.*w.*i));REL=real(Gw);IMA=imag(Gw);plot(REL,IMA)xlabel('imag');ylabel('real