发明名称 基于龙伯格滑模观测器的电子节气门自适应反步控制方法
摘要 本发明公开了一种基于龙伯格滑模观测器的电子节气门自适应反步控制方法,涉及汽车电子控制技术领域。首先,针对电子节气门开度变化不可测量问题,该方法基于电子节气门状态方程,采用龙伯格滑模观测器对其进行估计;然后,利用RBF神经网络的逼近特性,对非线性未知量——齿轮间隙扭矩进行逼近;最后,在李雅普诺夫稳定性理论的基础上,结合非线性反步控制方法,分别设计了控制律、RBF网络权值更新律及扰动自适应律。本发明能够较好地克服电子节气门控制中存在的非线性因素以及一些参数容易时变的难题,进一步提高了控制效果和动态响应性能。
申请公布号 CN104018944B 申请公布日期 2016.04.06
申请号 CN201410165505.3 申请日期 2014.04.23
申请人 重庆邮电大学 发明人 郑太雄;杨斌;李银国;王波
分类号 F02D11/10(2006.01)I;F02D41/02(2006.01)I 主分类号 F02D11/10(2006.01)I
代理机构 重庆市恒信知识产权代理有限公司 50102 代理人 刘小红
主权项 基于龙伯格滑模观测器的电子节气门自适应反步控制方法,其特征在于:首先采用龙伯格滑模观测器,根据节气门的实际开度θ和默认开度θ<sub>0</sub>对节气门开度的变化量进行估计,得到节气门开度变化量的估计值<img file="FDA00009043584000000114.GIF" wi="79" he="69" />其次采用RBF神经网络对齿轮间隙扭矩T<sub>g</sub>进行估计;最后结合李雅普诺夫稳定性理论获得自适应反步控制律,输出控制输入电压u控制电子节气门的直流电机对电子节气门的开度进行控制;所述龙伯格滑模观测器的表达式为:<math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mover><mover><mi>x</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>+</mo><msub><mi>l</mi><mn>1</mn></msub><mi>e</mi><mo>+</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mi>s</mi><mi>g</mi><mi>n</mi><mo>(</mo><mi>e</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mover><mover><mi>x</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>a</mi><mn>21</mn></msub><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo><mo>+</mo><msub><mi>a</mi><mn>22</mn></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>+</mo><mi>b</mi><mi>u</mi><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>sgn</mi><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><mi>sgn</mi><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>)</mo><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mi>T</mi><mi>g</mi></msub><mo>+</mo><mi>d</mi><mo>+</mo><msub><mi>l</mi><mn>2</mn></msub><mi>e</mi><mo>+</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mi>s</mi><mi>g</mi><mi>n</mi><mo>(</mo><mi>e</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000904358400000011.GIF" wi="1680" he="183" /></maths>式中,<math><![CDATA[<mrow><msub><mi>a</mi><mn>21</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo><msub><mi>a</mi><mn>22</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><msub><mi>K</mi><mi>t</mi></msub><msub><mi>K</mi><mi>v</mi></msub><mo>+</mo><msub><mi>K</mi><mi>f</mi></msub><msub><mi>R</mi><mi>a</mi></msub></mrow><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><msub><mi>JR</mi><mi>a</mi></msub></mrow></mfrac><mo>;</mo><mi>b</mi><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>t</mi></msub><msub><mi>K</mi><mrow><mi>c</mi><mi>h</mi></mrow></msub></mrow><mrow><msub><mi>K</mi><mi>l</mi></msub><msub><mi>JR</mi><mi>a</mi></msub></mrow></mfrac><mo>;</mo><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>t</mi><mi>f</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000904358400000012.GIF" wi="1717" he="150" /></maths><img file="FDA0000904358400000013.GIF" wi="262" he="134" />|d|≤D;<img file="FDA0000904358400000014.GIF" wi="230" he="71" />l<sub>1</sub>,l<sub>2</sub>,β<sub>1</sub>和β<sub>2</sub>分别为正常数,<img file="FDA0000904358400000015.GIF" wi="58" he="70" />为节气门开度估计值,x<sub>1</sub>=θ表示节气门实际开度,<img file="FDA0000904358400000016.GIF" wi="65" he="70" />为节气门开度变化量的估计值,T<sub>g</sub>为齿轮间隙扭矩,d为外部扰动,K<sub>sp</sub>为弹簧弹性系数,K<sub>l</sub>为齿轮传动比,J为折算到电机侧的系统总转动惯量,K<sub>t</sub>为电机扭矩常数,K<sub>v</sub>为电机反电动势常数,K<sub>f</sub>为滑动摩擦系数,K<sub>ch</sub>为斩波器增益,K<sub>pre</sub>为复位弹簧预紧力矩,K<sub>tf</sub>为库仑摩擦系数,D为d的取值界限;所述RBF神经网络包含输入层、隐含层和输出层,输入层包含两个输入节点,即节气门实际开度与期望开度的误差e<sub>1</sub>和误差变化<img file="FDA0000904358400000017.GIF" wi="70" he="62" />输出层包含一个输出节点,即齿轮间隙扭矩的估计值<img file="FDA0000904358400000018.GIF" wi="84" he="86" />该神经网络的隐含层及网络期望输出分别为:<math><![CDATA[<mrow><msub><mi>h</mi><mi>i</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><mo>|</mo><mi>j</mi><mo>-</mo><msub><mi>c</mi><mi>i</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>i</mi><mn>2</mn></msubsup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000904358400000019.GIF" wi="372" he="151" /></maths>T<sub>g</sub>=w<sup>*T</sup>h(j)+ε式中,c<sub>i</sub>和σ<sub>i</sub>分别为高斯基函数的中心值和宽度,<math><![CDATA[<mrow><mi>j</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>e</mi><mn>1</mn></msub></mtd><mtd><msub><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA00009043584000000110.GIF" wi="262" he="78" /></maths>j为网络输入,i为网络隐含层第i个节点,h=[h<sub>i</sub>]<sup>T</sup>为网络的高斯基函数输出,w<sup>*</sup>为网络的网络理想权值,ε为网络的逼近误差,ε≤ε<sub>N</sub>,ε<sub>N</sub>为逼近误差的最大值;所述自适应反步控制律的表达式为:<math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>u</mi><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mi>b</mi></mfrac><mo>(</mo><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><msub><mi>a</mi><mn>21</mn></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>a</mi><mn>22</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><mi>b</mi><mi>u</mi><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>sgn</mi><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mo>+</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>k</mi><mn>2</mn></msub><mi>sgn</mi><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mover><mi>T</mi><mo>^</mo></mover><mi>g</mi></msub><mo>+</mo><mover><mi>d</mi><mo>^</mo></mover><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mover><mi>e</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>d</mi></msub><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><msub><mi>e</mi><mn>2</mn></msub><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mi>sgn</mi><mo>(</mo><msub><mi>e</mi><mn>2</mn></msub><mo>)</mo><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00009043584000000111.GIF" wi="1078" he="216" /></maths>式中,<math><![CDATA[<mrow><msub><mi>a</mi><mn>21</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo><msub><mi>a</mi><mn>22</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><msub><mi>K</mi><mi>t</mi></msub><msub><mi>K</mi><mi>v</mi></msub><mo>+</mo><msub><mi>K</mi><mi>f</mi></msub><msub><mi>R</mi><mi>a</mi></msub></mrow><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><msub><mi>JR</mi><mi>a</mi></msub></mrow></mfrac><mo>;</mo><mi>b</mi><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mi>t</mi></msub><msub><mi>K</mi><mrow><mi>c</mi><mi>h</mi></mrow></msub></mrow><mrow><msub><mi>K</mi><mi>l</mi></msub><msub><mi>JR</mi><mi>a</mi></msub></mrow></mfrac><mo>;</mo><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>K</mi><mrow><mi>t</mi><mi>f</mi></mrow></msub><mrow><msubsup><mi>K</mi><mi>l</mi><mn>2</mn></msubsup><mi>J</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA00009043584000000112.GIF" wi="1718" he="151" /></maths><img file="FDA00009043584000000113.GIF" wi="262" he="143" />e<sub>2</sub>=x<sub>2</sub>‑x<sub>2d</sub>,c<sub>1</sub>,c<sub>2</sub>和c<sub>3</sub>分别为正常数;θ<sub>d</sub>为节气门期望开度;x<sub>1</sub>=θ表示节气门实际开度;x<sub>2</sub>为节气门开度的变化量;u为控制输入电压;R<sub>a</sub>为电机电阻;<img file="FDA0000904358400000021.GIF" wi="64" he="87" />为齿轮间隙扭矩的估计值;<img file="FDA0000904358400000022.GIF" wi="52" he="76" />为d的估计值;x<sub>2</sub>为虚拟控制律,x<sub>2d</sub>为理想虚拟控制律。
地址 400065 重庆市南岸区黄桷垭崇文路2号