发明名称 基于线性模型和滑模控制器的HCCI发动机燃烧正时控制方法
摘要 本发明请求保护一种基于线性模型和滑模控制器的HCCI发动机燃烧正时控制方法。根据离散非线性模型选择多元线性回归方法对状态方程的各个参数进行拟合,得到线性离散系统的状态方程。利用BP神经网络设计了黑箱模型对线性后的模型进行误差的修正。黑箱模型以HCCI发动机线性模型预测得到的燃烧正时<img file="DDA0000788529480000011.GIF" wi="96" he="57" />、进气歧管温度T<sub>man</sub>、进气歧管压力P<sub>man</sub>和燃油当量比φ为输入,修正后的燃烧正时θ<sub>CA50</sub>为输出。利用该黑箱模型和实际HCCI发动机的理想燃烧正时信号<img file="DDA0000788529480000012.GIF" wi="95" he="64" />,通过二分查找方法查找得到对应线性模型的理想燃烧正时信号<img file="DDA0000788529480000013.GIF" wi="94" he="60" />。并利用该理想燃烧正时信号<img file="DDA0000788529480000014.GIF" wi="96" he="60" />以及线性化后的模型,设计了带卡尔曼滤波器的离散滑模控制器。本发明保证线性模型预测燃烧正时的精度和燃烧正时跟踪理想输入的快速准确性。
申请公布号 CN105221282A 申请公布日期 2016.01.06
申请号 CN201510529105.0 申请日期 2015.08.25
申请人 重庆邮电大学 发明人 郑太雄;潘松;杨新琴;陶炬;徐浩;陈瑶;李玲
分类号 F02D45/00(2006.01)I 主分类号 F02D45/00(2006.01)I
代理机构 重庆市恒信知识产权代理有限公司 50102 代理人 刘小红;李金蓉
主权项 基于线性模型和滑模控制器的HCCI发动机燃烧正时控制方法,其特征在于:建立HCCI发动机离散非线性燃烧正时模型,并得到N组随θ<sub>IVC</sub>变化的T<sub>SOC</sub>、P<sub>SOC</sub>、θ<sub>CA50</sub>的数据;所述T<sub>SOC</sub>、P<sub>SOC</sub>、θ<sub>CA50</sub>和θ<sub>IVC</sub>的含义分别为,HCCI发动机起燃时刻缸内气体温度T<sub>SOC</sub>、起燃时刻缸内气体压力P<sub>SOC</sub>和燃烧正时曲轴角θ<sub>CA50</sub>,进气门关闭正时曲轴角θ<sub>IVC</sub>;T<sub>SOC</sub>、P<sub>SOC</sub>和θ<sub>CA50</sub>作为HCCI发动机线性模型的状态变量,θ<sub>IVC</sub>作为HCCI发动机线性模型的输出,则线性化后的状态方程表示为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>a</mi><mn>11</mn></msub></mtd><mtd><msub><mi>a</mi><mn>12</mn></msub></mtd><mtd><msub><mi>a</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>21</mn></msub></mtd><mtd><msub><mi>a</mi><mn>22</mn></msub></mtd><mtd><msub><mi>a</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>31</mn></msub></mtd><mtd><msub><mi>a</mi><mn>32</mn></msub></mtd><mtd><msub><mi>a</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mi>k</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>b</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><msub><mi>&theta;</mi><mrow><mi>I</mi><mi>V</mi><mi>C</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow>]]></math><img file="FDA0000788529450000011.GIF" wi="983" he="240" /></maths>根据已获得的θ<sub>IVC</sub>、T<sub>SOC</sub>、P<sub>SOC</sub>、θ<sub>CA50</sub>数据,利用多元线性回归的方法对参数a<sub>ij</sub>,b<sub>i</sub>,i=1,2,3,j=1,2,3进行估计,最终得到离散的状态方程x<sub>k+1</sub>=Ax<sub>k</sub>+Bu<sub>k</sub>其中x=[T<sub>SOC</sub> P<sub>SOC</sub> θ<sub>CA50</sub>]<sup>T</sup>u=[θ<sub>IVC</sub>]<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = '[' close = ']'><mtable><mtr><mtd><mi>A</mi></mtd><mtd><mi>B</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msubsup><mi>H</mi><mn>1</mn><mi>T</mi></msubsup><mi>X</mi><msup><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><mo>)</mo><mrow><mo>-</mo><mi>T</mi></mrow></msup></mtd></mtr><mtr><mtd><mrow><msubsup><mi>H</mi><mn>2</mn><mi>T</mi></msubsup><mi>X</mi><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><mo>)</mo></mrow><mrow><mo>-</mo><mi>T</mi></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>H</mi><mn>3</mn><mi>T</mi></msubsup><mi>X</mi><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><mo>)</mo></mrow><mrow><mo>-</mo><mi>T</mi></mrow></msup></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mi>H</mi><mi>X</mi><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><mo>)</mo></mrow><mrow><mo>-</mo><mi>T</mi></mrow></msup></mrow>]]></math><img file="FDA0000788529450000012.GIF" wi="908" he="327" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mn>3</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>N</mi></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>N</mi></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000788529450000013.GIF" wi="589" he="302" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>I</mi><mi>V</mi><mi>C</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>I</mi><mi>V</mi><mi>C</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>T</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mrow><mi>s</mi><mi>o</mi><mi>c</mi><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>C</mi><mi>A</mi><mn>50</mn><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>I</mi><mi>V</mi><mi>C</mi><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000788529450000014.GIF" wi="858" he="308" /></maths>其中A为系统矩阵,B为输入矩阵;将进气歧管温度T<sub>man</sub>、进气歧管压力P<sub>man</sub>、燃油当量比φ作为干扰因素加入所述HCCI发动机线性模型得到x<sub>k+1</sub>=Ax<sub>k</sub>+Bu<sub>k</sub>+Fw其中w=[T<sub>man</sub> P<sub>man</sub> φ]<sup>T</sup><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>f</mi><mn>11</mn></msub></mtd><mtd><msub><mi>f</mi><mn>12</mn></msub></mtd><mtd><msub><mi>f</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>21</mn></msub></mtd><mtd><msub><mi>f</mi><mn>22</mn></msub></mtd><mtd><msub><mi>f</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>31</mn></msub></mtd><mtd><msub><mi>f</mi><mn>32</mn></msub></mtd><mtd><msub><mi>f</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000788529450000015.GIF" wi="426" he="233" /></maths>F为干扰项矩阵,同理F矩阵中的f<sub>ij</sub>,i=1,2,3,j=1,2,3通过多元线性回归方法进行估计得到;利用BP神经网络建立一个黑箱模型对所述加入干扰因素的HCCI发动机线性模型误差进行修正,从而使预测的燃烧正时与实际的燃烧正时相近;该黑箱模型是以加入干扰因素的HCCI发动机线性模型预测得到的燃烧正时<img file="FDA0000788529450000021.GIF" wi="135" he="75" />进气歧管温度T<sub>man</sub>、进气歧管压力P<sub>man</sub>和燃油当量比φ为输入,修正后的燃烧正时θ<sub>CA50</sub>为输出的三层BP神经网络模型。
地址 400065 重庆市南岸区黄桷垭崇文路2号