一种盘磨系统的粉体粒度分布形状估计方法及其系统技术方案

本发明专利技术属于粉体粒度控制领域，具体涉及一种盘磨系统的粉体粒度分布形状估计方法及其系统，包括以下步骤：1)通过采集盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和粉体粒度分布形状的PDF；2)采用RBF神经网络对权值进行解耦计算；3)采用BP神经网络构建权值的非线性动态模型，由此获得权值向量的估计值；4)对粉体粒度分布形状的PDF进行估计。本发明专利技术应用于食品加工行业的粉体粒度分布形状的估计，同时也可用于冶金行业选矿过程的矿粉粒度、水泥生产过程的水泥粒度等复杂工业过程的粉体粒度分布形状的估计。

【技术实现步骤摘要】
一种盘磨系统的粉体粒度分布形状估计方法及其系统
本专利技术属于粉体粒度控制领域，具体涉及一种盘磨系统的粉体粒度分布形状估计方法及其系统。
技术介绍
在食品加工行业，面粉粒度常作为小麦粉质量等级的评价指标之一，必须符合一定的粒度要求，太大或太小都会影响生产指标要求和效率，通过对面粉粒度的控制，可以更好地把握面粉粉的品质变化，为后续的食品加工工序提供符合要求的原材料，从而提高生产效率及面粉质量。目前在实际生产中通常将粉体粒度的均值或者方差作为作为衡量面粉质量的重要指标，然而，均值或者方差并不能表征粉体粒度分布形状的全部信息，而粉体粒度分布形状包含了粉体的均值或者方差的所有信息。因此，将粉体粒度分布形状取代传统的均值或方差作为衡量面粉质量指标具有更为广泛的应用。在面粉生产过程中，影响盘磨系统的粉体粒度分布形状的因素有很多，如进料量、磨盘转速、磨盘间隙、磨盘压力等，而这些变量中，喂料量可以通过喂料器的振动频率进行调节，通过设置变频器频率可以有效控制带动动磨盘的主电机转速，磨盘间隙可采用机械或者液压调节装置完成等，因此，如何建立一个以进料量、磨盘转速、磨盘间隙、磨盘压力为输入变量，以粉体粒度分布形状为输出的盘磨系统的数学模型对粉体粒度分布形状的估计变得至关重要，和常规的建模方法不同的是，粉体粒度分布形状是一个空间分布，而非传统意义上的时序数据，因此，传统的建模方法(如机理建模或经验建模)将无法建立精确的刻画盘磨系统动态模型。因此，鉴于粉粒粒度分布形状在衡量面粉质量的重要性，如何实现对盘磨系统的粉体粒度分布形状估计，不但可以有效对小麦粉的品质变化进行估计，而且对降低能耗和提高生产效率提供一定的指导。
技术实现思路
针对上述存在的技术问题，本专利技术提供一种盘磨系统的粉体粒度分布形状估计方法，其特征在于，包括以下步骤：步骤1，采集盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和粉体粒度分布形状的PDF，即粉粒粒度分布形状的概率密度函数的数据；步骤2，采用具有高斯型激励函数的RBF神经网络，包括一组高斯基函数和相对应的权值，利用逆运算方法对不同时刻的与所述粉体粒度分布形状的PDF相对应的权值进行解耦计算；步骤3、采用BP神经网络构建权值的非线性动态模型，所述非线性动态模型的调节变量为当前时刻盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和所述权值，输出变量为下一时刻的所述权值，由此获得权值向量的估计值；步骤4、利用步骤3所述由BP神经网络构建的权值的非线性动态模型，结合步骤2所述RBF神经网络近似所述粉体粒度分布形状的PDF，得到粉体粒度分布形状的PDF的动态模型，用于对粉体粒度分布形状的PDF进行估计。所述步骤2包括以下步骤：步骤2-1，将所述进料量、磨盘压力、磨盘转速、磨盘间隙作为输入变量，所述粉体粒度分布的PDF作为输出变量；步骤2-2、采用具有高斯型激励函数的RBF神经网络近似所述粉体粒度分布形状的PDF，步骤2-3、利用逆运算方法对不同时刻粉体粒度分布形状的PDF相对应的权值进行解耦计算。所述步骤2.2中的近似方法为：设粉体粒度y为描述动态盘磨系统输出的一致有界随机过程变量，y∈[a,ξ]；在任一采样时刻k，u(k)为该时刻与所述输入变量对应的调节输入向量，u(k)∈Rm；在任一采样时刻k，粉体粒度分布形状可以通过其概率密度函数来描述，其定义式如下：式中，γ(y,u(k))为粉体粒度分布形状的PDF；P(a≤y＜ξ,u(k))表示盘磨系统在调节输入变量u(k)作用下的所述粉体粒度y落在区间[a,ξ]内的概率，本式表明粉体粒度分布形状的PDF由调节输入向量u(k)进行调节；设区间[a,b]为已知且连续且有界，采用如下具有高斯型激励函数的RBF神经网络来逼近输出所述PDF：首先，定义第l个网络节点的基函数Rl(y)如下式所示：式中，l表示第l个网络节点，n为网络节点总数，μl表示第l个网络节点的基函数的中心值，σl表示第l个网络节点的基函数的宽度；然后，利用RBF神经网络逼近原理，得出此时粉体粒度分布形状的PDF采用RBF神经网络的近似表示，即式中,wl(u(k))为与第l个网络节点的基函数相对应的权值；此外，由于粉体粒度分布形状的PDF需要满足如下自然隐含条件则此时粉体粒度分布形状的PDF近似表达式为γ(y,u(k))＝C(y)V(k)+L(y)式中V(k)＝[ω1(u(k)),ω2(u(k)),....,ωn-1(u(k))]T上式表明，一旦所述基函数Rl(y)确定，所述有界区间[a,b]已知，则在n个所述权值中有n－1个是相互独立的。所述步骤2.3中对所述权值进行解耦计算，具体方法为：将粉体粒度分布形状的PDF的近似表达式表示为如下向量形式对上式两边左乘[CT(y)Rn(y)]T，并对两边在所述区间[a,b]上进行积分，由此得到当矩阵非奇异时，上式所示权值可以通过矩阵求逆获得上式揭示了粉体粒度分布形状的PDF与所述权值之间关系，可以看出当所述基函数Rl(y)和粉体粒度分布形状的PDF数据可知时，便可获得相应的权值向量。所述步骤3中权值的非线性动态模型的构建方法为：在获得相应的权值向量之后，可以采用BP神经网络构建以盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙为输入变量，权值为输出变量的非线性动态模型，其动态可以采用如下表示V(k+1)＝f(V(k),u(k))式中，f(·)为盘磨系统的输入变量与权值之间的动态表达式，其采用BP神经网络进行建立。所述步骤4中所述动态模型为：当盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙获得后，可以采用所述动态模型对粉体粒度分布形状进行的在线连续估计。一种盘磨系统的粉体粒度分布形状估计方法的系统，所述系统采用上述的盘磨系统的粉体粒度分布形状估计方法；所述系统包括：数据采集单元、粉体粒度分布形状近似单元、权值解耦计算单元和粉体粒度分布形状估计单元；所述数据采集单元用于采集进料量、磨盘压力、磨盘转速、磨盘间隙、粉体粒度分布形状的PDF的数据；所述粉体粒度分布形状近似单元用于通过具有高斯型激励函数的RBF神经网络近似所述的粉体粒度分布的PDF；所述权值解耦计算单元用于利用逆运算对不同时刻粉体粒度分布形状的PDF相对应的权值进行解耦计算；所述粉体粒度分布形状估计单元用于将所述的当前时刻的进料量、磨盘压力、磨盘转速、磨盘间隙和权值作为输入变量，下一时刻的权值作为输出变量，利用BP神经网络构建权值的非线性动态模型，将BP神经网络构建的权值非线性动态模型结合RBF神经网络近似部分估计粉体粒度分布的PDF形状。所述的粉体粒度分布的PDF形状估计单元，包括粉体粒度分布形状的PDF近似模块和动态权值估计模块；所述粉体粒度分布形状的PDF近似模块用于通过高斯型激励函数的RBF神经网络近似粉体粒度分布的PDF形状；所述动态权值估计模块，用于将当前时刻的进料量、磨盘压力、磨盘转速、磨盘间隙和权值作为输入变量，下一时刻的权值作为输出变量，构建权值的BP神经网络的非线性动态模型。本专利技术的有益效果：本专利技术提出一种盘磨系统的粉体粒度分布形状估计方法及其系统，具体涉及到食品加工行业的粉体粒度分布形状的估计，同时也可用于冶金行业选矿过程的矿粉粒度、水泥生产过程的水泥粒度等复杂工业过程的粉体粒度分布形状的估计本文档来自技高网...

【技术保护点】
1.一种盘磨系统的粉体粒度分布形状估计方法，其特征在于，包括以下步骤：步骤1，采集盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和粉体粒度分布形状的PDF，即粉粒粒度分布形状的概率密度函数的数据；步骤2，采用具有高斯型激励函数的RBF神经网络，包括一组高斯基函数和相对应的权值，利用逆运算方法对不同时刻的与所述粉体粒度分布形状的PDF相对应的权值进行解耦计算；步骤3、采用BP神经网络构建权值的非线性动态模型，所述非线性动态模型的调节变量为当前时刻盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和所述权值，输出变量为下一时刻的所述权值，由此获得权值向量的估计值；步骤4、利用步骤3所述由BP神经网络构建的权值的非线性动态模型，结合步骤2所述RBF神经网络近似所述粉体粒度分布形状的PDF，得到粉体粒度分布形状的PDF的动态模型，用于对粉体粒度分布形状的PDF进行估计。

【技术特征摘要】
1.一种盘磨系统的粉体粒度分布形状估计方法，其特征在于，包括以下步骤：步骤1，采集盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和粉体粒度分布形状的PDF，即粉粒粒度分布形状的概率密度函数的数据；步骤2，采用具有高斯型激励函数的RBF神经网络，包括一组高斯基函数和相对应的权值，利用逆运算方法对不同时刻的与所述粉体粒度分布形状的PDF相对应的权值进行解耦计算；步骤3、采用BP神经网络构建权值的非线性动态模型，所述非线性动态模型的调节变量为当前时刻盘磨系统的进料量、磨盘压力、磨盘转速、磨盘间隙和所述权值，输出变量为下一时刻的所述权值，由此获得权值向量的估计值；步骤4、利用步骤3所述由BP神经网络构建的权值的非线性动态模型，结合步骤2所述RBF神经网络近似所述粉体粒度分布形状的PDF，得到粉体粒度分布形状的PDF的动态模型，用于对粉体粒度分布形状的PDF进行估计。2.根据权利要求1所述的盘磨系统的粉体粒度分布形状估计方法，其特征在于，所述步骤2包括以下步骤：步骤2-1，将所述进料量、磨盘压力、磨盘转速、磨盘间隙作为输入变量，所述粉体粒度分布的PDF作为输出变量；步骤2-2、采用具有高斯型激励函数的RBF神经网络近似所述粉体粒度分布形状的PDF，步骤2-3、利用逆运算方法对不同时刻粉体粒度分布形状的PDF相对应的权值进行解耦计算。3.根据权利要求2所述的盘磨系统的粉体粒度分布形状估计方法，其特征在于，所述步骤2.2中的近似方法为：设粉体粒度y为描述动态盘磨系统输出的一致有界随机过程变量，y∈[a,ξ]；在任一采样时刻k，u(k)为该时刻与所述输入变量对应的调节输入向量，u(k)∈Rm；在任一采样时刻k，粉体粒度分布形状可以通过其概率密度函数来描述，其定义式如下：式中，γ(y,u(k))为粉体粒度分布形状的PDF；P(a≤y＜ξ,u(k))表示盘磨系统在调节输入变量u(k)作用下的所述粉体粒度y落在区间[a,ξ]内的概率，本式表明粉体粒度分布形状的PDF由调节输入向量u(k)进行调节；设区间[a,b]为已知且连续且有界，采用如下具有高斯型激励函数的RBF神经网络来逼近输出所述PDF：首先，定义第l个网络节点的基函数Rl(y)如下式所示：式中，l表示第l个网络节点，n为网络节点总数，μl表示第l个网络节点的基函数的中心值，σl表示第l个网络节点的基函数的宽度；然后，利用RBF神经网络逼近原理，得出此时粉体粒度分布形状的PDF采用RBF神经网络的近似表示，即式中,wl(u(k))为与第l个网络节点的基函数相对应的权值；此外，由于粉体粒度分布形状的PDF需要满足如下自然隐含条件则此时粉体粒度分布形状的PDF近似表达式为γ(y,u(k))＝C(y)V(k)+L(y)式中V(k)＝[ω1(u(k)),ω2(u(k)),....,ωn-1(u(k))]T上式表明，一旦所述基函...

【专利技术属性】
技术研发人员：李明杰，周平，
申请(专利权)人：东北大学，
类型：发明
国别省市：辽宁,21