張 馳，高向東，蕭振林，陳曉輝

（1.廣東工業(yè)大學(xué)機(jī)電工程學(xué)院，廣東廣州510006；2.廣州番禺高勛染整設(shè)備制造有限公司，廣東番禺511400）

色噪聲下卡爾曼濾波磁光成像焊縫跟蹤算法

張 馳1，高向東1，蕭振林2，陳曉輝2

（1.廣東工業(yè)大學(xué)機(jī)電工程學(xué)院，廣東廣州510006；2.廣州番禺高勛染整設(shè)備制造有限公司，廣東番禺511400）

焊縫跟蹤是保證激光焊接質(zhì)量的前提，在激光焊接過程中，精確地識別和跟蹤焊縫非常重要。以緊密對接平板激光焊為試驗對象，研究一種色噪聲環(huán)境下應(yīng)用卡爾曼濾波最優(yōu)狀態(tài)估計對焊縫偏差進(jìn)行預(yù)測的方法。使用磁光傳感器攝取焊縫磁光圖像并提取焊縫位置參數(shù)，建立基于焊縫位置參數(shù)的系統(tǒng)測量方程和狀態(tài)方程。應(yīng)用擴(kuò)展?fàn)顟B(tài)向量的方法建立色噪聲環(huán)境下的卡爾曼濾波算法，對焊縫位置進(jìn)行最優(yōu)狀態(tài)估計，得到均方差極小值情況下的焊縫偏差預(yù)測最優(yōu)值，減小系統(tǒng)噪聲對焊縫位置識別的影響。試驗結(jié)果表明所提方法可以顯著提高焊縫跟蹤的準(zhǔn)確性。

焊縫跟蹤；磁光成像；色噪聲；卡爾曼濾波

0 前言

激光焊接是一種先進(jìn)的焊接工藝技術(shù)，在現(xiàn)代加工制造中發(fā)揮著越來越重要的作用。激光焊接能量密度極高，熱影響區(qū)域小，因此其對焊縫間隙有著嚴(yán)格的要求，為了保證激光焊接的質(zhì)量，要求焊縫間隙必須遠(yuǎn)小于光斑直徑。由于焊縫間隙小，使得激光束準(zhǔn)確對中焊縫成為實現(xiàn)高質(zhì)量焊接的重要因素，因此，焊接過程的精確跟蹤是激光焊接的關(guān)鍵所在[1]。目前，焊縫跟蹤方法主要有結(jié)構(gòu)光視覺傳感法、焊縫區(qū)直接圖像傳感法和激光頻閃攝像法。這三種方法對于焊縫寬度小于0.1 mm的微間隙焊縫都難以實現(xiàn)精確有效地自動檢測和跟蹤[2-3]，為此研究了一種基于磁光成像傳感的焊縫識別新方法[4]。

在激光焊接過程中，被焊材料的溫度、磁光傳感器與焊件的間距都會產(chǎn)生波動，并影響到焊縫磁光成像。因此焊縫磁光圖像必然包含很多噪聲干擾?？紤]到各種噪聲的干擾，探索一種卡爾曼濾波預(yù)測焊縫偏差的新方法。從磁光傳感器攝取的磁光圖像中提取焊縫位置參數(shù)并構(gòu)成狀態(tài)向量，建立系統(tǒng)狀態(tài)方程和測量方程。在此數(shù)學(xué)模型下應(yīng)用卡爾曼濾波對焊縫位置參數(shù)進(jìn)行最優(yōu)估計，減小系統(tǒng)噪聲與測量噪聲對焊縫位置檢測的影響，得到最小均方差條件下的焊縫偏差最優(yōu)預(yù)測值，為焊縫的糾偏提供實時有效的控制參量，提高焊縫跟蹤精度[5]。

1 磁光檢測微間隙焊縫試驗

根據(jù)法拉第旋光效應(yīng)，一束線偏振光通過旋光性介質(zhì)時，如果在介質(zhì)中沿光傳播方向施加外磁場B，則光通過介質(zhì)后，光的振動面轉(zhuǎn)過一定角度θ。轉(zhuǎn)角θ在材料確定的情況下，主要與光的波長、外磁場強(qiáng)度、光在介質(zhì)中傳播的路程有關(guān)。通過法拉第旋光效應(yīng)，構(gòu)建焊縫跟蹤試驗平臺，如圖1所示。試驗平臺包括YAG激光焊接系統(tǒng)、磁光傳感器、保護(hù)氣（氬）、工作臺、伺服控制電機(jī)和相應(yīng)夾具等。

由于磁場激勵器的影響，磁場在焊縫處會發(fā)生變化，磁光傳感器能檢測出這種磁場的變化，形成的磁光圖像則包含了焊縫位置信息。以低碳鋼平板緊密對接為試驗研究對象，焊縫間隙小于0.1 mm，實物如圖2a所示，圖2b為其對應(yīng)的磁光圖像。

圖1 焊縫跟蹤試驗裝置示意

圖2 焊縫實物圖與磁光圖像

2 焊縫位置參數(shù)和焊縫偏差

試驗材料選用兩塊厚1.5 mm的緊密對接的低碳鋼平板，焊件運動由工作臺的精密伺服電動機(jī)驅(qū)動，試驗條件如表1所示。為了充分反映實際生產(chǎn)環(huán)境中激光焊接的跟蹤情況，試驗采用偏斜焊縫路徑，激光頭行走路徑橫跨焊縫，依次經(jīng)歷左偏、對中、右偏三種情況，焊接路徑如圖3所示。磁光傳感器剛性連接在激光頭上，隨激光頭同步運動實時獲取焊縫信息。激光頭與磁光傳感器相對位置固定，因此根據(jù)磁光傳感器探測的焊縫偏差即可推算出激光束與焊縫之間的偏差信息。

表1 激光焊接磁光成像試驗條件

分析一組焊接動態(tài)序列磁光圖像，通過灰度拉伸、中值濾波、二值化、形態(tài)學(xué)變換和邊緣提取，得到焊縫中心，如圖2c所示。計算焊縫偏差，得到焊縫位置參數(shù)測量值，如圖4所示。焊接過程中焊縫偏差測量值與實際值較為接近，但由于系統(tǒng)噪聲的干擾，測量值有一定波動。

圖3 焊縫跟蹤路徑示意

圖4 焊縫偏差實際值與測量值

3 系統(tǒng)狀態(tài)方程與測量方程

建立精確的系統(tǒng)狀態(tài)轉(zhuǎn)移數(shù)學(xué)模型和測量模型是卡爾曼濾波預(yù)測焊縫位置的前提。焊接過程中，焊件相對于激光束在x-y平面上運動，定義焊縫位置在k時刻時為（x（k），y（k）），k+1時刻為（x（k+ 1），y（k+1）），k時刻，焊件x方向與激光束的相對運動速度為方向與激光束的相對運動速度為采樣時間為t，激光束在x和y兩個方向的加速度分別表示為ωx（k），ωy（k），則描述焊縫位置的方程可表述為[6]

由于激光束與焊件在x方向與y方向相互獨立地運動，而焊縫偏差主要產(chǎn)生在x方向，因此只考慮x方向即可。以焊縫當(dāng)前位置參數(shù)、當(dāng)前速度和前一時刻位置參數(shù)構(gòu)造狀態(tài)向量，即[6]

得到系統(tǒng)狀態(tài)方程為[6]

1.報社問題頻發(fā)，競爭激烈?！拔r代”環(huán)境下有的傳統(tǒng)媒體的生存狀況不容樂觀，同時，傳統(tǒng)媒體之間的競爭更加激烈，形成“適者生存，不適者被淘汰”的局面。據(jù)2014年4月23日《人民日報》報道，因虛假新聞問題，200余家報刊被查處，這種嚴(yán)峻情況反映了傳統(tǒng)媒體已經(jīng)出現(xiàn)了生存危機(jī)，紙質(zhì)媒體的銷路受阻，人們對紙媒態(tài)度的淡漠等問題凸現(xiàn)出來，一些新聞記者或編輯“鋌而走險”，結(jié)果導(dǎo)致報社受到嚴(yán)重懲罰。

4 色噪聲卡爾曼濾波模型

實際焊接過程中系統(tǒng)動態(tài)噪聲ωx（k）并非為白噪聲，而是焊接系統(tǒng)所產(chǎn)生的相關(guān)色噪聲。而在傳統(tǒng)Kalman濾波的推導(dǎo)中，系統(tǒng)噪聲和測量噪聲必須是相互獨立的高斯白噪聲，色噪聲的存在導(dǎo)致卡爾曼濾波的精度降低，因此在預(yù)測焊縫位置時無法直接用傳統(tǒng)卡爾曼濾波算法進(jìn)行處理，需要將有色噪聲歸入系統(tǒng)狀態(tài)向量，實現(xiàn)色噪聲的白化從而滿足卡爾曼濾波對噪聲的要求。

設(shè)系統(tǒng)動態(tài)噪聲服從自回歸模型，焊接過程色噪聲采用一階自回歸模型表示，對實際焊接來說能夠滿足要求[7]。設(shè)

式中e（k）為高斯白噪聲。

系統(tǒng)狀態(tài)方程及焊縫位置測量方程為

即

通過以上系統(tǒng)狀態(tài)向量維數(shù)擴(kuò)展處理將系統(tǒng)色噪聲白化，從而達(dá)到傳統(tǒng)卡爾曼濾波對噪聲的要求。應(yīng)用卡爾曼濾波對狀態(tài)方程進(jìn)行處理，其流程如圖5所示。

5 試驗結(jié)果與分析

基于卡爾曼濾波的焊縫跟蹤，是在觀測序列Z（k）已知的條件下，使得估計誤差達(dá)到最小方差的條件下，預(yù)測焊縫偏差。首先必須確定濾波的初始狀態(tài)，以起始相鄰兩幀圖像中的焊縫位置相關(guān)參數(shù)作為濾波的狀態(tài)初始值。在噪聲的影響下，相鄰時刻焊縫位置參數(shù)偏移量在0.02 mm內(nèi)隨機(jī)波動，取系統(tǒng)動態(tài)噪聲初值ωx（l）=0.01，得到濾波狀態(tài)的初始值為

圖5 卡爾曼濾波算法流程

計算k=2時刻的狀態(tài)一步預(yù)測值

k=2時刻的測量值為

計算測量的預(yù)測值為

可得到狀態(tài)向量最佳濾波值

完成一次濾波后，依次對每個時刻的焊縫位置參數(shù)進(jìn)行濾波處理，得到整個焊接過程的焊縫偏差，如圖6所示，可以看出，濾波后的焊縫偏差值比測量值更接近偏差的實際值。

圖6 焊縫偏差實際值、測量值與濾波值

為更好地比較濾波后的焊縫偏差與測量值的精確度，定義焊縫誤差其中α為焊縫偏差測量值或濾波值，β為焊縫偏差實際值。由此可知，焊縫誤差e越小，焊縫偏差越接近實際值，跟蹤精度越高。對試驗所得的數(shù)據(jù)進(jìn)行焊縫誤差的計算，得到曲線如圖7所示。濾波后的焊縫誤差與測量誤差對比如表2所示，由圖7和表2可知濾波后焊縫偏差更加接近實際值，而且濾波后焊縫偏差的波動更小和平滑。試驗驗證了過程色噪聲下卡爾曼濾波可以抑制系統(tǒng)噪聲的影響，得到更精確的焊縫偏差參數(shù)。

圖7 焊縫測量誤差和濾波誤差

表2 焊縫測量誤差與濾波誤差對比

6 結(jié)論

針對緊密對接的微間隙焊縫，磁光傳感器可以攝取到包含焊縫信息的磁光圖像。由磁光圖像中提取出焊縫的位置信息并構(gòu)成狀態(tài)向量，構(gòu)造系統(tǒng)狀態(tài)方程與測量方程，在過程色噪聲環(huán)境下應(yīng)用卡爾曼濾波實現(xiàn)對焊縫位置狀態(tài)的最優(yōu)估計，減小系統(tǒng)噪聲的干擾，為實時糾偏提供有效可靠的焊縫偏差信息，提高焊縫跟蹤精度。

[1]徐國建，王虹，唱麗麗，等.激光焊接的特性[J].電焊機(jī)，2010，40（11）：57-64.

[2]FANG Z J，XU D，TAN M.Visual seam tracking system for butt weld of thin plate[J].International Journal of Advanced Manufacturing Technology，2010（49）：519-526.

[3]房加強(qiáng)，何建萍，王付鑫，等.視覺傳感技術(shù)的應(yīng)用研究與發(fā)展趨勢[J].電焊機(jī)，2013，43（4）：46-50.

[4]GAO X D，LiuYH，YOUDY.Detection of micro-weld joint by magneto-optical imaging[J].Optics&Laser Technology，2014（62）：141-151.

[5]GAO X D，YOU D Y，SEIJI KATAYAMA.Seam Tracking MonitoringBasedonAdaptiveKalmanfilterEmbeddedElman Neural Network During High-power Fiber Laser Welding[J]. IEEE TRANSACTIONSONINDUSTRIALELECTRONICS，2012，59（11）：4315-4325.

[6]劉益，高向東，張弛.卡爾曼濾波磁光成像計盒維數(shù)焊縫跟蹤算法[J].電焊機(jī)，2015，45（1）：132-136.

[7]高向東，仲訓(xùn)杲，游德勇，KATAYAMA Seiji.色噪聲下卡爾曼濾波焊縫跟蹤算法與試驗研究[J].控制理論與應(yīng)用，2011，28（7）：931-935.

Kalman filtering algorithm with colored noise for seam tracking based on magneto-optical imaging

ZHANG Chi1，GAO Xiangdong1，XIAO Zhenlin2，CHEN Xiaohui2
（1.School of Electromechanical Engineering，Guangdong University of Technology，Guangzhou 510006，China；2. Guangzhou Panyu Gofront Dyeing&Finishing Machinery Manufacturer Ltd.，Panyu 511400，China）

Seam tracking is prerequisite for laser welding quality.It is very important to accurately detect weld position and track the weld path in laser welding process.For micro butt-joint laser welding of two low carbon steel plates，an approach is investigated to estimate the seam offset by using the Kalman filtering optimal state estimation under the colored noise environment.A magnetooptical sensor based on the principle of Faraday magneto effect is applied to acquire the magneto-optical images of welds.The weld position parameters can be extracted by using image processing techniques.The state-equation based on the weld position parameters and the measurement-equation for the weld position is established.A Kalman filtering algorithm under colored noise environment is developed by augmenting the state vectors to obtain the optimal seam offset prediction under the condition of minimum mean square error through the optimal state estimation of weld position.This algorithm can reduce the system noises interference influence. Experiment results show that the accuracy of seam tracking can be improved significantly by Kalman filtering algorithm.

seam tracking；magneto-optical imaging；colored noise；kalman filtering

TG409

A

1001-2303（2015）08-0081-05

10.7512/j.issn.1001-2303.2015.08.17

2014-05-16；

2015-08-03

收稿日期：國家自然科學(xué)基金（51175095）；廣東省自然科學(xué)基金(10251009001000001)；廣東省學(xué)科建設(shè)科技創(chuàng)新項目（2013KJCX0063）；廣東省協(xié)同創(chuàng)新與平臺環(huán)境建設(shè)專項資金項目（503149102077）；廣東省重大科技專項（2014B090921008）；廣州市科學(xué)研究專項（1563000554）；佛山市科技創(chuàng)新專項資金項目(2014AG10015)

張馳(1990—)，男，安徽人，在讀碩士，主要從事焊接自動控制的研究工作。