车辆主动悬架最优控制Word下载.docx

1、 （2） 式中： ； 汽车悬架可认为是一种连续线性的随机最优控制系统，由最优线性滤波器串接确定性调节器的最优反馈增益系数矩阵组成。这两部分参数可分别加以确定。对于控制要求的性能指标是二次函数积分型的调节器问题，外界干扰是高斯白噪声，综合性能指标为： （3） 此处认为汽车主动悬架的最优控制器为一个终端时间无限的线性调节器，问题仍是寻找最优控制，使目标函数J取极小。线性调节器的主要问题之一是如何选择Q、R阵以获得比较满意的控制过程动态响应，计算机仿真可以解决这个问题。 在悬架设计中，为提高汽车的操纵稳定性和行驶平顺性，应使簧载质量垂直加速度、悬架动扰度及轮胎动变形较小。此外，从实现控制的角度来看，

2、应使所需的控制能量较小。因此式（3）可写为 （4）或写为 （5）其中 这里，q1轮胎动变形加权系数 q2悬架动扰度加权系数式（3.9）中第一、二项为误差指标，表示在整个时间内系统实际状态与平衡之间的误差总和。这一积分越小，说明控制误差越小，性能越好。积分式中第三项为能量指标，表示在整个时间内支付能量的总和。系统状态转移是考控制u（t）来进行的，为要使系统误差很小，则需要支付很大的能量代价。最优反馈增益系数矩阵式可写成 （6）式中，增益值k1k4有明确的物理意义。k1可等效于一放置于簧载和非簧载质量间的弹簧，改变k1则影响簧载质量的固有频率；k2作用于簧载质量的绝对速度上，影响其悬挂阻尼；k3大

3、小涉及轮胎变形，对车轮的垂直弹跳频率产生影响；k4作用于非簧载质量的速度上，影响其非悬挂阻尼。二、 主动悬架系统的能控性，能观测性 能控性和能观测性是系统的一种特性，是现代控制理论中的两个基本概念。 状态完全能控的充分必要条件是能控矩阵满秩； 状态完全能观测的充分必要条件是能观测矩阵满秩。主动悬架系统参数的选取如下，即m1=36kg；m2=240kg；kt=160000N/m； 将参数值带入矩阵，利用Matlab中的函数Coctrb（A,B）求悬架系统的能控矩阵Co，利用函数rA=rank（Co）得矩阵的秩为rA=4,满秩，故系统是能控的。 利用Matlab中的函数Obobsv（A,C）求悬架

4、系统的能观测矩阵Ob，利用函数rBrank（Ob）得矩阵的秩为rB=4,满秩，故系统是能观测的。三、 主动悬架的频域仿真 为了求得主动悬架系统的最优控制u（t）,必须先求得反馈增益矩阵K，而K矩阵的求解决定于黎卡提代数方程的解P矩阵，这可以用计算机来实现求解。 程序用Matlab语言编写，给定一组矩阵A、B、Q、R的有关数据，经过计算，便可以最终得到相应的矩阵P、K的数值。下面取三组不同权系数q1，q2进行计算分析；1）取q1=3.35E5,q2=40.5E5时，由程序得k1=2012.5,k2=977.1,k3=1874.8,k4=31.3，并求得系统的传递函数及幅频特性，绘制系统的幅频特性

5、图%主动悬架q1=3.35e5;q2=40.5e5时的仿真程序：m1=36;m2=240;kt=160000;q1=3.35e5;q2=40.5e5;A=0 1 0 -1;0 0 0 0;0 0 0 -1;0 0 kt/m1 0;B=0;1/m2;0;-1/m1;D=0;1;0;C=0 0 0 0;1 0 0 0;0 0 1 0;E=1/m2;H=0;Q=q2 0 0 0;0 0 q1 0;0 0 0 0;R=1;K,P,F=lqr（A,B,Q,R）M=A-B*K;N=C-E*K;G=ss（M,D,N,H）;G1=tf（G）i=1;for s=0:0.1:80 s=s*2*pi*j;G11=（

6、7.811 *s3 + 580.4 *s2 + 3.727e004 *s + 1.422e-010）/（s4 + 4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004）;G12=（-4385 *s - 1.751e004）/（s4 + 4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004）;G13=（s3 + 4.942 *s2 + 64.29 *s - 2.145e-013）/（s4 + 4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004）;f（i）=abs（G11

7、）;h（i）=abs（G12）;g（i）=abs（G13）; i=i+1;ends=0:80;figureloglog（s,f,-,s,h,-.,s,g,:）legend（加速度,动扰度动载荷 图1. q13.35E5，q240.5E5的幅频特性图由图1可以看出主动悬架的车身加速度、悬架动扰度、轮胎动载荷幅频特性图同被动悬架相似，同样具有双峰，不同的是在低频固有频率附近，主动悬架的响应幅值明显减小，且变化平缓，主动悬架的减振性能较为突出；在高频固有频率附近，主动悬架的响应幅值变化较大。可知取该组权系数时，主动悬架的减振性能的改善程度不够理想；2）取q13.35E8，q240.5E8时，由程序得

8、k1=63640;k2=4863;k3=36146；k4904；及系统的传递函数和幅频特性，绘制幅频特性图%主动悬架q1=3.35e8;q2=40.5e8时的仿真程序:q1=3.35e8;q2=40.5e8; G11=（150.6*s3 + 1.673e004*s2 + 1.179e006*s + 1.653e-008）/（s4 + 45.36*s3 + 5473*s2 + 9.005e004*s + 1.179e006）; G12=（-3290*s - 7.332e004）/（s4 + 45.36*s3 + 5473*s2 + 9.005e004*s + 1.179e006）; G13=（s

9、3 + 45.36*s2 + 2033*s + 5.386e-012）/（s4 + 45.36*s3 + 5473*s2 + 9.005e004*s + 1.179e006）; f（i）=abs（G11）;图2. q13.35E8，q240.5E8的幅频特性图 由图2看出，主动悬架的低频共振频率明显地偏离了低频固有频率，与取前一组加权系数的主动悬架相比，悬架在高频附近幅值变化较大的现象得到很大改善，由于q1，q2主要为轮胎动变形和悬架动扰度的加权系数，可以看出相对于上一组加权系数，轮胎动变形和悬架动扰度的幅频特性得到了显著的改善，即车辆的平顺性和操纵稳定性得到显著提高。 3） 取q13.35E

10、9，q240.5E9时，由程序得k1201250，k27710，k361600，k42340，及系统的传递函数和幅频特性，绘制幅频特性图%主动悬架q1=3.35e9;q2=40.5e9时的仿真程序：q1=3.35e9;q2=40.5e9;Co=ctrb（A,B）; rA=rank（Co）;Ob=obsv（A,C） rB=rank（Ob）;G1=tf（G） i=1; G11=（256.7*s3 + 4.335e004*s2 + 3.727e006*s - 2.178e-008）/（s4 + 97.13*s3 + 9162*s2 + 1.427e005*s + 3.727e006）; G12=（-

11、2477*s - 9.938e004）/（s4 + 97.13*s3 + 9162*s2 + 1.427e005*s + 3.727e006）; G13= （s3 + 97.13*s2 + 6429*s + 1.635e-010）/（s4 + 97.13*s3 + 9162*s2 + 1.427e005*s + 3.727e006）; h（i）=abs（G12）; g（i）=abs（G13）;图3. q13.35E9，q240.5E9的幅频特性图 由图3可知主动悬架的低频共振频率同样明显地偏离了低频固有频率，在高频处，主动悬架的共振峰“几乎”已消失，知悬架在高频处对振动的抑制较为明显。 总之，权系数对悬架性能有较大影响。当q1、q2取得较大时，加速度值增大，动扰