发明名称 一种车辆行驶上下坡状态判定的方法和装置
摘要 本发明公开了一种利用加速度传感器和陀螺进行车辆行驶上下坡状态判定的方法和装置,所述方法包括步骤:S1、数据采集,获取车速、陀螺仪读值、加速度计读值,其中车速为V(t),陀螺仪测量值为gyro(t),加速度计读值为ax(t)、ay(t)、az(t),分别对应三轴加速度计的三个方向;S2、根据陀螺仪读值计算出转弯角速度;S3、测量安装角;S4、使用悬挂修正系数计算悬挂修正值;S5、计算转弯修正值;S6、计算实时坡度,平均坡度,判断坡度类型。通过本发明,可以更为精确地判断车辆当前所处的坡度状态。
申请公布号 CN103632062B 申请公布日期 2017.01.25
申请号 CN201310656432.3 申请日期 2013.12.06
申请人 北京乾图方园软件技术有限公司 发明人 马涛
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京庆峰财智知识产权代理事务所(普通合伙) 11417 代理人 刘元霞
主权项 一种利用加速度传感器和陀螺进行车辆行驶上下坡状态判定的方法,其特征在于包括步骤:S1、数据采集,获取车速、陀螺仪读值、加速度计读值,其中车速为V(t),陀螺仪测量值为gyro(t),加速度计读值为ax(t)、ay(t)、az(t),分别对应三轴加速度计的三个方向;S2、根据陀螺仪读值计算出转弯角速度,转弯角速度w=(gyro(t)‑mid)*gyro_scale_factor,其中mid为陀螺中值,gyro_scale_factor为转弯比例因子;S3、测量安装角,轴向安装角的计算公式:<img file="FDA0001033930350000011.GIF" wi="878" he="183" />其中xAngle为轴向安装角角度,ax(t),ay(t),az(t),分别对应加速度计三轴加速度计的三个方向的测量值;S4、使用悬挂修正系数计算悬挂修正值,悬挂修正值的计算公式为:change1(t)=av(t)*k1*cos(xAngle),其中k1为悬挂修正系数,通过对车辆进行直线加减速测试得出,获得方式为:通过手动调整使得加速度计在三个方向与车辆本身基本一致;车辆进行直线加减速行驶,记录加速度计测量出的车辆前进方向上的加速度值ax(t)和根据车速计算出的加速度值av(t);瞬间悬挂修正系数k1(t)的近似计算公式为:<maths num="0001"><math><![CDATA[<mrow><mi>k</mi><mn>1</mn><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>a</mi><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001033930350000012.GIF" wi="286" he="127" /></maths>S5、计算转弯修正值,转弯修正值的计算公式为:change2(t)=L*w(t)*w(t)*k2,其中L为车辆轴长,w(t)为转弯角速度,k2为转弯修正系数;转弯修正系数k2的获得方式为:a.通过手动调整使得加速度计在三个方向与车辆本身基本一致;b.车辆进行匀速转弯行驶,记录加速度计测量出的车辆前进方向上的加速度值ax(t)和根据陀螺仪测量出的转弯角速度w(t);c.通过下述公式计算<img file="FDA0001033930350000021.GIF" wi="461" he="127" />k2(t)为瞬间转弯修正系数;S6、计算实时坡度,平均坡度,判断坡度类型,实时坡度的计算公式为:ax_new(t)=ax(t)‑change1(t)+change2(t)ay_new(t)=ay(t)az_new(t)=az(t)‑av(t)*sin(xAngle)<maths num="0002"><math><![CDATA[<mrow><mi>A</mi><mi>n</mi><mi>g</mi><mi>l</mi><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><mi>a</mi><mi>x</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msqrt><mrow><mi>a</mi><mi>x</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>a</mi><mi>y</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>a</mi><mi>z</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>-</mo><mi>x</mi><mi>A</mi><mi>n</mi><mi>g</mi><mi>l</mi><mi>e</mi></mrow>]]></math><img file="FDA0001033930350000022.GIF" wi="1468" he="198" /></maths>平均坡度的计算公式:ax_new(t)<sub>n</sub>=ax(t)<sub>n</sub>‑change1(t)<sub>n</sub>+change2(t)<sub>n</sub>ay_new(t)<sub>n</sub>=ay(t)<sub>n</sub>az_new(t)<sub>n</sub>=az(t)<sub>n</sub>‑av(t)<sub>n</sub>*sin(xAngle)<maths num="0003"><math><![CDATA[<mrow><mi>A</mi><mi>n</mi><mi>g</mi><mi>l</mi><mi>e</mi><msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>n</mi></msub><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><mi>a</mi><mi>x</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>n</mi></msub></mrow><msqrt><mrow><mi>a</mi><mi>x</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>n</mi></msub><mn>2</mn></msup><mo>+</mo><mi>a</mi><mi>y</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>n</mi></msub><mn>2</mn></msup><mo>+</mo><mi>a</mi><mi>z</mi><mo>_</mo><mi>n</mi><mi>e</mi><mi>w</mi><msup><msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>n</mi></msub><mn>2</mn></msup></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>-</mo><mi>x</mi><mi>A</mi><mi>n</mi><mi>g</mi><mi>l</mi><mi>e</mi><mo>.</mo></mrow>]]></math><img file="FDA0001033930350000023.GIF" wi="1542" he="183" /></maths>
地址 100194 北京市海淀区东北旺北京中关村软件园孵化器1号楼B座1119室