发明名称 一种基于递归最小二乘的车辆侧倾角与俯仰角估计方法
摘要 发明公开了一种基于递归最小二乘的车辆侧倾角与俯仰角估计方法,本方法针对陆地行驶四轮车辆,建立符合其行驶特征的车辆动力学模型,进一步通过带遗忘因子的递归最小二乘方法实现对车辆侧倾角与俯仰角的实时、准确估计,且仅需低成本车载传感器,成本较低,所估计出的侧倾角与俯仰角信息可满足有关汽车组合导航与定位的需要。
申请公布号 CN103625475A 申请公布日期 2014.03.12
申请号 CN201310697908.8 申请日期 2013.12.18
申请人 东南大学 发明人 李旭;宋翔
分类号 B60W40/105(2012.01)I;B60W40/11(2012.01)I;B60W40/112(2012.01)I 主分类号 B60W40/105(2012.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 1.一种基于递归最小二乘的车辆侧倾角与俯仰角估计方法,其特征在于:本发明针对陆地行驶四轮车辆,建立符合其行驶特征的车辆动力学模型,进一步通过带遗忘因子的递归最小二乘(Recursive Least Squares,RLS)方法实现对车辆侧倾角与俯仰角的实时、准确估计,且仅需低成本车载传感器,成本较低;具体步骤包括:1)建立汽车行驶过程的动力学模型忽略地球旋转速度,假设车辆的俯仰角速度、侧倾角速度与垂向速度为零,则可建立车辆行驶过程的动力学方程为:<![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>v</mi><mo>.</mo></mover><mi>x</mi></msub><mo>=</mo><msub><mi>a</mi><mi>x</mi></msub><mo>+</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>y</mi></msub><mo>+</mo><mi>g</mi><mi>sin</mi><mi>&theta;</mi></mtd></mtr><mtr><mtd><msub><mover><mi>v</mi><mo>.</mo></mover><mi>y</mi></msub><mo>=</mo><msub><mi>a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><mi>g</mi><mi>sin</mi><mi></mi><mi>&phi;</mi><mi>cos</mi><mi>&theta;</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(1)中,v<sub>x</sub>,v<sub>y</sub>分别表示车辆的纵向速度和侧向速度,a<sub>x</sub>,a<sub>y</sub>分别表示车辆的纵向加速度和侧向加速度,ω<sub>z</sub>表示车辆的横摆角速度,上述定义都是针对车体坐标系,g表示重力加速度,φ,θ分别表示车辆的侧倾角与俯仰角,上标志“·”表示微分,如<img file="FDA0000440694710000014.GIF" wi="42" he="57" />表示对v<sub>x</sub>的微分;由(1)式可得<![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>&theta;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><msub><mover><mi>v</mi><mo>.</mo></mover><mi>x</mi></msub><mo>-</mo><msub><mi>a</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>y</mi></msub></mrow><mi>g</mi></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&phi;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mover><mi>v</mi><mo>.</mo></mover><mi>y</mi></msub></mrow><mrow><mi>g</mi><mi>cos</mi><mi>&theta;</mi></mrow></mfrac><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(2)中,纵向车速的微分可通过纵向车速对时间求导获取,考虑到车辆正常行驶中,v<sub>y</sub>与<img file="FDA0000440694710000015.GIF" wi="46" he="65" />数值较小因而可以忽略,同时,考虑到在大部分路面情况下,车辆的侧倾角与俯仰角通常是小角度,即有arcsin(·)≈·,则式(2)可简化为:<![CDATA[<math><mfenced open='' close='' separators=''><mtable><mtr><mtd><mi>&theta;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><msub><mover><mi>v</mi><mo>.</mo></mover><mi>x</mi></msub><mo>-</mo><msub><mi>a</mi><mi>x</mi></msub></mrow><mi>g</mi></mfrac><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><msub><mover><mi>v</mi><mo>.</mo></mover><mi>x</mi></msub><mo>-</mo><msub><mi>a</mi><mi>x</mi></msub></mrow><mi>g</mi></mfrac></mtd></mtr><mtr><mtd><mi>&phi;</mi><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mover><mi>v</mi><mo>.</mo></mover><mi>y</mi></msub></mrow><mrow><mi>g</mi><mi>cos</mi><mi>&theta;</mi></mrow></mfrac><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mrow><msub><mi>a</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub></mrow><mrow><mi>g</mi><mi>cos</mi><mi>&theta;</mi></mrow></mfrac></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mfenced></math>]]></maths>2)所需车载传感器安装由式(3)可知,只需测得车辆的纵向速度、纵向加速度、侧向加速度与横摆角速度,即可利用所建立并合理简化后的车辆行驶动力学方程,即式(3)来估计车辆的俯仰角和侧倾角;因此,仅需要两个低成本MEMS(Micro-Electro-Mechanic System,微机电系统)加速度传感器,一个低成本MEMS陀螺仪以及车速传感器即可满足测量要求;其中两个低成本MEMS加速度传感器安装于车辆质心位置附近,一个沿车体坐标系纵轴,用以测量纵向加速度,一个沿车体坐标系横轴,用以测量侧向加速度,低成本MEMS陀螺仪也安装于车辆质心位置附近,沿车体坐标系垂向轴安装,用以测量横摆角速度,车速传感器用于测量纵向车速,霍尔车速传感器或轮速传感器等传感器均可采用,在此不做限定,但要求车速测量精度误差小于0.05米/秒,在测得纵向车速信号后,将其对时间求导即可获取其微分;3)基于递归最小二乘的车辆侧倾角与俯仰角估计将式(3)表示为参数识别标准形式:<img file="FDA0000440694710000021.GIF" wi="1025" he="75" />式(4)中,k表示离散时刻,<![CDATA[<math><mrow><mi>&gamma;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mover><mi>&phi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>表示待估参数矩阵,其中,<img file="FDA0000440694710000023.GIF" wi="42" he="69" />与<img file="FDA0000440694710000024.GIF" wi="39" he="78" />分别表示待估计的车辆俯仰角与侧倾角;<![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>v</mi><mo>.</mo></mover><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub><mo>-</mo><msub><mi>a</mi><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mi>y</mi><mo>_</mo><mi>m</mi></mrow></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>z</mi><mo>_</mo><mi>m</mi></mrow></msub><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>表示系统输出矩阵,a<sub>x_m</sub>、a<sub>y_m</sub>与ω<sub>z_m</sub>分别表示利用低成本MEMS传感器所测得的纵向加速度、侧向加速度与横摆角速度,v<sub>x_m</sub>表示通过车速传感器所获取的车辆纵向速度,<img file="FDA0000440694710000026.GIF" wi="97" he="75" />表示v<sub>x_m</sub>的微分,由车速传感器所测得的纵向速度信号对时间求导来获取,即在每一离散时刻k,有:<![CDATA[<math><mrow><msub><mover><mi>v</mi><mo>.</mo></mover><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>v</mi><mrow><mi>x</mi><mo>_</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(5)中,dt表示采样时间间隔,本发明中,dt=0.01(秒);<img file="FDA0000440694710000028.GIF" wi="570" he="167" />表示输入回归矩阵,本发明中上角标<sup>T</sup>表示对矩阵转置;则利用带遗忘因子的递归最小二乘(Recursive Least Squares,RLS)算法实时估计车辆侧倾角与横摆角的估计步骤如下:(1)计算系统输出矩阵y(k),并计算输入回归矩阵<img file="FDA0000440694710000029.GIF" wi="170" he="82" />(2)计算增益矩阵K(k);<img file="FDA00004406947100000210.GIF" wi="679" he="140" />其中,方差矩阵<img file="FDA00004406947100000211.GIF" wi="698" he="126" />参数λ为遗忘因子,能有效减少不再与模型相关的旧数据的影响,并防止协方差发散,通常取值范围在[0.9,1],本发明取0.975;(3)计算待估参数矩阵γ(k);<img file="FDA0000440694710000031.GIF" wi="883" he="75" />其中I为2×2单位矩阵,至此,可实时估计出车辆侧倾角与横摆角。
地址 221700 江苏省徐州市丰县中阳大道18号