CN106059412B - 基于信度规则库推理的他励直流电动机转速控制方法 - Google Patents
基于信度规则库推理的他励直流电动机转速控制方法 Download PDFInfo
- Publication number
- CN106059412B CN106059412B CN201610321324.4A CN201610321324A CN106059412B CN 106059412 B CN106059412 B CN 106059412B CN 201610321324 A CN201610321324 A CN 201610321324A CN 106059412 B CN106059412 B CN 106059412B
- Authority
- CN
- China
- Prior art keywords
- value
- motor
- formula
- sample
- rule
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P7/00—Arrangements for regulating or controlling the speed or torque of electric DC motors
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于信度规则库推理的他励直流电动机转速控制方法。所建立的置信规则库,能够描述闭环控制系统偏差量和偏差量积分与控制量(输出量)之间存在的复杂非线性关系,其中置信规则的前项属性为输入量的参考值,后项为输出量的置信结构。对于被输入激活的规则,通过证据推理(ER)算法将被激活规则后项中的置信结构进行融合,从融合结果中推理出控制量估计值。给出非线性优化学习模型,解决专家给定的BRB初始参数不精确的问题。训练后BRB输出的控制量估计值,可以作为被控对象的输入,产生控制作用。
Description
技术领域
本发明涉及一种基于信度规则库推理的他励直流电动机智能控制方法,属于智能控制领域。
背景技术
在工业控制中,PID控制是应用最为广泛的一种控制规律,它具有结构简单、可靠性高和控制器参数易于调节等优点。但是,它是一种线性控制器,适用的被控对象是用传递函数描述的线性系统或近似线性系统。随着工业系统的结构与功能日益复杂,其中的被控对象的系统参数通常是随时间变化的、且这种变化有时是不确定的,这样就难以用准确的数学模型对其进行建模。在这种情况下,就必须通过有经验的工程师或是专家实时调整比例、积分和微分项的相关变量,以便取得较好控制效果。此时,偏差量、偏差量积分和偏差量微分与PID控制器输出的控制量之间就构成了非线性的映射关系,当系统参数的不确定性变化越频繁,则这种映射关系的非线性程度就越高。若能够基于这些专家经验和控制数据,建模这种复杂的非线性关系,从而设计出可以在计算机上运行的算法,则能大大提高控制效率,提升控制系统的智能化水平。
发明内容
本发明针对被控对象系统参数不确定性变化的情况,设计一种基于置信规则库(BRB)的他励直流电动机控制器。所建立的置信规则库,能够描述闭环控制系统偏差量和偏差量积分与控制量(输出量)之间存在的复杂非线性关系,其中置信规则的前项属性为输入量的参考值,后项为输出量的置信结构。对于被输入激活的规则,通过证据推理(ER)算法将被激活规则后项中的置信结构进行融合,从融合结果中推理出控制量估计值。同时给出非线性优化学习模型,解决专家给定的BRB初始参数不精确的问题。训练后BRB输出的控制量估计值,可以作为被控对象的输入,产生控制作用。
本发明提出一种基于信度规则库推理的他励直流电动机转速控制方法,该方法包括以下各步骤:
步骤(1)确定他励直流电动机传递函数,当被控对象输入为电枢电压UC(单位:V),输出为电动机的转速y(单位:r/min),在不加负载的情况下,被控对象的传递函数为:
上式中传递函数增益系数Ku=1/Ce,Ce(单位:V·s/rad)是电动机的电势常数,Ce∈[0,1],电动机的电磁时间常数Ta(单位:秒s)为:
这里La(单位:H)和Ra(单位:Ω)分别为电枢回路电感和电阻,La∈[0,1],Ra∈[0,10],电机时间常数Tm(单位:s)为:
这里J(单位:N·m·s2)为转动部分折算到电机轴上的总转动惯量,J∈[0,1],Cm(单位:N·m/A)是转矩常数,且Cm=(30/π)Ce;
步骤(2)确定变参数情况下他励直流电动机闭环控制系统的控制样本数据集其中,S中的每一个元素是一个三维向量 ep和分别表示闭环控制系统的偏差量和偏差量积分(输入量),UC P表示控制量,ANum表示采集总样本个数,100<ANum<+∞,样本数据集可以通过以下两种方式采集:
方式一:搭建基于PI控制器的他励直流电动机闭环控制系统,在电机运行的第1-4个月,20℃-25℃温度范围下,随机采集Num1个数据;在第5-8个月,80℃-100℃温度范围下,随机采集Num2个数据;在第9-12个月,-10℃-0℃温度范围下,随机采集Num3个数据,其中Num1、Num2和Num3都属于[100,+∞];此三组不同温度下样本数据总和组成集合S;当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%],其中上升时间tr是指响应曲线从稳态值的10%上升到90%所需的时间,调节时间ts是指在响应曲线的稳态线上,用稳态值的百分数(这里取2%)作一个允许误差范围,响应曲线达到并永远保持在这一允许误差范围内所需时间,σp是指最大超调量,其计算公式为:
式中c(tp)为阶跃响应的最大值,c(∞)为系统响应的稳态值;
方式二:利用Simulink仿真搭建基于PI控制器的他励直流电动机闭环控制系统,对被控对象传递函数中的电磁时间常数Ta和电势常数Ce施加扰动,模拟Ta和Ce的不确定性变化,两参数施加扰动后变为Ta(1+σai)和Ce(1+σei),其中,σai和σei分别是Ta和Ce的扰动量,σai∈[-5%,5%]和σei∈[-5%,5%],i=1,2,…,m,表示第i次扰动,共计扰动了m次,0<m<100,第i次扰动下对应的被控对象传递函数为:
式中π表示圆周率,取值保留四位有效数字;
在每种扰动情况下,获得Num个样本,100≤Num≤500,此时得到样本总数ANum=m×Num;同样当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%];
步骤(3)建立置信规则库(BRB),它能够描述闭环控制系统偏差量e和偏差量积分eI(输入量)与控制量UC(输出量)之间存在的复杂非线性关系,令e0,1=e,e0,2=eI,该BRB由如下形式的信度规则组成:
Rk:若则[(D1,β1,k),(D2,β2,k),…,(DN,βN,k)] (6)
式(6)中,为规则Rk中前项属性e0,i的参考值,且有Qi为的取值集合空间,其中的元素满足 表示采集到样本集S中前向属性e0,i最小值,表示采集到样本集S中前向属性e0,i最大值,mi表示参考值取值的个数mi≥1;从集合Q1,Q2中分别抽取m1和m2元素作为e0,1,e0,2参考值,共计可以产生L=m1×m2条规则,L≥1,k=1,2,3…,L为规则的编号;
Rk的规则权重为θk,满足0≤θk≤1,k=1,2,…,L,系统的偏差量e和偏差量积分eI作为规则Rk的前项属性,它们相应的属性权重分别为δ1,δ2,并有0≤δi≤1,i=1,2;
式(6)中,Rk后项属性分别为D1,D2,…,DN,并有Lu≤D1<D2<…<DN≤Ru,Lu表示采集到样本集S中UC最小值,表示采集到样本集S中UC最大值,N表示后项属性参考值个数,N≥2;β1,k,β2,k,…,βN,k分别为D1,D2,…,DN的信度值,并有0≤βa,k≤1,a∈{1,2…,N};
步骤(4)将S中某个样本中的ep,作为规则输入变量e,eI的取值,带入到BRB的每一条规则Rk中,得到输出结果为:
这里,是激活权重与后项初始置信度经信度融合后得到的信度值,计算公式为:
其中,
式(8)与式(9)中,wk为第k个输入ek,的第k条规则的激活权重,这里:
其中,wk∈[0,1];为相对属性权重,表达式为:
式(10)中,为第k条规则中第i个输入e0,i与相应参考值Ek,i的匹配度;这里Ek,i∈Ai,Ai={Ai,j|j=1,2,...,Ji}和Ai,1<Ai,2<…<Ai,Ji,的计算公式如下:
由式(7)最后算出控制量估计值为:
步骤(5)建立BRB的优化模型,将S中的ep,UC P样本作为训练样本,其中ep,作为步骤(4)中所构建BRB的输入,产生控制量估计值V表示BRB模型构成的参数向量V=(βa,k,θk,δi),k=1,2,…,L,a=1,2,…,N,i=1,2,则可以定义目标函数为:
式中,Q=ANum为步骤(2)中采集样本个数,UC,m=UC P,m=p=1,2,…,ANum;基于训练样本的优化过程,就是要通过调整BRB参数的取值,使得目标函数的值达到极小,由此可得模型的最优参数值;优化可通过Matlab中的非线性优化函数Fmincon实现;在式(15)-(18)给出的约束条件下,找到ξ(V)取最小值时,V中指标参数的最优值;
0≤βa,k≤1 (15)
0≤θk≤1 (17)
0≤δi≤1 (18)
具有最优指标参数取值的规则库为优化后的信度规则库;
有益效果:利用置信规则库(BRB)描述闭环控制系统偏差量和偏差量积分与控制量(输出量)之间存在的复杂非线性关系,给出非线性优化学习模型,解决专家给定的BRB初始参数不精确的问题。训练后BRB输出的控制量估计值,可以作为被控对象的输入,产生控制作用。
附图说明
图1是PI闭环控制系统。
图2是置信规则库控制器设计流程图。
图3是BRB闭环控制系统模型结构框图。
具体实施方式
本发明提出一种基于信度规则库推理的他励直流电动机转速控制方法,其特征在于该方法包括以下各步骤:
1.确定他励直流电动机传递函数
当被控对象输入为电枢电压UC(单位:V),输出为电动机的转速y(单位:r/min),在不加负载的情况下,被控对象的传递函数为:
上式中传递函数增益系数Ku=1/Ce,Ce(单位:V·s/rad)是电动机的电势常数,Ce∈[0,1]电动机的电磁时间常数Ta(单位:秒s)为:
这里La(单位:H)和Ra(单位:Ω)分别为电枢回路电感和电阻,La∈[0,1],Ra∈[0,10],电机时间常数Tm(单位:s)为:
这里J(单位:N·m·s2)为转动部分折算到电机轴上的总转动惯量,J∈[0,1],Cm(单位:N·m/A)是转矩常数,且Cm=(30/π)Ce。
为便于理解,这里举例说明,假设已知他励直流电动机初始额定电枢电压CH=220V,额定电枢电流IH=55A,额定转速nH=1000r/min,电枢电阻为R=2(Ω),转动惯量J=0.5N·m·s2,电磁时间常数Ta=0.017s和电势常数Ce=0.192V·s/rad,由式(2)得Tm=0.0473s,又Ku=1/Ce=5.2083,带入式(1)可得他励直流电机传递函数G(S)=5.2083/(0.000 8041s^2+0.0473s+1)。
2.确定变参数情况下他励直流电动机闭环控制系统的控制样本数据集
其中,S中的每一个元素是一个三维向量 ep和分别表示闭环控制系统的偏差量和偏差量积分(输入量),UC P表示控制量,ANum表示采集总样本个数,100<ANum<+∞,样本数据集可以通过以下两种方式采集:
方式一:搭建基于PI控制器的他励直流电动机闭环控制系统,在电机运行的第1-4个月,20℃-25℃温度范围下,随机采集Num1个数据;在第5-8个月,80℃-100℃温度范围下,随机采集Num2个数据;在第9-12个月,-10℃-0℃温度范围下,随机采集Num3个数据,其中Num1∈[100,+∞]、Num2∈[100,+∞]和Num3∈[100,+∞];此三组不同温度下样本数据总和组成集合S。当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%],其中上升时间tr是指响应曲线从稳态值的10%上升到90%所需的时间,调节时间ts是指在响应曲线的稳态线上,用稳态值的百分数(这里取2%)作一个允许误差范围,响应曲线达到并永远保持在这一允许误差范围内所需时间,σp是指最大超调量,其计算公式为:
式中c(tp)为阶跃响应的最大值,c(∞)为系统响应的稳态值。
方式二:利用Simulink仿真搭建基于PI控制器的他励直流电动机闭环控制系统,对被控对象传递函数中的电磁时间常数Ta和电势常数Ce施加扰动,模拟Ta和Ce的不确定性变化,两参数施加扰动后变为Ta(1+σai)和Ce(1+σei),其中,σai和σei分别是Ta和Ce的扰动量,σai∈[-5%,5%]和σei∈[-5%,5%],i=1,2,…,m,表示第i次扰动,共计扰动了m次,0<m<100,第i次扰动下对应的被控对象传递函数为:
式中π表示圆周率,取值保留四位有效数字。
在每种扰动情况下,获得Num个样本,100≤Num≤500,此时得到样本总数ANum=m×Num。同样当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%]。
为便于理解,这里举例说明,在电机参数不变情况下,传递函数还是以上文中给出的G1(S)=5.2083/(0.000 8041s^2+0.0473s+1)为例。但是由于他励直流电动机电枢电阻和电感存在容差,并且长时间运行会引起电机自身温度变化,以及电机负载变化引起的转动惯量的变化,这些因素都会导致电机中的电磁时间常数Ta和电势常数Ce会随时间发生不确定性变化。此例给出Ta和Ce分别在σai∈[-5%,5%]和σei∈[-5%,5%]范围内随机变化时,假设它随机变化值分别为σa=+5%,σe=+5%它的传递函数G2(S)=4.96/(0.000 7665s^2+0.04294s+1)。
将正常情况下传递函数G1(S)和参数加扰动后传递函数G2(S)分别放入图1被控对象中,调节Kp和KI使得产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%];每种情况下,获得Num=117个样本,此时得到样本总数ANum=m×Num=2×117=234。部分样本如表1所示。
表1部分采集样本
如图3所示,3.建立置信规则库(BRB)
能够描述闭环控制系统偏差量e和偏差量积分eI(输入量)与控制量UC(输出量)之间存在的复杂非线性关系,令e0,1=e,e0,2=eI,该BRB由如下形式的信度规则组成:
Rk:若则[(D1,β1,k),(D2,β2,k),…,(DN,βN,k)] (6)
式(6)中,为规则Rk中前项属性e0,i的参考值,且有为的取值集合空间,其中的元素满足 表示采集到样本集S中前向属性e0,i最小值,表示采集到样本集S中前向属性e0,i最大值,mi表示参考值取值的个数mi≥1;从集合Q1,Q2中分别抽取m1和m2元素作为e0,1,e0,2参考值,共计可以产生L=m1×m2条规则,L≥1,k=1,2,3…,L为规则的编号;
并有,Rk的规则权重为θk,满足0≤θk≤1,k=1,2,…,L,系统的偏差量e和偏差量积分eI作为规则Rk的前项属性,它们相应的属性权重分别为δ1,δ2,并有0≤δi≤1,i=1,2;
式(6)中,Rk后项属性分别为D1,D2,…,DN,并有Lu≤D1<D2<…<DN≤Ru,Lu表示采集到样本集S中UC最小值,表示采集到样本集S中UC最大值,N表示后项属性参考值个数,N≥2;β1,k,β2,k,…,βN,k分别为D1,D2,…,DN的信度值,并有0≤βa,k≤1,a∈{1,2…,N};
为便于理解,还以上文中给出为例,举例说明如何建立初始置信规则库,由步骤(2)确定变参数情况下他励直流电动机闭环控制系统的控制样本数据集对所获输入样本S的分析可以确定中每个分量的变化范围,输入-0.1888≤e≤10,0≤eI≤0.6255;输出为74≤UC≤204.1426。将偏差量设置m1=6个参考点,偏差量积分设置m2=7个参考点;输出控制量设置4个参考点,具体设置输入变量e、eI和输出变量UC的参考值(语义值)如表2所示。
表2输入(e,eI)与输出UC语义值与参考值
以上的语义值中P、N、L、M和S分别代表“正”、“负”、“大”、“中”和“小”。
进而,可以给出BRB系统中的置信规则为:
其中,Ek,i∈Ai,则共计可以生成L=m1×m2=42条规则,表3中给出了全部规则,其中的信度赋值βk,l为初始值,θk,δI,i=1,2初值均取1。
表3置信规则库的全部规则
4.将S中某个样本中的ep,作为规则输入变量e,eI的取值,带入到BRB的每一条规则Rk中,经推理得到输出结果为:
这里,是激活权重与后项初始置信度经信度融合后得到的信度值,计算公式为:
其中,
式(8)与式(9)中,wk为第k个输入ek,的第k条规则的激活权重,这里:
其中,wk∈[0,1];为相对属性权重,表达式为:
式(10)中,为第k条规则中第i个输入e0,i与相应参考值Ek,i的匹配度;这里Ek,i∈Ai,Ai={Ai,j|j=1,2,...,Ji}和Ai,1<Ai,2<…<Ai,Ji,的计算公式如下:
由式(7)最后算出控制量估计值为:
为便于理解,还以上文中给出为例,假设初始BRB的输入为表1第2个样本,即当S=(e,eI)=(9.9792,0.0300)时,偏差量e∈[7,10]=[PM2,PL],偏差量积分eI∈[0,0.1860]=[Z,PS1]可知由式(12)可知激活了初始BRB四条规则分别是表3第29条规则PM2AND Z、第30条规则PM2AND PS1、第36条规则PL AND Z和第37条规则PL AND PS1。由式(10)可得规则权重分别为w29=0.0058,w30=0.0011,w36=0.8331,w37=0.1600,直观上可知样本点最接近第36规则,所以第36条激活的权重最大。再由式(8)-(9)的到ER推理后的输出信度结构:β1=0.9633,β2=0.0349,β3=0.0019,β4=0。最后由式(14)可得控制量的估计值:
按上述步骤计算整体误差
5.建立BRB的优化模型,将S中的ep,样本作为训练样本,其中ep,作为步骤(4)中所构建BRB的输入,产生控制量估计值V表示BRB模型构成的参数向量V=(βa,k,θk,δi),k=1,2,…,L,a=1,2,…,N,i=1,2,则可以定义目标函数为:
式中,Q=ANum为步骤(2)中采集样本个数,ANum;基于训练样本的优化过程,就是要通过调整BRB参数的取值,使得目标函数的值达到极小,由此可得模型的最优参数值;优化可通过Matlab中的非线性优化函数Fmincon实现;在式(15)-(18)给出的约束条件下,找到ξ(V)取最小值时,V中指标参数的最优值;
0≤βa,k≤1 (15)
0≤θk≤1 (17)
0≤δi≤1 (18)
具有最优指标参数取值的规则库为优化后的信度规则库;
为便于理解,还以上文中给出为例,用UC,m=UC p,m=1,2,…,ANum表示第m个训练样本输入相对应的实际输出;表示第m个输入训练样本进行BRB模型推理的控制量估计值,目标函数为式(14),约束条件为(15)-(18)。通过Matlab中的非线性优化函数Fmincon实现得到表4优化后置信规则库,优化后的属性权重δ1=1,δ2=0.8143。
表4优化后置信规则库
以下结合附图,详细介绍本发明方法的实施例:
本发明方法的流程图如图2所示,核心部分是:利用置信规则库(BRB)描述闭环控制系统偏差量和偏差量积分与控制量(输出量)之间存在的复杂非线性关系,给出非线性优化学习模型,解决专家给定的BRB初始参数不精确的问题。训练后BRB输出的控制量估计值,可以作为被控对象的输入,产生控制作用。
以下结合模型他励直流电机为例。详细介绍本发明方法的各个步骤,并通过实验结果验证置信规则库(BRB)控制器相比传统PID控制器优越性。
1.确定他励直流电动机传递函数
已知他励直流电动机初始额定电枢电压CH=220V,额定电枢电流IH=55A,额定转速nH=1000r/min,电枢电阻为R=2(Ω),转动惯量J=0.5N·m·s2,电磁时间常数Ta=0.017s和电势常数Ce=0.192V·s/rad当被控对象输入为电枢电压UC(单位:V),输出为电动机的转速y(单位:r/min),在不加负载的情况下,被控对象的传递函数为:
2.利用Simulink仿真搭建图1所示基于PI控制器的他励直流电动机闭环控制系统,对被控对象传递函数中的电磁时间常数Ta和电势常数Ce施加扰动,模拟Ta和Ce的不确定性变化,两参数施加扰动后变为Ta(1+σai)和Ce(1+σei),其中,σai和σei分别是Ta和Ce的扰动量,σai∈[-5%,5%]和σei∈[-5%,5%],共计扰动了9次对应的被控对象传递函数如表5所示。
表5参数不确定变化的他激直流电动机传递函数
将表5中9组传递函数(包括正常和加扰动)分别放入图1被控对象中,调节Kp和KI使得产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0,2%];每种情况下,获得Num=117个样本,此时得到样本总数ANum=m×Num=9×117=1053。部分样本如表6所示。
表6(a)部分1-234采集样本
表6(b)部分235-468采集样本
表6(c)部分469-702采集样本
表6(d)部分703-936采集样本
表6(e)部分937-1053采集样本
3.建立置信规则库(BRB)
由步骤(2)确定变参数情况下他励直流电动机闭环控制系统的控制样本数据集对所获输入样本的分析S可以确定中每个分量的变化范围,输入-0.1920≤e≤10,0≤eI≤0.6728;输出为74≤UC≤204.1426。将偏差量设置m1=8个参考点,偏差量积分设置m2=9个参考点;输出控制量设置5个参考点,具体设置输入变量e、eI和输出变量UC的参考值(语义值)如表7所示。
表7输入(e,eI)与输出UC语义值与参考值
以上的语义值中P、N、L、M和S分别代表“正”、“负”、“大”、“中”和“小”。
进而,可以给出BRB系统中的置信规则为:
其中,Ek,i∈Ai,则共计可以生成L=m1×m2=72条规则,表8中给出了全部规则,其中的信度赋值βk,l为初始值,θk,δI,i=1,2初值均取1。
表8置信规则库的全部规则
4.将S中各个样本中的ep,作为规则输入变量e,eI的取值,带入到BRB的每一条规则Rk中,经推理得到控制量估计输出结果如表9所示。
表9(a)由ER推理得到1-468部分控制量估计值
表9(b)由ER推理得到469-936部分控制量估计值
表9(c)由ER推理得到937-1053部分控制量估计值
5.建立BRB的优化模型
由式(14)计算控制量观测值与估计值均方误差:
此时控制量作用于被控对象不能满足控制效果。因此利用步骤(5)非线性优化模型优化BRB参数,使得目标函数的值达到极小,由此可得模型的最优参数值。优化可通过Matlab中的非线性优化函数Fmincon实现,得到表10优化后置信规则库,优化后的属性权重δ1=1,δ2=0.7719。
表10优化后置信规则库
Claims (3)
1.基于信度规则库推理的他励直流电动机转速控制方法,其特征在于,该方法包括以下各步骤:
步骤(1)确定他励直流电动机传递函数,当被控对象输入为电枢电压UC,输出为电动机的转速y,在不加负载的情况下,被控对象的传递函数为:
上式中传递函数增益系数Ku=1/Ce,Ce是电动机的电势常数,电动机的电磁时间常数Ta为:
这里La和Ra分别为电枢回路电感和电阻,La∈[0,1],Ra∈[0,10],电机时间常数Tm为:
这里J为转动部分折算到电机轴上的总转动惯量,Cm是转矩常数,且Cm=(30/π)Ce;
步骤(2)确定变参数情况下他励直流电动机闭环控制系统的控制样本数据集其中,S中的每一个元素是一个三维向量 ep和分别表示闭环控制系统的偏差量和偏差量积分,UC P表示控制量,ANum表示采集总样本个数,100<ANum<+∞;
步骤(3)建立置信规则库BRB,它能够描述闭环控制系统偏差量e和偏差量积分eI与控制量UC之间存在的复杂非线性关系,令e0,1=e,e0,2=eI,该BRB由如下形式的信度规则组成:
Rk:若且则[(D1,β1,k),(D2,β2,k),…,(DN,βN,k)] (6)
式(6)中,为规则Rk中前项属性e0,i的参考值,且有Qi为的取值集合空间,其中的元素满足 表示采集到样本集S中前项属性e0,i最小值,表示采集到样本集S中前项属性e0,i最大值,mi表示参考值取值的个数mi≥1;从集合Q1,Q2中分别抽取m1和m2元素作为e0,1,e0,2参考值,共计可以产生L=m1×m2条规则,L≥1,k=1,2,3…,L为规则的编号,i=1,2;
Rk的规则权重为θk,满足0≤θk≤1,k=1,2,…,L,系统的偏差量e和偏差量积分eI作为规则Rk的前项属性,它们相应的属性权重分别为δ1,δ2,并有0≤δ1≤1,0≤δ2≤1;
式(6)中,Rk后项属性分别为D1,D2,…,DN,并有Lu≤D1<D2<…<DN≤Ru,Lu表示采集到样本集S中UC最小值,Ru表示采集到样本集S中UC最大值,N表示后项属性参考值个数,N≥2;β1,k,β2,k,…,βN,k分别为D1,D2,…,DN的信度值,并有0≤βa,k≤1,a∈{1,2…,N};
步骤(4)将S中某个样本中的ep,作为规则输入变量e,eI的取值,带入到BRB的每一条规则Rk中,得到输出结果为:
这里,是激活权重与给定初始置信度经信度融合后得到的信度值,计算公式为:
其中,
式(8)与式(9)中,wk为第k个输入ek,的第k条规则的激活权重,这里:
其中,wk∈[0,1];为相对属性权重,表达式为:
式(10)中,为第k条规则中第i个输入e0,i与相应参考值Ek,i的匹配度;这里Ek,i∈Ai,Ai={Ai,j|j=1,2,...,Ji}和的计算公式如下:
由式(7)最后算出控制量估计值为:
步骤(5)建立BRB的优化模型,将S中的ep,UC P样本作为训练样本,其中ep,作为步骤(4)中所构建BRB的输入,产生控制量估计值V表示BRB模型构成的参数向量V=(βa,k,θk,δi),k=1,2,…,L,a=1,2,…,N,则可以定义目标函数为:
式中,Q=ANum为步骤(2)中采集样本个数,m=p=1,2,…,ANum;基于训练样本的优化过程,就是要通过调整BRB参数的取值,使得目标函数的值达到极小,由此可得模型的最优参数值;优化可通过Matlab中的非线性优化函数Fmincon实现;在式(15)-(18)给出的约束条件下,找到ξ(V)取最小值时,V中指标参数的最优值;
0≤βa,k≤1 (15)
0≤θk≤1 (17)
0≤δi≤1 (18)
具有最优指标参数取值的规则库为优化后的信度规则库。
2.根据权利要求1所述的基于信度规则库推理的他励直流电动机转速控制方法,其特征在于,所述的样本数据集采集为:
搭建基于PI控制器的他励直流电动机闭环控制系统,在电机运行的第1-4个月,20℃-25℃温度范围下,随机采集Num1个数据;在第5-8个月,80℃-100℃温度范围下,随机采集Num2个数据;在第9-12个月,-10℃-0℃温度范围下,随机采集Num3个数据,其中Num1、Num2和Num3都属于[100,+∞];此三组不同温度下样本数据总和组成集合S;当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%],其中上升时间tr是指响应曲线从稳态值的10%上升到90%所需的时间,调节时间ts是指在响应曲线的稳态线上,用稳态值的百分数2%作一个允许误差范围,响应曲线达到并永远保持在这一允许误差范围内所需时间,σp是指最大超调量,其计算公式为:
式中c(tp)为阶跃响应的最大值,c(∞)为系统响应的稳态值。
3.根据权利要求1所述的基于信度规则库推理的他励直流电动机转速控制方法,其特征在于,所述的样本数据集采集为:
利用Simulink仿真搭建基于PI控制器的他励直流电动机闭环控制系统,对被控对象传递函数中的电磁时间常数Ta和电势常数Ce施加扰动,模拟Ta和Ce的不确定性变化,两参数施加扰动后变为Ta(1+σai)和Ce(1+σei),其中,σai和σei分别是Ta和Ce的扰动量,σai∈[-5%,5%]和σei∈[-5%,5%],i=1,2,…,m,表示第i次扰动,共计扰动了m次,0<m<100,第i次扰动下对应的被控对象传递函数为:
式中π表示圆周率,取值保留四位有效数字;
在每种扰动情况下,获得Num个样本,100≤Num≤500,此时得到样本总数ANum=m×Num;同样当这些样本中的UC P作用于被控对象后,产生的输出要满足以下控制系统的稳态和动态性能指标取值:tr∈[0,0.1],ts∈[0,0.15]和σp∈[0%,2%];tr上升时间,是指响应曲线从稳态值的10%上升到90%所需的时间,ts指在响应曲线的稳态线上,用稳态值的百分数2%作一个允许误差范围,响应曲线达到并永远保持在这一允许误差范围内所需时间;σp是指最大超调量。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2016100599337 | 2016-01-28 | ||
CN201610059933 | 2016-01-28 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106059412A CN106059412A (zh) | 2016-10-26 |
CN106059412B true CN106059412B (zh) | 2018-11-13 |
Family
ID=57177670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610321324.4A Active CN106059412B (zh) | 2016-01-28 | 2016-05-12 | 基于信度规则库推理的他励直流电动机转速控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106059412B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443766A (zh) * | 2018-09-10 | 2019-03-08 | 中国人民解放军火箭军工程大学 | 一种重载车辆变速箱齿轮安全性分析方法 |
CN109507876B (zh) * | 2019-01-25 | 2021-12-28 | 杭州电子科技大学 | 一种基于信度推理的电推船舶电机pid参数整定方法 |
CN110084518A (zh) * | 2019-04-29 | 2019-08-02 | 杭州电子科技大学 | 一种基于混合输入置信规则推理模型的边坡稳定性评判方法 |
CN110850206A (zh) * | 2019-11-13 | 2020-02-28 | 武汉理工大学 | 基于置信规则推理的逆变器电容老化故障诊断方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955580A (zh) * | 2014-05-04 | 2014-07-30 | 杭州电子科技大学 | 基于信度规则库推理的集成电路参数成品率估计方法 |
CN105139086A (zh) * | 2015-08-13 | 2015-12-09 | 杭州电子科技大学 | 基于优化置信规则推理的轨道高低不平顺幅值估计方法 |
CN105243393A (zh) * | 2015-10-27 | 2016-01-13 | 长春工业大学 | 一种基于特征的复杂机电系统故障预报方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101927246B1 (ko) * | 2012-12-12 | 2019-03-12 | 한국전자통신연구원 | 전동기 위치 검출부 및 브러쉬리스 전동기 시스템 |
-
2016
- 2016-05-12 CN CN201610321324.4A patent/CN106059412B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955580A (zh) * | 2014-05-04 | 2014-07-30 | 杭州电子科技大学 | 基于信度规则库推理的集成电路参数成品率估计方法 |
CN105139086A (zh) * | 2015-08-13 | 2015-12-09 | 杭州电子科技大学 | 基于优化置信规则推理的轨道高低不平顺幅值估计方法 |
CN105243393A (zh) * | 2015-10-27 | 2016-01-13 | 长春工业大学 | 一种基于特征的复杂机电系统故障预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106059412A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106059412B (zh) | 基于信度规则库推理的他励直流电动机转速控制方法 | |
Mohideen et al. | Real-coded Genetic Algorithm for system identification and tuning of a modified Model Reference Adaptive Controller for a hybrid tank system | |
CN104378038B (zh) | 基于人工神经网络的永磁同步电机参数辨识方法 | |
US10466659B2 (en) | Method for designing a non-linear controller for non-linear processes | |
Emhemed et al. | Modelling and simulation for industrial DC motor using intelligent control | |
CN106026822A (zh) | 伺服电机驱动系统的惯量在线辨识方法及系统 | |
CN106597840B (zh) | 一种基于产生式规则推理的pid参数整定方法 | |
CN108696210A (zh) | 基于参数辨识的直流电机电流环控制器参数自整定方法 | |
CN110716430B (zh) | 一种采用等效扰动补偿的伺服系统快速吸引重复控制方法 | |
CN104270055A (zh) | 一种交流伺服系统位置环ip控制器的参数自校正方法 | |
CN106227964B (zh) | 基于扩张状态观测器的迟滞非线性系统参数辨识方法 | |
CN108762083A (zh) | 一种基于加速度观测器的自动控制系统 | |
CN109167546B (zh) | 一种基于数据生成模型的异步电机参数在线辨识方法 | |
CN105867138A (zh) | 一种基于pid控制器的稳定平台控制方法及装置 | |
Liu et al. | $ H_\infty $ Control for a Class of Discrete-Time Systems via Data-Based Policy Iteration With Application to Wind Turbine Control | |
Prusty et al. | A novel fuzzy based adaptive control of the four tank system | |
CN102707617A (zh) | 一种模糊PID算法ActiveX控件的实现方法 | |
CN109743002B (zh) | 伺服系统控制器、前馈控制信号确定方法、惯量辨识方法 | |
Maraba et al. | Speed control of an asynchronous motor using pid neural network | |
CN112436166A (zh) | 质子交换膜燃料电池空气供给系统的滑模控制方法 | |
CN106371321A (zh) | 一种焦化炉炉膛压力系统模糊网络优化pid控制方法 | |
CN108983618A (zh) | 基于凸多胞形的pmsm鲁棒h∞输出反馈控制设计方法 | |
Perdukova et al. | A model-based fuzzy control of an induction motor | |
CN113283006A (zh) | 基于rbf神经网络控制燃料电池电堆多参数响应方法 | |
CN106444389A (zh) | 一种废塑料裂解温度系统的模糊rbf网络优化pi控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210506 Address after: 210001 room 7032, No.11 Zhanyuan Road, Qinhuai District, Nanjing City, Jiangsu Province Patentee after: Nanjing Smart Water Transportation Technology Co.,Ltd. Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang Patentee before: HANGZHOU DIANZI University |
|
TR01 | Transfer of patent right |