主权项 |
基于吸引律的二分之一周期重复控制方法,包括以下步骤:(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>ρ</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><msubsup><mi>βe</mi><mi>k</mi><mi>α</mi></msubsup></mrow>]]></math><img file="FDA0000782679740000011.GIF" wi="440" he="83" /></maths>参数ρ,β分别为可调整参数,其取值范围满足β>0,0<ρ<1;吸引幂次α=p/q,0<p<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>[</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>]</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><msub><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>±</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>′</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>±</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>ρ</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msub><mrow><mi>β</mi><mo>(</mo><mi>e</mi></mrow><mi>k</mi><mi></mi></msub><msup><mo>)</mo><mi>α</mi></msup><mo>-</mo><msubsup><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>]</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>′</mo></msup></msub><mo>=</mo><msup><mrow><mo>[</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>]</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><msub><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>±</mo><msub><mi>y</mi><mrow><msup><mi>k</mi><mo>′</mo></msup><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><msup><mi>A</mi><mo>′</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>±</mo><msub><mi>y</mi><msup><mi>k</mi><mo>′</mo></msup></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>ρ</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>βe</mi><mi>k</mi><mi>α</mi></msubsup><mo>-</mo><msubsup><mi>d</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>]</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>ρ</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>k</mi></msub><mo>-</mo><mi>β</mi><msubsup><mi>e</mi><mi>k</mi><mi>α</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>。 |