CN101017373A - 工业过程多分辨率软测量仪表及方法 - Google Patents

工业过程多分辨率软测量仪表及方法 Download PDF

Info

Publication number
CN101017373A
CN101017373A CNA2006101555563A CN200610155556A CN101017373A CN 101017373 A CN101017373 A CN 101017373A CN A2006101555563 A CNA2006101555563 A CN A2006101555563A CN 200610155556 A CN200610155556 A CN 200610155556A CN 101017373 A CN101017373 A CN 101017373A
Authority
CN
China
Prior art keywords
data
formula
multiresolution
module
sigma
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.)
Pending
Application number
CNA2006101555563A
Other languages
English (en)
Inventor
刘兴高
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CNA2006101555563A priority Critical patent/CN101017373A/zh
Publication of CN101017373A publication Critical patent/CN101017373A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Complex Calculations (AREA)

Abstract

一种工业过程多分辨率软测量仪表,包括与工业对象连接的现场智能仪表、用于存放历史数据的数据存储装置及上位机,智能仪表、数据存储装置及上位机依次相连,所述的上位机为软测量智能控制器,所述的软测量智能控制器包括标准化处理模块、偏最小二乘模块、多分辨率分解模块、径向基函数神经网络建模模块、多分辨率重构模块、信号采集模块以及软测量模块。以及提出了一种软测量方法。本发明提供一种同时考虑复共线性和非线性特性、多分辨率特性、能够得到良好的软测量效果的工业过程多分辨率软测量系统及方法。

Description

工业过程多分辨率软测量仪表及方法
(一)技术领域
本发明涉及工业过程软测量领域,特别地,涉及一种工业过程多分辨率软测量仪表和方法。
(二)背景技术
在现代流程工业中,大量关键性过程状态、产品质量等参数缺乏在线直接测量手段。这已成为制约生产安全、产品质量、产量及生产效益进一步提高的瓶颈。软测量技术正式解决此类问题的有效途径。
利用工业实测数据,采用统计的方法建立工业过程的软测量模型,避开了复杂的机理分析,模型对观测数据的拟和程度高,求解相对方便,是软测量建模的热点。常用的统计建模方法有系统辨识、回归分析、人工神经元网络、统计学习理论等。
但是目前常用的软测量仪表及方法,往往只考虑了数据的复共线性和非线性特性,而没有考虑到数据的多分辨率特性,对于复杂工业过程往往难以得到较好的软测量结果。而工业生产过程从本质上来说是具有多分辨率特性的:工业实测数据既有高频信号,又有低频信号;微观、宏观反应机理的多分辨率特性;过程变量总是在不同采样率下得到的;相应的控制或操作也是在不同分辨率下发生的。多分辨率分析可以充分利用数据中所包含的信息,在此基础上建立的模型更贴近实际工业过程,可以得到更好的软测量效果。
(三)发明内容
为了克服已有的工业过程软测量系统的没有考虑数据的多分辨率特性、不能得到良好的软测量效果的不足,本发明提供一种同时考虑复共线性和非线性特性、多分辨率特性、能够得到良好的软测量效果的工业过程多分辨率软测量系统及方法。
本发明解决其技术问题所采用的技术方案是:
一种工业过程多分辨率软测量仪表,包括与工业过程对象连接的现场智能仪表、用于存放历史数据的数据存储装置及上位机,智能仪表、数据存储装置及上位机依次相连,所述的上位机为软测量智能处理器,所述的软测量智能处理器包括:
标准化处理模块,用于对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i - - - ( 1 )
2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX I - TX ‾ ) - - - ( 2 )
3)标准化: X = TX - TX ‾ σ X - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A20061015555600094
为训练样本的均值;
偏最小二乘模块,用于采用NIPALS算法提取偏最小二乘主成分,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j = - - - ( 5 )
u j ^ = t j b j , j = 1,2 . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u ^ j , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u ^ j q j T , j = 1,2 . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式(11)、(12)计算:
w j T = u ^ j T X u ^ j T u j ^ , = 1,2 , . . . , a - - - ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
多分辨率分解模块,对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间VO可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
Wj={ψj,k(t)ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解,其算式为(15):
Figure A20061015555600105
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
Figure A20061015555600106
D j f ( t ) = Σ k ∈ Z d j , k ψ j , k ( t ) - - - ( 17 )
径向基函数神经网络建模模块,用于建立软测量模型,采用如下过程:
1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α;
2)用最小二乘学习算法确定RBF中心矢量Ci
3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18):
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn(1≤i≤N)为RBF中心;(ω)i(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
多分辨率重构模块,用于将各分辨率下的软测量值相加得到总的软测量值;信号采集模块,用于依照设定的每次采样的时间间隔,从数据库中采集数据;软测量模块,用于对待检测数据VX用训练时得到的
Figure A20061015555600113
和δx 2进行标准化处理,并将标准化处理后的数据依次进过训练得到的偏最小二乘模块和多分辨率分解模块处理后作为径向基函数神经网络建模模块的输入,将输入代入训练得到的径向基函数神经网络模型,再经过多分辨率重构模块得到软测量函数值。
作为优选的一种方案:所述的软测量智能处理器还包括:模型更新模块,用于定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
作为优选的再一种方案:所述的软测量仪表还包括DCS系统,所述的DCS系统由数据接口、控制站和历史数掘库构成,所述的数据存储装置为DCS系统的历史数据库,所述的软测量智能处理器还包括:结果显示模块,用于将软测量结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。
作为优选的另一种方案:所述的现场智能仪表、DCS系统、软测量智能处理器通过现场总线依次连接。
一种用所述的工业过程多分辨率软测量仪表实现的软测量方法,包括以下步骤:
(1)、确定软测量所用的关键变量,从历史数据库中采集系统正常时所述变量的数据作为训练样本TX;
(2)、设置偏最小二乘模块的提取隐变量个数、多分辨率分解模块的分解层数、径向基函数神经网络建模模块的核参数参数,并设定采样周期;
(3)、训练样本TX在软测量智能处理器中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
3.1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i , - - - ( 1 )
3.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) , - - - ( 2 )
3.3)标准化: X = TX - TX ‾ σ x , - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A20061015555600124
为训练样本的均值;
(4)、采用NIPALS算法提取偏最小二乘主成分,在不损失过多信息的同时,降低系统维数,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u j ^ = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u ^ j , j = 1,2 , . . . , a ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u ^ j q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式计算,其算式为(11)、(12):
w j T = u ^ j T X u ^ j T u ^ j , j = 1,2 , . . . , a - - - ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
(5)、对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度。空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解为(15):
Figure A20061015555600141
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
D j f ( t ) = Σ K ∈ Z d j , k ψ j , k ( t ) - - - ( 17 ) ;
(6)、在各个分辨率上分别建立RBF模型,采用如下过程:
6、1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α;
6、2)用最小二乘学习算法确定RBF中心矢量Ci
6、3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18):
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
(7)、将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待测量数据VX;对VX用训练时得到的
Figure A20061015555600145
和δx 2进行标准化处理,并将标准化处理后的数据依次用训练时得到的偏最小二乘模块和多分辨率分解模块参数进行处理得到不同分辨率上的分解信号,作为对应各分辨率上径向基函数神经网络模型的输入,将输入代入训练得到的RBF模型,得到各分辨率上的软测量函数值;
(8)、将各分辨率下的软测量值相加得到总的软测量值;
作为优选的一种方案:所述的软测量方法还包括:(9)、定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
作为优选的再一种方案:所述的软测量仪表还包括DCS系统,所述的DCS系统由数据接口、控制站和历史数据库构成,所述的数据存储装置为DCS系统的历史数据库,在所述的(8)中计算得到软测量值,将结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。
本发明的技术构思为:本发明分别针对工业过程数据的复共线性、多分辨率特性和非线性特性,将偏最小二乘法、多分辨率分析和径向基函数神经网络三种方法相结合,建立工业过程软测量模型。
针对工业过程数据中的复共线性、多分辨率特性和非线性特性,充分利用了偏最小二乘法,多分辨率理论和神经网络各自的特点,将偏最小二乘法的解相关性能力,多分辨率分析对信息很强的分解与重构能力以及神经网络的多变量非线性映射能力很好地结合了起来,发挥了各自的优势。由此建立的工业过程软测量模型更符合过程的实际特性,可以得到更好的预测和泛化效果。
本发明的有益效果主要表现在:1、同时考虑了复共线性、多分辨率特性和非线性特性;2、建立的模型更贴近实际工业过程,可以得到更好的软测量效果。
(四)附图说明
图1是本发明所提出的软测量系统的硬件结构图;
图2是本发明所提出的软测量智能处理器的功能模块图;
(五)具体实施方式
下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
实施例1
参照图1、图2,一种工业过程多分辨率软测量仪表,包括与工业过程对象1连接的现场智能仪表2、用于存放历史数据的数据存储装置5及上位机6,智能仪表1、数据存储装置5及上位机6依次相连,所述的上位机6包括:
标准化处理模块7,用于对数据库中采集系统正常时的数据进行标准化处理,各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成:
1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i , - - - ( 1 )
2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) , - - - ( 2 )
3)标准化: X = TX - TX ‾ σ x , - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A20061015555600164
为训练样本的均值;
偏最小二乘模块8,用于采用NIPALS算法提取偏最小二乘主成分,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u j ^ = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u j ^ , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u j ^ q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式(11)、
(12)计算:
w j T = u ^ j T X u ^ j T u j ^ , j = 1,2 , . . . , a ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
多分辨率分解模块9,对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解,其算式为(15):
Figure A20061015555600172
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
D j f ( t ) = Σ k ∈ Z d j , k ψ j , k ( t ) - - - ( 17 )
径向基函数神经网络(RBF)建模模块10,用于建立软测量模型,采用如下过程:
1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α;
2)用最小二乘学习算法确定RBF中心矢量Ci
3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18);
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn
(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
多分辨率重构模块11,用于将各分辨率下的软测量值相加得到总的软测量值;
信号采集模块12,用于依照设定的每次采样的时间间隔,从数据库中采集数据;
软测量模块13,用于对待检测数据VX用训练时得到的 和δx 2进行标准化处理,并将标准化处理后的数据依次进过训练得到的偏最小二乘模块和多分辨率分解模块处理后作为径向基函数神经网络建模模块的输入,将输入代入训练得到的径向基函数神经网络模型,再经过多分辨率重构模块得到软测量函数值。
所述的软测量智能处理器6还包括:模型更新模块14,用于定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
所述软测量仪表还包括DCS系统,所述的DCS系统由数据接口3、控制站4、数据库5构成;智能仪表2、DCS系统、上位机6通过现场总线依次相连;软测量智能处理器6还包括结果显示模块15,用于将软测量结果传给DCS系统,并在DCS的控制站显示过程状态,同时通过DCS系统和现场总线将过程状态信息传递到现场操作站进行显示。
所述智能处理器6的硬件部分包括:I/O元件,用于数据的采集和信息的传递;数据存储器,存储运行所需的数据样本和运行参数等;程序存储器,存储实现功能模块的软件程序;运算器,执行程序,实现指定的功能;显示模块,显示设置的参数和运行结果。
当软测量仪表待检测过程已配有DCS系统时,样本实时动态数据的检测、存储利用DCS系统的实时和历史数据库,软测量功能主要在上位机上完成。
当软测量仪表待检测过程没有配备DCS系统时,采用数据存储器来替代DCS系统的实时和历史数据库的数据存储功能,并将软测量仪表制造成包括I/O元件、数据存储器、程序存储器、运算器、显示模块几大构件的不依赖于DCS系统的一个独立的完整的片上系统,在不管检测过程是否配备DCS的情况下,都能够独立使用,更有益于推广使用。
本实施例的工业过程多分辨率软测量仪表,包括与工业过程对象1连接的现场智能仪表2、DCS系统以及软测量智能处理器6,所述的DCS系统由数据接口3、控制站4、数据库5构成;智能仪表2、DCS系统、软测量智能处理器6通过现场总线依次相连,所述的软测量智能处理器6包括:
标准化处理模块7,用于对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i , - - - ( 1 )
2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) , - - - ( 2 )
3)标准化: X = TX - TX ‾ σ x , - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A20061015555600194
为训练样本的均值
偏最小二乘(PLS)模块8,用于采用NIPALS算法提取偏最小二乘主成分,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u ^ j = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u j ^ , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u j ^ q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式(11)、(12)计算:
w j T = u ^ j T X u ^ J T u j ^ , j = 1,2 , . . . , a ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
多分辨率分解(MR)模块9,对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
VJ={J,k(t)|J,k(t)=2-J/(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解,其算式为(15):
Figure A20061015555600206
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
Figure A20061015555600211
D j f ( t ) = Σ k ∈ Z d j , k ψ j , k ( t ) - - - ( 17 )
径向基函数神经网络建模(RBF)模块10,用于用于建立软测量模型,采用如下过程:
1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α;
2)用最小二乘学习算法确定RBF中心矢量Ci
3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18):
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
多分辨率重构模块11,将各分辨率下的软测量值相加得到总的软测量值;
信号采集模块12,用于设定每次采样的时间间隔,从数据库中采集数据;
软测量模块13,用于对待检测数据VX用训练时得到的 和δx 2进行标准化处理,并将标准化处理后的数据依次进过训练得到的PLS和MR模块处理后作为RBF建模模块的输入,将输入代入训练得到的RBF模型,再经过多分辨率重构模块得到软测量函数值。
所述的软测量智能处理器还包括:模型更新模块14,用于定期将离线检测的实际数据加到训练集中,以更新RBF模型。
所述的软测量智能处理器还包括:结果显示模块15,用于将软测量结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。
实施例2
参照图1、图2,一种工业过程多分辨率软测量仪表实现的软测量方法,所述的软测量方法包括以下步骤:
(1)、确定软测量所用的关键变量,从历史数据库中采集系统正常时所述变量的数据作为训练样本TX;
(2)、设置偏最小二乘模块的提取隐变量个数、多分辨率分解模块的分解层数、径向基函数神经网络建模模块的核参数参数,并设定采样周期;
(3)、训练样本TX在软测量智能处理器中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
3.1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i , - - - ( 1 )
3.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) , - - - ( 2 )
3.3)标准化: X = TX - TX ‾ σ x , - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A20061015555600224
为训练样本的均值;
(4)、采用NIPALS算法提取偏最小二乘主成分,在不损失过多信息的同时,降低系统维数,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u j ^ = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u j ^ , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u j ^ q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式计算,其算式为(11)、(12):
w j T = u ^ j T X u ^ j T u j ^ , j = 1,2 , . . . , a - - - ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
(5)、对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度。空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解为(15):
Figure A20061015555600241
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
Figure A20061015555600242
D j f ( t ) = Σ K ∈ Z d j , k ψ j , k ( t ) - - - ( 17 ) ;
(6)、在各个分辨率上分别建立RBF模型,采用如下过程:
6、1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定
形状参数α;
6、2)用最小二乘学习算法确定RBF中心矢量Ci
6、3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,
其算式为(18):
f ( X ) = ω 0 + Σ i = 0 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
(7)、将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待测量数据VX;对VX用训练时得到的
Figure A20061015555600245
和δx 2进行标准化处理,并将标准化处理后的数据依次用训练时得到的偏最小二乘模块和多分辨率分解模块参数进行处理得到不同分辨率上的分解信号,作为对应各分辨率上径向基函数神经网络模型的输入,将输入代入训练得到的RBF模型,得到各分辨率上的软测量函数值;
(8)、将各分辨率下的软测量值相加得到总的软测量值;
所述的软测量方法还包括:(9)、定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
所述的数据存储装置5为DCS系统的历史数据库,所述的DCS系统由数据接口3、控制站4和历史数据库5构成,智能仪表2、DCS系统、软测量智能处理器6通过现场总线依次相连;在所述的(8)中计算得到软测量值,将结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。

Claims (7)

1、一种工业过程多分辨率软测量仪表,包括与工业过程对象连接的现场智能仪表、用于存放历史数据的数据存储装置及上位机,智能仪表、数据存储装置及上位机依次相连,其特征在于:所述的上位机为软测量智能处理器,所述的软测量智能处理器包括:
标准化处理模块,用于对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
1)计算均值: TX ‾ = 1 N Σ i = 1 N T X i - - - ( 1 )
2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) - - - ( 2 )
3)标准化: X = TX - TX ‾ σ x - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A2006101555560002C4
为训练样本的均值;
偏最小二乘模块,用于采用NIPALS算法提取偏最小二乘主成分,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 1 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u j ^ = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
p j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u j T ^ Y u j T ^ u j ^ , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 , . . . , a - - - ( 9 )
F = Y - u j ^ q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式(11)、(12)计算:
w j T = u ^ j T X u ^ j T u ^ j , j = 1,2 , . . . a - - - ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
多分辨率分解模块,对数据进行多分辨率分解;其具体过程如下:
原始信号空间V0可以分解为一系列逼近空VJ与细节空间Wj,其中J是最粗的尺度,空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
     VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
     Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解,其算式为(15):
Figure A2006101555560003C4
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
Figure A2006101555560003C5
D j f ( t ) = Σ k ∈ Z d j , k ψ j , k ( t ) - - - ( 17 )
径向基函数神经网络建模模块,用于建立软测量模型,采用如下过程:
1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α:
2)用最小二乘学习算法确定RBF中心矢量Ci
3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18):
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;Ci∈Rn(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
多分辨率重构模块,用于将各分辨率下的软测量值相加得到总的软测量值;
信号采集模块,用于依照设定的每次采样的时间间隔,从数据库中采集数据;
软测量模块,用于对待检测数据VX用训练时得到的 和δx 2进行标准化处理,并将标准化处理后的数据依次进过训练得到的偏最小二乘模块和多分辨率分解模块处理后作为径向基函数神经网络建模模块的输入,将输入代入训练得到的径向基函数神经网络模型,再经过多分辨率重构模块得到软测量函数值。
2、如权利要求1所述的工业过程多分辨率软测量仪表,其特征在于:所述的软测量智能处理器还包括:模型更新模块,用于定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
3、如权利要求1或2所述的工业过程多分辨率软测量仪表,其特征在于:所述的软测量仪表还包括DCS系统,所述的DCS系统由数据接口、控制站和历史数据库构成,所述的数据存储装置为DCS系统的历史数据库,所述的软测量智能处理器还包括:结果显示模块,用于将软测量结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。
4、如权利要求3所述的工业过程多分辨率软测量仪表,其特征在于:所述的现场智能仪表、DCS系统、软测量智能处理器通过现场总线依次连接。
5、一种用如权利要求1所述的工业过程多分辨率软测量仪表实现的软测量方法,其特征在于,所述的软测量方法包括以下步骤:
(1)、确定软测量所用的关键变量,从历史数据库中采集系统正常时所述变量的数据作为训练样本TX;
(2)、设置偏最小二乘模块的提取隐变量个数、多分辨率分解模块的分解层数、径向基函数神经网络建模模块的核参数参数,并设定采样周期;
(3)、训练样本TX在软测量智能处理器中,对数据进行标准化处理,使得各变量的均值为0,方差为1,得到输入矩阵X,采用以下过程来完成,其算式为(1)、(2)、(3):
3.1)计算均值: TX ‾ = 1 N Σ i = 1 N TX i , - - - ( 1 )
3.2)计算方差: σ x 2 = 1 N - 1 Σ i = 1 N ( TX i - TX ‾ ) , - - - ( 2 )
3.3)标准化: X = TX - TX ‾ σ x , - - - ( 3 )
其中,TX为训练样本,N为训练样本数,
Figure A2006101555560005C2
为训练样本的均值;(4)、采用NIPALS算法提取偏最小二乘主成分,在不损失过多信息的同时,降低系统维数,其具体过程如下:其算式为:(4)、(5)、(6),
t j = Σ p = 1 P x p w pj = Xw j - - - ( 4 )
u j = Σ q = 1 Q y q c kj = Yc j - - - ( 5 )
u j ^ = t j b j , j = 1,2 , . . . a - - - ( 6 )
载荷向量根据标准最小二乘回归计算,其算式为(7)、(8):
P j T = t j T X t j T t j , j = 1,2 , . . . , a - - - ( 7 )
q j T = u ^ j T Y u ^ j T u ^ j , j = 1,2 , . . . , a - - - ( 8 )
迭代的最后一步为计算残差矩阵,其算式为(9)、(10):
E = X - t j p j T , j = 1,2 . . . , a - - - ( 9 )
F = Y - u j ^ q j T , j = 1,2 , . . . , a - - - ( 10 )
列向量w和c分别为自变量和因变量提取主成分的权重向量,通过下式计算,其算式为(11)、(12):
w j T = u ^ j T X u ^ j T u ^ j , j = 1,2 , . . . , a - - - ( 11 )
c j T = t j T Y t j T t j , j = 1,2 , . . . , a - - - ( 12 )
若主成分得分、权向量、载荷向量以及内模型回归系数均以矩阵形式表示为T、W、C、P、Q、B,同时定义:
                R=W(PTW)-1
T=XR
Y ^ = U ^ Q T = TBQ T = XRBQ T ;
(5)、对数据进行多分辨率分解,充分利用数据的多分辨率特性,以获取更多信息。其具体过程如下:
原始信号空间V0可以分解为一系列逼近空间VJ与细节空间Wj,其中J是最粗的尺度,也称为分解的尺度。空间VJ由尺度函数{J,k(t),k∈Z}张成,其算式为(13):
      VJ={J,k(t)|J,k(t)=2-J/2(2-Jt-k)}    (13)
而Wj由小波函数{ψj,k(t),j=1,...,J,k∈Z}张成,其算式为(14):
      Wj={ψj,k(t)|ψj,k(t)=2-j/2ψ(2-Jt-k)}    (14)
其中j是尺度因子,k是平移因子;
根据上述描述,函数f(t)∈L2(R)可以分解为(15):
其中,第一项代表近似信息,第二项为细节信息,逼近因子aJ,k与细节因子dj,k通过Mallat算法计算;
近似信号AJf(t)与细节信号Djf(t)(j=1,2,...,J)的定义如下式(16)、(17):
Figure A2006101555560006C3
D j f ( t ) = Σ k ∈ Z d j , k ψ j , k ( t ) - - - ( 17 ) ;
(6)、在各个分辨率上分别建立RBF模型,采用如下过程:
6、1)选用高斯函数Φ(ν)=exp(-ν22)作为网络的激活函数,给定形状参数α;
6、2)用最小二乘学习算法确定RBF中心矢量Ci
6、3)正交优选确定最佳隐含层数目、网络输出权值得到软测量模型,其算式为(18):
f ( X ) = ω 0 + Σ i = 1 N ω i Φ ( | | X - C i | | ) - - - ( 18 )
其中,X∈Rn是输入向量;Φ(·)为从R+→R的一个非线性函数;C1∈Rn(1≤i≤N)为RBF中心;ωi(1≤i≤N)为连接权值,ω0为偏置量;N为隐含层的神经元数;‖·‖是欧氏范数;
(7)、将采集的数据传送到DCS实时数据库中,在每个定时周期从DCS数据库的实时数据库中,得到最新的变量数据作为待测量数据VX;对VX用训练时得到的
Figure A2006101555560007C1
和δx 2进行标准化处理,并将标准化处理后的数据依次用训练时得到的偏最小二乘模块和多分辨率分解模块参数进行处理得到不同分辨率上的分解信号,作为对应各分辨率上径向基函数神经网络模型的输入,将输入代入训练得到的RBF模型,得到各分辨率上的软测量函数值;
(8)、将各分辨率下的软测量值相加得到总的软测量值;
6、如权利要求5所述的工业过程多分辨率软测量方法,其特征在于:所述的软测量方法还包括:(9)、定期将离线检测的实际数据加到训练集中,以更新径向基函数神经网络模型。
7、如权利要求5或6所述的工业过程多分辨率软测量方法,其特征在于:所述的数据存储装置为DCS系统的历史数据库,所述的DCS系统由数据接口、控制站和历史数据库构成,在所述的(8)中计算得到软测量值,将结果传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示。
CNA2006101555563A 2006-12-28 2006-12-28 工业过程多分辨率软测量仪表及方法 Pending CN101017373A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2006101555563A CN101017373A (zh) 2006-12-28 2006-12-28 工业过程多分辨率软测量仪表及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2006101555563A CN101017373A (zh) 2006-12-28 2006-12-28 工业过程多分辨率软测量仪表及方法

Publications (1)

Publication Number Publication Date
CN101017373A true CN101017373A (zh) 2007-08-15

Family

ID=38726426

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2006101555563A Pending CN101017373A (zh) 2006-12-28 2006-12-28 工业过程多分辨率软测量仪表及方法

Country Status (1)

Country Link
CN (1) CN101017373A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102601951A (zh) * 2012-03-12 2012-07-25 浙江大学 基于超声信号的塑料注射过程模腔压力检测方法
CN101520453B (zh) * 2009-02-10 2012-12-26 浙江大学 一种催化裂化粗汽油干点软测量仪表
CN103605287A (zh) * 2013-08-05 2014-02-26 浙江大学 循环流化床锅炉床温预测系统及方法
CN103675006A (zh) * 2013-09-22 2014-03-26 浙江大学 最小二乘的工业熔融指数软测量仪表和方法
CN103674778A (zh) * 2013-09-22 2014-03-26 浙江大学 Rbf粒子群优化的工业熔融指数软测量仪表和方法
CN105160396A (zh) * 2015-07-06 2015-12-16 东南大学 一种利用现场数据建立神经网络模型的方法
CN105793789A (zh) * 2013-12-05 2016-07-20 拜耳技术服务有限责任公司 用于过程单元中的全部过程区段的自动的监视和状态确定的计算机实现的方法和系统
CN107057749A (zh) * 2016-08-04 2017-08-18 石化盈科信息技术有限责任公司 汽油稳定塔分离精度的在线软测量系统及测量方法
CN107064054A (zh) * 2017-02-28 2017-08-18 浙江大学 一种基于cc‑pls‑rbfnn优化模型的近红外光谱分析方法
CN108492026A (zh) * 2018-03-06 2018-09-04 宁波大学 一种基于集成正交成分最优化回归分析的软测量方法
CN110197031A (zh) * 2019-05-29 2019-09-03 北京化工大学 一种用于高密度聚乙烯串级聚合反应的软测量方法
CN118430092A (zh) * 2024-07-01 2024-08-02 深圳灿态信息技术有限公司 一种基于mcc系统的数据通用采集方法

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101520453B (zh) * 2009-02-10 2012-12-26 浙江大学 一种催化裂化粗汽油干点软测量仪表
CN102601951B (zh) * 2012-03-12 2014-05-07 浙江大学 基于超声信号的塑料注射过程模腔压力检测方法
CN102601951A (zh) * 2012-03-12 2012-07-25 浙江大学 基于超声信号的塑料注射过程模腔压力检测方法
CN103605287B (zh) * 2013-08-05 2016-03-30 浙江大学 循环流化床锅炉床温预测系统及方法
CN103605287A (zh) * 2013-08-05 2014-02-26 浙江大学 循环流化床锅炉床温预测系统及方法
CN103675006A (zh) * 2013-09-22 2014-03-26 浙江大学 最小二乘的工业熔融指数软测量仪表和方法
CN103674778A (zh) * 2013-09-22 2014-03-26 浙江大学 Rbf粒子群优化的工业熔融指数软测量仪表和方法
CN103674778B (zh) * 2013-09-22 2015-08-26 浙江大学 Rbf粒子群优化的工业熔融指数软测量仪表和方法
CN103675006B (zh) * 2013-09-22 2015-09-09 浙江大学 最小二乘的工业熔融指数软测量仪表和方法
CN105793789A (zh) * 2013-12-05 2016-07-20 拜耳技术服务有限责任公司 用于过程单元中的全部过程区段的自动的监视和状态确定的计算机实现的方法和系统
CN105160396A (zh) * 2015-07-06 2015-12-16 东南大学 一种利用现场数据建立神经网络模型的方法
CN105160396B (zh) * 2015-07-06 2018-04-24 东南大学 一种利用现场数据建立神经网络模型的方法
CN107057749A (zh) * 2016-08-04 2017-08-18 石化盈科信息技术有限责任公司 汽油稳定塔分离精度的在线软测量系统及测量方法
CN107057749B (zh) * 2016-08-04 2018-12-07 石化盈科信息技术有限责任公司 汽油稳定塔分离精度的在线软测量系统及测量方法
CN107064054A (zh) * 2017-02-28 2017-08-18 浙江大学 一种基于cc‑pls‑rbfnn优化模型的近红外光谱分析方法
CN107064054B (zh) * 2017-02-28 2019-08-02 浙江大学 一种基于cc-pls-rbfnn优化模型的近红外光谱分析方法
CN108492026A (zh) * 2018-03-06 2018-09-04 宁波大学 一种基于集成正交成分最优化回归分析的软测量方法
CN108492026B (zh) * 2018-03-06 2021-05-11 宁波大学 一种基于集成正交成分最优化回归分析的软测量方法
CN110197031A (zh) * 2019-05-29 2019-09-03 北京化工大学 一种用于高密度聚乙烯串级聚合反应的软测量方法
CN110197031B (zh) * 2019-05-29 2021-03-26 北京化工大学 一种用于高密度聚乙烯串级聚合反应的软测量方法
CN118430092A (zh) * 2024-07-01 2024-08-02 深圳灿态信息技术有限公司 一种基于mcc系统的数据通用采集方法
CN118430092B (zh) * 2024-07-01 2024-09-20 深圳灿态信息技术有限公司 一种基于mcc系统的数据通用采集方法

Similar Documents

Publication Publication Date Title
CN101017373A (zh) 工业过程多分辨率软测量仪表及方法
CN106197999B (zh) 一种行星齿轮故障诊断方法
CN100480926C (zh) 一种基于小波分析的工业生产过程故障诊断系统及方法
CN105092892B (zh) 一种车辆加速度数据的获取方法和装置
CN108399248A (zh) 一种时序数据预测方法、装置及设备
CN105629958B (zh) 一种基于子时段mpca‑svm的间歇过程故障诊断方法
JP2002525757A (ja) 統計的プロセス管理のための多次元的方法及びシステム
CN105115594A (zh) 基于小波熵和信息融合的齿轮箱振动信号故障特征提取方法
CN113919540B (zh) 一种生产过程运行状态的监测方法及相关设备
CN107180367A (zh) 一种基于机器学习的车辆估价的方法、存储介质及装置
CN108549908A (zh) 基于多采样概率核主成分模型的化工过程故障检测方法
CN100470418C (zh) 基于盲源信号分析的聚丙烯熔融指数软测量仪表及方法
CN116108371B (zh) 基于级联异常生成网络的云服务异常诊断方法与系统
CN111461187A (zh) 一种建筑物沉降智能检测系统
CN114065919A (zh) 一种基于生成对抗网络的缺失值补全方法及介质
CN111680725A (zh) 基于重构贡献的气体传感器阵列多故障隔离算法
CN112200048A (zh) 一种基于回归模型的旋转设备故障预测方法、系统及可读存储介质
CN100470427C (zh) 基于仿生智能的工业软测量仪表及软测量方法
CN112906288A (zh) 基于正常样本数据的开关柜状态评估方法
CN110209145B (zh) 一种基于核矩阵近似的二氧化碳吸收塔故障诊断方法
Xue et al. A novel framework for motor bearing fault diagnosis based on multi-transformation domain and multi-source data
CN114755010A (zh) 一种旋转机械振动故障诊断方法及其系统
CN114580525A (zh) 一种面向含缺失数据的手势动作分类方法
CN117454113A (zh) 一种基于CNN-BiLSTM及RF模型的电力变压器故障预警方法
CN113688773B (zh) 一种基于深度学习的储罐穹顶位移数据修复方法及其装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication