现代控制理论:直流电机模型.doc
- 文档编号:30814057
- 上传时间:2024-01-30
- 格式:DOC
- 页数:23
- 大小:1.22MB
现代控制理论:直流电机模型.doc
《现代控制理论:直流电机模型.doc》由会员分享,可在线阅读,更多相关《现代控制理论:直流电机模型.doc(23页珍藏版)》请在冰豆网上搜索。
目录
1.直流电机…………………………………………………………3
2.状态空间表达式…………………………………………………6
3.对角标准型及相关分析…………………………………………7
4.系统状态空间表达式求解………………………………………8
5.系统能控性和能观性……………………………………………8
6.系统输入输出传递函数…………………………………………9
7.两种方法判断开环稳定性………………………………………9
8.闭环极点配置……………………………………………………10
9.全维状态观测器设计……………………………………………13
10.带状态观测器的状态反馈控制系统的相关跟踪图…………17
10.带状态观测器的闭环状态反馈系统相关分析………………21
11.结束语…………………………………………………………22
现代控制理论基础
结课作业
选题:
直流电机模型
姓名:
班级:
测控1003
学号:
201002030313
第I条1直流电动机的介绍
节1.011.1研究的意义
直流电机是现今工业上应用最广的电机之一,直流电机具有良好的调速特性、较大的启动转矩、功率大及响应快等优点。
在伺服系统中应用的直流电机称为直流伺服电机,小功率的直流伺服电机往往应用在磁盘驱动器的驱动及打印机等计算机相关的设备中,大功率的伺服电机则往往应用在工业机器人系统和CNC铣床等大型工具上。
[1]
节1.021.2直流电动机的基本结构
直流电动机具有良好的启动、制动和调速特性,可以方便地在宽范围内实现无级调速,故多采用在对电动机的调速性能要求较高的生产设备中。
直流伺服电机的电枢控制:
直流伺服电机一般包含3个组成部分:
图1.1
①磁极:
电机的定子部分,由磁极N—S级组成,可以是永久磁铁(此类称为永磁式直流伺服电机),也可以是绕在磁极上的激励线圈构成。
②电枢:
电机的转子部分,为表面上绕有线圈的圆形铁芯,线圈与换向片焊接在一起。
③电刷:
电机定子的一部分,当电枢转动时,电刷交替地与换向片接触在一起。
直流电动机的启动
电动机从静止状态过渡到稳速的过程叫启动过程。
电机的启动性能有以下几点要求:
1)启动时电磁转矩要大,以利于克服启动时的阻转矩。
2)启动时电枢电流要尽可能的小。
3)电动机有较小的转动惯量和在加速过程中保持足够大的电磁转矩,以利于缩短启动时间。
直流电动机调速可以有:
(1)改变电枢电源电压;
(2)在电枢回路中串调节电阻;
(3)改变磁通,即改变励磁回路的调节电阻Rf以改变励磁电流。
本文章所介绍的直流伺服电机,其中励磁电流保持常数,而有电枢电流进行控制。
这种利用电枢电流对直流伺服电机的输出速度的控制称为直流伺服电机的电枢控制。
如图1.2
图1.2
——定义为电枢电压(伏特)。
——定义为电枢电流(安培)。
——定义为电枢电阻(欧姆)。
——定义为电枢电感(亨利)。
——定义为反电动势(伏特)。
——定义为励磁电流(安培)。
——定义为电机产生的转矩(牛顿•米)
——定义为电机和反射到电机轴上的负载的等效粘带摩擦系数(牛顿•米∕度•)
—定义为电机和反射到电机轴上的负载的等效转动惯量(千克•)。
节1.031.3建立数学模型
电机所产生的转矩,正比于电枢电流I与气隙磁通Φ的乘积,即:
(1-1)
而气隙磁通Φ又正比于激励电流,故式(1-1)改写为
(1-2)
对于激磁电流为常数,合并为一个常数K,称为电机力矩常数。
电枢电流I的正负即代表电机的正反转。
当电枢转动时,在电枢中感应出与电机转轴角速度成正比的电压,称为反电动势,即
(1-3)
其中称为反电动势常数。
电机的速度是由电枢电压E控制,应用基尔霍夫电压定律导出电枢电流I的微分方程式为:
(1-4)
电枢电流I产生力矩,用来克服系统含负载的惯性和摩擦,可得
(1-5)
由式(1-3)与式(1-4)合并移项后可得:
(1-6)
式(1-5)移项后可得:
(1-7)
将式(1-6)与式(1-7)以状态方程式来表示如下:
(1-8)
令R=1、L=0.2、1、=0.1、=5、K=0.5,,代入式(1-8)可得:
A=、B=
、
设,,则
1-9
1、系统状态空间表达式
B=
…………………………………………………………………………………………
MATLAB相关源程序
>>G=ss(A,B,C,D)
a=
x1x2
x1-5-5
x20.1-0.02
b=
u1
x15
x20
c=
x1x2
y101
d=
u1
y10
…………………………………………………………………………………………
2、化为对角标准型并分析
系统特征方程:
==0
系统特征根:
=-4.8975=-0.1225
特征向量:
其逆矩阵为:
则=
变换后状态空间表达式:
由于线性变换矩阵P是非奇异的,因此,状态空间表达式中的系统矩阵A与是相似矩阵,具有相同的基本特征,行列式相同、秩相同、迹相同、特征多项式相同、特征值相同。
…………………………………………………………………………………………
MATLAB相关源程序
>>[P,d]=eig(A)求A的特征值和特征向量
P=
-0.99980.7158
0.0205-0.6983
d=
-4.89750
0-0.1225
>>inv(P)*A*P化A为对角线标准型
ans=
-4.8975-0.0000
0.0000-0.1225
>>P*d*inv(P)
ans=
-5.0000-5.0000
0.1000-0.0200
…………………………………………………………………………………………
3、系统状态空间表达式的求解
在第2个处理中已将系统矩阵A转换为对角线标准型,且矩阵A的特征值互异,则状态转移矩阵为:
==
=
设初始时间,则
线性定常非齐次状态方程的解为:
x(t)=x(0)+d
4、系统的能控性和能观性
A矩阵是2*2矩阵,即n=2
(1)能控性:
若系统能控性矩阵U=的秩为n,则系统状态完全能控。
U
满秩,故系统可控。
(2)能观性:
若系统能观测性矩阵V的秩为n,则系统状态能观。
满秩,故系统可观。
5、系统的输入输出传递函数
………………………………………………………………………………………………………
MATLAB相关源程序
>>[num,den]=ss2tf(A,B,C,D)求系统传递函数
num=
000.5000
den=
1.00005.02000.6000
………………………………………………………………………………………………………6、系统开环稳定性分析
(1)特征根方法
在经典控制理论中,对系统稳定性的分析基于特征方程的所有根是否分布在根平面的左半部分。
所有特征根都分布在左半平面则系统稳定;如果至少有一个特征根分布在右半平面则系统不稳定;如果没有右半平面的根,但在虚轴上有根(即有纯虚根),则系统是临界稳定的。
在以上处理过程中已求出系统特征根为=-4.8975=-0.1225这两个特征根均分布在根平面的左半部分,故系统稳定。
(2)Lyapunov第二法
研究系统的稳定性时,可令u=0。
显然|A|,故原点是系统的平衡状态。
由=得系统的状态方程为
选取李氏函数为(正定)
则沿任意轨迹对时间的导数:
(负定)
又由于当,,故根据相关定理知,平衡点是大范围内渐进稳定的。
8、闭环极点配置
(1)闭环极点配置的目的
控制系统的性能主要取决于系统极点在根平面上的分布,所以通过极点配置改变极点的分布,使得闭环系统阶跃响应的上升时间比开环系统阶跃响应的上升时间缩短,从而获得所希望的动态性能。
(2)闭环极点配置的方法
通过选择状态反馈矩阵K,通过状态反馈使闭环系统的极点恰好配置在根平面上所期望的位置。
(3)充要条件:
系统完全能控。
(4)分析过程
原控制系统
(5)配置过程
反馈控制期望闭环极点:
(原始极点:
)
加状态反馈后特征方程:
………………①
期望的闭环特征方程:
………………②
使以上①②两式多项式对应项的系数相等,得到2个代数方程,即可求出状态反馈阵=
(6)闭环极点配置后系统的传递函数
原系统开环阶跃响应
>>step(G)
极点配置后阶跃响应
由图可看出,极点配置后阶跃响应上升时间比原系统上升时间缩短了三倍左右,故极点配置达到预期效果。
………………………………………………………………………………………………………
MATLAB相关源程序
>>T=[-15-0.4]期望闭环极点
T=
-15.0000-0.4000
>>K=acker(A,B,T)求状态反馈阵
K=
2.076010.3848
>>eig(A-B*K)极点配置验证
ans=
-15.0000
-0.4000
>>[num,den]=ss2tf(A-B*K,B,C,D)求闭环传递函数
num=
0-0.00000.5000
den=
1.000015.40006.0000
………………………………………………………………………………………………………
8、全维状态观测器设计
(1)目的
在现代控制理论中,按各种最优原则建立起来的最优控制系统、解耦系统都离不开状态反馈,然而系统的状态变量不是都能用物理方法量测得到的,有些根本无法量测,因而给状态反馈的物理实现造成了困难,故通过设计一个全维状态观测器实现状态重构,从而获知系统的状态变量。
(2)原则
a、观测器以被观测系统的输入和输出为其输入量,其输出量即为原系统的一个状态渐进估计。
b、被观测系统应是完全能控的。
(3)分析过程
原系统状态方程
观测器状态方程
(4)设计过程
观测器期望极点:
求反馈阵L:
观测器特征多项式:
………………①
期望的特征多项式:
………………②
使以上①②两式多项式对应项的系数相等,得到2个代数方程,即可求出反馈阵
(5)验证配置结果
观测器特征多项式:
求得特征根故,阵求取正确
………………………………………………………………………………………………………
MATLAB相关源程序:
Po=[-20,-10]
Po=
-20-10
>>L=acker(A',C',Po)'
L=
745.0000
24.9800
>>eig(A-L*C)验证观测器极点是否配置在期望极点
ans=
-10
-20
………………………………………………………………………………………………………
(6)形成系统状态估计器
>>est=estim(G,L)
a=
x1x2
x1-5-750
x20.1-25
b=
u1
x1745
x224.98
c=
x1x2
y101
y210
y301
d=
u1
y10
y20
y30
Inputgroups:
NameChannels
Measurement1
Outputgroups:
NameChannels
OutputEstimate1
StateEstimate2,3
Continuous-timemodel.
(7)阶跃响应
>>step(est)
(8)分别得到输出和观测状态的传递函数
>>tf(est)
Transferfunctionfrominputtooutput...
24.98s+199.4
#1:
----------------
s^2+30s+200
745s-110
#2:
----------------
s^2+30s+200
24.98s+199.4
#3:
----------------
s^2+30s+200
Inputgroups:
NameChannels
Measurement1
Outputgroups:
NameChannels
OutputEstimate1
StateEstimate2,3
9、带状态观测器的状态反馈控制系统的状态变量图
跟踪正弦波
>>T=1:
300;plot(T/50,y(1:
300,:
));ylabel('y');legend('系统','观测');
>>T=1:
300;plot(T/50,x1(1:
300,:
));ylabel('x1');legend('系统','观测');
>>T=1:
300;plot(T/50,x2(1:
300,:
));ylabel('x2');legend('系统','观测');
跟踪阶跃信号
>>T=1:
500;plot(T/50,y(1:
500,:
));ylabel('y');legend('系统','观测');
>>T=1:
500;plot(T/50,x1(1:
500,:
));ylabel('x1');legend('系统','观测');
>>T=1:
500;plot(T/50,x2(1:
500,:
));ylabel('x2');legend('系统','观测');
10、带状态观测器的闭环状态反馈控制系统的相关分析
带状态观测器的状态反馈系统由原系统、观测器和控制器三部分组成。
(1)、状态空间表达式
设能控能观测的受控系统为……………………①
状态反馈控制律为……………………②
状态观测器方程为……………………③
由以上三式得整个闭环系统的状态空间表达式为:
即
(2)、闭环传递函数
带观测器状态反馈闭环系统的传递函数阵与是否采用观测器反馈无关,即等于直接状态反馈闭环系统的传递函数阵,即,观测器渐进给出并不影响组合系统的特性。
故,闭环传递函数
(3)、闭环稳定性分析
由观测器构成状态反馈的闭环系统,其特征多项式等于状态反馈部分的特征多项式和观测部分的特征多项式的乘积,两者相互独立。
故可知其特征值即为状态反馈部分和观测部分的特征根,即均位于根平面的左半平面,故由观测器构成状态反馈的闭环系统是稳定的,
五、结束语
(1)所作工作和内容总结
本次大作业我选择的是直流电机模型,通过对直流电机的结构与原理分析,从而构造数学模型,把实际动态问题抽象出来,建立线性动态系统在状态空间的模型,然后运用现代控制理论研究该系统在输入作用下状态运动过程的规律和改变这些规律的可能性与措施;建立和揭示系统的结构性质、动态行为和性能之间的关系。
主要的处理内容有①列出状态空间表达式,化为对角标准型状态空间表达式并进行分析。
对于线性定常系统,若系统的特征值互异,则必存在非奇异变换阵P,经过的变换,可将状态方程化为对角线标准型。
②系统状态空间表达式的求解。
求解一个系统的状态空间表达式首先得求状态转移矩阵,线性定常系统在状态空间中任意时刻的状态是通过状态转移矩阵由初始状态在某一时刻内的转移。
③系统的能控性和能观性判断。
状态方程描述了输入引起状态的变化过程,输出方程则描述了由状态变化引起的输出变化,能控性和能观性分别分析输入对状态的控制能力以及输出对状态的反映能力。
④系统的输入输出传递函数;分别用特征根方法和Lyapunov第二法分析系统的开环稳定性。
特征根法是基于特征方程的根是否分布在根平面左半部分来判断的;Lyapunov第二法判稳时要先解出平衡点然后选取李氏函数进而判稳。
⑤极点配置。
采用状态反馈将系统的闭环极点配置到合适的值,使得闭环统阶跃响应的上升时间比开环系统阶跃响应的上升时间缩短3倍左右。
⑥设计全维状态观测器。
对于一个能观测的系统,他它的状态变量尽管不能直接量测,但是通过其输入和输出以及它们的导数,可以把它重构出来,故可设计一个状态观测器。
⑦带状态观测器的状态反馈控制系统的状态变量图;状态空间表达式和闭环传递函数、闭环系统的稳定性分析。
(2)遇到的问题与分析解决
状态空间表达式的求解。
本系统是一个线性非齐次状态方程,所以在求解时应注意到它的求解公式。
观测器设计。
在观测器设计时,观测器的期望极点最先选的原系统极点的3到5倍,分别为-14和-0.5,但通过MATLAB仿真后,其跟踪效果很不理想,故通过不断更改期望极点,最终选择极点为-20,-10,使L阵达到较合适的值,从而使观测器的状态变量能较快的趋向原系统的状态变量。
在用MATLAB进行相关运算时,由于对MATLAB的使用不熟练,有些程序和使用不会,在图书馆借了一本《基于MATLAB的现代控制》书后相关编程问题得到解决。
在用simulink进行仿真时,出现图形不对、显示不完整的问题,在与同学交流讨论后进行了更改并解决。
(3)体会
通过这次大作业,我收获了很多,这些都是无法通过考试得到的,在上课期间,对于这门课一直感觉很抽象,自己看书也没看懂,它跟实际问题脱轨,所以一直没学明白,只是靠死记硬背,并不真正理解,通过这次作业,我对于现代控制理论这门课有了一个全面的掌握,通过对一个实际问题的数学建模并运用现代控制理论分析系统的状态空间描述、能控性、能观测性、稳定性分析、状态反馈、状态观测器补偿和设计,对所学内容有一个全面、连贯的应用,真正的理解了这些知识点的含义、原理和方法。
另一方面,我进一步掌握了MATLAB的使用,这个软件对于工科学生是一个很有用处的软件,在今后的的学习工作中都会用到,熟悉掌握它的应用是很有必要的。
其次,这次作业更是加深我对word软件的掌握能力,尤其是数学公式的编写,整个过程更是锻炼了我的耐心和细心。
-23-
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 现代 控制 理论 直流电机 模型