现代控制理论大作业Word文档下载推荐.docx
- 文档编号:13073279
- 上传时间:2022-10-04
- 格式:DOCX
- 页数:21
- 大小:172.09KB
现代控制理论大作业Word文档下载推荐.docx
《现代控制理论大作业Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《现代控制理论大作业Word文档下载推荐.docx(21页珍藏版)》请在冰豆网上搜索。
图1.1
①磁极:
电机的定子部分,由磁极N—S级组成,可以是永久磁铁(此类称为永磁式直流伺服电机),也可以是绕在磁极上的激励线圈构成。
②电枢:
电机的转子部分,为表面上绕有线圈的圆形铁芯,线圈与换向片焊接在一起。
③电刷:
电机定子的一部分,当电枢转动时,电刷交替地与换向片接触在一起。
直流电动机的启动
电动机从静止状态过渡到稳速的过程叫启动过程。
电机的启动性能有以下几点要求:
1)启动时电磁转矩要大,以利于克服启动时的阻转矩。
2)启动时电枢电流要尽可能的小。
3)电动机有较小的转动惯量和在加速过程中保持足够大的电磁转矩,以利于缩短启动时间。
直流电动机调速可以有:
(1)改变电枢电源电压;
(2)在电枢回路中串调节电阻;
(3)改变磁通,即改变励磁回路的调节电阻Rf以改变励磁电流。
本文章所介绍的直流伺服电机,其中励磁电流保持常数,而有电枢电流进行控制。
这种利用电枢电流对直流伺服电机的输出速度的控制称为直流伺服电机的电枢控制。
如图1.2
图1.2
Ea——定义为电枢电压(伏特)。
Ia——定义为电枢电流(安培)。
Ra——定义为电枢电阻(欧姆)。
La——定义为电枢电感(亨利)。
Eb——定义为反电动势(伏特)。
If——定义为励磁电流(安培)。
Tm——定义为电机产生的转矩(牛顿•米)
Bm——定义为电机和反射到电机轴上的负载的等效粘带摩擦系数(牛顿•米∕度•秒-1)
Jm—定义为电机和反射到电机轴上的负载的等效转动惯量(千克•米2)。
1.3建立数学模型
电机所产生的转矩Tm,正比于电枢电流I与气隙磁通Φ的乘积,即:
Tm=K1nIaΦ(1-1)
而气隙磁通Φ又正比于激励电流If,故式(1-1)改写为
Tm=K1K2IfIa=KIa(1-2)
对于激磁电流If为常数,K1K2If合并为一个常数K,称为电机力矩常数。
电枢电流I的正负即代表电机的正反转。
当电枢转动时,在电枢中感应出与电机转轴角速度成正比的电压,称为反电动势,即
Eb=Kbωm=Kbdθmdt(1-3)
其中Kb称为反电动势常数。
电机的速度是由电枢电压E控制,应用基尔霍夫电压定律导出电枢电流I的微分方程式为:
LadIadt+RIa+Eb=Ea(1-4)
电枢电流I产生力矩,用来克服系统含负载的惯性和摩擦,可得
Jmd2θmdt2+Bmdθmdt=T=KIa(1-5)
由式(1-3)与式(1-4)合并移项后可得:
dIadt=-RaLaIa-KbLaωm+1LaEa(1-6)
式(1-5)移项后可得:
dωmdt=KJmIa-BmJmωm(1-7)
将式(1-6)与式(1-7)以状态方程式来表示如下:
ddtIaωm=-RaLa-KbLaKJm-BmJmIaωm+1La0Ea
yt=01Iωm+0Ea(1-8)
令R=1、L=0.2、Kb=1、Bm=0.1、Jm=5、K=0.5,[1]p229,代入式(1-8)可得:
A=-RaLa-KbLaKJm-BmJm=-5-50.1-0.02、B=1La0=50
C=[01]、D=[0]
设x1=Ia,x2=ωm,则
x=-5-50.1-0.02x+50u
y=[01]x1-9
2.1对所建的模型进行分析
A=-5-50.1-0.02;
B=01;
C=01]
2.2求矩阵的特征值和特征向量
(1)特征值
对于线性定常系统x=Ax+Buy=Cx
则λI-A=det(λI-A)=λn+a1λn-1+…+an-1λ+an
称为系统的特征多项式,令其等于零,即得到系统的特征方程
λI-A=λn+a1λn-1+…+an-1λ+an=0
式中的A为n*n的系统矩阵。
特征方程的根λii=1,2,…,n称为系统的特征值。
因为上述系统为线性定常系统,则λI-A=0所求的根为系统的特征值。
解得λ1=-4.8975;
λ2=-0.1125
得到的系统特征根都为负,系统稳定。
(2)特征向量
设λi是系统一个特征值,若存在一个n维非零向量pi,满足
Api=λipi
或
λiI-APi=0
则称Pi为系统相对应于特征值λi的特征向量。
2.3将状态方程化为对角标准型
对于线性定常系统,若系统的特征值λ1,λ2,…,λn互异,必存在非奇异变换矩阵P,经过x=Px或x=p-1x的变换,可将状态方程化为对角线标准型,即
X=λ1⋯0⋮⋱⋮0⋯λnx+bu
λ1和λ2互异,必存在非奇异变换矩阵P,经过x=Px的变换,将状态方程化对角为标准型。
由APi=λiPi
求出矩阵P1=-0.99980.0205P2=0.7158-0.6983
P=-0.99980.02050.7158-0.6983
P-1=-1.0217-0.0300-1.0473-1.4628
X=P-1APX+P-1BU
Y=CPX
X=-4.897500-0.1225X+-1.0473-1.4628U
Y=[0.0205-0.6983]X
得到新的矩阵:
A’=-4.897500-0.1225;
B’=-1.0473-1.4628;
C’=0.0205-0.6983
2.4从状态空间表达式求取传递函数阵
线性定常系统的状态空间表达式为:
x=Ax+Buy=Cx
对上式取拉氏变换,可得
sX(s)-X(0)=AX(s)+BU(s)
Y(s)=CX(s)+DU(s)
设初始条件X(0)=0,则有
sX(s)=AX(s)+BU(s)
X(s)=C(SI-A)-1BU(s)
Y(s)=[C(SI-A)-1B]U(s)
=GsU(s)
得到传递函数:
Gs=C(sI-A)-1B
代入数值:
得到Gs=0.5s2+5.02s+0.6
根据传递函数可以写出新的能控标准Ι型的状态空间:
x=01-0.6-5.02x+01u
y=0.50
还可以写出标准ΙΙ型(能观型实现)
x=0-0.61-5.02x+0.50u
y=01
2.5系统状态空间表达式的求解
设线性定常系统的齐次状态方程为
x=Ax(2-1)
在这里初始值为t=t0,初始状态为x(t0)。
系统齐次状态方程在初始状态x(t0)激励下的解x(t)(其中t≥t0),称为系统的自由解或零输入解。
设齐次状态方程的解x(t)为t的向量幂级数形式
即
x(t)=b0+b1t+b2t2+…++biti+…(2-2)
式(2-2)代入式(2-1),得
b1+2b2+…+ibiti-1+…
=A(b0+b1t+b2t2+…++biti+…)(2-3)
由于式(2-2)是式(2-1)的解,所以式(2-3)对所有时间t均成立,故式(2-3)等号两边t的同次幂级数应相等,即
b1=Ab0
b2=12Ab1=12A2b0
…
bi=1iAbi-1=1i!
Aib0
对式(2-2),若令t=0,可得
b0=x(0)(2-4)
x(0)为状态向量x(t)的初始值,即定常系统的初始状态。
将bi(i=1,2,…)及b0代入式(2-2),得到
x(t)=b0+Ab0t+12iA2b0t2+1i!
Aib0ti+…
=(I+At+12iA2t2+…+1i!
Aiti+…)x(0)(2-5)
仿照标量指数函数e-at展开成幂级数形式
e-at=1+at+12iat2+…+1i!
aiti+…(2-6)
将式(2-5)括号内n*n矩阵的无穷项级数和称为矩阵指数函数,记为eAt,即
eAt=I+At+12iA2t2+…+1i!
Aiti+…(2-7)
则齐次状态方程的解可表示为x(t)=eAtx(0)
这里
求eAt这里t=1s
根据前面的分析,我们求出了系统的特征值和特征向量
eAt=Pe∧tP-1
eAt=-0.0114-0.91860.01840.9035
2.6Lyapunov第二法分析系统的开环稳定性
线性定常连续系统
x=Ax
在平衡状态Xe=0处,渐进稳定的充要条件是:
对给定的一个正定对称矩阵Q,存在一个正定的对称矩阵P,且满足矩阵方程:
AΤP+PA=-Q
而标量函数v(x)=xTPx是这个系统的一个二次型形式的李雅普诺夫函数。
我们这里我们选定对称矩阵Q时,常取Q=I,于是得
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 现代 控制 理论 作业
![提示](https://static.bdocx.com/images/bang_tan.gif)