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

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

2、进一步加深对MATLAB软件的认识和理解，学会使用MATLAB语言来进行系统建模、系统校正中的性能参数的求解、以及系统仿真与调试。4、此外，还可以通过本次设计来锻炼独立思考和动手解决控制系统实际问题的能力。2设计任务和要求2.1初始条件磁盘驱动系统的物理结构如图1所示，结构框图如图2所示。图1 磁盘驱动系统物理示意图R(s)-C(s)AmplifierCoilLoad传感器N(s)-图2 磁盘驱动系统的结构框图图2中，Amplifier模块代表控制器，此处用比例控制器；coil模块代表执行机构，此处为电机，该模块的传递函数简化为,load环节代表荷载，该模块的传递函数为。根据对原系统的分析，选

3、择适当的校正装置，设计满足给定性能指标的补偿环节。2.2设计要求（1）分析原系统的稳定性和性能指标；（2）设计一个补偿环节，使系统满足，相角裕度 的性能指标；（3）画出系统在（1）校正前后的奈奎斯特曲线和波特图；（4）用Matlab画出上述每种情况的斜坡响应曲线；（5）用Matlab画出校正前后系统的根轨迹；（6）对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析计算的过程，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。3原系统性能分析3.1校正前系统性能指标分析由初始条件可得：其开环传递函数为：此时可得 闭环传递函数为：在单位阶跃响应中，由单位阶

4、跃图可得上升时间超调时间调节时间超调量在matlab中仿真单位阶跃响应程序如下：num=400000;den=1 1020 20000 400000;sys=tf(num,den);%建立闭环传递函数模型p=roots(den);%计算系统特征根判断系统稳定性t=0:0.01:5;%设定仿真时间为5sfigure(1)step(sys,t);grid%求取系统单位阶跃响应xlabel(t);ylabel(c(t);title(step response);单位阶跃响应曲线图如下：图3 单位阶跃响应图在matlab中的单位斜坡响应程序如下：figure(2)u=t;%定义输入为斜坡信号lsim(

5、sys,u,t,0);grid%求取系统的单位斜坡响应xlabel(t);ylabel(c(t);title(ramp response);其斜坡响应图如下：图4 单位斜坡响应曲线3.2校正前系统稳定性能分析由程序运算可得，该闭环传递函数的特征根如下：p =1.0e+03 * -1.0004 + 0.0000i -0.0098 + 0.0174i -0.0098 - 0.0174i所以可得,由于系统的特征根都具有负实部或者说都位于复平面的左半部，故此时的系统是一个稳定系统。4设计校正系统4.1校正前系统的频域性能指标由系统给定的开环传递函数可得：其原函数的一阶斜坡信号输入得到的速度误差系数为同

6、时，其原函数的相角裕度为校正前的伯德图在matlab中的程序如下：G=tf(400000,1,1020,20000,0);%建立开环系统模型figure(1)margin(G);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其伯德图如下： 图5 校正前的伯德图由伯德图可得：原传递函数的截至频率为：原传递函数的穿越频率为：原传递函数的相角裕度为：原传递函数的幅值裕度为：4.2校正系统方案选择在进行串联校正时，校正环节的结构通常采用超前校正，滞后校正，滞后-超前校正三种类型，也就是工程上常用的PID调节器。在此次校正中，采用分析法来确定装置的结构和参数。但必须先明确校正方案。对于超

7、前校正，对频率在之间的输入信号有微分作用，在该频率范围内，超前校正具有超前相角。超前校正的基本原理就是利用超前相角补偿系统的滞后相角，改善系统的动态性能，如增加相角裕度，提高稳定性等。并且在超前校正中，一般取，即用超前校正补偿的相角一般不超过。而滞后校正要求使滞后校正的零，极点靠的很近，使之产生的滞后相角很小，还要求滞后校正的零极点靠近原点，尽量不影响中频段。滞后校正的原理是滞后校正的高频段是负增益，因此，滞后校正对系统中高频噪音有削弱作用，增强抗干扰能力。可以利用滞后校正的这一低通滤波所造成的高频衰减特性，降低系统的截至频率，提高系统的相角裕度，以改善系统的暂态性能。因此，可以看出超前校正利

8、用超前网络的超前特性，但滞后并不是利用相位的滞后特性，而是利用其高频衰减特性，这是不符合本次设计的要求的。在本次设计中要求改善系统的动态性能，使之不会因外界干扰而发生系统不稳定现象，否则会使磁盘驱动读取系统无法正常读取数据。同时，需要补偿的相角较小，在超前校正环节中容易实现。因此，采用超前校正的方法对系统进行校正。4.3校正步骤（1）根据给定的稳态系统系能指标，确定系统的开环增益K.因为超前校正不改变系统的稳态指标，所以第一步仍然是先调整放大器，使系统满足稳态性能指标。因此，由题目给定的已知参数可得：（2） 利用（1）求得的K，绘制未校正前系统的伯德图。如下图所示：图6 校正前的伯德图（3）

9、在伯德图上量取未校正系统的相位裕度和幅值裕度，并计算为使相位裕度达到给定指标所需补偿的超前相角。其中，为给定的相位裕度指标；为未校正系统的相位裕度；为附加的角度。加一个正相角，即 取其中，为校正后的截至频率。所以，。（4） 取，并由求出。即所需补偿的相角由超前校正装置来提供。可得：（5） 为使超前校正装置的最大超前相角出现在校正后系统的截至频率上，即，取未校正系统的幅值为时的频率作为校正后系统的截至频率。可得：作直线与未校正系统对数幅频特性曲线相交于，取。其在伯德图上取点如下图所示：。图7 校正前的伯德图可得：超前校正传递函数为（6） 校验指标：绘制校正后的伯德图，检验是否满足给定的性能指标。

10、当系统仍不满足要求时，则增大值，从（3）开始重新计算。校正后伯德图的程序为：num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);figure(1)margin(G2);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其生成的伯德图如下：图8 校正后的伯德图由伯德图可得：校正后的穿越频率校正后的截至频率校正后的相角裕度校正后的幅值裕度由校正前后的相角裕度比较可得：该校正满足系统要求的校正指标。所以其校正后的开环传递函数为：4.4校正前后比较（1） 校正前后伯德图比较在matlab中校正

11、前后伯德图的程序如下：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);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其校正前后的伯德图如下：图9 校正前后的伯德图从校正前后的伯德图可以看出：该系统的动态稳定性增强，磁盘驱动读取系统时更加稳定，受外界干扰造成的影响减弱。（2） 校正前后系统的根轨迹图比较在matlab中其校正前的根轨迹程序如下：G1=tf(4000

12、00,1,1020,20000,0);%建立开环系统模型figure(2)rlocus(G1);%绘制根轨迹其校正前的根轨迹图如下：图10 校正前的根轨迹图 此时开环传递函数有0个零点，三个极点，分别为0，-20，-1000。其校正后的根轨迹程序如下：num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);figure(2)rlocus(G2)；%绘制根轨迹其校正后的根轨迹图如下：图11 校正后的根轨迹图校正后，系统的零点增加，极点也增加。此时有一个零点，为-0.068。四个极点，为0，-0.04,-

13、20,-1000。由校正前后根轨迹图可得：增加的零极点隔虚轴越近，其影响越小。由于校正后其零点接近坐标原点，从而使得闭环系统的阻尼比减小，使得超调量增大。（3） 校正前后奈奎斯特图比较系统校正后的奈奎斯特程序如下：num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);figure(3)nyquist(G2);%绘制奈奎斯特图其校正后的奈奎斯特图如下：图12 校正后的奈奎斯特图由图可判断此时的系统是个稳定系统。原系统的奈奎斯特源程序如下：G=tf(400000,1,1020,20000,0);%建立开

14、环系统模型figure(1)nyquist(G);%绘制奈奎斯特图其生成的奈奎斯特图如下：图13 校正前的奈奎斯特图由奈奎斯特可以看出：校正前的系统是一个稳定系统。由两图比较可得：两图都是稳定的，并且增加的零极点对其性能稳定性没多大影响。（4） 校正前后斜坡响应曲线比较校正后的斜坡响应程序如下所示：num=27200,400000;den=conv(1,0,conv(1,1000,conv(1,20,0.04,1);G2=tf(num,den);t=0:0.01:100; %设计响应时间figure(4)u=t; %定义输入为斜坡信号lsim(G2,u,t,0);gird %求取系统的单位斜坡

15、响应xlabel(t);ylabel(c(t);title(ramp response);其生成的校正后的斜坡响应曲线如下：图14 校正后的斜坡响应曲线此时可求出其误差为： 其校正前的斜坡输入响应曲线程序如下：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; %定义输入为斜坡信号lsim(G1,u,t,0);gird %求取系统的单位斜坡响应xlabel(t);ylabel(c(t);title(ramp response);其生成的斜坡响应曲线如下

16、：图15 校正前的斜坡响应曲线校正前的误差为：由上述两图分析可得，校正后系统的响应图形变化不大。并且求得的响应误差相等。综合上面的校正前后伯德图，奈奎斯特图，根轨迹图，斜坡响应曲线图可得：系统的动态稳定性增强，利用其相位超前特性来增大系统的相位裕度改变系统的开环频率特性,使系统的瞬态响应变快。更加符合磁盘驱动读取系统的要求。心得体会通过本次自动控制原理课程设计，更深入的掌握了在自动控制原理理论课上学习到的各种基本控制理论，同时训练了其他各方面的能力；以下是我本次课程设计的具体体会。（1）理论联系实际的真实设计，由于所设计的知识涉及很广，这迫使我查阅许多相关资料文献，学习相关知识，这就训练了我综

17、合运用经典控制理论和相关课程知识的能力。（2）掌握了自动控制原理中各种校正装置的作用及用法，学会根据不同的系统性能指标要求进行合理的系统设计，并调试满足系统的指标；对于串联超前校正，通过校正前后系统Bode图的比较，我知道校正提高了系统相角裕度，降低了截止频率，并且引入了开环零极点。（3）进一步加深对MATLAB软件的认识和理解，学会了使用MATLAB语言来进行基本的数学运算，这让我感觉到MATLAB这个软件的优越性，用它来解复杂的数学方程式避免了传统复杂的解答过程，这都给我们提供了极大的便捷。对于本次设计，我更加深入的了解了MATLAB中各种系统语言，学会了用它来进行系统建模、参数求解、系统仿真与调试、以及各种自动控制原理中的各种图形的绘画。（4）我