网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于神经网络的荧光油膜厚度与灰度研究  PDF

  • 钱泓江
  • 董秀成
  • 徐椰烃
  • 蒋金洋
  • 陈桂芳
西华大学电气与电子信息学院,成都 610039

中图分类号: V19

最近更新:2023-02-17

DOI:10.16356/j.1005-2615.2021.03.020

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

在全局摩阻测量中,薄油膜技术可以很好地表征表面摩阻的分布情况。用特定波长的紫外线照射添加了荧光显色分子的不同厚度的油膜,油膜将发出不同的亮度。利用该原理通过检测受激发的荧光油膜灰度值可解算出相应油膜的厚度。本次采用BP神经网络及极限学习机(Extreme learning machine, ELM)神经网络搭建模型完成了荧光油膜厚度与灰度关系的预测,运用Hopfield神经网络完成了相应参数的辨识。实验表明,ELM神经网络模型、BP神经网络模型及插值法模型的预测误差分别为5.150%、5.485%和5.935%。通过Hopfield神经网络辨识,光源功率、光距和曝光系数等影响因素的参数误差率控制在1%左右,达到实际工程运用的要求。与传统插值法相比,通过神经网络可获得更高的精度,为荧光油膜灰度与厚度研究提供了一种可行的方法。

利用荧光油膜在激光照射下的显色反应,由荧光油膜灰度表征出油膜厚度这一关键因素,进而通过得到油膜厚度解算该处的摩阻分布。早在1977年,由Tanner[

1]提出通过油膜的激光干涉来测量表面摩阻,为摩阻的简化测量提供了坚实的基础。1993年,Monson[2]对油膜干涉法进行完善,油膜干涉理论发展成熟。而1999年,Brown[3]精确建立了在流体表面摩擦应力作用下的油膜演化方程,形象地表征出了油膜厚度随时间的变化率与其自身体积力、摩阻等因素的关系。对使用油膜进行全局摩阻的直接测量提供了可能性。在2008年,由Liu[4]提出用荧光油膜的灰度成像来表征油膜厚度这一新方法,建立了荧光油膜灰度与厚度的线性模型,该模型提供了理论基础,但还无法直接用于工程测量中。2011年,李鹏[5]将荧光油膜灰度与厚度线性方程引入相对表面剪切力模型,得出相对剪切力与成像平面灰度分布间的映射模型,为进行全局剪切力直接测量方法的研究和建立求解模型提供了有利的基础。2012年,李[6]提供了较为简便的数据采集方法,他通过用盖玻片和载玻片等高透光光学玻璃构成斜面进行荧光油膜灰度、厚度标定数据的采集,并初步验证了该线性关系的正确性。

在标定过程中,一般选择插值法或数据拟合的方式来处理数据使其平滑。传统的插值法或者数据拟合的方式在一般厚度(如毫米级)的条件下比较实用,但是在面对高精度(如微米级)时,则需要进行算法改进来达到精度要求,其改进的过程较为繁琐,也存在着一定的困难,而两者都有各自的缺点,比如数据拟合方式,由于该方法不一定通过所有数据点,在精度方面也存在一定的局限性;传统插值法在靠近数据两端点插值时,会存在“龙格现象[

7],为了避免该问题,则采用分段插值形式。

本次采用神经网络对荧光油膜建立模型,针对神经网络对单输入、单输出模型有较好预测的特点,与传统插值法结果进行比对,验证了通过神经网络建立预测模型,有精度更高、预测更准确的优势。

1 实验理论分析

1.1 插值法原理

本次采用一元全区插值法进行荧光油膜灰度与厚度关系的模型建立,其数学模型表示为

z=i=kk+7yij=kjik+7[(t-xj)/(xi-xj)] (1)

式中:xi从属于n个结点,xj为不同于xi的其他结点。通过7次拉格朗日插值计算求得插值点t处的函数近似值z,且插值点t在该n个结点中(n=8)选取满足

xk<xk+1<xk+2<xk+3<t<xk+4<xk+5<xk+6<xk+7 (2)

当插值点t位于区间外时,则属特殊情况,将取该端的n/2个结点进行插值计算。

1.2 BP神经网络

BP神经网络为多层前馈网络,遵循误差反向传播训练机制,围绕梯度下降法中心思想展开,采用梯度搜索方[

8],最终目标是使实际输出值与期望输出值的误差均方值达到最小。

BP网络结构如图12所示。图中k为迭代次数,u(k)y(k)为网络的输入。y(k)为受控对象实际输出,yn(k)为网络的输出。参考文献[

9]的方法,将受控对象实际输出y(k)及输入u(k)的值作为整个网络的输入,并将系统输出与网络输出的误差目标函数作为神经网络调参标准,从而搭建网络。

图1  BP神经网络逼近的结构

Fig.1  Structure of BP neural network approximation

图2  用于逼近的BP网络

Fig.2  BP network for approximation

前向传播是网络的输出运算过程,隐层神经元的输入过程实际为所有输入的加权函数求和过程,令隐层神经个数为q,则有

xj=i=1qwijxi (3)

式中:wij为网络层间连接权值;xj为限层结点;xj为输入层结点。

yn(k)=xk,则网络的输出与预想输出的误差为

e(k)=y(k)-yn(k) (4)

误差性能指标函数为

E=12e(k)2 (5)

逆向传播通过δ学习算法,对各层间的权值进行调整。输出层与隐层的连接权值计算为

Δwj2=-ηEwj2=ηe(k)xkwj2=ηe(k)xj' (6)

式中xj'为激发函数。k+1时刻网络的权值为

wj2(k+1)=wj2(k)+Δwj2 (7)

隐层及输入层连接权值计算为

Δwij=-ηEwij=ηe(k)×ynwij (8)

式(8)中的权值梯度为

ynwij=ynxj'×xj'xj×xjwij=wj2×xj'xjxi=wj2×xj'(1-xj')×xi (9)

k+1时刻网络权值为

wij(k+1)=wij(k)+Δwij (10)

当上一次的权值对本次权值产生了影响,就需加入约束因子α,此时的权值有

wj2(k+1)=wj2(k)+Δwj2+α(wj2(k)-wj2(k-1)) (11)
wij(k+1)=wij(k)+Δwij+α(wij(k)-wij(k-1)) (12)

式中:η为学习速率,α为约束因子,η0,1α0,1。雅可比矩阵(输入输出灵敏度信息)算法表达式为

ykukynkuk=ynkxj'×xj'xj×xjx1=jwj2xj'1-xj'w1j (13)

1.3 ELM机器学习算法辨识原理

极限学习机(Extreme learning machine, ELM)神经网络较BP神经网络具有更高的精度和运算速[

10]。令有N个任意样本(Xi,Ti),其中Xi=[xi1,xi2,,xin]TRnTi=[ti1,ti2,,tim]TRm,则对于有P个隐层结点的单隐层神经网络为

i=1Pβig(WiXj+bi)=oj      j=1,,N (14)

式中:g(x)为激活函数,Wiβi分别为输入、输出权重,bi为偏置。该网络目标函数为

j=1Noj-tj=0 (15)

则可得矩阵形式

Hβ=K (16)

式中:H为隐层输出,β为输出权重,K为期望值。训练过程中得到

H(Wi,bi)βi-K=minH(Wi,bi)βi-K (17)

可得损失目标函数

E=j=1N(i=1Pβig(WiXj+bi)-tj)2 (18)

1.4 Hopfield神经网络参数辨识原理

假定用于参数辨识的约束函数与Hopfield网络约束函数相等,对其状态方程矩阵进行分解变换可得网络的权值矩阵W和神经元外输入矩阵I,将其代入网络状态方程式进行运算,收敛状态可得到较稳定的辨识参数结果。令一个系统的状态方程为

x˙=Ax+Bu (19)

式中:AB为待辨识的参数矩阵,取P=A11A12A21A22B21B22Tx=x1x2T,令x为状态矢量,u为控制输入,则其参数辨识过程就是对向量p的求解过程。设可调系统为

x˙p=θ1x+θ2u (20)

式中:θ1= a11a12a21a22 θ2= b1b2 ,取V=a11a12a21a22b21b22T,则可得到状态偏差方程

e˙=(A-θ1)x+(B-θ2)u (21)

定义基于状态偏差变化率的系统参数辨识函数为

E=12e˙Te˙ (22)
e˙Te˙=[x˙-θ1x-θ2u]T[x˙-θ1x-θ2u] (23)
E=12[x˙-θ1x-θ2u]T[x˙-θ1x-θ2u] (24)

Hopfield网络约束函数趋于收敛的过程,就是估计矩阵θ1θ2接近实际矩阵AB的过[

11]。网络第i个神经元的动态微分方程满足

Ciduidt=j=1nwijvj-uiRi+Iivi=g(ui) (25)

式中:g(ui)=ρ1-e-λu1+e-λuρ>0Ci为调整系数;Ri为网络收敛因子。

Ri,并取Ci=1,则状态方程可以转换为

duidt=i=1nwijvj+Ii (26)

Hopfield网络的约束函数为

EN=-12i=1nj=1nwijvivj+i=1n1Ri0vigi-1vdv-i=1nIivi (27)

由于Ri,取Hopfield网络的输出等于待辨识参数,则

EN=-12i=1nj=1nwijvivj-i=1nIivi (28)

式(26)式(28)可以看出,有

ENvi=-i=1nwijvj-Ii (29)

取参数辨识的约束函数与Hopfield网络约束函数相[

12⁃13],即E=EN,则由式(29)可得

Evi=-i=1nwijvj-Ii (30)

通过E=E1+E2,可以换算为

dEdV=dE1dV+dE2dV=-WV-I (31)
dEdV=E1v1E1v6+E2v1E2v6 (32)

通过式(32)解出权值WI

W=-x12x1x200x1u0x2x1x2200x2u000x12x1x20x1u00x2x1x220x2uux1ux200u2000ux1ux20u2 (33)
I=x1x˙1x2x˙1x1x˙2x2x˙2ux˙1ux˙2T (34)

将式(3334)的WI代入式(26),可得到稳定的ui,最终辨识结果的输出为

F=g(u1)      G=g(u2) (35)

2 实验结果与讨论

2.1 荧光油膜标定过程

荧光油膜的标定过程是一个数据实时采集的过程,本次通过搭建载玻片斜面的方法进行油膜灰度与厚度相关信息的采[

6],为了方便像素点的读取和像素单位的转换,对相机采集的荧光油膜灰度图进行了网格化处理,如图34所示。图4中间的红色线条则表示本次的荧光油膜灰度与厚度数据采集部分,且该网格xy轴像素间隔均为16像素。

图3  荧光油膜灰度图

Fig.3  Grayscale of fluorescent oil film

图4  油膜网格图

Fig.4  Oil film grid

由盖破片和载玻片构成的斜面所采集的数据并不是直接为荧光油膜灰度与厚度的关系数据,而是灰度与横坐标x长度的关系,还需进行三角变换处理,如图5所示。可以得到

h=xtank (36)

图5  斜面三角变换

Fig.5  Transformation of slope trigonometric

采集的初始数据经几何变换处理,如表12所示。

表1  荧光油膜灰度及厚度初始数据
Table 1  Initial grayscale and thickness of fluorescence oil film ( mm )
序号12345678910
油膜灰度值 32 34 36 40 41 44 47 52 54 57
油膜厚度值 0.000 00 0.005 53 0.006 64 0.007 14 0.010 51 0.014 29 0.021 43 0.026 55 0.028 57 0.033 19
序号 11 12 13 14 15 16 17 18 19 20
油膜灰度值 59 62 65 67 71 72 74 77 80 82
油膜厚度值 0.038 16 0.042 68 0.049 23 0.050 01 0.062 50 0.063 05 0.071 44 0.078 58 0.096 79 0.100 01
序号 21 22 23 24 25 26 27 28 29 30
油膜灰度值 84 86 88 90 92 93 95 97 99 100
油膜厚度值 0.114 3 0.128 59 0.148 78 0.157 16 0.171 46 0.185 73 0.200 02 0.274 89 0.285 74 0.300 03
表2  荧光油膜灰度⁃厚度测试数据
Table 2  Test data of grayscale and of thickness of fluorescence oil film
序号12345678910
油膜灰度值 33 35 43 58 68 79 85 87 96 101
油膜厚度值 0.003 87 0.006 08 0.012 72 0.035 95 0.057 15 0.085 72 0.121 44 0.142 87 0.228 60 0.347 35

2.2 神经网络预测

2.2.1 BP神经网络预测

将标定实验得到荧光油膜灰度与厚度初始数据(表1)作为神经网络的训练数据来源,神经网络设置的结构为1⁃4⁃4⁃1⁃1,网络的初始权值取随机值,其训练过程如图6所示。网络的训练次数为500次,训练的神经网络、结果误差图分别如图78所示,实验中可以通过不断调整训练次数、最小误差以及权值等参数来使神经网络模型达到最优值。

图6  训练过程图

Fig.6  Training process diagram

图7  神经网络训练

Fig.7  Neural network training

图8  误差曲线

Fig.8  Error curves

将训练好的模型通过表2中的测试数据来进行预测,其预测结果如图9所示。

图9  BP神经网络预测结果

Fig.9  Prediction results of BP neural network

2.2.2 ELM神经网络预测

ELM网络较BP网络具有运算速度快的优势,且精度较高,训练模型的预测结果如图10所示。

图10  ELM神经网络预测结果

Fig.10  Prediction results of ELM neural network

2.2.3 对比分析

绘制BP神经网络、ELM网络及传统插值法预测曲线,如图11所示,然后将预测值与真实荧光油膜厚度值进行比较便可以得到预测误差,分析神经网络的预测精度。各算法预测数据结果如表3所示。

图11  数据预测对比

Fig.11  Comparison of data prediction

表3  神经网络与插值法预测数据及其误差比较
Table 3  Comparison of prediction data and error of neural network and interpolation method
测试灰度

真实

厚度/mm

BP预测

厚度/mm

ELM预测

厚度/mm

插值预测

厚度/mm

BP误差/%ELM误差/%插值误差/%
33 0.003 87 0.003 91 0.003 62 0.003 39 1.033 6.460 12.40
35 0.006 08 0.004 83 0.006 63 0.006 75 20.56 9.046 11.02
43 0.012 72 0.011 90 0.012 12 0.013 81 6.447 4.171 8.569
58 0.035 95 0.034 52 0.035 54 0.035 67 3.978 1.140 0.779
68 0.057 15 0.054 71 0.053 29 0.053 59 4.270 6.750 6.229
79 0.085 72 0.088 16 0.090 43 0.092 02 2.846 5.495 7.350
85 0.121 44 0.117 21 0.118 99 0.120 90 3.483 2.070 0.445
87 0.142 87 0.140 96 0.131 86 0.139 26 1.337 7.706 2.527
96 0.228 60 0.240 36 0.233 49 0.233 67 2.783 2.139 2.218
101 0.347 35 0.319 17 0.324 68 0.320 23 8.113 6.527 7.808
平均误差/% 5.485 5.150 5.935

表3中的数据可以直观看出,ELM神经网络的误差低达5.150%,而传统插值法的误差最大,为5.935%,故通过神经网络建立预测模型具有一定的优势,且无需分析模型理论结构;与此同时,BP神经网络的误差为5.485%,相较ELM网络的误差较高,故ELM神经网络在模型预测方面比BP神经网络占据一定优势,且运算速度快;经实验分析,通过神经网络建立预测模型,较传统插值法的精度有一定的提高。

2.3 Hopfield神经网络参数辨识

Hopfield神经网络进行模型的参数辨识,需要根据经验公式或者手册给出初定的参数[

13]。引用Liu[4]提出的数学模型,有

I(x,y)=aIex(x,y)h(x,y) (37)

式中:I(x,y)为荧光亮度(灰度值);Iex(x,y)为初始光强;h(x,y)为荧光油膜的厚度;a为荧光效率等参数。因为影响油膜灰度值的因素不止油膜本身的厚度因素,还有周围环境因素以及油膜的自身特性,如浓度,所以对该线性模型做变换,可以得到

h=(n1+n2+n3+)I+(m1+m2+) (38)

因影响因素可能为非线性系统和常数,所以nimj系数定义为影响因素的权重系数,根据经验式,可以得到假定影响因素参数,如表4所示。

表4  影响因素的假定权重参数值
Table 4  Assumed weight parameter values of influencing factors
影响因素符号类型假定权重初始值
功率 W 非线性 0.13
光距 R 非线性 0.21
曝光系数 T 常数 33.3

通过Hopfield神经网络参数辨识,可得以下辨识结果,如图12~14所示。

图12  紫外光源功率参数辨识

Fig.12  Power parameter identification of UV light source

图13  紫外光源距离参数辨识

Fig.13  Distance parameter identification of UV light source

图14  曝光系数参数辨识

Fig.14  Identification of exposure coefficient parameters

通过Hopfield神经网络所辨识的结果参数如表5所示。

表5  辨识后的参数值
Table 5  Parameter value after identification
影响因素初始权重辨识后的权重误差百分比/%
功率 0.13 0.128 9 0.85
光距 0.21 0.207 8 1.05
曝光系数 33.3 33.129 6 0.51

通过Hopfield神经网络进行参数辨识,能够很好地找出功率、光距和曝光系数等参数对荧光油膜灰度成像的影响程度,且辨识结果误差范围在1%左右,达到实际工程应用的标准。

3 结 论

经过实验分析论证,ELM神经网络模型、BP神经网络模型及插值法模型的预测误差分别为5.150%、5.485%和5.935%。通过神经网络搭建模型较传统插值法精度有所提高;而ELM神经网络具备比BP神经网络运算速度快、精度更高的优势。利用Hopfield神经网络进行参数辨识,其辨识得到参数的误差均控制在1%左右,后续可以继续调整Hopfield神经网络的参数,使误差率达到期望的效果。与传统插值法进行对比,利用神经网络搭建模型可以获得更高的精度,为后续荧光油膜灰度与厚度研究提供了一种可行的方法。

References

1

Tanner L. A skin friction meter, using the viscosity balance principle,suitable for use with flatrocurved metal surfaces[J]. Phys E:Scientific Instr197710278-284. [百度学术] 

2

Monson D JMatree G G. Boundary-layer transition and global skin friction measurement with an oil-fringe imaging technique[J]. Society of Automotive Engineers19931021829-1843. [百度学术] 

3

Brown J LNaughton J W. The thin oil film equation[M]. Washington DCNASA_AMES Research Center1999. [百度学术] 

4

Liu TianshuMontefort J. Global luminescent oil-friction meter[J]. AIAA Journal2008462):476-485. [百度学术] 

5

李鹏明晓.风力机叶片的全局表面摩擦力测量的荧光油膜法[J].南京航空航天大学学报2011435):581-585. [百度学术] 

Li PengMing Xiao. Fluorescent oil film method for global surface friction measurement of wind turbine blades[J]. Journal of Nanjing University of Aeronautics & Astronautics2011435): 581-585. [百度学术] 

6

李鹏.全局表面摩擦应力直接测量技术研究[D].南京南京航空航天大学2012. [百度学术] 

Li Peng. Research on direct measurement technology of global surface friction stress[D]. NanjingNanjing University of Aeronautics and Astronautics2012. [百度学术] 

7

张雨浓李名鸣陈锦浩. 龙格现象难题破解之系数与阶次双确定方法[J]. 计算机工程与应用2013493): 44-49. [百度学术] 

Zhang YunongLi MingmingChen Jinhaoet al. Double determination of coefficient and order for solving the Runge problem[J]. Computer Science and Technology2013493): 44-49. [百度学术] 

8

刘冰心王宁张冬. 基于BP神经网络的智能电网配电系统改进算法的研究[J]. 现代电子技术20123521): 143-144,148. [百度学术] 

Liu BingxinWANG NingZHANG Dong. Research on improved algorithm of smart grid distribution system based on BP neural network[J]. Modern Electronic Technology20123521): 143-144,148. [百度学术] 

9

霍冠良宁志华.基于BP神经网络的含褶皱复合材料强度预测[J].南京航空航天大学学报2020523): 460-467. [百度学术] 

Huo GuanliangNing Zhihua. BP neural network for strength prediction of folded composites[J]. Journal of Nanjing University of Aeronautics & Astronautics2020523): 460-467. [百度学术] 

10

Zhai JShao QWang X. Architecture selection of ELM networks based on sensitivity of hidden nodes[J]. Neural Processing Letters2016442): 471-489. [百度学术] 

11

卓新建章祥荪徐宗本. 关于Hopfield-型神经网络收敛性和求能量函数极小解的可靠性的研究[J]. 系统科学与数学20012): 152-162. [百度学术] 

Zhuo Xinjian. ZHANG Xiangsun, XU Zongben. Research on convergence of Hopfield-type neural networks and reliability of minimum solutions of energy functions[J]. Systems Science and Mathematics20012): 152-162. [百度学术] 

12

章毅周明天王平安. 关于神经网络的能量函数 [J]. 计算机研究与发展1999367): 794-799. [百度学术] 

Zhang YiZHOU MingtianWANG Pingan. On the energy function of neural networks[J]. Computer Research and Development1999367): 794-799. [百度学术] 

13

郎世城. 基于Hopfield神经网络的船舶电网故障诊断研究[D].重庆重庆交通大学2019. [百度学术] 

Lang Shicheng.Research on ship power grid fault diagnosis based on Hopfield neural network[D]. ChongqingChongqing Jiaotong University2019. [百度学术] 

14

孙跃谭晶晶唐春森.IGBT机理建模及其基于神经网络的参数辨识方法[J].西南交通大学学报2015506): 1143-1149,1163. [百度学术] 

Sun YueTan JingjingTang Chunsen. IGBT mechanism modeling and parameter identification me⁃thod based on neural network[J]. Journal of Southwest Jiaotong University2015506): 1143-1149, 1163. [百度学术] 

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!