主权项 |
1、一种引擎曲轴转角与转速的估测方法,主要应用于一引擎控制单元内或该引擎控制单元的周边,该引擎控制单元可检知曲轴转角感知器的曲轴齿讯及一计时器的讯号输出;包括使用:一讯号处理单元,自引擎启动运转时令计时器开始计时,以读取感知器检知的曲轴齿讯,及各齿讯之间由计时器所读取的时间差,并根据下列二方程式计算取得曲轴转速ω<sub>k</sub>与加速度α<sub>k</sub>,表示如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>ω</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>θ</mi><mi>k</mi></msub><mo>-</mo><msub><mi>θ</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><msub><mi>Δt</mi><mi>k</mi></msub></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>α</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>ω</mi><mi>k</mi></msub><mo>-</mo><msub><mi>ω</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>Δt</mi><mi>k</mi></msub><mo>-</mo><msub><mi>Δt</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac></mrow></math>]]></maths>式中:ω<sub>k</sub>为第k个齿讯的曲轴转速,θ<sub>k</sub>为第k个齿讯的曲轴转角,θ<sub>k-1</sub>为第k-1个齿讯的曲轴转角,△t<sub>k</sub>为第k-1个齿讯到第k个齿讯之间时差,△t<sub>k-1</sub>为第k-2个齿讯到第k-1个齿讯之间时差,α<sub>k</sub>为第k个齿讯的曲轴加速度,ω<sub>k-1</sub>为第k-1个齿讯的曲轴转速;在上述计时器归零并重新计时后,根据曲轴转速ω<sub>k</sub>及曲轴加速度α<sub>k</sub>,计算取得动态曲轴取样转角θ<sub>m</sub>、转速ω<sub>m</sub>及加速度α<sub>m</sub>,以下列方程式表示:<maths num="0003"><![CDATA[<math><mrow><msub><mi>θ</mi><mi>m</mi></msub><mo>=</mo><msub><mi>θ</mi><mi>k</mi></msub><mo>+</mo><msub><mi>ω</mi><mi>k</mi></msub><mi>Δt</mi><mo>+</mo><mfrac><mrow><msub><mi>α</mi><mi>k</mi></msub><msup><mi>Δt</mi><mn>2</mn></msup></mrow><mn>2</mn></mfrac></mrow></math>]]></maths>ω<sub>m</sub>=ω<sub>k</sub>+Δtα<sub>k</sub>α<sub>m</sub>=α<sub>k</sub>式中:△t为第k个齿讯到取样需求的时间差;及一由离散动态模型设计而成的估测单元,此离散动态模型以下列状态方程式表示:x<sub>k+1</sub>=Fx<sub>k</sub>+Gu<sub>k</sub>+Γw<sub>k</sub>式中:<maths num="0004"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mrow></math>]]></maths>x<sub>k</sub>=[θ<sub>k</sub> ω<sub>k</sub>]<sup>T</sup>,<maths num="0005"><![CDATA[<math><mrow><mi>G</mi><mo>=</mo><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>T</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></mrow></math>]]></maths>T为观测系统取样时间,u<sub>k</sub>=α<sub>m</sub>,Γ=G为外界对观测系统的干扰矩阵,w<sub>k</sub>为外界对观测系统的干扰;估测单元的输出向量则以下列方程式表示:y<sub>k</sub>=Hx<sub>k</sub>+v<sub>k</sub>式中:由θ<sub>m</sub>可知H=[10],v<sub>k</sub>为感知器量测误差;且估测单元包含由一估测器配合另一回馈增益值矩阵,进行闭回路状态估测,该估测器实质上为一观测系统<img file="C200610080522C00031.GIF" wi="103" he="78" />其与回馈增益值矩阵L一同以下列方程式表示:<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>F</mi><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><mi>G</mi><msub><mi>u</mi><mi>k</mi></msub><mo>+</mo><mi>L</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><mi>H</mi><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>式中:<maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msup><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mover><mi>θ</mi><mo>^</mo></mover><mi>k</mi></msub></mtd><mtd><msub><mover><mi>ω</mi><mo>^</mo></mover><mi>k</mi></msub></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup></mrow></math>]]></maths>为估测器的观测系统状态;L=[L<sub>1</sub> L<sub>2</sub>]<sup>T</sup>为回馈增益值矩阵,L<sub>1</sub>及L<sub>2</sub>为回馈增益值矩阵L的值,用以滤除杂讯干扰,并估测取得曲轴转角与转速。 |