鲁棒控制作业.docx
- 文档编号:9217623
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:12
- 大小:196.17KB
鲁棒控制作业.docx
《鲁棒控制作业.docx》由会员分享,可在线阅读,更多相关《鲁棒控制作业.docx(12页珍藏版)》请在冰豆网上搜索。
鲁棒控制作业
HarbinInstituteofTechnology
《鲁棒控制》课程作业
课程名称:
鲁棒控制
院系:
航天学院
班级:
姓名:
学号:
教师:
马杰、贺风华
哈尔滨工业大学
作业1
(1)为什么1927年Black在贝尔实验室利用高增益抑制真空管特性变化对放大器精度的影响?
由系统灵敏度的公式可知:
S=
所以增大系统的增益,G增大,所以S减小,所以在低频段处,高增益会使灵敏度降低。
而灵敏度在数值上等于从干扰到输出的传递函数,因此灵敏度降低能抑制干扰对放大器精度的影响,包括真空管特性。
但是,如果一味地提高增益会使Nyquist图更靠近(-1,j0)点,影响了系统的稳定性,,因此也不能过度提高增益去抑制干扰对精度的影响。
(2)Bode稳定性判定定理是否具有鲁棒性?
请说明原因。
Bode稳定判据是指在Bode图上,根据开环频率特性的Bode图在幅频特性大于0dB的频段内相频特性正负穿越-180°线的次数之差来判断系统是否稳定。
若直接利用Bode稳定判据对系统的标称对象进行稳定性判别,只能判断标称系统的稳定性,而不能判断标称性能的内稳定性,也不能判断在有摄动情况下的系统的鲁棒稳定性。
综上所述,此稳定判据不具有鲁棒性。
作业2
考虑如下形式的零点不确定性,该形式适合表示零点从左半平面穿越到右半平面的情况。
给定zp的一个范围,绘制该不确定性的10个随机采样bode图。
由题目要求可得:
,其中
则被控对象可以表示为:
选定
给定的
范围[2,8]
matlab语句:
zp=ureal('zp',5,'range',[2,8]);
Gp=tf([1,zp],[11024]);
bode(usample(Gp,10))
绘图结果:
作业2(第三堂课作业,参数不确定性,matlab绘图)
考虑实际被控对象模型
该被控对象中含有时间滞后,其中
。
求取标准化的加权函数,表示为成型不确定性形式,并画出实际被控对象以及
的幅频特性图。
解答:
选取标称模型为
,
则相对不确定
寻找
,使其幅频特性能够覆盖
。
从幅频特性调试寻找,
是满足条件的加权函数
matlab程序语句:
s=tf('s');
fort=0:
0.01:
0.1
s=tf('s');
g=exp(-t*s)-1;
bode(usample(g,10))
holdon
end
w=0.11*s+0.5;
bode(w);
绘图结果:
经过放大(如下图)可以看出,设计成功。
实际被控对象Bode图如下:
作业三SISO鲁棒控制系统设计
题目:
结合课程学习内容,请查阅资料,完成一个SISO鲁棒控制系统设计过程,包括不确定性模型的建立,加权函数的选择,控制器设计,并给出仿真结果(包括Matlab仿真代码)。
一、设计简介
设计名称:
双手协调机器人单个关节控制系统的设计
某双手协调机器人单个关节的控制系统为单位负反馈系统,被控对象为机械臂,控制模型用传递函数可表示为:
其中,k为考虑到执行部件放大作用而产生的比例系数,经过查阅资料,其标称值与不确定性范围如下表所示:
参数
k
a
标称值
4
0.5
范围
3-5
0.25-0.75
由此可得其标称函数为:
为了确定其不确定性,对式
(1)进行不确定性采样,Matlab程序如下:
k=ureal('k',4,'range',[3,5]);
a=ureal('a',0.7,'range',[0.25,0.75]);
G0=tf(4,[10.50]);
G=tf(k,[1a0]);
Gw=G-G0;
bode(usample(Gw,20))
由于原系统中不确定性参数较多,为了简化分析和设计处理,决定将系统的多个参数的不确定性描述为动态不确定性。
若选择加性不确定性,则系统可描述为:
其中
为系统标称传递函数,
表示系统摄动界函数,
为标准化后的摄动函数。
分析系统的特点,对于此系统,可选取其摄动界函数为:
经过反复试凑,K=2时,可以覆盖系统的全部摄动。
故选取摄动界函数为:
Matlab程序如下:
k=ureal('k',4,'range',[3,5]);
a=ureal('a',0.7,'range',[0.25,0.75]);
G0=tf(4,[10.50]);
G=tf(k,[1a0]);
Gw=G-G0;
Gt=tf(2,[10.50]);
bode(usample(Gw,20),usample(Gt,20))
放大后的图像:
如上页图所示,其中实线表示摄动界函数,如下图放大后,可见它在低频段紧贴不确定性的采样边界,而在高频段留出一定的裕量防止系统因为高频谐振模态及增大的参数不确定性,因此这是一个较为合适的加权函数。
2.设计系统要求性能
由资料分析及推算,本系统工作频段带宽频率大致为
。
选取典型的权函数形式为
其中
为最小带宽频率,A为低频段灵敏度上界,与系统类型还有最大稳态跟踪误差有关,M为高频段灵敏度峰值指标。
为了满足性能要求:
选取最小带宽频率为
最大稳态跟踪误差
选择为0.01,由于
可得A=100。
对于S的最大峰值,选取典型值M=2。
则可以求得加权灵敏度函数
画出
的Bode图如下:
3.控制器的设计
设计控制器时,在这里选用Matlab中的mixsyn函数求解。
由前式(5)(7)可得选取的
和
分别为:
并选取
。
利用Matlab编写程序如下:
k=ureal('k',4,'range',[3,5]);
a=ureal('a',0.7,'range',[0.25,0.75]);
G=tf(k,[0.7-a]);
Wt=tf(2,[0.00174131]);
Ws=tf([12000],[2200000]);
[Q,W,E,R]=mixsyn(G,Wt,Ws,1);
[num,den]=ss2tf(Q.a,Q.b,Q.c,Q.d);
C=tf(num,den)
由上述程序可以得到控制器如下(其中参数较大的问题在最后进行分析):
再通过Matlab绘制最终系统的开环频率特性Bode图,程序如下:
T=C*G
bode(usample(T,10))
则可得到Bode图如下:
由图可见,最终系统稳定。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 鲁棒控制 作业
![提示](https://static.bdocx.com/images/bang_tan.gif)