专利名称：：基于泡沫尺寸统计分布的浮选过程故障诊断方法基于泡沫尺寸统计分布的浮选过程故障诊断方法本发明涉及选矿领域，图像处理技术，概率统计等领域，具体为浮选泡沫图像尺寸统计分布的表征方法，密度估计方法，带时滞的非线性随机动态系统的模型分析，药剂添加量故障的检测与诊断。浮选过程是矿物加工中广泛应用的选矿方法，涉及到极其复杂的物理化学过程，为一类典型复杂工业过程，表现出多变量、非线性、大时滞等特点。随着计算机应用技术、图像处理技术及智能控制等技术的飞速发展，研发了浮选泡沫图像视觉监控系统，对浮选过程进行在线实时监控，解决现场因人工操作主观判断造成的浮选生产过程不稳定、指标波动问题，为厂矿人员分析处理生产指标监控问题提供强有力工具。浮选泡沫特征作为判断浮选效果好坏的一个重要的依据，它包含有大量与操作变量和产品质量有关的信息。通常在实际浮选过程中，有经验的操作工人主要通过观察浮选槽表面泡沫的视觉信息来完成浮选过程的操作，但由于对泡沫结构好坏的判断并没有一个统一的标准。近年来，基于机器视觉和图像处理技术的支撑，常规上对获取的泡沫图像进行分割后，采用简单的均值和方差特征来表征泡沫结构，从而较为粗糙的建立泡沫结构与浮选操作状态的关系。传统方法假设泡沫尺寸分布是基于正态高斯分布，而事实上，浮选泡沫尺寸分布有着很大的特殊性，典型的浮选泡沫尺寸分布呈现非高斯分布，其概率密度分布曲线存在较高的峰值和较大的偏斜度，具有左偏斜(粗选槽泡沫)或者右偏斜(扫选槽泡沫)、长拖尾等特点。所以简单的研究高斯分布的平均值和方差会造成对泡沫现状特征描述的较大误差，会严重影响浮选过程的故障检测和诊断的准确性。本发明的目的在于解决浮选泡沫尺寸统计分布的准确描述和浮选过程药剂添加故障诊断的问题，提供了一种全面表征泡沫结构特征的方法，设计的非参数估计算子逼近尺寸分布的概率密度曲线，将泡沬尺寸PDF转化为核函数的动态权系数，建立的带时滞的非线性权动态模型，基于该模型，运用线性矩阵不等式(LMI)得到合理的优化故障检测方法和诊断方法。6综合利用传感技术、图像处理技术和数据分析技术，设计自适应滤波器，运用随机分布故障检测理论，提供保证检测性能的方法。本发明采用摄像机、光源、图像采集卡、计算机及其附属部件构成系统硬件平台，由此获取浮选槽泡沫图像，并应用形态学和分水岭方法分割泡沫图像，系统泡沫图像分割软件采用C+十编程语言开发，采用simulink建立带时滞的非线性不确定权动态模型。本发明主要内容如下首先通过浮选泡沫图像采集平台，用工业摄像机获得泡沬图像，并对图像进行形态学操作和分水岭分割，对分割后图像进行泡沫结构特征的概率密度统计。利用设计的非参数估计算子逼近出输出概率密度(PDF)曲线，将泡沬尺寸PDF转化为动态权系数，构造出基于PDF的非线性动态随机系统模型。当发生故障时，利用故障前后PDF曲线的动态变化来检测和诊断故障。基于泡沫尺寸的浮选过程故障诊断策略的应用，能实时获得的泡沫形状密度，诊断出浮选泡沫的健康状况，对供矿断料、调整剂Na2C03的加入量失常等故障能给出准确的预报警，能及时给出控制决策，指导加药量调节、优化生产。图1泡沫图像分割和气泡尺寸统计分布的软件界面图2理想泡沬图像和故障泡沬图像图3故障发生前后输出PDF的3Dmesh图图4带时滞非线性随机系统故障分析模型图5逼近输出PDF的非参数估计方法比较图6故障发生时残差信号范数的响应下面结合附图对本发明作进一步的详细说明。为准确全面的描述泡沬形状分布，对由摄像机获取的泡沫图像进行形态学操作和分水岭分割，设计的非参数估计算子逼近尺寸分布的概率密度曲线，非参数估计算子包括核函数算子，直方图估计，小波估计算子，B样条估计算子。将泡沫尺寸PDF转化为动态权系数，建立的带时滞的非线性权动态模型，基于该模型，运用线性矩阵不等式得到合理的优化故障检测方法和诊断方法。通过搭建的浮选泡沫图像采集平台获取泡沬图像，运用开发的系统泡沬图像分割软件，对获取泡沬图像进行形态学操作和分水岭分割，如图1的软件界面为采集的粗选槽泡沫图像分水岭分割后的结果，红色柱状图为泡沫尺寸统计的直方图显示。浮选流程中不同的浮选槽(如粗选槽、粗扫槽和精选槽)，其泡沫尺寸统计分布各有特点。目前存在的传统方法假设泡沫尺寸分布是基于正态高斯分布，只考虑均值和方差等特征。而事实上，典型的粗选槽泡沫尺寸分布呈现非高斯分布，其概率密度分布曲线存在较高的峰值和较大的偏斜度，通常为左偏斜、长拖尾等特点。图2所示为粗选槽采集的理想泡沫图像和药剂故障发生时的泡沫图像。这里的药剂故障指的是调整剂Na2C03的加入量过高，导致pH值超过正常范围(9.49.8)。从分水岭分割后的图像结果来看，当pH值超过正常范围时，泡沫尺寸有增大的趋势，如图3。为准确地判断该故障，提出泡沫尺寸统计分布的浮选过程故障诊断方法具体如下-1.采用核密度估计逼近泡沫尺寸概率密度函数(PDF)设动态随机系统的输入为M(O，输出为XOe[",6],则输出;^)在[。,0范围的概率为S^《)=f/kerO，")血(1)式中/to(^")为输出概率密度函数，其对应的物理含义为浮选泡沫图像分割后泡沫尺寸的概率密度函数(PDF)分布，"(0是控制输入，即浮选系统中的调整剂Na2C03的加入量。根据函数逼近原则，可用下述设计的核密度估计算子来逼近/ke,0c,"):/kerW=^>,(^)(2)其中，/kCTW:逼近的概率密度函数PDF，:第i个核函数的权系数，《("'):第i个核函数，:第i个核函数的x轴中点，&:核函数的窗宽。根据泡沫的尺寸分布，在逼近曲线时选用了30个核函数，核函数以Epanechnikov函数为原型，构建符合浮选过程系统的核函数如式(3):《(^")^d—(^)2)化K[《—(3)4A为了保证输出PDF的积分和为1，核函数应满足>(^^)=1，考虑实际泡沬尺寸分布，需选取合适的h值，其中《=100*/，i=l,2，...,30，；^=200。各个核函数都固定不变，可获得各个核函数相对应的权系数w,，用来表征输出PDF。图5对比了不同非参数估计算子对概率密度的逼近，包括设计的核方法估计，直方图估计法，B样条估计算子。2.构造输出PDF模型为避免由反馈控制所产生的权值出现负值，引入带逼近误差的输出PDF平方根模型V/kerO，w，F)=t(w，F)《0)+w，F)(4),'=i式中，尺,(2)(/=1,2，.』)是定义在[，6]上的选定的基函数，w(z，",F)是逼近PDF曲线带来的误差。w,(w)(/"，2,."")是与"(/)有关的权函数。记『(z)=[Wl(",F),W2(W，F)，…d(",(5)A-〖《(z)《。(z)电A2=〖尺o(z)&(z)&,A3=#0.从￡九々,")^=1可知，权函数向量只有"-l个是相对独立的，则式(4)可改写为=《r(z),+(z)+哞，",)(6)式中K(z)-K。(z)-^^(z),))=:^^1,A。-A,A广A入。/z(『(0)为A3八3第w个核函数、(z)对应的权系数。对浮选泡沫尺寸的概率密度函数的描述，就转化成一组动态权系数模型。3.基于输出PDFs的故障检测滤波因量测信息为输出概率密度分布，为了检测故障，须设计如下故障检测滤波器二爿雄)+^雄一o++-+丄s(0;,=蹄)；(7)雄)是估计状态向量，力，4/，h，A，￡为参数矩阵，"为时滞，丄eir"是待定的滤波器增益，残差信号柳由量测PDFs和估计PDFs之差的积分来确定，其中c7(2)e"是一个定义在["，6]上的给定的权向量。事实上，残差信号W)可以认为是定义在故障前后两幅图像泡沬尺寸PDF之间的距离。记雄)-;c(0-j^)为误差状态向量,力(,)=柳-^)，将式(6)代入式(7)，则可得到误差系统为欢0=(^—)对0+-力—ir2[/2(Ex:0))--丄A(O(8)其中，r,=j[ct(z)A^(z)^^,r2=￡o"(z)A:(z>fe,A(0=fcr(z)6)(z,M，i^fec残差的表达式为f(0=r,雄)+r2(/z(蹄))-力(蹄)))+△W(9)因逼近误差有界，设h(z，",尸)^5，那么可以得到r~cr(z)w(z，m,F)S5，^=^丄cr(z)dz(10)对误差系统进行Lyapimov稳定性分析，并建立线性矩阵不等式求解，可以得到系统的稳定性条件如式(11):雄)||2a=maxjsup||x(0||，77—《-"/so从而可得到残差的范数可以用下式(12)表示(ii)(12)在浮选过程中，根据式(11)求得判断故障的临界误差状态向量的范数值"，将a代入式(12)求得残差范数的阈值P来检测故障。如果残差范数大于这个阈值;9就可以判断系统发生故障，如果小于^则可以判定系统只是发生较小的波动，借此系统就可以采取不同的调节方案来处理。如图6所示虚线为求解得到的阈值，实线为残差范数的响应。沿时间轴方向，残差范数没有超过设定的阈值时，检测滤波器提示系统运行正常，约22s后，当残差范数超过设定的阈值时，则提示系统发生故障。检测滤波器能及时指导实际操作，采取有效的药剂调节。4.基于输出PDFs的故障诊断滤波为了量测故障，设计了如下故障诊断滤波器<formula>formulaseeoriginaldocumentpage11</formula>与检测滤波的设计不同，诊断滤波是在根据式(12)判断出系统发生故障后设计，增加了故障项，其中卢(,)就是对故障F的估计。n和^是诊断滤波的经验参数，可根据稳定性分析得到。设#(0=,(/)-》(0，则估计的故障误差系统为<formula>formulaseeoriginaldocumentpage11</formula>考虑实际中故障的特点，假设故障是有界的。对故障误差系统进行Lyapunov稳定性分析，建立线性矩阵不等式求解，就可以得到滤波器的增益和经验参数，使滤波器可以很好的跟踪发生故障的大小，从而可以对故障进行精确地诊断和调节以尽快消除故障。权利要求1.一种基于泡沫尺寸统计分布的浮选过程故障诊断方法，其特征在于首先通过浮选泡沫图像采集平台，用工业摄像机获得泡沫图像，并对图像进行形态学操作和分水岭分割，对分割后图像进行泡沫结构特征的概率密度统计，利用设计的非参数估计算子逼近出输出概率密度PDF曲线，将泡沫尺寸PDF转化为动态权系数，构造出基于PDF的非线性动态随机系统模型，当发生故障时，利用故障前后PDF曲线的动态变化来检测和诊断故障，具体包括以下步骤①采用核密度估计逼近泡沫尺寸概率密度函数设动态随机系统的输入为u(t)，输出为y(t)∈[a，b]，则输出y(t)在[a，ξ)范围的概率为<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mtext>a≤y</mtext><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>≤</mo><mi>ξ</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>ξ</mi></msubsup><msub><mi>f</mi><mi>ker</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mi>dx</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中fker(x，u)为输出概率密度函数，其对应的物理含义为浮选泡沫图像分割后泡沫尺寸的概率密度函数分布，u(t)是控制输入，即浮选系统中的调整剂Na2CO3的加入量，根据函数逼近原则，用下述设计的核密度估计算子来逼近fker(x，u)<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>^</mo></mover><mi>ker</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><mi>K</mi><mrow><mo>(</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math></maths>其中，id="icf0003"file="A2009100435120002C3.tif"wi="13"he="6"top="183"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>逼近的概率密度函数PDF，wi第i个核函数的权系数，id="icf0004"file="A2009100435120002C4.tif"wi="20"he="8"top="203"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>第i个核函数，Xi第i个核函数的x轴中点，h核函数的窗宽；根据泡沫的尺寸分布，在逼近曲线时选用了30个核函数，核函数以Epanechnikov函数为原型，构建符合浮选过程系统的核函数如式(3)<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mi>K</mi><mrow><mo>(</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mo>=</mo><mfrac><mn>3</mn><mn>4</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>/</mo><mi>h</mi><mo>,</mo></mrow>]]></math>id="icf0005"file="A2009100435120003C1.tif"wi="63"he="11"top="27"left="52"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>x∈[Xi-h，Xi+h](3)为了保证输出PDF的积分和为1，核函数应满足<mathsid="math0004"num="0004"><math><![CDATA[<mrow><munderover><mo>∫</mo><mi>a</mi><mi>b</mi></munderover><mi>K</mi><mrow><mo>(</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub></mrow><mi>h</mi></mfrac><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>,</mo></mrow>]]></math>id="icf0006"file="A2009100435120003C2.tif"wi="29"he="12"top="43"left="137"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>考虑实际泡沫尺寸分布，需选取合适的h值，其中Xi＝100*i，i＝1，2，...，30，h＝200，各个核函数都固定不变，可获得各个核函数相对应的权系数wi，用来表征输出PDF；②构造输出PDF模型为避免由反馈控制所产生的权值出现负值，引入带逼近误差的输出PDF平方根模型<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msqrt><msub><mover><mi>f</mi><mo>^</mo></mover><mi>ker</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow></msqrt><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow><msub><mi>K</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>+</mo><mi>ω</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中，Ki(z)(i＝1，2，...n)是定义在[a，b]上的选定的基函数，ω(z，u，F)是逼近PDF曲线带来的误差，wi(u)(i＝1，2，...n)是与u(t)有关的权函数，记K0(z)＝[k1(z)，k2(z)，...kn-1(z)]TW(z)＝[w1(u，F)，w2(u，F)，...wn-1(u，F)]T(5)<mathsid="math0006"num="0006"><math><![CDATA[<mrow><msub><mi>Λ</mi><mn>1</mn></msub><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><msubsup><mi>K</mi><mn>0</mn><mi>T</mi></msubsup><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>K</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>dz</mi><mo>,</mo><msub><mi>Λ</mi><mn>2</mn></msub><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><msub><mi>K</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>k</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>dz</mi><mo>,</mo><msub><mi>Λ</mi><mn>3</mn></msub><mo>=</mo><msubsup><mrow><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>k</mi></mrow><mi>n</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>dz</mi><mo>≠</mo><mn>0</mn><mo>.</mo></mrow>]]></math></maths>从<mathsid="math0007"num="0007"><math><![CDATA[<mrow><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><msub><mover><mi>f</mi><mo>^</mo></mover><mi>ker</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mi>dz</mi><mo>=</mo><mn>1</mn></mrow>]]></math>id="icf0009"file="A2009100435120003C5.tif"wi="30"he="8"top="190"left="39"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>可知，权函数向量只有n-1个是相对独立的，则式(4)可改写为<mathsid="math0008"num="0008"><math><![CDATA[<mrow><msqrt><msub><mover><mi>f</mi><mo>^</mo></mover><mi>ker</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow></msqrt><mo>=</mo><msup><mi>K</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mrow><mo>(</mo><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>k</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>+</mo><mi>ω</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中<mathsid="math0009"num="0009"><math><![CDATA[<mrow><mi>K</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>K</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mfrac><msub><mi>Λ</mi><mn>2</mn></msub><msub><mi>Λ</mi><mn>3</mn></msub></mfrac><msub><mi>k</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math>id="icf0011"file="A2009100435120003C7.tif"wi="43"he="10"top="232"left="33"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0010"num="0010"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><msqrt><msub><mi>Λ</mi><mn>3</mn></msub><mo>-</mo><msup><mi>W</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>Λ</mi><mn>0</mn></msub><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msqrt><msub><mi>Λ</mi><mn>3</mn></msub></mfrac><mo>,</mo></mrow>]]></math>id="icf0012"file="A2009100435120003C8.tif"wi="52"he="12"top="230"left="79"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0011"num="0011"><math><![CDATA[<mrow><msub><mi>Λ</mi><mn>0</mn></msub><mo>=</mo><msub><mi>Λ</mi><mn>1</mn></msub><msub><mi>Λ</mi><mn>3</mn></msub><mo>-</mo><msubsup><mi>Λ</mi><mn>2</mn><mi>T</mi></msubsup><msub><mi>Λ</mi><mn>2</mn></msub><mo>,</mo></mrow>]]></math>id="icf0013"file="A2009100435120003C9.tif"wi="33"he="4"top="234"left="134"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>h(W(t))为第n个核函数kn(z)对应的权系数，对浮选泡沫尺寸的概率密度函数的描述，就转化成一组动态权系数模型；③基于输出PDFs的故障检测滤波因量测信息为输出概率密度分布，为了检测故障，设计如下故障检测滤波器<mathsid="math0012"num="0012"><math><![CDATA[<mrow><mfencedopen='{'close=''><mtable><mtr><mtd><mover><mover><mi>x</mi><mo>^</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mi>d</mi></msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>Hu</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>H</mi><mi>d</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>Lϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>;</mo></mtd></mtr><mtr><mtd><mi>ϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mrow><mo>(</mo><msqrt><msub><mi>f</mi><mi>ker</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>F</mi><mo>)</mo></mrow></msqrt><mo>-</mo><msqrt><msub><mover><mi>f</mi><mo>^</mo></mover><mi>ker</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></msqrt><mo>)</mo></mrow><mi>dz</mi><mo>;</mo></mtd></mtr><mtr><mtd><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>E</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>;</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math></maths>id="icf0015"file="A2009100435120004C2.tif"wi="6"he="4"top="79"left="30"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>是估计状态向量，A，Ad，H，Hd，E为参数矩阵，d为时滞，L∈Rm×p是待定的滤波器增益，残差信号ε(t)由量测PDFs和估计PDFs之差的积分来确定，其中σ(z)∈Rn×1是一个定义在[a，b]上的给定的权向量；记<mathsid="math0013"num="0013"><math><![CDATA[<mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math>id="icf0016"file="A2009100435120004C3.tif"wi="29"he="4"top="103"left="37"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>为误差状态向量，<mathsid="math0014"num="0014"><math><![CDATA[<mrow><mover><mover><mi>x</mi><mo>~</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>x</mi><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mover><mover><mi>x</mi><mo>^</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math>id="icf0017"file="A2009100435120004C4.tif"wi="29"he="5"top="103"left="107"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>将式(6)代入式(7)，则可得到误差系统为<mathsid="math0015"num="0015"><math><![CDATA[<mrow><mover><mover><mi>x</mi><mo>~</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><mi>L</mi><msub><mi>Γ</mi><mn>1</mn></msub><mo>)</mo></mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mi>d</mi></msub><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>-</mo><mi>L</mi><msub><mi>Γ</mi><mn>2</mn></msub><mo>[</mo><mi>h</mi><mrow><mo>(</mo><mi>Ex</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>E</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>LΔ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math></maths>其中，<mathsid="math0016"num="0016"><math><![CDATA[<mrow><msub><mi>Γ</mi><mn>1</mn></msub><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><msup><mi>K</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>Edz</mi><mo>,</mo><msub><mi>Γ</mi><mn>2</mn></msub><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>k</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>dz</mi><mo>,</mo><mi>Δ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>ω</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow><mi>dz</mi><mo>.</mo></mrow>]]></math>id="icf0019"file="A2009100435120004C6.tif"wi="142"he="9"top="134"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>残差的表达式为<mathsid="math0017"num="0017"><math><![CDATA[<mrow><mi>ϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Γ</mi><mn>1</mn></msub><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Γ</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>h</mi><mrow><mo>(</mo><mi>Ex</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>E</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>Δ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math></maths>因逼近误差有界，设|ω(z，u，F)|≤δ，那么可以得到<mathsid="math0018"num="0018"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><mi>Δ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>=</mo><mo>|</mo><mo>|</mo><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>ω</mi><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>u</mi><mo>,</mo><mi>F</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>≤</mo><mover><mi>δ</mi><mo>~</mo></mover><mo>,</mo></mrow>]]></math>id="icf0021"file="A2009100435120004C8.tif"wi="66"he="11"top="180"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0019"num="0019"><math><![CDATA[<mrow><mover><mi>δ</mi><mo>~</mo></mover><mo>=</mo><mi>δ</mi><msubsup><mo>∫</mo><mi>a</mi><mi>b</mi></msubsup><mi>σ</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>dz</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math>id="icf0022"file="A2009100435120004C9.tif"wi="67"he="9"top="181"left="116"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>对误差系统进行Lyapunov稳定性分析，并建立线性矩阵不等式求解，得到系统的稳定性条件如式(11)<mathsid="math0020"num="0020"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>≤</mo><mi>α</mi><mo>=</mo><mi>max</mi><mo>{</mo><munder><mi>sup</mi><mrow><mo>-</mo><mi>d</mi><mo>≤</mo><mi>t</mi><mo>≤</mo><mn>0</mn></mrow></munder><mo></mo><mo>|</mo><mo>|</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>,</mo><msup><mi>η</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>δ</mi><mo>~</mo></mover><mo>|</mo><mo>|</mo><mi>R</mi><mo>|</mo><mo>|</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math></maths>从而可得到残差的范数可以用下式(12)表示<mathsid="math0021"num="0021"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><mi>ϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>></mo><mi>β</mi><mo>:</mo><mo>=</mo><mi>α</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>Γ</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><mo>+</mo><mo>|</mo><mo>|</mo><msub><mi>Γ</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo><mo>|</mo><mo>|</mo><msub><mi>U</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo><mo>|</mo><mo>|</mo><mi>E</mi><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>+</mo><mover><mi>δ</mi><mo>~</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math></maths>在浮选过程中，根据式(11)求得判断故障的临界误差状态向量的范数值α，将α代入式(12)求得残差范数的阈值β来检测故障，如果残差范数大于这个阈值β就可以判断系统发生故障，如果小于β则可以判定系统只是发生较小的波动，借此系统就可以采取不同的调节方案来处理，当残差范数超过设定的阈值时，则提示系统发生故障，检测滤波器能及时指导实际操作，采取有效的药剂调节；④基于输出PDFs的故障诊断滤波为了量测故障，设计了如下故障诊断滤波器<mathsid="math0022"num="0022"><math><![CDATA[<mrow><mfencedopen='{'close=''><mtable><mtr><mtd><mover><mover><mi>x</mi><mo>^</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mi>d</mi></msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>Hu</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>H</mi><mi>d</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>Lϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>J</mi><mover><mi>F</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mover><mover><mi>F</mi><mo>^</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>γ</mi><mn>1</mn></msub><mover><mi>F</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>γ</mi><mn>2</mn></msub><mi>ϵ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math></maths>与检测滤波的设计不同，诊断滤波是在根据式(12)判断出系统发生故障后设计，增加了故障项，其中id="icf0026"file="A2009100435120005C2.tif"wi="7"he="5"top="118"left="89"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>就是对故障F的估计，则估计的故障误差系统为<mathsid="math0023"num="0023"><math><![CDATA[<mrow><mover><mover><mi>x</mi><mo>~</mo></mover><mo>·</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><msub><mi>LΓ</mi><mn>1</mn></msub><mo>)</mo></mrow><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mi>d</mi></msub><mover><mi>x</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>J</mi><mover><mi>F</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>L</mi><msub><mi>Γ</mi><mn>2</mn></msub><mo>[</mo><mi>h</mi><mrow><mo>(</mo><mi>Ex</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>E</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mi>LΔ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math></maths>考虑实际中故障的特点，假设故障是有界的，对故障误差系统进行Lyapunov稳定性分析，建立线性矩阵不等式求解，得到滤波器的增益和经验参数，使滤波器可以很好的跟踪发生故障的大小，从而对故障进行精确地诊断和调节以尽快消除故障。全文摘要一种基于泡沫尺寸统计分布的浮选过程故障诊断方法，本发明对获取的泡沫图像进行分水岭分割后，采用概率密度函数(PDF)准确描述泡沫尺寸统计特性，设计适于泡沫尺寸分布的估计算子来逼近概率密度分布，从而将泡沫尺寸PDF转化为非参数估计的动态权系数，进而建立带时滞的非线性权动态模型，基于线性矩阵不等式得到可行的优化故障检测方法和诊断方法。本发明可用于浮选流程的连续生产过程中药剂故障的检测和诊断，通过跟踪输出概率密度，分析估计曲线的权系数，设计滤波器检测出故障，并对故障进行准确的跟踪，有利于操作人员及时发现故障，稳定生产过程，实现浮选生产操作的优化。文档编号G06T7/00GK101566848SQ20091004351公开日2009年10月28日申请日期2009年5月27日优先权日2009年5月27日发明者唐朝晖,朱红求,杜建江,桂卫华,许灿辉,谢永芳,阳春华申请人:中南大学