发明名称 一种基于模糊算法的Smith预估器参数估计方法
摘要 本发明阐述了一种基于模糊算法的Smith预估器参数估计方法,该方法可用于实时估计Smith预估器参数,协助Smith预估器准确、快速地提高时滞系统的PID控制性能。该方法利用时滞系统的输出与输入实时地估计Smith预估器的比例系数K<sub>p</sub>(t)以及滞后环节系数τ(t),利用模糊算法实时估计惯性环节系数T(t)。该发明的优点和特点是:能够有效地实时估计出Smith预估器参数K<sub>p</sub>(t)和τ(t),大大提高了实时估计的准确性;惯性环节系数T(t)模糊算法模型的建立,不仅提高了T(t)实时估计的准确性与适应性,而且消除了外部干扰信号对惯性环节系数T(t)实时估计的干扰。该方法提高了基于PID控制的时滞系统的控制质量和准确性,并且保证了PID控制的实时性与稳定性。
申请公布号 CN104950669A 申请公布日期 2015.09.30
申请号 CN201510260549.9 申请日期 2015.05.20
申请人 江苏华豪航海电器有限公司 发明人 赵琳;齐兵;王通达;王兴元;刘峰;倪方忠
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 周全
主权项 一种基于模糊算法的Smith预估器参数估计方法,其特征在于,包括以下步骤:S1,设计时滞系统比例参数K<sub>p</sub>(t)实时估计算法,比例参数K<sub>p</sub>(t)的估计表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>K</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mi>N</mi><mo>&times;</mo><msub><mi>T</mi><mi>B</mi></msub><mrow><mo>(</mo><mtext>t</mtext><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000720773710000011.GIF" wi="378" he="211" /></maths>式中,T<sub>i</sub>(t)为t时刻时滞系统的实时输出,T<sub>B</sub>(t)为t时刻时滞系统的实时输入,<img file="FDA0000720773710000012.GIF" wi="165" he="140" />为N次实验下时滞系统的实时输出总和,N×T<sub>B</sub>(t)为N次实验下时滞系统的实时输入总和;S2,设计时滞系统滞后环节参数τ(t)实时估计算法,滞后环节参数τ(t)的估计表达式为:τ(t)=τ(t‑1)+Δτ(t)式中,τ(t‑1)为t‑1时刻滞后环节参数,Δτ(t)为t时刻滞后环节参数修正值;S3,建立惯性环节系数T(t)模糊算法模型,惯性环节系数T(t)模糊算法模型先在全范围内初步确定惯性环节系数T(t)所处的子区间,然后利用模糊算法模型对惯性环节系数T(t)进行高精度辨识;包括:①将被控对象的惯性环节参数T(t)在全范围内分为若干子区间,初步确定其所在的区间;②选取各子区间最小值作为Smith预估器的惯性环节参数T(t)的参考值,并以对应的系统输出作为系统辨识基准;令惯性环节参数T(t)在各子区间内逐渐增加,并记录下系统输出及其与系统辨识基准间的误差值;③建立模糊算法模型;a.选取系统辨识误差为输入语言变量,惯性环节参数为输出语言变量,其语言值均为5个,即{ZE,PSS,PS,PM,PL},令系统辨识误差和惯性环节参数T(t)语言变量的论域分别为X、Y,即X=Y={0,1,2,3,4,5,6,7,8};b.选取量化因子k<sub>e</sub>将输入变量转换为对应语言变量模糊集的论域,已知系统输入参考值为100,R<sub>s</sub>为系统实际输入值,量化因子k<sub>e</sub>的表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>k</mi><mi>e</mi></msub><mo>=</mo><mfrac><mn>100</mn><msub><mi>R</mi><mi>s</mi></msub></mfrac></mrow>]]></math><img file="FDA0000720773710000021.GIF" wi="189" he="141" /></maths>由此可得到输入、输出语言变量赋值表如下所示:<img file="FDA0000720773710000022.GIF" wi="1638" he="652" />c.根据输入、输出语言变量赋值表,令E为输入语言变量,T为输出语言变量,则模糊算法模型的控制规则如下:R<sub>1</sub>:if E=ZE then T=ZER<sub>2</sub>:if E=PSS then T=PSSR<sub>3</sub>:if E=PS then T=PSR<sub>4</sub>:if E=PM then T=PMR<sub>5</sub>:if E=PL then T=PL控制规则对应的模糊关系矩阵采用与乘运算,如下所示:R<sub>1</sub>=ZE(E)×ZE(T)=ZE<sup>T</sup>οZER<sub>2</sub>=PSS(E)×PSS(T)=PSS<sup>T</sup>οPSSR<sub>3</sub>=PS(E)×PS(T)=PS<sup>T</sup>οPSR<sub>4</sub>=PM(E)×PM(T)=PM<sup>T</sup>οPMR<sub>5</sub>=PL(E)×PL(T)=PL<sup>T</sup>οPL则模糊算法模型的总模糊关系矩阵为:R=R<sub>1</sub>∪R<sub>2</sub>∪R<sub>3</sub>∪R<sub>4</sub>∪R<sub>5</sub>d.模糊推理:模糊集合T<sub>0</sub>、T<sub>2</sub>、T<sub>4</sub>、T<sub>6</sub>和T<sub>8</sub>为输出语言变量{0,2,4,6,8}的模糊推理具体结果:T<sub>0</sub>=ZEοRT<sub>2</sub>=PSSοRT<sub>4</sub>=PSοRT<sub>6</sub>=PMοRT<sub>8</sub>=PLοR模糊集合T<sub>1</sub>、T<sub>3</sub>、T<sub>5</sub>和T<sub>7</sub>为输出语言变量{1,3,5,7}的模糊推理具体结果:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mn>0</mn></msub><mo>+</mo><msub><mi>T</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000720773710000032.GIF" wi="245" he="123" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>T</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mn>2</mn></msub><mo>+</mo><msub><mi>T</mi><mn>4</mn></msub></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000720773710000033.GIF" wi="251" he="125" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>T</mi><mn>5</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mn>4</mn></msub><mo>+</mo><msub><mi>T</mi><mn>6</mn></msub></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000720773710000034.GIF" wi="250" he="128" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>T</mi><mn>7</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mn>6</mn></msub><mo>+</mo><msub><mi>T</mi><mn>8</mn></msub></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000720773710000035.GIF" wi="246" he="124" /></maths>e.模糊判决:采用重心法作为模糊判决算法,输出模糊集合的反模糊化公式为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mn>8</mn></munderover><msub><mi>T</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mn>8</mn></munderover><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000720773710000036.GIF" wi="427" he="277" /></maths>式中,T<sub>k</sub>(t)为惯性环节参数的模糊推理结果,μ(v<sub>k</sub>)为输出语言变量量化变量,模糊判决结果即为所求t时刻的惯性环节参数;S4,通过S1、S2和S3能够实时估计出时滞系统的各参数,并且利用所估计的时滞系统的参数修正Smith预估器参数,Smith预估器利用已更新的内部参数实时计算时滞系统的PID控制修正量,从而实现对时滞系统进行最优化控制。
地址 225211 江苏省扬州市江都区大桥工业园区