单闭环磁盘驱动读取系统的校正器设计

1、目录1设计目的与意义12设计任务和要求12.1 初始条件1.2.2 设计要求2.3原系统性能分析33.1 校正前系统性能指标分析33.2 校正前系统稳定性能分析54设计校正系统64.1 校正前系统的频域性能指标64.2 校正系统方案选择74.3 校正步骤74.4 校正前后比较11心得体会18参考文献19武汉理工大学自动控制原理课程设计说明书单闭环磁盘驱动读取系统的校正器设计1设计目的与意义1、通过本次课程设计，能够培养自己理论联系实际的设计思想，训练自己综合运用经典控制理论和相关课程知识的能力。2、掌握自动控制原理中各种校正装置的作用及用法，根据不同的系统性能指标要求进行合理的系统设计，并调试

2、满足系统的指标。3、进一步加深对MATLA歆件的认识和理解，学会使用MATLA引言来进行系统建模、系统校正中的性能参数的求解、以及系统仿真与调试。4、止匕外，还可以通过本次设计来锻炼独立思考和动手解决控制系统实际问题的能力。2设计任务和要求2.1 初始条件磁盘驱动系统的物理结构如图1所示，结构框图如图2所示。图2磁盘驱动系统的结构框图图2中，Amplifier模块代表控制器，此处用比例控制器Ka=80；coil模块代表执武汉理工大学自动控制原理课程设计说明书行机构，此处为电机，该模块的传递函数简化为GKs)=5000,load环节代表荷载，该s1000模块的传递函数为G2(s)=1。s(s20

3、)根据对原系统的分析，选择适当的校正装置，设计满足给定性能指标的补偿环节。2.2 设计要求(1)分析原系统的稳定性和性能指标；(2)设计一个补偿环节，使系统满足Kv=20,相角裕度PM增加5-10二的性能指标;(3)画出系统在(1)校正前后的奈奎斯特曲线和波特图；(4)用Matlab画出上述每种情况的斜坡响应曲线；(5)用Matlab画出校正前后系统的根轨迹；(6)对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析计算的过程，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。武汉理工大学自动控制原理课程设计说明书3原系统性能分析3.1 校正前系统性能指

4、标分析由初始条件可得:其开环传递函数为:此时可得闭环传递函数为：G(s)=20sss(1)(100020K=201)20sJs-+1；+/i+2011000人20)在单位阶跃响应中，由单位阶跃图可得上升时间tr-0.12s超调时间t0.182sp调节时间ts=0.496s超调量：.,%=1.171100%=17%1在matlab中仿真单位阶跃响应程序如下：num=400000;den=1102020000400000;sys=tf(num,den);%建立闭环传递函数模型p=roots(den);%计算系统特征根判断系统稳定性t=0:0.01:5;%设定仿真时间为5sfigure(1)step

5、(sys,t);grid%t取系统单位阶跃响应武汉理工大学自动控制原理课程设计说明书xlabel('t');ylabel('c(t)');title('stepresponse');单位阶跃响应曲线图如下:图3单位阶跃响应图在matlab中的单位斜坡响应程序如下：figure(2)u=t;%定义输入为斜坡信号lsim(sys,u,t,0);grid%求取系统的单位斜坡响应xlabel('t');ylabel('c(t)');title('rampresponse');其斜坡响应图如下：武汉理工大学自

6、动控制原理课程设计说明书图4单位斜坡响应曲线3.2 校正前系统稳定性能分析由程序运算可得，该闭环传递函数的特征根如下：p=1.0e+03*- 1.0004+0.0000i- 0.0098+0.0174i- 0.0098-0.0174i所以可得，由于系统的特征根都具有负实部或者说都位于复平面的左半部，故此时的系统是一一个稳定系统。武汉理工大学自动控制原理课程设计说明书4设计校正系统4.1校正前系统的频域性能指标由系统给定的开环传递函数可得:其原函数的一阶斜坡信号输入得到的速度误差系数为kv=limsG(s)H(s)s0二20400000lims-s,0s(s1000)(s20)同时，其原函数的相

7、角裕度为1=1800+/G(j0)校正前的伯德图在matlab中的程序如下：G=tf(400000,1,1020,20000,0);%建立开环系统模型figure(1)margin(G);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其伯德图如下：BodeDiagramFrequency(rad/s)JBOVeauLHQaMIOWHaeS0HP图5校正前的伯德图武汉理工大学自动控制原理课程设计说明书由伯德图可得:原传递函数的截至频率为：6c=141rad/s原传递函数的穿越频率为：&=15.5rad/sy原传递函数的相角裕度为：=50.9dey原传递函数的幅值裕度为：G

8、(j8c)|=34.2db4.2 校正系统方案选择在进行串联校正时，校正环节的结构通常采用超前校正，滞后校正，滞后-超前校正三种类型，也就是工程上常用的PID调节器。在此次校正中，采用分析法来确定装置的结构和参数。但必须先明确校正方案。对于超前校正，对频率在工1之间的输入信号有微aTT分作用，在该频率范围内，超前校正具有超前相角。超前校正的基本原理就是利用超前相角补偿系统的滞后相角，改善系统的动态性能，如增加相角裕度，提高稳定性等。并且在超前校正中，一般取a=510,即用超前校正补偿的相角一般不超过65。而滞后校正要求使滞后校正的零，极点靠的很近，使之产生的滞后相角很小，还要求滞后校正的零极点

9、靠近原点，尽量不影响中频段。滞后校正的原理是滞后校正的高频段是负增益，因此，滞后校正对系统中高频噪音有削弱作用，增强抗干扰能力。可以利用滞后校正的这一低通滤波所造成的高频衰减特性，降低系统的截至频率，提高系统的相角裕度，以改善系统的暂态性能。因此，可以看出超前校正利用超前网络的超前特性，但滞后并不是利用相位的滞后特性，而是利用其高频衰减特性，这是不符合本次设计的要求的。在本次设计中要求改善系统的动态性能，使之不会因外界干扰而发生系统不稳定现象，否则会使磁盘驱动读取系统无法正常读取数据。同时，需要补偿的相角较小，在超前校正环节中容易实现。因此,采用超前校正的方法对系统进行校正。4.3 校正步骤(

10、1)根据给定的稳态系统系能指标，确定系统的开环增益K.因为超前校正不改变系统的稳态指标，所以第一步仍然是先调整放大器，使系统满足稳态性能指标。因此，由题目给定的已知参数可得：武汉理工大学自动控制原理课程设计说明书G(s)=20s(S,1)(S1)100020K=20(2)利用(1)求得的K,绘制未校正前系统的伯德图。如下图所示:BodeDiagramFrequency(rad/s)IDtf/veduFnowaMLyeQ.vesanp图6校正前的伯德图(3)在伯德图上量取未校正系统的相位裕度和幅值裕度，并计算为使相位裕度达到给定指标所需补偿的超前相角A<P=h¥0+w0。其中，丫

11、为给定的相位裕度指标；¥0为未校正系统的相位裕度；名为附加的角度。加一个正相角巴即ZG0(js-H(j0C)>lZG0(j®c)H(J®c)sC-c取G0(j«c)H(jQc)-G0(j«c)H(Jc)一.,.一4二其中，0c为校正后的截至频率。所以，邛=15，(4)取中m=中,并由武汉理工大学自动控制原理课程设计说明书1sinma二1-sinm求出a。即所需补偿的相角由超前校正装置来提供。可得：a=1.7,(5)为使超前校正装置的最大超前相角出现在校正后系统的截至频率0c上，即'8m=取未校正系统的幅值为-10lga(db)时的

12、频率作为校正后系统的截至频率8c。IIIc可得：作-10lg1.7db直线与未校正系统对数幅频特性曲线相交于e=18.8rad's,取6c=18.8rads。其在伯德图上取点如下图所示：loedvesahpBodeDiagramGm=34.2dB(at141rad/s),Pm=50.9deg(at15.7rad/s)Frequency(rad/s)图7校正刖的伯德图可得：T=0.04s超前校正传递函数为Gc(s)J0.068sc10.04s(6)校验指标：绘制校正后的伯德图，检验是否满足给定的性能指标。当系统仍不满足要求时，则增大名值，从(3)开始重新计算。校正后伯德图的程序为:武汉理

13、工大学自动控制原理课程设计说明书num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);figure(1)margin(G2);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其生成的伯德图如下：BodeDiagramFrequency(rad/s)LMtfcednwIrhQaMluedvesanp图8校正后的伯德图由伯德图可得：校正后的穿越频率a=18.9rad；sg校正后的截至频率.c=174rad.s校正后的相角裕度=60.6deg校正后的幅值裕度G'(j切3=33.2db

14、由校正前后的相角裕度比较可得：该校正满足系统要求的校正指标10武汉理工大学自动控制原理课程设计说明书所以其校正后的开环传递函数为:，10.064s20Gs)10.04ss(上1)(£1)1000204.4校正前后比较(1)校正前后伯德图比较在matlab中校正前后伯彳惠图的程序如下:G1=tf(400000,1,1020,20000,0);%建立开环系统模型num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);figure(1)bode(G1,G2);%绘制伯德图，计算幅值和相角裕度，以及

15、对应的穿越频率和截止频率其校正前后的伯德图如下:JBOVeauLHQaM-90一BodeDiagram0-100-200-300-270-135-180-22510210610410Frequency(rad/s)0810图9校正前后的伯德图11武汉理工大学自动控制原理课程设计说明书从校正前后的伯德图可以看出：该系统的动态稳定性增强，磁盘驱动读取系统时更加稳定，受外界干扰造成的影响减弱。(2)校正前后系统的根轨迹图比较在matlab中其校正前的根轨迹程序如下：G1=tf(400000,1,1020,20000,0);%建立开环系统模型figure(2)rlocus(G1);法制根轨迹其校正前的

16、根轨迹图如下：0.3410.161e+035000.340.16I-4-RootLocus30000.760.640.50.860.940.985.5e+033e+032.5e+032e+031.5e+030.9850.940.860.760.640.5sanoceororxAWTanrQam20001000臀-1000-2000-3000-4000-3000-2000-10000100020001、RealAxis(seconds)图10校正前的根轨迹图此时开环传递函数有0个零点，三个极点，分别为0,-20,-1000。其校正后的根轨迹程序如下：num=27200,400000;den=co

17、nv(1,0,conv(1,1000,conv(1,20,0.04,1);12武汉理工大学自动控制原理课程设计说明书G2=tf(num,den);figure(2)rlocus(G2)；修制根轨迹其校正后的根轨迹图如下:sanoceQTPXAWTanTQan图11校正后的根轨迹图校正后，系统的零点增加，极点也增加。此时有一个零点，为-0.068。四个极点，为0,-0.04,-20,-1000。由校正前后根轨迹图可得：增加的零极点隔虚轴越近，其影响越小。由于校正后其零点接近坐标原点，从而使得闭环系统的阻尼比减小，使得超调量增大。(3)校正前后奈奎斯特图比较系统校正后的奈奎斯特程序如下：num=2

18、7200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);13武汉理工大学自动控制原理课程设计说明书figure(3)nyquist(G2);修制奈奎斯特图其校正后的奈奎斯特图如下:NyquistDiagramPXAVran-QamlSystem:G2Real:-0.156Imag:0.0699Frequency(rad/s):-60.5-0.5-0.4-0.3-0.2-0.10RealAxis图12校正后的奈奎斯特图由图可判断此时的系统是个稳定系统。原系统的奈奎斯特源程序如下：G=tf(400000,1,102

19、0,20000,0);%S立开环系统模型figurenyquist(G);法制奈奎斯特图其生成的奈奎斯特图如下：14-1.2-1-0.8-0.6-0.4-0.2RealAxisNyquistDiagram5oo511武汉理工大学自动控制原理课程设计说明书orxAvyanraaE-2020图13校正前的奈奎斯特图由奈奎斯特可以看出：校正前的系统是一个稳定系统。由两图比较可得：两图都是稳定的，并且增加的零极点对其性能稳定性没多大影响。(4)校正前后斜坡响应曲线比较校正后的斜坡响应程序如下所示：num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0

20、.04,1);G2=tf(num,den);t=0:0.01:100;%设计响应时间figure(4)u=t;%定义输入为斜坡信号lsim(G2,u,t,0);gird批取系统的单位斜坡响应xlabel('t');ylabel('c(t)');title('rampresponse');15武汉理工大学自动控制原理课程设计说明书其生成的校正后的斜坡响应曲线如下:4X10109876543eanrrpmA0102030405060708090100Time(seconds)图14校正后的斜坡响应曲线此时可求出其误差为:=limsGsHs=limss

21、-0s0400000ss1000s20=20_1ess=,kv=1=0.0250LinearSimulationResults其校正前的斜坡输入响应曲线程序如下:num1=0,400000;den1=conv(1,0,conv(1,1000,1,20);G1=tf(num1,den1);t=0:0.01:100;流计响应时间figure(4)u=t;%定义输入为斜坡信号16武汉理工大学自动控制原理课程设计说明书lsim(G1,u,t,0);gird就取系统的单位斜坡响应xlabel('t');ylabel('c(t)');title('rampresponse');其生成的斜坡响应曲线如下：图15校正前的斜坡响应曲线校正前的误差为:kv=limsGsHs=limss-0s-010.064s40000010.04sss1000s20-20120=0.051ess一，kv由上述两图分析可得，校正后系统的响应图形变化不大。并且求得的响应误差相等。综合上面的校正前后伯德图，奈奎斯特图，根轨迹图，斜坡响应曲线图可得：系统的动态稳定性增强，利用其相位超前特性来增大系统的相位裕度改变系统的开环频率特性，使系统的瞬态响应变快。更加符合磁盘驱动读取系统的要求。17武汉理工大学自动控制原理课程设计说明书心得体会通过本次自动控制原理