发明名称 一种基于ARX模型和PI预测函数的无刷直流电动机控制方法
摘要 本发明公开了一种基于ARX模型和PI预测函数的无刷直流电动机控制方法,通过PI与预测函数相结合的控制方法对预测函数进行改进。控制量增加了两个可调参数,使得控制调节更加灵活,通过调节两个参数可改善无刷直流电动机的转速响应性能。本发明在相同的实验条件下,PI预测函数控制比预测函数控值上升快,快速到达稳定状态,且静态误差小。PI预测函数控制可根据实际情况调节<i>K<sub>p</sub></i>、<i>K<sub>i</sub></i>两个参数,可达到良好的控制效果。与包含多个基函数的预测函数控制的计算量相比较,PI预测函数算法计算量小,控制灵活方便,无需重新修改控制程序,只需调节两个参数,同时解决了预测函数模型不匹配导致的控制效果变差的问题。
申请公布号 CN103633912B 申请公布日期 2016.01.13
申请号 CN201310582158.X 申请日期 2013.11.18
申请人 南京信息工程大学 发明人 郭伟;陈一帆;徐金成;周丽
分类号 H02P21/00(2006.01)I 主分类号 H02P21/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 一种基于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信号施加到电动机的功率驱动电路上,控制三相逆变桥的功率管的开断状态,从而控制无刷直流电动机的转速。
地址 215101 江苏省苏州市吴中区木渎镇中山东路70号吴中科技创业园2号楼2310室