华科Matlab大作业_铣床控制器设计

得到单位斜坡响应如图5所示：
图5原系统单位斜坡响应
ess=u(100)-u1(100)
ess =
2.5151
因为原系统为1型系统
计算结果与仿真结果基本一致。
阶跃响应如图6：
step(Gf)
图6校正前系统的阶跃响应
由以上两图知，原系统的单位斜坡响应的稳态误差为2.5，超调量为3.75%。
2
对于超前校正器，它需要两个参数α、T来描述，为了题目的稳态误差的要求，增加一个参数K，他的传递函数为：
G2=tf(num2,den2);
Gf2=feedback(G2,1)
Transfer function:
626 s + 40
--------------------------------------------
162 s^4 + 972.8 s^3 + 815.9 s^2 + 631 s + 40
2010级
《MATAB语言与控制系统仿真》课程
铣床控制器设计
一、
铣床是以各类电动机为动力的传动装置与系统的对象以实现生产过程自动化的技术装置。电气系统是其中的主干部分，在国民经济各行业中的许多部门得到广泛应用。自动铣床具有工作平稳可靠，操作维护方便，运转费用低的特点，已成为现代生产中的主要设备。自动铣床控制系统的设计是一个很传统的课题，现在随着各种先进精确的诸多控制仪器的出现，铣床控制的设计方案也越来越先进，越来越趋于完美。我们之前学习了自动控制原理，要学会把学习的知识用在现实生活的具体设计之中，才能更深的体会和理解这些知识。这次我们通过Matlab仿真工具Simulink，完成对该控制系统仿真简单建模，进行简单仿真分析和研究。
四、
1
图3Simulink对校正前系统仿真
图4Simulink仿真原系统的单位斜坡响应
在scope中只能看到比较粗略的斜坡响应，Matlab没有求斜坡响应的函数，为了计算其稳态误差，仍考虑用step函数。根据闭环传递函数的定义
对于单位阶跃信号有
则
对于单位斜坡信号有
由此可以看出，求系统的单位斜坡响应只需将闭环传递函数的分母多项式乘以拉氏算子“s”，再使用step函数即可。
在该处幅值为 20.3dB。
校正环节不能影响后面的相角特性，零点转折频率设置在新增益穿越频率的1/10处，即
可得
于是校正器的传递函数可以写为
从而可得进行滞后校正后，系统的开环传递函数为
滞后校正器设计完成，校正后单位斜坡响应如图12：
num2=[40*15.65 40];
den2=conv([1 6 5 0],[161.97 1]);
从图中可以看出校正前系统的增益穿越频率为2.57rad/s，相角裕度PM=-6.02 。设计要求超调量限定在20%，则要求的阻尼比大约为0.45。算出要求的相位超调角
在图7中找到对应增益为-10.32dB的点
从而得到
于是超前校正网络的传递函数便求得为
为了抵消衰减，总的开环增益提高三倍，故校正后的开环传递函数为
二、
图1铣床的基本模型
现有一个控制器的传递函数为的磨床，要对它进行校正，使它达到以下的要求：
1、对于斜坡输入R(s)= ，输出地稳态误差为 。
2、对于阶跃输入，输出的超调量要小于20%。
图2系统的抽象框图
图2系统的开环传递函数G(s)：
三、
控制系统的动态性能、稳定小，但是也将减小系统的阻尼比，是系统的超调量和振荡增加；同样，增加开环积分环节可以提高系统类型，是系统跟踪输入信号的能力加强，却有可能导致系统动态性能恶化，甚至不稳定。
num3=[626 40];
den3=[162 972.8 815.9 631 40 0];
G3=tf(num3,den3);
t=[0:0.1:1000];
u=t;
u2=step(G3,t);
plot(t,u);
hold on
plot(t,u2)
图12校正后单位斜坡响应
ess=u(10000)-u2(10000)
最后，感谢老师对我们的指导！
六、
[1]飞思科技产品研发中心. MATLAB7辅助控制系统设计与仿真．北京：电子工业出版社，2005．
[2]Morris Driels.Linear Control Systems Engineering(线性控制系统工程).北京：清华大学出版社，2005
可得
于是校正器的传递函数可以写为
从而可得进行滞后校正后，系统的开环传递函数为
由于重新设计的过程中没有改变K的值，所以稳态误差不会改变，只需要验证阶跃响应的超调量。重新校正后系统的阶跃响应如图15：
num2=[40*21.645 40];den2=conv([1 6 5 0],[339.18 1]);
在图7中找到对应增益为-21.16dB的点
从而得到
于是超前校正网络的传递函数便求得为
为了抵消衰减，总的开环增益提高三倍，故校正后的开环传递函数为
用Matlab重新作出阶跃响应如图10：
num2=[40*1.592 40];
den2=conv([1 6 5 0],[0.0122 1]);
G2=tf(num2,den2);
图8校正后单位斜坡响应
求稳态误差：
u2=step(G3,t);
ess=u(100)-u2(100)
ess =
0.1250
满足第一个设计要求。
校正后系统的bode图如图8图9：
图9校正后阶跃响应
因为47.6%的超调量太大，故将进行下一步的滤波器设计。这次设法获得80 的超前相位，重复上述计算过程，可得：
num=2;
den=conv([1 0],conv([1 1],[1 5]));
G=tf(num,den);
Gf=feedback(G,1);
d=[1 6 5 2 0];
Gs=tf(num,d);
t=[0:0.1:10];
u=t;
u1=step(Gs,t);
plot(t,u);
hold on
plot(t,u1)
因此，为了使控制系统同时具有满意的动态性能和稳态性能，就需要对系统加入一些校正环节。校正的基本思想就是：低频段具有足够高的增益值，以保证稳态误差足够小；中频段具有-20dB/dec的斜率，保证系统的稳定性；高频段具有足够负的斜率，保证抑制噪声干扰能力。
本次实验采用超前校正环节改善系统性能，观察它的阶跃响应和斜坡响应，确定需要采用哪几个参数来描述这个校正器，再根据题目的各个条件，分别算出各个参数需要满足的范围，最后用Matlab对经过校正后的系统进行仿真，如果发现波形还有些不够满足题目要求，再对参数进行微小的调整，最终使得它满足题目的要求。
为了满足第一个设计要求， =1/8，所以 ，此时开环传递函数：
即超前校正器的增益为K=20。
用matlab做出校正前系统的伯德图，如图7
num1=[40];
den1=[1 6 5 0];
G1=tf(num1,den1);
Gf1=feedback(G1,1);
margin(G1)
图7校正前系统bode图
校正后单位斜坡响应如图8：
num2=[40*0.749 40];
den2=conv([1 6 5 0],[0.0696 1]);
G2=tf(num2,den2);
Gf2=feedback(G2,1)
Transfer function:
29.96 s + 40
-------------------------------------------------
本次实验是第一次需要用Matlab作为工具进行设计的实验，但是也是按部就班的按照公式来就可以了，难度不是很大。实验设计中，我先是采用了超前校正，但是经过几次尝试之后发现超前校正无法提供足够的相位超调角，导致阶跃响应超调量过大，然后我果断的选择了滞后校正器完成本次设计。使用Matlab不断地调整参数，没有得过且过，最终满足了所有要求。温故而知新，通过查找资料我对以前学过的知识有了新的认识。
另外，自己提高比较大的部分就是Matlab软件的应用。由于设计性试验不少计算和仿真都需要用Matlab来实现，通过课程学习和这次作业，我从对Matlab一窍不通，到基本熟悉Matlab在控制系统中的各种应用。看书也发现这个软件是多么的强大，自己学会的只是冰山一角。这次实验验激发了我对这个软件的兴趣，在以后的学习中要多多了解。
ess =
0.1250
step(Gf2)
图13校正后阶跃响应
稳态误差满足要求，阶跃响应如图13，可以看出超调量为27.8%，不满足第二个设计要求，需要对滞后校正器进行适当修改。取相位超调角为60 ，重新进行上述计算：
图14校正前系统bode图
由图14知，频率 =0.462rad/s时满足相位裕度要求， 。
Gf2=feedback(G2,1);
step(Gf2)
图10重新校正后阶跃响应
超调量依然不能满足设计要求，经过多次尝试发现由超前校正器去校正系统无法无法提供足够的相位超调量，所以下面选择滞后校正器。
3
为了方便阅读，重新将校正前系统bode图粘贴如下：
图11校正前系统bode图
系统单位斜坡响应稳态误差只与K的大小有关，因此K的取值依然为20不变（详细见上面相关推导）。设计要求超调量限定在20%，则要求的阻尼比大约为0.45。由上面的波特图知，频率 =0.639rad/s时满足相位裕度要求。算出要求的相位超调角
G2=tf(num2,den2);Gf2=feedback(G2,1);
step(Gf2)
图15重新校正后阶跃响应
超调量只有17.1%，满足要求，整个控制器设计完成，系统框图如图16：
图16校正后系统框图