[发明专利]基于流模型的机器组负载预测控制方法

权利要求书

1、基于流模型的机器组负载预测控制方法，其特征在于，所述方法是在机器组负载预测控制计算机上依次按以下步骤实现的：

步骤(1)：初始化，设定以下参变量

采样时间间隔，每隔时间间隔T给出前道工序中各机器组任务输出率，所述机器组由加工能力相似的多台机器组成，而采样时刻则用k表示；

机器组加工能力，是在单位时间内机器组所能完成的加工任务的加工时间总和，前道工序中机器组i的加工能力用u_i表示，i＝1，…，m，用矩阵表示为U＝[u₁，u₂，…，u_m]^T，后道工序中机器组j的加工能力用v_j表示，j＝1，…，n，用矩阵表示为V＝[v₁，v₂，…，v_n]^T，m和n分别为前、后道工序中机器组的数目；

前道工序机器组任务输出率，前道工序机器组i基于生产工艺约束在k采样时刻加工完成的任务被安排到后道工序由机器组j加工的比例，用c_ij(k)表示，用矩阵表示为：

C(k)n×m=[C1(k),C2(k),···,Cn(k)]T]]>

=[C1(k),C2(k),···,Cm(k)]]]>

=c11(k)c21(k)···cm1(k)c12(k)c22(k)···cm2(k)············c1n(k)c2n(k)···cmn(k)]]>

其中0≤c_ij≤1且Σj=1ncij=1,]]>i＝1，…，m；j＝1，…，n；

机器组负载，工序中某个机器组前等待加工任务的加工时间总和，后道工序机器组j在k时刻的负载表示为y_j(k)，用矩阵表示为Y(k)^n×1＝[y₁(k)，y₂(k)，…，y_n(k)]^T；

后道工序中机器组j在k时刻的机器负载期望值表示为y_j^r(k)，用矩阵表示为Yr(k)n×1=[y1r(k),y2r(k),···,ynr(k)]T;]]>

中间工序的加工延迟时间用d表示，其是指任务从前道工序加工完毕开始，通过中间工序的加工，到达后道工序所用的平均时间单位，含等待时间；

给定的控制周期用T_all表示；

控制周期T_all内后道工序机器组j的总负载用Load_j表示；

步骤(2)：用机器组负载信息采集装置采集机器组负载实时信息，机器组负载信息采集装置由PLC采集装置、嵌入式系统采集装置、DCS系统采集装置中的一种或它们的组合构成；

步骤(3)：所述的机器组负载预测控制计算机从所述采集装置中读取机器组负载实时信息，依次按以下步骤进行机器组负载预测控制：

步骤(3.1)：按下式确定各采样间隔时间内后道工序各机器组负载期望值y_j^r(k)，

yjr(k)=LoadjTall/T,]]>j＝1，…，n

步骤(3.2)：按下式建立后道工序机器组负载预测控制问题：

J(k+d+1)=12[Y(k+d+1)-Yr(k+d+1)]T·][Y(k+d+1)-Yr(k+d+1)]]]>

+12Σi=1mCi(k)T·Ci(k)]]>

求控制律C(k)，使得min{J(k+d+1)}，其中对所有i满足ACⁱ(k)＝1，A＝[1，1，…，1]^1×m；

步骤(3.3)：按下述步骤建立后道工序机器组负载预测控制模型

步骤(3.3.1)：建立后道工序机器组负载1阶预测控制模型

y_j(k+1)＝max{y_j(k)-v_jT+[c_1j(k-d)u₁+…+c_mj(k-d)u_m]T，0}

＝max{y_j(k)-Tv_j+TC_j(k-d)U，0}，j＝1，…，n    ，

步骤(3.3.2)：采用具有L条模糊规则的自适应神经模糊推理系统ANFIS建立后道工序机器组负载d阶预测控制模型

即：

Y(k+d+1)n×1=Xn×1+Σq=0dNqn×n·[C1(k-d+q),···,Cm(k-d+q)]n×m·Um×1]]>

=Xn×1+Σq=0dNqn×n·Σi=1muiCi(k-d+q)]]>

其中：X^n×1＝[E₁，E₂，…，E_n]^T，Ej=Σl=1Lhl(k)·[αl(yj(k)-Tvj)+βl]]]>

q＝0，1，…，d，Fjq=T·Σl=1Lhl(k)γql,]]>

h_l(k)＝f_l(y_j(k)-Tv_j，TC_j(k-d)U，…，TC_j(k-1)U)，是ANFIS输入为

y_j(k)-Tv_j，TC_j(k-d)U，…，TC_j(k-1)U时第l条模糊规则的激发值，f_l为模糊规则激发值计算函数，其为第l条模糊规则中ANFIS各输入变量所对应模糊数的隶属度函数的乘积，其中，上述模糊数的隶属度函数采用钟形函数，钟形函数中的参数即为ANFIS中待定的前件参数，而ANFIS中待定的后件参数为α^l，γ₀^l，γ₁^l，…，γ_d^l，β^l；

步骤(3.3.3)：采用如下步骤确定ANFIS中的前件参数和后件参数；

步骤(3.3.3.1)：对所有j＝1，2，…，n，随机产生C_j(k)的值，并按步骤(3.3.1)所述公式计算y_j(k+1)值，从而产生若干用于训练ANFIS的输入输出训练数据集，其中，输入数据为Y(k)，C(k-d)，…，C(k-1)，输出数据为Y(k+d+1)；

步骤(3.3.3.2)：采用步骤(3.3.2.1)生成的训练数据集及ANFIS经典学习算法，确定ANFIS中待定的前件参数和后件参数；

步骤(3.4)：根据步骤(3.3)得到的后道工序机器组负载d阶预测控制模型，采用拉格朗日松弛方法按下式计算最佳控制率C(k)，其中Cⁱ(k)为：

Ci(k)=ui(ATAAAT-I)·[I-(Σi=1mui2)·N‾(ATAAAT-I)]-1·[X‾+(Σi=1mui)·N‾ATAAT]+ATAAT,]]>

其中：I为单位矩阵，X‾=NdT[X+Σq=0d-1(Σi=1muiCi(k-d+q))-Yr(k+d+1)],]]>N‾=NdTNd.]]>