CN1606235A - He-Ne激光器拍波噪声滤除装置及其噪声滤除方法 - Google Patents

He-Ne激光器拍波噪声滤除装置及其噪声滤除方法 Download PDF

Info

Publication number
CN1606235A
CN1606235A CN 200410081242 CN200410081242A CN1606235A CN 1606235 A CN1606235 A CN 1606235A CN 200410081242 CN200410081242 CN 200410081242 CN 200410081242 A CN200410081242 A CN 200410081242A CN 1606235 A CN1606235 A CN 1606235A
Authority
CN
China
Prior art keywords
signal
noise
filter
noise filtering
convergence
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
CN 200410081242
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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN 200410081242 priority Critical patent/CN1606235A/zh
Publication of CN1606235A publication Critical patent/CN1606235A/zh
Pending legal-status Critical Current

Links

Images

Abstract

He-Ne激光器拍波噪声滤除装置及其噪声滤除方法,其特点是将用双纵模双频激光干涉仪以拍波干涉的原理对高精度大尺寸零件作无导轨测量装置中节点判别光路输出的大尺寸测量信号和参考信号分别通过高频放大器(1)连接功分器(2),功分器连接混频器(4),其中一路功分器与混频器间接有延迟线(3),得到两路有固定相位差的低频信号。利用此两路低频信号中拍波演变出的信号高度相关,而所含噪声不相关的特点,将其依次输入基于数字信号处理器(DSP)(7)、硬件的变步长自适应滤波器(5)和改进型小波滤波器(6),最终得到大大滤除了激光器噪声的高信噪比拍波信号。

Description

He-Ne激光器拍波噪声滤除装置及其噪声滤除方法
一、技术领域
本发明涉及大尺寸零件无导轨测量的He-Ne激光器拍波噪声滤除装置及其噪声滤除方法,属于激光精密测试领域。
二、背景技术
在用双纵模双频激光干涉仪以拍波干涉的原理对高精度大尺寸零件作无导轨测量的技术中,He-Ne激光器拍波本身的噪声是限制测试精度的主要因素。由于该噪声的广谱性及随机性,一般的滤波方法的效果是很有限的。为此提出了各种自适应滤波器算法和小波变换去噪算法,应该说这些方法对消除拍波噪声有一定效果。但这些自适应滤波器均假定有主通道和参考通道两个信号输入通道,主通道接收到信号s和不相关噪声n0,即s+n0形成对消除器的原输入;参考通道接收到一个与信号不相关的、而与噪声n0相关的噪声n1,形成对消除器的参考输入,参考输入的噪声n1被滤波器产生一个输出y,它尽可能地去翻版n0,该输出从原输入s+n0减去,产生系统的输出为e=s+n0-y。实际系统中很难找到完全符合要求的参考通道输入,因而满足要求的假定条件经常是不现实的。而且这些自适应滤波器多采用最小均方误差算法(LMS),它具有结构简单和性能稳定等优点,在自适应信号处理领域得到了广泛的应用。但是LMS算法的稳态误差与步长成正比,收敛速度与步长成反比,这一内在的矛盾使得LMS算法的收敛速度和稳态误差特性不能同时满足,必须在收敛速度和稳态误差这两个性能指标之间进行权衡,而这通常依靠大量的实验与经验,要达到较好的效果需要大量的时间、精力,并对使用者有较高的素质要求。而小波变换去噪一般采用阈值去噪法。这类方法由D.L.Dohono在1995年首先提出,此方法根据信号与噪声的幅值和频率不同,对含噪信号进行小波变换得到小波系数序列,通过阈值去除小波系数中的噪声元素,用经过处理的小波系数进行小波反空域相关滤波法。此方法根据Mallat等的研究结果,信号成分的奇异性具有正的Lipschitz指数,其小波系数的模值将随着尺度的增大而增大。随机噪声的Lipschitz指数则为负数,其模值将随着尺度的增大而迅速减小。将相邻尺度的小波变换值的相关量进行归一化处理并与小波变换值比较来判断噪声与信号,以噪声在各尺度上的方差作为终止迭代的标准。此方法能够较好地保留信号中的细节成分。但此法在迭代收敛时,如果剔除的小波系数的方差和本尺度上的噪声方差之比恰好等于1,那么我们认为剔除的小波系数都是噪声系数;如果小于1,将有很多噪声被当成信号保存下来,从而降低滤噪效果。在达到收敛时,这个比值总是小于1,因此会引入随机噪声。
在用双纵模双频激光干涉仪以拍波干涉的原理对高精度大尺寸零件作无导轨测量的技术中,节点的确定有两路大尺寸拍波信号:一路为被测量调制相位的测量信号;一路为相位不变的参考信号。通过对它们比相找零得到测量结果的。由于拍波信号频率很高(与He-Ne激光器长短有关,可达500MHz~1000MHz),很难对其有效滤波,因而直接比相找零测量的精度受限于较低的信噪比,通常很难得到较高的测量精度。
三、发明内容
本发明的目的是针对现有技术的不足而提供一种He-Ne激光器拍波噪声滤除装置及其噪声滤除方法,其特点是采用混频器将拍波信号频率在保留被测相位信息的条件下降至低频,形成两路有固定相位差的低频信号,将其作为自适应滤波器的两路输入,以变步长LMS自适应算法实现自适应滤波,再以能保证迭代收敛时剔除的小波系数的方差和本尺度上的噪声方差之比恰好等于1的改进型小波去噪算法滤波,经以上环节滤除激光器拍波噪声后进行比相找零测量,从而得到高的测量精度。
本发明的目的由以下技术措施实现的。
测量信号和参考信号分别通过高频放大器连接功分器,功分器连接混频器,其中一路功分器与混频器间接有延迟线,得到两路有固定相位差的低频信号。
将上述两路有固定相位差的低频信号输入自适应滤波器。两路信号X0和X1均含有背景噪声,即X0=S0+N0,X1=S1+N1,由于拍波演变出的信号S0与S1高度相关,而噪声N0与N1不相关,自适应滤波器的作用是通过调节滤波器的系数,使输出Yn逼近期望的信号S0。可以证明:
E [ e n 2 ] = E [ ( S 0 - Y n ) 2 ] + E [ N 0 2 ]
使均方误差取最小值,则:
min E [ e n 2 ] = min E [ ( S 0 - Y n ) 2 ] + E [ N 0 2 ]
上式表明:当误差en的平均功率为最小时,输出Yn则为S0最好的最小均方估计。
为同时满足收敛速度和稳态误差这两个性能指标之间相互矛盾的要求,采用变步长自适应滤波算法,其迭代公式为:
Y(n)=X(n)W(n)
e(n)=d(n)-Y(n)
W(n+1)=W(n)+2μ(n)e(n)X(n)
μ ( n + 1 ) = α μ ( n ) + γ e ( n ) 2
式中Y(n)为实际输出,X(n)为输入信号,W(n)为权系数,X(n)=[x(n),x(n-1),x(n-2),...x(n-L+1)]T表示时刻n的输入信号矢量,T为矢量标志,由最近L个信号采样值构成,W(n)=[ω0(n),ω1(n),ω2(n),...ωL-1(n)]T是时刻n自适应滤波器的权系数矢量。L是自适应滤波器的阶数,d(n)是期望输出值,在这里为S0,e(n)是误差,μ(n)是控制稳定性和收敛速度的参量,称之为步长因子,变步长因子μ(n)迭代公式中,0<α<1,γ>0。
这种变步长自适应滤波算法在初始收敛阶段或未知系统参数发生变化时,步长应比较大,以便有较快的收敛速度或对时变系统的跟踪速度;而在算法收敛后,不管主输入端干扰信号有多大,都应保持很小的调整步长以达到很小的稳态失调噪声。这种算法能同时满足收敛速率、跟踪速率及权失调噪声之间相互矛盾的要求。
为满足信号处理实时性的要求,该滤波器通过数字信号处理器(DSP)高速运算实现。基于DSP硬件的变步长自适应滤波器大大提高了滤波器的时效性,可以实时处理0~几MHz的信号,有效滤除激光器拍波信号中的噪声。
将经过变步长自适应滤波器处理输出的信号输入其后的改进型小波滤波器,该滤波器首先利用空域相关法剔除噪声系数,在噪声系数λ>1时反复迭代,直到λ<1。此时停止迭代将引人噪声。为此在这一步引人调整系数c,寻找最佳的调整系数c,使λ=1,我们注意到,当λ>1时,lnλ>0;当λ<1时,lnλ<0。利用此性质由下式自适应寻找c的最佳估计:
c→c+μlnλ
在此式中,μ是调整步长的一个正的常数。设置c的初值为1,当λ<1时,lnλ<0,迭代后c值将减小,用此减小的c值去乘相关系数,这样做抽取的细节点数将减少,而剔除的小波系数的方差将增加,λ也将增加。如果λ的值仍旧小于1,则λ代人上式继续再次迭代,剔除的小波系数的方差将再次增加,λ也将再次增加,如果迭代后入的值大于1,lnλ>0,迭代后c值将增大,用此增大的c值去乘相关系数,这样做抽取的细节点数将增加,而剔除的小波系数的方差将减小,λ也将减小。可见,无论λ>1或λ<1,本文的方法都能使λ趋向1的方向靠近。这样一直迭代下去,直到入λn+1n≤σm×10-3。这种迭代的结果使剔除的小波系数的方差等于噪声的方差,可以达到滤波的最佳效果。常数μ的作用是控制稳定性和收敛速度的参量,可以叫做步长因子,μ值较大时,c可以较快地收敛到平稳值,同时出现较大的随机抖动。μ较小时,c值可以较为平稳和精确但又较为缓慢地收敛到平稳值。
上述改进型小波滤波器采用自适应调整的方法,在小波变换值的相关量里加入自适应微调系数,使在迭代收敛时剔除的小波系数的方差和本尺度上的噪声方差之比恰好等于1,避免引入随机噪声,从而提高了空域相关法的滤波效果。同样为满足信号处理实时性的要求,该改进型小波滤波器也通过数字信号处理器(DSP)高速运算实现。
本发明具有以下优点:
1.将大尺寸拍波测量信号和参考信号分别通过高频放大器连接功分器,功分器连接混频器,其中一路功分器与混频器间接有延迟线,使拍波信号在保留被测相位信息的条件下降至低频,形成两路有固定相位差的低频信号,易于后续数字信号电路处理。
2.可将上述两路有固定相位差的低频信号X0和X1直接输入变步长自适应滤波器,避免了一般自适应滤波器均需要有参考输入通道,且参考输入通道只含有与主通道信号不相关、而与主通道噪声相关的噪声的苛刻要求。X0和X1均含有背景噪声,即X0=S0+N0,X1=S1+N1,由于拍波演变出的信号S0与S1高度相关,而噪声N0与N1不相关,自适应滤波器通过调节滤波器的系数,能使输出Yn逼近期望的信号S0。具体采用通过数字信号处理器(DSP)高速运算实现的变步长自适应滤波算法,该算法在初始收敛阶段或未知系统参数发生变化时,步长应比较大,以便有较快的收敛速度或对时变系统的跟踪速度;而在算法收敛后,不管主输入端干扰信号有多大,都应保持很小的调整步长以达到很小的稳态失调噪声。这种算法同时满足了收敛速率、跟踪速率及权失调噪声之间相互矛盾的要求,可以有效滤除激光器拍波信号中的噪声。
3.将经过变步长自适应滤波器处理输出的信号输入其后的改进型小波滤波器,再次滤波。该滤波器利用空域相关法剔除噪声系数,为避免λ<1时停止迭代引人随机噪声,引人调整系数c,利用c→c+μlnλ自适应寻找最佳的调整系数c,使λ=1,提高了空域相关法的滤波效果。
4.通过低频化处理及双次滤波,可以得到大大滤除了激光器噪声的高信噪比拍波信号进而获得更高的比相测量精度。
四、附图说明
图1为He-Ne激光器拍波噪声滤除装置框图。
1高频放大器,2功分器,3延迟线,4混频器,5变步长自适应滤波器,6改进型小波滤波器,7数字信号处理器(DSP)。
图2为变步长自适应滤波器框图。
8自适应滤波器,9合成器
五、具体实施方式
下面通过实施例对本发明进行具体描述,有必要在此指出的是以下实施例只用于对本发明进行进一步说明,不能理解为对本发明保护范围的限制,该领域的技术熟练人员可以根据上述本发明的内容对本发明做出一些非本质的改进和调整。
实施例:
在用双纵模双频激光干涉仪以拍波干涉的原理对高精度大尺寸零件作无导轨测量的装置中,节点判别光路输出两路大尺寸拍波信号,如图1所示:一路为被测量调制相位的测量信号,一路为相位不变的参考信号,拍波信号频率很高(与He-Ne激光器长短有关,可达500MHz~1000MHz)。将测量信号和参考信号分别通过高频放大器1连接功分器2,功分器连接混频器4,其中一路功分器与混频器间接有延迟线3,得到两路有固定相位差的低频信号X0和X1,将此两路信号X0和X1输入变步长自适应滤波器5,该滤波器通过数字信号处理器(DSP)7高速运算实现。此两路有固定相位差的低频信号X0和X1均含有背景噪声,即X0=S0+N0,X1=S1+N1,由于拍波演变出的信号S0与S1高度相关,而噪声N0与N1不相关,自适应滤波器通过调节滤波器的系数,能使输出Yn逼近期望的信号S0。具体采用变步长自适应滤波算法,该算法在初始收敛阶段或未知系统参数发生变化时,步长应比较大,以便有较快的收敛速度或对时变系统的跟踪速度;而在算法收敛后,不管主输入端干扰信号有多大,都应保持很小的调整步长以达到很小的稳态失调噪声。这种算法同时满足了收敛速率、跟踪速率及权失调噪声之间相互矛盾的要求,可以有效滤除激光器拍波信号中的噪声。再将变步长自适应滤波器输出信号Yn输入同样用数字信号处理器(DSP)高速运算实现的改进型小波滤波器6,该滤波器采用自适应调整的方法,在小波变换值的相关量里加入自适应微调系数,使在迭代收敛时剔除的小波系数的方差和本尺度上的噪声方差之比恰好等于1,避免引入随机噪声,从而提高了空域相关法的滤波效果。最终得到高信噪比拍波信号进而获得更高的比相测量找零精度,即无导轨大尺寸测量中至关重要的节点判定精度。

Claims (3)

1.He-Ne激光器拍波噪声滤除装置,其特征在于该装置由测量信号和参考信号分别通过高频放大器(1)连接功分器(2),功分器(2)连接混频器(4),其中一路功分器与混频器间接有延迟线(3),得到两路有固定相位差的低频信号,利用此两路低频信号中拍波演变出的信号高度相关,而所含噪声不相关的特点,将其依次输入基于数字信号处理器DSP(7)硬件的变步长自适应滤波器(5)和改进型小波滤波器(6),最终得到大大滤除了激光器噪声的高信噪比拍波信号。
2.按照权利要求1所述He-Ne激光器拍波噪声滤除装置的噪声滤除方法,其特征在于变步长自适应滤波器通过数字信号处理器DSP硬件高速运算,其迭代公式为:
Y(n)=X(n)W(n)
e(n)=d(n)-Y(n)
W(n+1)=W(n)+2μ(n)e(n)X(n)
μ ( n + 1 ) = α μ ( n ) + γ e ( n ) 2
式中Y(n)为实际输出,X(n)为输入信号,W(n)为权系数,X(n)=[x(n),x(n-1),x(n-2),...x(n-L+1)]T表示时刻n的输入信号矢量,T为矢量标志,由最近L个信号采样值构成,W(n)=[ω0(n),ω1(n),ω2(n),...ωL-1(n)]T是时刻n自适应滤波器的权系数矢量,L是自适应滤波器的阶数,d(n)是期望输出值,e(n)是误差,μ(n)是控制稳定性和收敛速度的参量,称之为步长因子,变步长因子μ(n)迭代公式中,0<α<1,γ>0,使变步长自适应滤波算法在初始收敛阶段或未知系统参数发生变化时,步长比较大,有较快的收敛速度或对时变系统的跟踪速度;而在算法收敛后,不管主输入端干扰信号有多大,都保持很小的调整步长以达到很小的稳态失调噪声,这种算法能同时满足收敛速率、跟踪速率及权失调噪声之间相互矛盾的要求。
3.按照权利要求1或2所述He-Ne激光器拍波噪声滤除装置的噪声滤除方法,其特征在于改进型小波滤波器采用自适应调整的方法,在小波变换值的相关量里加入自适应微调系数c,由下式自适应寻找c的最佳估计:
c→c+μlnλ
在此式中,μ是调整步长的一个正的常数,使在迭代收敛时剔除的小波系数的方差和本尺度上的噪声方差之比恰好等于1,避免引入随机噪声,从而提高了空域相关法的滤波效果。
CN 200410081242 2004-11-16 2004-11-16 He-Ne激光器拍波噪声滤除装置及其噪声滤除方法 Pending CN1606235A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200410081242 CN1606235A (zh) 2004-11-16 2004-11-16 He-Ne激光器拍波噪声滤除装置及其噪声滤除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200410081242 CN1606235A (zh) 2004-11-16 2004-11-16 He-Ne激光器拍波噪声滤除装置及其噪声滤除方法

Publications (1)

Publication Number Publication Date
CN1606235A true CN1606235A (zh) 2005-04-13

Family

ID=34765650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200410081242 Pending CN1606235A (zh) 2004-11-16 2004-11-16 He-Ne激光器拍波噪声滤除装置及其噪声滤除方法

Country Status (1)

Country Link
CN (1) CN1606235A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101968658A (zh) * 2010-09-27 2011-02-09 中国计量学院 基于外差式激光干涉仪的纳米静态锁定定位方法
CN102832908A (zh) * 2012-09-20 2012-12-19 西安科技大学 基于小波变换与变步长lms自适应滤波的信号降噪方法
CN103427789A (zh) * 2013-07-23 2013-12-04 四川大学 一种基于分数阶计算方程的图书馆图文信息去噪滤波器
CN105929238A (zh) * 2016-04-13 2016-09-07 西安理工大学 一种气固两相流交流电信号最小均方误差自适应滤波方法
CN106301289A (zh) * 2016-08-03 2017-01-04 广东工业大学 利用自适应滤波算法消除泥浆脉冲信号中的泵冲噪声的方法
CN106559722A (zh) * 2015-09-25 2017-04-05 中国科学院声学研究所 基于人耳听觉特性的音频重放系统均衡方法
CN111219572A (zh) * 2019-11-29 2020-06-02 华中科技大学 一种混频器支架
CN111638494A (zh) * 2020-05-31 2020-09-08 西南电子技术研究所(中国电子科技集团公司第十研究所) 数字波束合成系统多通道幅相校正方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101968658A (zh) * 2010-09-27 2011-02-09 中国计量学院 基于外差式激光干涉仪的纳米静态锁定定位方法
CN101968658B (zh) * 2010-09-27 2012-07-11 中国计量学院 基于外差式激光干涉仪的纳米静态锁定定位方法
CN102832908A (zh) * 2012-09-20 2012-12-19 西安科技大学 基于小波变换与变步长lms自适应滤波的信号降噪方法
CN102832908B (zh) * 2012-09-20 2015-06-17 西安科技大学 基于小波变换与变步长lms自适应滤波的信号降噪方法
CN103427789A (zh) * 2013-07-23 2013-12-04 四川大学 一种基于分数阶计算方程的图书馆图文信息去噪滤波器
CN103427789B (zh) * 2013-07-23 2016-06-08 四川大学 一种基于分数阶计算方程的图书馆图文信息去噪滤波器
CN106559722B (zh) * 2015-09-25 2018-12-04 中国科学院声学研究所 基于人耳听觉特性的音频重放系统均衡方法
CN106559722A (zh) * 2015-09-25 2017-04-05 中国科学院声学研究所 基于人耳听觉特性的音频重放系统均衡方法
CN105929238A (zh) * 2016-04-13 2016-09-07 西安理工大学 一种气固两相流交流电信号最小均方误差自适应滤波方法
CN106301289A (zh) * 2016-08-03 2017-01-04 广东工业大学 利用自适应滤波算法消除泥浆脉冲信号中的泵冲噪声的方法
CN111219572A (zh) * 2019-11-29 2020-06-02 华中科技大学 一种混频器支架
CN111219572B (zh) * 2019-11-29 2021-01-05 华中科技大学 一种混频器支架
CN111638494A (zh) * 2020-05-31 2020-09-08 西南电子技术研究所(中国电子科技集团公司第十研究所) 数字波束合成系统多通道幅相校正方法

Similar Documents

Publication Publication Date Title
Yegnanarayana et al. Significance of group delay functions in spectrum estimation
US8093911B2 (en) Time-of-flight measurement based on transfer function and simulated exponential stimulus
CN1606235A (zh) He-Ne激光器拍波噪声滤除装置及其噪声滤除方法
EP1540844B1 (en) Methods for transmitting a waveform having a controllable attenuation and propagation velocity
Duru et al. On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation
CN106383107A (zh) 一种基于平滑滤波联合vs‑lms的光谱信号降噪方法
CN108445539A (zh) 一种消除地震子波旁瓣干扰的方法、设备及系统
CN102289715A (zh) 基于前向线性预测的自适应小波神经网络去噪建模方法
Sun et al. Wavelet transform based de-noising method for self mixing interferometry signals
CN110808929A (zh) 相减策略的实复转换式信噪比估计算法
Ghogho et al. Fast computation of the exact FIM for deterministic signals in colored noise
Yurasova et al. Dynamic measurement errors correction adaptive to noises of a sensor
Jiang et al. Temporal acceleration of time‐domain integral‐equation solvers for electromagnetic scattering from objects residing in lossy media
Sims et al. Localization in One Dimensional Random Media:¶ A Scattering Theoretic Approach
Vandersteen et al. System identification for data acquisition characterization
CN112883787A (zh) 一种基于频谱匹配的短样本低频正弦信号参数估计方法
Biswal et al. Time frequency analysis and FPGA implementation of modified S-transform for de-noising
Medvedev et al. Fast computation of the direct scattering transform by fourth order conservative multi-exponential scheme
US20230393184A1 (en) Device and methods for phase noise measurement
Kleene et al. Concept of a complex envelope faber polynomial approach for the solution of Maxwell’s equations
Chen et al. A Denoising Method for Scanning Lidar Mie Scattering Echo Signal Based on MSW
Wang et al. Z-ADALINE based high-precision wide-frequency signal measurement algorithm for power electronic power grid
CN113721201B (zh) 一种线性调频信号调频率的估计方法
Zhang et al. The design and calibration of a low-cost underwater sound velocity profiler
Sun et al. FPGA based filter design for self-mixing interferometry signals

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication