发明名称 基于吸引律的二分之一周期重复控制方法
摘要 基于吸引律的二分之一周期重复控制方法,给定具有半周期对称性的参考信号,依据检测获得的输出信号,比较模块产生跟踪误差信号;提供具有有限时间吸引性质的一种幂次吸引律,形成理想误差动态;依据其限定的动态特性构造e/v信号变换模块;构造半周期反馈环节,形成重复控制信号。将由重复控制器提供的控制信号作为被控伺服对象的输入,通过完全消除半周期对称干扰信号,实现位置伺服系统跟随参考信号变化。具体的控制器参数整定工作可依据表征系统收敛性能的指标进行,且提供了表征跟踪误差收敛过程的单调减区域、绝对吸引层和稳态误差带边界。
申请公布号 CN103197556B 申请公布日期 2015.11.25
申请号 CN201310092161.3 申请日期 2013.03.21
申请人 浙江工业大学 发明人 孙明轩;胡轶;吴星;许利达;何海港;张杰;邬玲伟
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 基于吸引律的二分之一周期重复控制方法,包括以下步骤:(1)建立位置伺服系统动态特性的差分方程数学模型:A(q<sup>‑1</sup>)y<sub>k</sub>=q<sup>‑d</sup>B(q<sup>‑1</sup>)u<sub>k</sub>+w<sub>k</sub>其中,d表示延迟,u<sub>k</sub>和y<sub>k</sub>分别表示k时刻的输入和输出信号,w<sub>k</sub>为k时刻的干扰信号,A(q<sup>‑1</sup>)和B(q<sup>‑1</sup>)为关于q<sup>‑1</sup>的多项式,A(q<sup>‑1</sup>)=1+a<sub>1</sub>q<sup>‑1</sup>+…+a<sub>n</sub>q<sup>‑n</sup>B(q<sup>‑1</sup>)=b<sub>0</sub>+b<sub>1</sub>q<sup>‑1</sup>+…+b<sub>m</sub>q<sup>‑m</sup>其中,q<sup>‑1</sup>是一步延迟算子,n为A(q<sup>‑1</sup>)的阶数,m为B(q<sup>‑1</sup>)的阶数,a<sub>1</sub>,…,a<sub>n</sub>,b<sub>0</sub>,…,b<sub>m</sub>为系统参数且b<sub>0</sub>≠0;d为整数,且d≥1;(2)给定参考信号r<sub>k</sub>,该参考信号具有二分之一周期对称特性:P1.r<sub>k</sub>=±r<sub>k‑N/2</sub>或P2.r<sub>k</sub>=±r<sub>k′</sub>其中,k′=(ceil(2k/N)‑1)N‑k,k≥N/2,r<sub>k‑N/2</sub>,r<sub>k′</sub>分别表示k‑N/2,k′时刻的参考信号;(3)根据参考信号二分之一周期对称特性,构造等效扰动d<sub>k</sub>对于P1<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>k</mi></msub><mo>=</mo><msub><mi>w</mi><mi>k</mi></msub><mover><mo>+</mo><mo>-</mo></mover><msub><mi>w</mi><mrow><mi>k</mi><mo>-</mo><mi>N</mi><mo>/</mo><mn>2</mn></mrow></msub></mrow>]]></math><img file="FDA0000782679740000012.GIF" wi="317" he="68" /></maths>对于P2<img file="FDA0000782679740000013.GIF" wi="291" he="71" />k'=(ceil(2k/N)‑1)N‑k,k≥N/2其中,w<sub>k‑N/2</sub>,w<sub>k'</sub>分别表示k‑N/2,k'时刻位置伺服系统所受的干扰信号;(4)构造带干扰抑制作用的幂次吸引律,提供的离散形式的幂次吸引律为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>e</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&rho;</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><msubsup><mi>&beta;e</mi><mi>k</mi><mi>&alpha;</mi></msubsup></mrow>]]></math><img file="FDA0000782679740000011.GIF" wi="440" he="83" /></maths>参数ρ,β分别为可调整参数,其取值范围满足β&gt;0,0&lt;ρ&lt;1;吸引幂次α=p/q,0&lt;p&lt;q且p,q为奇数;(5)设计二分之一周期重复控制器对<img file="FDA0000782679740000028.GIF" wi="329" he="83" />情形,二分之一周期重复控制器表达为u<sub>k</sub>=±u<sub>k‑N/2</sub>+v<sub>k</sub>其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>v</mi><mi>k</mi></msub><mo>=</mo><msup><mrow><mo>&lsqb;</mo><msup><mi>q</mi><mrow><mo>-</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msup><mi>B</mi><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&lsqb;</mo><msub><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&PlusMinus;</mo><msub><mi>y</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>-</mo><mi>N</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>+</mo><msup><mi>A</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>&PlusMinus;</mo><msub><mi>y</mi><mrow><mi>k</mi><mo>-</mo><mi>N</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&rho;</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msub><mrow><mi>&beta;</mi><mo>(</mo><mi>e</mi></mrow><mi>k</mi><mi></mi></msub><msup><mo>)</mo><mi>&alpha;</mi></msup><mo>-</mo><msubsup><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000782679740000029.GIF" wi="1777" he="122" /></maths>对<img file="FDA00007826797400000210.GIF" wi="268" he="68" />情形,二分之一周期重复控制器表达为u<sub>k</sub>=±u<sub>k'</sub>+v<sub>k'</sub>其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>v</mi><msup><mi>k</mi><mo>&prime;</mo></msup></msub><mo>=</mo><msup><mrow><mo>&lsqb;</mo><msup><mi>q</mi><mrow><mo>-</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msup><mi>B</mi><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&lsqb;</mo><msub><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&PlusMinus;</mo><msub><mi>y</mi><mrow><msup><mi>k</mi><mo>&prime;</mo></msup><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><msup><mi>A</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>&PlusMinus;</mo><msub><mi>y</mi><msup><mi>k</mi><mo>&prime;</mo></msup></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&rho;</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>&beta;e</mi><mi>k</mi><mi>&alpha;</mi></msubsup><mo>-</mo><msubsup><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000782679740000021.GIF" wi="1557" he="122" /></maths>(6)建立理想误差动态特性建立具有干扰抑制项的误差动态方程;将u<sub>k</sub>作为位置伺服系统的控制输入信号,测量获得位置伺服系统输出信号y<sub>k</sub>,能够跟随参考信号r<sub>k</sub>变化;从上述二分之一周期重复控制器,得理想误差动态方程:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>e</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&rho;</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><mi>&beta;</mi><msubsup><mi>e</mi><mi>k</mi><mi>&alpha;</mi></msubsup><mo>+</mo><msub><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msubsup><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000782679740000022.GIF" wi="1067" he="108" /></maths>公式1式中,<img file="FDA0000782679740000023.GIF" wi="63" he="73" />用于补偿等效扰动d<sub>k</sub>,<img file="FDA0000782679740000024.GIF" wi="56" he="68" />取扰动d<sub>k</sub>的平均值或者扰动d<sub>k</sub>在上一时刻值d<sub>k‑1</sub>。
地址 310014 浙江省杭州市下城区潮王路18号