发明名称 融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法
摘要 本发明涉及一种融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法,其包括以下步骤:1)设置一个质心侧偏角观测系统,包括车速传感器、纵向加速度传感器、横向加速度传感器、横摆角速度传感器及控制器,分别得到原始信号:纵向加速度、横向加速度、横摆角速度;2)对、、分别进行卡尔曼滤波处理,得到处理后的估计值:纵向加速度、横向加速度、横摆角速度;3)分别采用基于卡尔曼滤波和基于信号积分的方式进行质心侧偏角观测:4)对步骤3)两种方法的结果进行加权处理,得到质心侧偏角观测值。基于卡尔曼滤波和加速度积分两种方法进行加权处理,对质心侧偏角进行观测。不仅具有较广的适用范围,而且能够在较低成本下得到较准确的质心侧偏角观测结果。
申请公布号 CN102009653B 申请公布日期 2013.03.27
申请号 CN201010540903.0 申请日期 2010.11.10
申请人 清华大学 发明人 罗禹贡;江青云;褚文博;李克强;连小珉;刘力;杨殿阁;郑四发;王建强
分类号 B60W40/10(2012.01)I;B60W40/12(2012.01)I 主分类号 B60W40/10(2012.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;关畅
主权项 1.一种融合卡尔曼滤波和加速度积分的车轮质心侧偏角观测方法,其特征在于包括以下步骤:1)设置一个质心侧偏角观测系统,包括车速传感器、纵向加速度传感器、横向加速度传感器、横摆角速度传感器及控制器,分别得到原始信号:纵向加速度<img file="FSB00001004124500011.GIF" wi="42" he="51" />、横向加速度<img file="FSB00001004124500012.GIF" wi="46" he="59" />、横摆角速度<img file="FSB00001004124500013.GIF" wi="31" he="52" />;2)对纵向加速度<img file="FSB00001004124500014.GIF" wi="43" he="53" />、横向加速度<img file="FSB00001004124500015.GIF" wi="44" he="59" />、横摆角速度<img file="FSB00001004124500016.GIF" wi="32" he="51" />信号分别进行卡尔曼滤波处理,具体方法如下:假设需要做滤波处理的原始信号为<img file="FSB00001004124500017.GIF" wi="29" he="39" />,实际信号为x,处理后的估计值为<img file="FSB00001004124500018.GIF" wi="26" he="41" />,x(k)表示第k时刻x的瞬时值,ΔT是采样时间间隔,将该信号在时间域上第k时刻进行泰勒展开得到:<img file="FSB00001004124500019.GIF" wi="1499" he="109" />忽略二阶以上的高阶项,选取x(k)及其一阶导数x′(k)为状态变量[x(k) x′(k)]<sup>T</sup>,则式(1)写成如下所示的状态方程:<img file="FSB000010041245000110.GIF" wi="1319" he="210" />其中,w=x″(k),视为随机白噪声误差,假设传感器测量信号为<img file="FSB000010041245000111.GIF" wi="331" he="51" />,其中n<sub>sm</sub>是真实信号所叠加的噪声,视为随机白噪声误差,则有:<img file="FSB000010041245000112.GIF" wi="1229" he="137" />式(2)、(3)构成了完整的卡尔曼滤波器的空间方程结构:X(k)=A·X(k-1)+W(k)(4)<img file="FSB000010041245000113.GIF" wi="481" he="63" />其中<img file="FSB000010041245000114.GIF" wi="302" he="138" />,<img file="FSB000010041245000115.GIF" wi="263" he="138" />,<img file="FSB000010041245000116.GIF" wi="337" he="187" />,<img file="FSB000010041245000117.GIF" wi="249" he="61" />,H=[1 0],N(k)=n<sub>sm</sub>;若假设Q、R分别为W(k)、N(k)的协方差矩阵,构建的卡尔曼滤波器为: <img file="FSB00001004124500021.GIF" wi="1065" he="79" /><img file="FSB00001004124500022.GIF" wi="1590" he="143" />P(k)=(1-Kg(k)·H)·(A·P(k-1)·A<sup>T</sup>+Q)通过式(5),得到处理后的估计值:纵向加速度<img file="FSB00001004124500023.GIF" wi="43" he="53" />、横向加速度<img file="FSB00001004124500024.GIF" wi="45" he="60" />、横摆角速度<img file="FSB00001004124500025.GIF" wi="27" he="51" />;3)分别采用基于卡尔曼滤波和基于信号积分的方式进行质心侧偏角观测,具体方法如下:①基于卡尔曼滤波的质心侧偏角观测质心侧偏角β的定义式为:<img file="FSB00001004124500026.GIF" wi="1056" he="114" />其中,v<sub>x</sub>,v<sub>y</sub>分别为纵向车速和横向车速,对式(6)微分得到:<img file="FSB00001004124500027.GIF" wi="1676" he="129" />同时,根据车辆运动学,质心侧偏角和横摆角速度之间具有如下关系:<img file="FSB00001004124500028.GIF" wi="1058" he="121" />对式(7)、(8)离散化后得到质心侧偏角扩展卡尔曼滤波器空间方程结构如下:<img file="FSB00001004124500029.GIF" wi="1770" he="145" /><img file="FSB000010041245000210.GIF" wi="76" he="46" /><img file="FSB000010041245000211.GIF" wi="712" he="131" />其中,N<sub>β</sub>(k)为真实信号β(k)所叠加的噪声,则将(9)式改为扩展卡尔曼滤波器空间方程结构的标准形式如下:X<sub>β</sub>(k)=A<sub>β</sub>·X<sub>β</sub>(k-1)+W<sub>β</sub>(k)(10)<img file="FSB000010041245000212.GIF" wi="568" he="71" />其中,X<sub>β</sub>(k)=β(k),<img file="FSB000010041245000213.GIF" wi="435" he="125" />,<img file="FSB000010041245000214.GIF" wi="450" he="129" />,H<sub>β</sub>=1,<img file="FSB00001004124500031.GIF" wi="1211" he="139" />;此处,W<sub>β</sub>(k)为随机白噪声误差,若假设Q<sub>β</sub>、R<sub>β</sub>分别为W<sub>β</sub>(k)、N<sub>β</sub>(k)的协方差矩阵,构建的卡尔曼滤波器为:<img file="FSB00001004124500032.GIF" wi="1108" he="82" /><img file="FSB00001004124500033.GIF" wi="1432" he="159" />P<sub>β</sub>(k)=(1-Kg<sub>β</sub>(k)·H<sub>β</sub>)·(A<sub>β</sub>·P<sub>β</sub>(k-1)·A<sub>β</sub><sup>T</sup>+Q<sub>β</sub>)根据式(11)得出质心侧偏角的观测值;②基于信号积分的质心侧偏角观测当式(9)的观测方程中出现v<sub>x</sub>(k)·γ(k)=0的情况时,卡尔曼滤波器将不再适用,对应的工况是车辆纯侧偏,此时,式(7)变化为如下形式:<img file="FSB00001004124500034.GIF" wi="1169" he="116" />将式(12)进行离散化,得到基于信号积分的观测结果如下:<img file="FSB00001004124500035.GIF" wi="1537" he="165" />4)对步骤3)两种方法的结果进行加权处理,得到质心侧偏角观测值,具体方法如下:设利用扩展的卡尔曼滤波方式得到的观测结果为<img file="FSB00001004124500036.GIF" wi="70" he="65" />,利用信号积分方式得到的观测结果为<img file="FSB00001004124500037.GIF" wi="62" he="67" />,质心侧偏角观测器的最终观测结果为<img file="FSB00001004124500038.GIF" wi="35" he="63" />,则有如下关系式:<img file="FSB00001004124500039.GIF" wi="1490" he="67" />当<img file="FSB000010041245000310.GIF" wi="174" he="67" />时,认为未出现明显的横摆,此时K=1,完全采取信号积分的观测方式,即<img file="FSB000010041245000311.GIF" wi="146" he="67" />;当<img file="FSB000010041245000312.GIF" wi="148" he="66" />时,认为此时出现明显的横摆,此时K=0,将由扩展的卡尔曼滤波观测方式,即<img file="FSB000010041245000313.GIF" wi="156" he="66" />;当<img file="FSB000010041245000314.GIF" wi="279" he="66" />时,将两种质心侧偏角观测结果进行加权融合,得到质心侧偏角观测器的最终观测结果<img file="FSB000010041245000315.GIF" wi="34" he="65" />。
地址 100084 北京市海淀区清华园1号清华大学汽车系