1.一种基于多变量互信息优选的质量预测与监测方法,其特征在于,包括以下步骤:首先,离线建模阶段包括如下所示步骤(1)至步骤(11);
步骤(1):从生产过程对象的历史数据库中找出能反映产品质量的指标所对应的数据组成输出矩阵Y∈Rn×k,与输出Y相对应的采样数据组成输入矩阵X∈Rn×m,其中,n为训练样本数,m为过程测量变量数,k为质量指标数,R为实数集,Rn×m表示n×m维的实数矩阵;
步骤(2):计算输出矩阵Y中各列向量的均值μ1,μ2,…,μk与标准差δ1,δ2,…,δk后,按照公式 对Y中各行向量实施标准化处理得到标准化后的输出矩阵 其中行向量y与 分别表示矩阵Y与 中的任意一个行向量,输出均值向量μ=[μ1,μ2,…,μk]、输出标准差对角矩阵 中对角线上的元素为δ1,δ2,…,δk;
步骤(3):对矩阵X实施标准化处理,得到标准化后的输入矩阵 初始化c=1并根据如下所示公式初始化一条长度为m的二进制数b:
步骤(4):根据二进制数b中元素1的位置,对应地从矩阵 中选取相应的列向量组成矩阵Z;
步骤(5):利用如下所示的多变量互信息计算过程,计算矩阵Z与矩阵 之间的相关性φc;
①将矩阵Z与矩阵 合并成一个矩阵 设置距离近邻个数r=3并初始化p=1;
②计算矩阵A中第p行向量与A中第q=1,2,…,n行向量之间的距离,然后按数值大小做升序排列,并记录第r+1个数值θp;
③计算矩阵Z中第p行向量与 中第q=1,2,…,n行向量之间的距离 然后记录距离大小满足条件 的个数βp;
④计算矩阵 中第p行向量与 中第q=1,2,…,n行向量之间的距离 然后记录距离大小满足条件 的个数δp;
⑤利用记录的βp与δp计算其相应的Digamma函数ψ(βp)与ψ(δp),其中Digamma函数满足条件ψ(c+1)=ψ(x)+1/x,x为函数自变量,当x=1时,ψ(1)=-0.57721;
⑥判断是否满足条件:p<n?若是,则置p=p+1后返回至②;若否,则分别计算得到ψ(β1),ψ(β2),…,ψ(βn)与ψ(δ1),ψ(δ2),…,ψ(δn)的平均值,记为 与ξ;
⑦按照下式计算矩阵Z与 之间的互信息大小:
那么,矩阵Z与矩阵 之间的相关性φc为
步骤(6):判断是否满足条件:c<2m-1+2m-2+…+20?若是,则置c=c+1与b=b+1后返回至步骤(4);若否,则得到C=2m-1+2m-2+…+20个相关性数值φ1,φ2,…,φC;
步骤(7):将φ1,φ2,…,φC中的最大值所对应的下标号c转换成m位的二进制数后,记录该二进制数中元素1所在位置为集合Φ,则Φ中优选出的变量即为与输出质量指标 密切相关的输出变量;
步骤(8):根据集合Φ中存储的变量位置从矩阵 中选出对应的列向量组成输出相关矩阵 而 中其余的列向量则组成输出不相关矩阵步骤(9):建立输出相关矩阵 与输出矩阵 之间的软测量模型: 其中表示软测量模型拟合出的函数关系,E为预测误差,根据应用对象的不同,可用算法有偏最小二乘回归、神经网络、支持向量回归;
步骤(10):根据公式 计算出软测量模型的输出估计值 并计算 的协方差矩阵 上标号T表示矩阵或向量的转置;
步骤(11):计算矩阵 的协方差矩阵 并根据公式 与计算出控制上限Dlim与Qlim;
其次,离线建模阶段完成后,当过程对象测量到新的样本数据时,即开始实施在线质量指标预测与监测,包括如下所示步骤(12)至步骤(16);
步骤(12):采集最新采样时刻的样本数据x∈R1×m,并对其实施步骤(3)中与矩阵X相同的标准化处理得到向量步骤(13):根据集合Φ从行向量 中选出相应的列组成行向量 而 中其余元素则组成行向量步骤(14):根据公式 计算出当前采样时刻的软测量模型的输出 并根据公式计算当前采样时刻的质量指标数据y;
步骤(15):根据如下所示公式计算与质量指标相关的监测统计量D,和与质量不相关的监测统计量Q:步骤(16):若D≤Dlim且Q≤Qlim,则当前工况运行正常,返回步骤(12)继续实施对下一个采样时刻的质量指标预测与监测;若D>Dlim,则过程对象运行出现了与质量相关的故障;若Q>Qlim,则过程对象运行出现了与质量不相关的故障。