现代控制理论大作业 北科

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

现代控制理论大作业分析对象:汽车悬架系统

指导老师:周晓敏

专业:机械工程

姓名:白国星

学号:S2*******

1、建模

悬架就是车轮或车桥与汽车承载部分之间具有弹性的连接装置的总称,具有传递载荷、缓与冲击、衰减振动以及调节汽车行驶中的车身位置等作用。传统汽车悬驾系统就是被动悬驾,其参数不能改变,无法控制其对不同路面激励的响应,因此对不同路面的适应性较差。为提高汽车的行驶平顺性、操纵稳定性与制动性等性能,人们开始用主动悬架系统来代替传统的被动悬架系统。主动悬架系统能根据路面的情况通过一个动力装置改变悬挂架的参数,改善汽车的各方面性能。

对悬驾系统进行仿真计算首先要建立悬驾系统动力学模型,随后对所建立的模型进行仿真分析。为了简化模型,取汽车的一个车轮的悬驾系统进行研究,该模型可简化为一维二自由度的弹簧阻尼质量系统,图1所示为该模型的模拟图。

图1 悬架系统模型的模拟图

其中u 为动力装置的作用力,w 为路面位移,x1为车身位移,x2为悬驾位移,用车身位移来度量车身的振动情况,并视为系统的输出。路面状况以w 为尺度,并视为系统的一个干扰输入。当汽车从平面落入坑时,w 可用一个阶跃信号来模拟。u 为主动悬架的作用力,它就是系统的控制量。

进行受力分析,由牛顿第二规律可得车身悬架系统的动力学方程为:

()()()()()

1121212212122s s t m x K x x b x x u

m x K x x b x x u K w x ⎧=-+-+⎪⎨

=-+--+-⎪⎩ 设系统状态变量为

:

[]1

2

12x x x x x =

则上面系统动力学方程可改写为状态空间表达式:

x Ax Bu

y Cx Du

=+⎧⎨

=+⎩ 其中:

()1

1

1

1222

200

100001s

s

s t s K K b b A m m m m K K K b b m m m m ⎡⎤⎢⎥⎢⎥⎢⎥--=⎢⎥⎢⎥-+⎢⎥-⎢⎥⎣

12

200

001

01t B m K m m ⎡⎤

⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦

[]1000C = []00D = u u w ⎡⎤=⎢⎥⎣⎦

Matlab 系统模型程序代码: m1=800;m2=320;ks=10000;b=30000; kt=10*ks;

a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1; ks/m2 -(ks+kt)/m2 b/m2 -b/m2]; b=[0 0;0 0;1/m1 0;-1/m2 kt/m2]; c=[1 0 0 0]; d=[0 0]; G=ss(a,b,c,d)

[num,den]=ss2tf(a,b,c,d,1);

tfsys=tf(num,den) eig(a)

figure(1)

impulse(G)

figure(2)

step(G)

系统状态空间模型::

系统传递函数:

系统矩阵特征值:

系统矩阵特征值都在左半平面,故系统稳定。系统脉冲响应:

系统阶跃响应:

2、系统能控性与能观性分析

Matlab系统能控性与能观性分析程序代码: m1=800;m2=320;ks=10000;b=30000;

kt=10*ks;

a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1;

ks/m2 -(ks+kt)/m2 b/m2 -b/m2];

b=[0 0;0 0;1/m1 0;-1/m2 kt/m2];

c=[1 0 0 0];

d=[0 0];

N=size(a);

n=N(1)

ud=ctrb(a,b);%能控矩阵%

ob=obsv(a,c);%能观矩阵%

if rank(ud)==n%判断能控性%

disp('系统能控')

elseif rank(ud)

disp('系统不能控')

end

if rank(ob)==n %判断能观性%

disp('系统能观')

elseif rank(ob)

disp('系统不能观')

end

运行结果:

3、系统稳定性分析

Matlab系统稳定性分析程序代码:

m1=800;m2=320;ks=10000;b=30000;

kt=10*ks;

a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1;

ks/m2 -(ks+kt)/m2 b/m2 -b/m2];

Q=eye(4,4);

P=lyap(a,Q);%求解矩阵P

flag=0;

n=length(a);

for i=1:n

det(P(1:i,1:i))

if(det(P(1:i,1:i))<=0)

flag=1;

end

end

if flag==1

disp('系统不稳定');

else

disp('系统稳定');

end

运行结果:

4、系统状态反馈控制器设计

若系统期望性能指标为超调量小于5%,峰值时间小于0、5秒,即:

(

)

2

2

exp 1505

1%

.p n t w σξπ

ξξ

=--≤=

≤-

则:

0707

9

.n w ξ≥≥

为计算方便取

0707

10

.n w ξ==

则主导极点为:

2121707707..n n w jw j λλξξ==-±-=-±

另两个极点应选择成使其与原点距离远大于主导极点与原点的距离,

相关文档
最新文档