主权项 |
一种基于ARX模型和PI预测函数的无刷直流电动机控制方法,其特征在于,具体包括以下步骤:步骤1、使用ARX模型作为无刷直流电动机被辨识系统的数学模型,ARX模型表示为:A(z<sup>‑1</sup>)y(k)=B(z<sup>‑1</sup>)u(k‑1)+ξ(k)其中,A(z<sup>‑1</sup>)=1+a<sub>1</sub>z<sup>‑1</sup>+…+a<sub>na</sub>z<sup>‑na</sup>,B(z<sup>‑1</sup>)=1+b<sub>1</sub>z<sup>‑1</sup>+…+b<sub>nb</sub>z<sup>‑nb</sup>,A(z<sup>‑1</sup>)、B(z<sup>‑1</sup>)为关于平移算子z<sup>‑1</sup>的多项式,a,b分别为A(z<sup>‑1</sup>)、B(z<sup>‑1</sup>)关于平移算子z<sup>‑1</sup>的系数,下标na、nb分别为输入、输出阶次,ξ(k)为零均值、方差为σ<sup>2</sup>的系统噪声干扰,u(k‑1)为k‑1时刻的输入电压,y(k)为k时刻的输出转速;通过递推最小二乘法离线辨识ARX模型中的未知参数A(z<sup>‑1</sup>)、B(z<sup>‑1</sup>),得到ARX预测模型;步骤2、根据步骤1所辨识得到的ARX预测模型,通过下式计算控制量u(k):<img file="FDA0000804173030000011.GIF" wi="317" he="78" />其中:<img file="FDA0000804173030000012.GIF" wi="352" he="86" /><img file="FDA0000804173030000013.GIF" wi="1353" he="86" /><img file="FDA0000804173030000014.GIF" wi="1325" he="89" /><b>f</b><sub>n</sub>=[f<sub>1</sub>(j),f<sub>2</sub>(j),…,f<sub>N</sub>(j)],j=0,1,…,P‑1<b>f</b><sub>n</sub>(0)=[f<sub>1</sub>(0),f<sub>2</sub>(0),…,f<sub>N</sub>(0)]<sup>T</sup><b>d</b>=[d(k+1),d(k+2),…,d(k+P)]<sup>T</sup>Δ<b>d</b>=[Δd(k+1),Δd(k+2),…,Δd(k+P)]<sup>T</sup><b>G</b><sub>j</sub>=[g<sub>0</sub>,g<sub>1</sub>,…,g<sub>j</sub>]<sup>T</sup>;j=0,1,…,P‑1Δ<b>G</b><sub>j</sub>=[g<sub>0</sub>,g<sub>1</sub>‑g<sub>0</sub>…,g<sub>j</sub>‑g<sub>j‑1</sub>]<sup>T</sup>;j=0,1,…,P‑1<img file="FDA0000804173030000021.GIF" wi="566" he="146" />d(k+j)=(1‑β<sup>j</sup>)[c(k)‑y(k)]+[1‑S<sub>j</sub>(z<sup>‑1</sup>)]y<sub>m</sub>(k)‑H<sub>j</sub>(z<sup>‑1</sup>)u(k‑1) (1) 其中,u(k)是系统第k时刻的控制量输出;K<sub>p</sub>、K<sub>i</sub>分别为比例项系数、积分项系数;f<sub>n</sub>(j)为基函数在t=jT<sub>s</sub>时的值,T<sub>s</sub>为采样周期;<b>f</b><sub>n</sub>为基函数的值构成的矩阵,N为基函数的个数,n为基函数个数的索引,n=1,2,…,N;Q和R分别表示误差加权矩阵和控制加权矩阵;c(k)为k时刻无刷直流电动机的转速设定值;y(k)为k时刻无刷直流电动机的实际转速;y<sub>m</sub>(k)为k时刻无刷直流电动机的模型输出转速;d(k+j)是由式(1)表示的多项式;<b>d</b>(k+j)为d(k+j)的向量表示形式,简写为<b>d</b>;Δ<b>d</b>(k+j)为向量<b>d</b>(k+j)与<b>d</b>(k+j‑1)的向量差,简写为Δ<b>d</b>;<b>G</b><sub>j</sub>是由G<sub>j</sub>(z<sup>‑1</sup>)多项式中的系数组成的向量,Δ<b>G</b><sub>j</sub>是<b>G</b><sub>j</sub>与<b>G</b><sub>j‑1</sub>的向量差;T<sub>s</sub>是采样时间;T<sub>r</sub>是参考轨迹的期望响应时间;β为参考轨迹的衰减因子;P为预测步长;j为第j步预测时刻;S<sub>j</sub>(z<sup>‑1</sup>)、H<sub>j</sub>(z<sup>‑1</sup>)是各系数的递推解,是由A(z<sup>‑1</sup>)、B(z<sup>‑1</sup>)和第j步预测时刻唯一确定的多项式,表示如下:S<sub>j</sub>(z<sup>‑1</sup>)=s<sub>j,0</sub>+s<sub>j,1</sub>z<sup>‑1</sup>+…+s<sub>j,na‑1</sub>z<sup>‑(na‑1)</sup>H<sub>j</sub>(z<sup>‑1</sup>)=h<sub>j,0</sub>+h<sub>j,1</sub>z<sup>‑1</sup>+…+h<sub>j,nb‑1</sub>z<sup>‑(nb‑1)</sup>其中,s、h分别为S<sub>j</sub>(z<sup>‑1</sup>)、H<sub>j</sub>(z<sup>‑1</sup>)关于平移算子z<sup>‑1</sup>的系数;步骤3、无刷直流电动机控制系统包括速度环和电流环,所述速度环的控制为PI预测函数控制,其控制参数为步骤2计算所得的控制量u(k),所述电流环的控制为比例控制;PI预测函数中优化目标函数取下式:<img file="FDA0000804173030000022.GIF" wi="1205" he="154" /><b>e</b>(k+j)=<b>d</b>(k+j)‑<b>G</b><sub>j</sub>u(k),Δ<b>e</b>(k+j)=Δ<b>d</b>(k+j)‑Δ<b>G</b><sub>j</sub>u(k)其中,K<sub>p</sub>为比例增益,K<sub>i</sub>为积分系数,P为预测时域,M为控制时域,Q、 R为控制量加权因子;将控制量u(k)以可执行文件的形式加载到DSP的RAM中,DSP的捕获单元读取位置信号,计算得出电动机的转速,将转速与速度参考值比较得到速度偏差值,经过速度环PI预测函数控制器后得到电流参考值;将电流参考值与实际电流反馈值进行比较得到电流偏差值,再经过比例调节得到占空比可变的PWM信号,将PWM信号施加到电动机的功率驱动电路上,控制三相逆变桥的功率管的开断状态,从而控制无刷直流电动机的转速。 |