发明名称 强脉冲噪声环境下基于方差比盲分离的复杂结构材料高精度透视检测方法
摘要 本发明属于激光波数调频干涉方法技术领域,尤其是涉及一种基于方差比盲分离的复杂结构材料高精度透视检测方法,适用于工业环境下,树脂基复合材料、多层半导体器件等内部结构的高精度检测。本发明提出强脉冲噪声环境下基于方差比盲分离方法的复杂结构材料高精度透视检测,能够在强脉冲噪声环境下,通过对干涉光强序列进行斯皮尔曼相关分析与傅里叶变换相结合,检测出多层复杂结构材料的内部结构信息,克服了传统的傅里叶变换方法不能对复杂结构材料的内部结构信息进行盲分离检测的缺陷。
申请公布号 CN105928877B 申请公布日期 2017.03.15
申请号 CN201610269852.X 申请日期 2016.04.27
申请人 广东工业大学 发明人 谢胜利;吕梓亮;何昭水;周延周;周郭许
分类号 G01N21/00(2006.01)I;G06F17/14(2006.01)I 主分类号 G01N21/00(2006.01)I
代理机构 广州市红荔专利代理有限公司 44214 代理人 张文
主权项 强脉冲噪声环境下基于方差比盲分离的复杂结构材料高精度透视检测方法,其特征在于,基于方差比盲分离方法在强脉冲噪声环境下对复杂结构材料进行高精度透视检测,具体步骤为:1)调频激光器波数输出,使激光照射复杂结构材料上,与此同时,CCD相机拍摄N<sub>s</sub>张干涉图像;拍摄第n<sub>s</sub>张干涉图像坐标(x,y)处的光强表达式为:<maths num="0001"><math><![CDATA[<mrow><mi>I</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msqrt><mrow><msub><mi>I</mi><mi>p</mi></msub><mo>&CenterDot;</mo><msub><mi>I</mi><mi>q</mi></msub></mrow></msqrt><mi>c</mi><mi>o</mi><mi>s</mi><mo>&lsqb;</mo><mn>2</mn><mo>&CenterDot;</mo><msub><mi>&Lambda;</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>&CenterDot;</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150483870000011.GIF" wi="1622" he="134" /></maths>其中,I<sub>p</sub>,I<sub>q</sub>为复杂结构材料中界面p,q的反射光强,Λ<sub>pq</sub>为界面p,q的光程差,M包括复杂结构材料和干涉参考面的总层数,k(n<sub>s</sub>)为拍摄第n<sub>s</sub>张干涉图像时对应的激光波数输出;2)对干涉图像每个像素点构造光强向量I<sub>1</sub>=[I(1),I(2),…,I(N<sub>s</sub>‑1)]<sup>T</sup>,I<sub>2</sub>=[I(2),I(3),…,I(N<sub>s</sub>)]<sup>T</sup>,T表示矩阵转置;利用斯皮尔曼相关系数,对I<sub>1</sub>构造自相关矩阵C<sub>xx</sub>,对I<sub>1</sub>,I<sub>2</sub>构造互相关矩阵C<sub>xy</sub>;3)对C<sub>xx</sub>求解关于C<sub>xy</sub>的广义特征值:<maths num="0002"><math><![CDATA[<mrow><msub><mi>&eta;</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&pi;</mi><mo>&CenterDot;</mo><mfrac><mrow><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>k</mi></mrow><mrow><msub><mi>N</mi><mi>s</mi></msub><mo>-</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mi>&xi;</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&pi;</mi><mo>&CenterDot;</mo><mfrac><mrow><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>k</mi></mrow><mrow><msub><mi>N</mi><mi>s</mi></msub><mo>-</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150483870000012.GIF" wi="1622" he="143" /></maths>上式中,j为虚数单位,Δk为激光波数扫频范围,f<sub>pq</sub>为界面p,q的干涉信号频率;4)通过式(2)计算界面p,q的干涉信号频率f<sub>pq</sub>为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>N</mi><mi>s</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo><mo>&CenterDot;</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>(</mo><msub><mi>&eta;</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn><mi>&pi;</mi><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>k</mi></mrow></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150483870000013.GIF" wi="1622" he="126" /></maths>5)对单点像素干涉光强信号作关于波数的傅里叶变换并求解界面p,q的干涉相位,<maths num="0004"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>&lsqb;</mo><mfrac><mrow><mi>Im</mi><mrow><mo>(</mo><mover><mi>I</mi><mo>~</mo></mover><mo>(</mo><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mi>Re</mi><mrow><mo>(</mo><mover><mi>I</mi><mo>~</mo></mover><mo>(</mo><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0001150483870000014.GIF" wi="541" he="166" /></maths>上式中,<maths num="0005"><math><![CDATA[<mrow><mover><mi>I</mi><mo>~</mo></mover><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>n</mi><mi>s</mi></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>s</mi></msub></munderover><mi>I</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mi>j</mi><mo>&CenterDot;</mo><mn>2</mn><mi>&pi;</mi><mo>&CenterDot;</mo><msub><mi>f</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>&CenterDot;</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&CenterDot;</mo><mi>W</mi><mo>&lsqb;</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150483870000015.GIF" wi="1638" he="149" /></maths>上式中,W为窗函数;6)对CCD相机拍摄的干涉图像每个像素点重复进行步骤2)~5),得到界面p,q的全幅面干涉相位,相位解卷绕后即可在强脉冲噪声环境下测量出材料的复杂结构信息。
地址 510006 广东省广州市番禺区广州大学城外环西路100号