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

目录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设计目的与意义1、通过本次课程设计，能够培养自己理论联系实际的设计思想，训练自己综合运用经典控制理论和相关课程知识的能力。2、掌握自动控制原理中各种校正装置的作用及用法，根据不同的系统性能指标要求进行合理的系统设计，并调试满足系统的指标。3、进一步加深对MATLAB软件的认识和理解，学会使用MATLAB语言来进行系统建模、系统校正中的性能参数的求解、以及系统仿真与调试。4、此外，还可以通过本次设计来锻炼独立思考和动手解决控制系统实际问题的能力。2设计任务和要求2.1初始条件磁盘驱动系统的物理结构如图1所示，结构框图如图2所示。图1磁盘驱动系统物理示意图图2磁盘驱动系统的结构框图图2中，Amplifier模块代表控制器，此处用比例控制器80aK；coil模块代表执R(s)-C(s)AmplifierRLmsKCoilLoad）（1J1ss1)(sH传感器N(s)-aK武汉理工大学《自动控制原理》课程设计说明书2行机构，此处为电机，该模块的传递函数简化为10005000(s)G1s,load环节代表荷载，该模块的传递函数为)20(1(s)G2ss。根据对原系统的分析，选择适当的校正装置，设计满足给定性能指标的补偿环节。2.2设计要求（1）分析原系统的稳定性和性能指标；（2）设计一个补偿环节，使系统满足20vK，相角裕度01-5增加PM的性能指标；（3）画出系统在（1）校正前后的奈奎斯特曲线和波特图；（4）用Matlab画出上述每种情况的斜坡响应曲线；（5）用Matlab画出校正前后系统的根轨迹；（6）对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析计算的过程，并包含Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。武汉理工大学《自动控制原理》课程设计说明书33原系统性能分析3.1校正前系统性能指标分析由初始条件可得：其开环传递函数为：)120()11000(20G(s)sss此时可得20K闭环传递函数为：201201100020ssss在单位阶跃响应中，由单位阶跃图可得上升时间s12.0tr超调时间s182.0tp调节时间s496.0ts超调量%17%1001117.1%在matlab中仿真单位阶跃响应程序如下：num=[400000];den=[1102020000400000];sys=tf(num,den);%建立闭环传递函数模型p=roots(den);%计算系统特征根判断系统稳定性t=0:0.01:5;%设定仿真时间为5sfigure(1)step(sys,t);grid%求取系统单位阶跃响应武汉理工大学《自动控制原理》课程设计说明书4xlabel('t');ylabel('c(t)');title('stepresponse');单位阶跃响应曲线图如下：stepresponset(seconds)c(t)00.10.20.30.40.50.60.70.80.9100.20.40.60.81System:sysTime(seconds):0.182Amplitude:1.17System:sysTime(seconds):0.12Amplitude:0.995System:sysTime(seconds):0.496Amplitude:1图3单位阶跃响应图在matlab中的单位斜坡响应程序如下：figure(2)u=t;%定义输入为斜坡信号lsim(sys,u,t,0);grid%求取系统的单位斜坡响应xlabel('t');ylabel('c(t)');title('rampresponse');其斜坡响应图如下：武汉理工大学《自动控制原理》课程设计说明书500.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.70.80.91rampresponset(seconds)c(t)图4单位斜坡响应曲线3.2校正前系统稳定性能分析由程序运算可得，该闭环传递函数的特征根如下：p=1.0e+03*-1.0004+0.0000i-0.0098+0.0174i-0.0098-0.0174i所以可得,由于系统的特征根都具有负实部或者说都位于复平面的左半部，故此时的系统是一个稳定系统。武汉理工大学《自动控制原理》课程设计说明书64设计校正系统4.1校正前系统的频域性能指标由系统给定的开环传递函数可得：其原函数的一阶斜坡信号输入得到的速度误差系数为20)20)(1000(400000lim0)()(lim0ssssssHssGskv同时，其原函数的相角裕度为）（jG180校正前的伯德图在matlab中的程序如下：G=tf(400000,[1,1020,20000,0]);%建立开环系统模型figure(1)margin(G);%绘制伯德图，计算幅值和相角裕度，以及对应的穿越频率和截止频率其伯德图如下：100101102103104105-270-225-180-135-90Phase(deg)BodeDiagramGm=34.2dB(at141rad/s),Pm=50.9deg(at15.7rad/s)Frequency(rad/s)-200-150-100-50050System:GFrequency(rad/s):15.5Magnitude(dB):0.118System:GFrequency(rad/s):141Magnitude(dB):-34.1Magnitude(dB)图5校正前的伯德图武汉理工大学《自动控制原理》课程设计说明书7由伯德图可得：原传递函数的截至频率为：sradc141原传递函数的穿越频率为：sradg5.15原传递函数的相角裕度为：deg9.50原传递函数的幅值裕度为：dbc2.34jG）（4.2校正系统方案选择在进行串联校正时，校正环节的结构通常采用超前校正，滞后校正，滞后-超前校正三种类型，也就是工程上常用的PID调节器。在此次校正中，采用分析法来确定装置的结构和参数。但必须先明确校正方案。对于超前校正，对频率在T1~a1T之间的输入信号有微分作用，在该频率范围内，超前校正具有超前相角。超前校正的基本原理就是利用超前相角补偿系统的滞后相角，改善系统的动态性能，如增加相角裕度，提高稳定性等。并且在超前校正中，一般取10~5a，即用超前校正补偿的相角一般不超过65。而滞后校正要求使滞后校正的零，极点靠的很近，使之产生的滞后相角很小，还要求滞后校正的零极点靠近原点，尽量不影响中频段。滞后校正的原理是滞后校正的高频段是负增益，因此，滞后校正对系统中高频噪音有削弱作用，增强抗干扰能力。可以利用滞后校正的这一低通滤波所造成的高频衰减特性，降低系统的截至频率，提高系统的相角裕度，以改善系统的暂态性能。因此，可以看出超前校正利用超前网络的超前特性，但滞后并不是利用相位的滞后特性，而是利用其高频衰减特性，这是不符合本次设计的要求的。在本次设计中要求改善系统的动态性能，使之不会因外界干扰而发生系统不稳定现象，否则会使磁盘驱动读取系统无法正常读取数据。同时，需要补偿的相角较小，在超前校正环节中容易实现。因此，采用超前校正的方法对系统进行校正。4.3校正步骤（1）根据给定的稳态系统系能指标，确定系统的开环增益K.因为超前校正不改变系统的稳态指标，所以第一步仍然是先调整放大器，使系统满足稳态性能指标。因此，由题目给定的已知参数可得：武汉理工大学《自动控制原理》课程设计说明书8）（）（）（120s11000ss20sG20K（2）利用（1）求得的K，绘制未校正前系统的伯德图。如下图所示：100101102103104105-270-225-180-135-90System:GFrequency(rad/s):15.5Phase(deg):-129System:GFrequency(rad/s):142Phase(deg):-180Phase(deg)BodeDiagramGm=34.2dB(at141rad/s),Pm=50.9deg(at15.7rad/s)Frequency(rad/s)-200-150-100-50050System:GFrequency(rad/s):15.5Magnitude(dB):0.104System:GFrequency(rad/s):138Magnitude(dB):-33.7Magnitude(dB)图6校正前的伯德图（3）在伯德图上量取未校正系统的相位裕度和幅值裕度，并计算为使相位裕度达到给定指标所需补偿的超前相角00-。其中，为给定的相位裕度指标；0为未校正系统的相位裕度；为附加的角度。加一个正相角，即)()(0)'('cjG0cjHcjGcjH）（c'c取)()(0)'()'(0cjHcjGcjHcjG其中，'c为校正后的截至频率。所以，15。（4）取m，并由武汉理工大学《自动控制原理》课程设计说明书9mmasin1sin1求出a。即所需补偿的相角由超前校正装置来提供。可得：7.1a（5）为使超前校正装置的最大超前相角出现在校正后系统的截至频率'c上，即'cm，取未校正系统的幅值为)(lg10-