发明名称 车辆前方目标物的探测方法
摘要 车辆前方目标物的探测方法属于车辆行车信息感知与处理技术领域。其特征在于,它采用四阶卡尔曼滤波方法对车辆采集到的目标物的测量数据进行滤波处理,能够得到更平滑精确的目标物距离和相对速度信息,并进一步提出了目标的测量连续性判断方法,将雷达丢失的数据和来自另一目标物的数据进行相应的处理。
申请公布号 CN1299125C 申请公布日期 2007.02.07
申请号 CN200410042538.5 申请日期 2004.05.21
申请人 清华大学 发明人 李克强;侯德藻;罗禹贡;杨殿阁;连小珉
分类号 G01S13/93(2006.01);B60T17/18(2006.01) 主分类号 G01S13/93(2006.01)
代理机构 代理人
主权项 1、车辆前方目标物的探测方法,其特征在于,它依次含有由车辆处理器控制执行的以下步骤:1)初始化:给定雷达测量系统距离测量噪声方差σ<sub>r</sub><sup>2</sup>、角度测量噪声方差σ<sub>θ</sub><sup>2</sup>、系统噪声方差Q<sub>u</sub>、系统状态向量初值X<sub>0e</sub>、系统状态估计误差协方差阵初值P<sub>0e</sub>;2)测量目标物距离D<sub>n</sub>和方位角A<sub>n</sub>,其中,下标n表示测量和计算的次数,其初值为1;现时测量值以相应量加下标n表示;3)计算目标物在本车行驶方向,即x方向上距离的现时观测值Z<sub>n</sub>:Z<sub>n</sub>=D<sub>n</sub>cosA<sub>n</sub>4)计算本次x向距离测量的噪声方差:<maths num="001"><![CDATA[ <math><mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup><msup><mrow><mo>&CenterDot;</mo><mi>cos</mi></mrow><mn>2</mn></msup><msub><mi>A</mi><mi>n</mi></msub><mo>+</mo><msup><msub><mi>D</mi><mi>n</mi></msub><mn>2</mn></msup><mo>&CenterDot;</mo><msubsup><mi>&sigma;</mi><mi>&theta;</mi><mn>2</mn></msubsup><msup><mi>sin</mi><mn>2</mn></msup><msub><mi>A</mi><mi>n</mi></msub></mrow></math>]]></maths>5)按下述四阶卡尔曼滤波递推公式进行滤波计算:<maths num="002"><![CDATA[ <math><mrow><msub><mi>X</mi><mi>ny</mi></msub><mo>=</mo><mi>F</mi><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>e</mi></mrow></msub></mrow></math>]]></maths><maths num="003"><![CDATA[ <math><mrow><msub><mi>P</mi><mi>ny</mi></msub><mo>=</mo><mi>F</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>e</mi></mrow></msub><mo>&CenterDot;</mo><msup><mi>F</mi><mo>&prime;</mo></msup><mo>+</mo><mi>G</mi><mo>&CenterDot;</mo><msub><mi>Q</mi><mi>u</mi></msub><mo>&CenterDot;</mo><msup><mi>G</mi><mo>&prime;</mo></msup></mrow></math>]]></maths><maths num="004"><![CDATA[ <math><mrow><msub><mi>K</mi><mi>n</mi></msub><mo>=</mo><msub><mi>P</mi><mi>ny</mi></msub><mo>&CenterDot;</mo><msup><mi>H</mi><mo>&prime;</mo></msup><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mi>H</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mi>ny</mi></msub><mo>&CenterDot;</mo><msup><mi>H</mi><mo>&prime;</mo></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>P<sub>ne</sub>=(I-K<sub>n</sub>·H)·P<sub>ny</sub>X<sub>ne</sub>=X<sub>ny</sub>+K<sub>n</sub>·(Z<sub>n</sub>-H·X<sub>ny</sub>)其中,<maths num="005"><![CDATA[ <math><mrow><msubsup><mi>X</mi><mi>ne</mi><mo>&prime;</mo></msubsup><mo>=</mo><mo>[</mo><msub><mi>x</mi><mi>n</mi></msub><mo>,</mo><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msub><mo>,</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>n</mi></msub><mo>,</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>n</mi></msub><mo>]</mo></mrow></math>]]></maths>表示当前目标物信息的动态最优估计,是矩阵转置符号,x<sub>n</sub>表示目标物x向距离值的最优估计,<img file="C2004100425380002C6.GIF" wi="38" he="49" />表示车辆与目标物间相对速度的最优估计;X<sub>ny</sub>为目标物信息的一步预报估计,F为系统矩阵:<maths num="006"><![CDATA[ <math><mrow><mi>F</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd><mtd><msup><mi>T</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mtd><mtd><msup><mi>T</mi><mn>3</mn></msup><mo>/</mo><mn>6</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd><mtd><msup><mi>T</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>T为系统采样时间;P<sub>ne</sub>是当前系统状态估计误差协方差阵,G是系统误差影响矩阵,<maths num="007"><![CDATA[ <math><mrow><mi>G</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>Q<sub>u</sub>是系统噪声方差;P<sub>ny</sub>为系统状态一步预报估计误差协方差阵;K<sub>n</sub>为卡尔曼滤波增益矩阵;H为观测矩阵,且H=[1 0 0 0];I为单位阵;6)储存P<sub>ne</sub>、X<sub>ne</sub>的值,用于继续的递推计算和输出;储存P<sub>ny</sub>、X<sub>ny</sub>的值,用于测量连续性判断计算;7)测量连续性的判断:7.1)计算观测值的一步预报估计Z<sub>ny</sub>:Z<sub>ny</sub>=H·X<sub>ny</sub>7.2)与X<sub>ny</sub>相应的方差矩阵S:<maths num="008"><![CDATA[ <math><mrow><mi>S</mi><mo>=</mo><mi>H</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mi>ny</mi></msub><mo>&CenterDot;</mo><msup><mi>H</mi><mo>&prime;</mo></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup></mrow></math>]]></maths>7.3)定义观测因子:d<sup>2</sup>=v′·S<sup>-1</sup>·v其中:v=Z<sub>n</sub>-Z<sub>ny</sub>7.4)依据d<sup>2</sup>值判断最新测量值与上次测量值来自同一目标的概率是否等于或超过阈值,如等于或超过,则认为通过滤波估计获得的目标物信息是准确的,则输出本次距离、相对速度和方位角,并返回第2)步进行下一次测量;如小于该阈值,则表明雷达对目标物测量信息丢失,或测量值来自另一个目标,需要进一步判断小于上述阈值的次数是否连续超过预定值e:若未超过预定值e,则假设在目标物测量信息丢失的这一段时间内,目标物运动状态维持不变,按照卡尔曼滤波假设,得到目标物状态的计算方法如下式:<maths num="009"><![CDATA[ <math><mrow><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><msub><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mi>T</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><msup><mi>T</mi><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>6</mn></mfrac><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><msup><mi>T</mi><mn>3</mn></msup></mrow></math>]]></maths><maths num="010"><![CDATA[ <math><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mi>n</mi></msub><mo>=</mo><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mi>T</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><msup><mi>T</mi><mn>2</mn></msup></mrow></math>]]></maths><maths num="011"><![CDATA[ <math><mrow><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>n</mi></msub><mo>=</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub><mi>T</mi></mrow></math>]]></maths><maths num="012"><![CDATA[ <math><mrow><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>n</mi></msub><mo>=</mo><msub><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msub></mrow></math>]]></maths>输出当前目标物在x方向的距离的最优估计x<sub>n</sub>,相对速度<img file="C2004100425380003C6.GIF" wi="36" he="47" />的最优估计,以及方位角A<sub>n</sub>;返回第2)步进行下一次测量;若超过预定值e,则认为有新的目标物出现,将最新的测量目标物设为跟踪目标物,其测量信息作为目标物的状态初值,即最新目标物x向距离的测量值作为X<sub>ne</sub>中x<sub>n</sub>的值,X<sub>ne</sub>中其他状态值设为0,将P<sub>ne</sub>设为单位阵,返回第2)步开始对新的目标物进行跟踪测量。
地址 100084北京市100084-82信箱