华科电气MATLAB大作业文档格式.docx
- 文档编号:13633070
- 上传时间:2022-10-12
- 格式:DOCX
- 页数:23
- 大小:253.45KB
华科电气MATLAB大作业文档格式.docx
《华科电气MATLAB大作业文档格式.docx》由会员分享,可在线阅读,更多相关《华科电气MATLAB大作业文档格式.docx(23页珍藏版)》请在冰豆网上搜索。
图1中,D(s)为外部扰动,N(s)为测量噪声干扰。
铣床的传递函数为:
首先,在没有控制器
的情况下,看看系统的输出结果。
1.系统的稳态误差理论分析
系统的稳态误差为:
式中,
,
。
利用终值定理,系统斜坡响应的稳态误差为:
显然速度误差
远大于
,误差太大,不满足要求。
2.系统稳态误差仿真分析
用Simulink画出校正前的斜坡输入仿真图如图2所示。
图2.校正前的Simulink斜坡输入仿真
设输入斜坡为
,利用Simulink仿真,在同一示波器中记录了输入
和输出
的波形图。
为方便观察,把示波器中的曲线用MATLAB命令画出如图3所示。
其中,画图代码如下:
curve=plot(ry(:
1),ry(:
2),'
-g'
ry(:
3),'
-r'
)
set(curve
(1),'
linewidth'
2)%设置曲线r(t)的粗细为2
set(curve
(2),'
2)%设置曲线y(t)的粗细为2
legend('
y(t)'
'
r(t)'
)%设置曲线名称人r(t),y(t)
xlabel('
仿真时间(s)'
)%X坐标轴名称标注
ylabel('
幅值'
)%Y轴坐标轴标注
title('
稳态误差'
)%所画图的名称
gridon%添加网格
axis([025025]);
%坐标范围控制
set(gca,'
xtick'
[012...25]);
ytick'
图3.校正前斜坡输入仿真结果
图3中红色为斜坡输入
,绿色为输出
,为方便观察系统的稳态误差,把坐标控制在[0,25]区间,可以看出,在25S时,已经基本达到稳态,此时的误差可从图中得到约为2.5左右,可见与理论相符。
3.阶跃响应仿真分析
图4.校正前系统的Simulink阶跃输入仿真
用同样的方法可以得到阶跃响应的仿真曲线如图5所示。
图5.校正前的阶跃输入仿真结果
从图5中可以大致的看出该系统的阶跃响应的各个性能指标。
为了更加方便准确的获得阶跃响应的各个性能指标,下面利用传递函数和MATLAB代码仿真。
代码如下:
图6.无控制器时的阶跃输入仿真结果
num=[2];
den=conv(conv([1,0],[1,1]),[1,5])
G=tf(num,den)
Gf=feedback(G,1);
t=[0:
0.01:
20];
u=1;
plot(t,u,'
);
holdon
step(Gf,'
-b'
)
den=
1650
G=
2
-----------------
s^3+6s^2+5s
Continuous-timetransferfunction.
从图6中可以看出,在没有控制器时,系统的阶跃响应的超调量为3.75%,符合系统要求。
但是,总体来说,由于系统的速度稳态误差太大,此系统需要改进。
由系统的传递函数
可知,该系统为1型系统,根据系统的型别与稳态误差的关系可知,1型系统的单位斜坡响应稳态误差为
,其中
为速度误差系数。
所以,为减小系统的斜坡响应稳态误差,需要适当的增大
的值。
根据要求,要使稳态误差小于
,也就是要求
,则根轨迹增益要求
下面画出根轨迹增益为1的开环传递函数
的根轨迹如图7所示。
图7.开环传递函数的根轨迹
画根轨迹代码如下:
num=[1];
den=conv(conv([1,0],[1,1]),[1,5]);
G=tf(num,den);
rlocus(G)
由图7可以看出,根轨迹过虚轴时的根轨迹增益约为30.7,所以,当根轨迹增益大于30.7时,系统不稳定,这与上面分析的
矛盾。
所以,只靠增大系统增益,并不能满足系统的性能要求。
根据上面分析,为使系统满足性能要求,必须要更进一步改进。
首先,我们看到,为满足斜坡响应稳态误差的要求,则需要
,即根轨迹增益
,于是,我们不妨取
来看看怎样改进这个新的系统。
四.根轨迹法设计相位滞后环节
为了使系统满足超调量小于20%的要求,在
的前提下,采用根轨迹的方法来设计滞后环节,使其满足要求。
1.相位滞后环节设计
滞后环节的传递函数为
于是有
即
将
代入易得
下面画出开环传递函数
的根轨迹如图所示。
rlocus(G);
gridon
图8.校正前系统根轨迹
根据超调量低于20%,则要求阻尼比大于0.45。
为留足够裕度,于是取阻尼比为0.6。
由图8可知,当阻尼比为0.601时,系统的增益为1.35,即
于是根据
可以得到
取
,则
于是滞后环节的传递函数为
于是,整个系统的传递函数为
2.加入相位滞后环节的仿真分析
利用MATLAB程序可以获得新系统的阶跃响应如图9所示。
图9.加入滞后环节后的阶跃响应
仿真代码如下:
num=2.7*[1,0.01];
den=conv(conv(conv([1,0],[1,1]),[1,5]),[1,0.00054]);
由图9可见,系统的阶跃响应超调量为11.3%,小于20%,满足要求。
五.超前校正设计
1.超前校正器设计
取根轨迹增益
后,我们画出新的传递函数的波特图如图10所示。
图10.
时的波特图
画波特图的代码如下:
num=[50];
bode(G)
gridon;
从图10中可以看出,当幅频响应中增益为0时,此时对应的相位为
明显,
,所以系统不稳定。
所以,为使
时系统稳定,我们可以引入超前校正环节,来改变系统零极点的分布,从而改变根轨迹,增大幅频响应为0时的相位。
校正环节的传递函数为
下面,主要是确定
的大小。
根据要求,系统的阶跃响应超调量小于20%,所以,根据相位裕度
可知,相位裕度为
,为留足够裕度,则系统相位需提前大约
根据
于是可以得到
幅频曲线上移大小为
于是在原波特图图11中可以看到,当
时的频率为
,所以可认为新系统的穿越频率为
,可以得到
于是,校正环节的传递函数为
所以,整个系统的传递函数为
图11.
2.超前校正仿真分析
画出引入超前校正后系统的波特图和根轨迹如图12和图13所示。
num=50*[1.4482,1];
den=conv(conv(conv([1,0],[1,1]),[1,5]),[0.0114,1]);
D=zpk(G);
figure
(1)
figure
(2)
D
D=
6351.8(s+0.6905)
-----------------------
s(s+87.72)(s+5)(s+1)
Continuous-timezero/pole/gainmodel.
图12.超前校正后系统的波特图
图13.超前校正后系统的根轨迹
从图12中可以看出,校正后系统的穿越频率约为
,对应的相位约为
,相位裕度为
,系统稳定。
同时,从图13中可以看出,校正环节引入后,系统分别增加了一个零点和一个极点,从而改变了系统的跟轨迹。
可以看到,根轨迹通过虚轴时的根轨迹增益约为
,所以,在要求的根轨迹增益
下,系统是稳定的,满足要求。
下面看看引入校正环节后的系统的斜坡响应和阶跃响应。
利用前面的方法,仿真结果如下:
图14.超前校正后系统的Simulink斜坡输入仿真
图15.超前校正后系统斜坡输入仿真
图16.超前校正后系统阶跃输入仿真
从图15中可以看出,以为没有改变之前设定的增益
值,斜坡响应的稳态误差为0.1,也就是
,与理论相符;
但是从图6中可以看出,系统的阶跃响应虽然在校正后能够稳定,但是超调量达到40.2%,远大于要求的20%。
所以该系统任然不可以使用,需要进一步改进。
下面,不妨采用滞后校正来试试能否满足要求。
六.滞后校正设计
1.滞后校正器设计
滞后校正环节的传递函数为
根据前面分析,希望得到系统的阻尼比为0.45,于是需要相位裕度
,为留足裕度,取需要引起的相位滞后量为
图17.校正前系统的波特图
从图17中可以看出,当相位为
时,频率为0.457rad/s,所以,新系统的穿越频率大约为
,对应的幅频响应的增益为26dB。
,从而
同时,
,可得
于是,滞后校正环节的传递函数为
于是,整个系统带的传递函数为
2.仿真分析
画出引入滞后校正器后的波特图和阶跃响应图如图18和图19所示。
画图代码如下:
num=50*[21.88,1];
den=conv(conv(conv([1,0],[1,1]),[1,5]),[436.5,1]);
图18.滞后校正后系统的波特图
图19.滞后校正后系统阶跃输入仿真
引入滞后校正后,因为没有改变
的值,所以,斜坡输入稳态误差依然为0.1,此处不再讨论。
从图18中可以看出,引入滞后校正器后,系统的穿越频率为
,相应的相位裕度为
,可见系统有足够的稳定裕度。
另外,从图197中可以看出,系统阶跃响应的超调量为16.9%,小于20%,满足系统的性能要求。
于是可见,引入滞后校正器后,系统的斜坡输入稳态误差和阶跃输入响应的超调量两个指标都满足题目对其要求,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电气 MATLAB 作业