[发明专利]一种球磨机负荷参数软测量方法

权利要求书

1.一种球磨机负荷参数软测量方法，其特征在于，包括以下步骤：

S100、将球磨机负荷参数定义为标签，记录已知标签、采集已知标签工况下的辅助变量、未知标签工况下的辅助变量，已知标签和已知标签工况下的辅助变量形成的建模数据所在领域作为源域，未知标签和未知标签工况下的辅助变量形成的实时数据所在领域作为目标域；

S200、对采集的已知标签工况下的辅助变量和未知标签工况下的辅助变量进行快速傅里叶变换预处理，分别得到源域数据Xs和目标域数据Xt，将源域数据Xs和目标域数据Xt进行组合，得到输入数据X＝[Xs,Xt]，对源域中标签Ys进行单位化处理得到标签Ysa；

S300、对X和Ysa用联合分布适配方法降维处理得到新的特征表示Z，将Z划分为源域特征Zs和目标域特征Zt，对源域中的标签Ys进行归一化处理得到标签Ysy，用Zs和Ysy建立多任务-最小二乘支持向量机模型；

S400、对目标域中标签Yt进行归一化处理得到归一化参数opt1，用多任务-最小二乘支持向量机模型对Zt预测，将预测值用opt1反归一化并与Yt比较得到误差。

2.根据权利要求1所述的一种球磨机负荷参数软测量方法,其特征在于:所述标签为料球比、磨矿浓度和充填率的一种或多种。

3.根据权利要求1所述的一种球磨机负荷参数软测量方法,其特征在于:所述辅助变量为球磨机筒壁的振动信号、球磨机轴承上的振动信号和球磨机滚筒周围的振声信号的一种或多种。

4.根据权利要求1所述的一种球磨机负荷参数软测量方法,其特征在于:所述步骤S300包括：

S310、对所述输入数据X和所述标签Ysa用联合分布适配方法进行降维处理，求解以下公式的k个最小的特征向量，即得到新的特征表示Z＝ATX，且k＜m，其中A是正交变换矩阵，

( X Σ c = 0 C M c X T + η I ) A = XHX T A Φ ]]>

其中，n＝ns+nt，ns表示源域样本数量，nt表示目标域样本数量；是中心矩阵，其中，1是n×n的全1矩阵；η是保证最优化问题明确定义的正则化参数；是A的k个最大特征；

S320、特征表示Z区分为源域特征Zs和目标域特征Zt，其中Zs＝Z(1:k,1:ns)，表示取Z中的前ns列，Zt＝Z(1:k,ns+1:ns+nt)，表示取Z中的后nt列，n＝ns+nt；

S330、对源域数据标签Ys1,Ys2,Ys3分别进行归一化处理得到Ysy1,Ysy2,Ysy3；

S340、以Zs作为以下回归函数方程的输入数据X，Ysy1,Ysy2,Ysy3作为多任务标签，建立多任务-最小二乘支持向量机模型，并优化参数λ和γ，

f i ( x ) = Σ l = 1 m Σ j = 1 n l α l , j K ( x l , j , x ) + m λ Σ j = 1 n i α l , j K ( x l , j , x ) + b j , ]]>

其中λ、γ为正则化因子。