張 檣,周西峰,郭前崗
(南京郵電大學(xué) 自動(dòng)化學(xué)院, 江蘇 南京 210046)
基于改進(jìn)的DE算法的超聲回波參數(shù)估計(jì)
張 檣,周西峰,郭前崗
(南京郵電大學(xué) 自動(dòng)化學(xué)院, 江蘇 南京 210046)
利用最小二乘法將超聲回波參數(shù)估計(jì)問(wèn)題轉(zhuǎn)化為優(yōu)化問(wèn)題,采用改進(jìn)的差分進(jìn)化算法(DE算法)對(duì)待優(yōu)化函數(shù)進(jìn)行優(yōu)化處理。仿真研究說(shuō)明,該方法在噪聲環(huán)境下依然有效,并且不依賴(lài)于初始值選擇,不需計(jì)算梯度,可以在全局范圍內(nèi)搜索。
DE算法;高斯牛頓法;超聲回波信號(hào);信號(hào)處理
超聲無(wú)損檢測(cè)一直是無(wú)損檢測(cè)的研究熱點(diǎn),超聲信號(hào)波形中包含有被檢測(cè)物體的幾何形狀、大小、缺陷故障等一系列有價(jià)值的信息。模型化的超聲信號(hào)分析方法能有效提取信號(hào)中有用信息,其中每個(gè)回波模型都是一個(gè)非線性函數(shù)包含有一系列的參數(shù):回波帶寬、中心頻率、相位、幅度、超聲回波的渡越時(shí)間。這些參數(shù)實(shí)時(shí)反映超聲信號(hào)的波形,并且與反射器的物理特性以及傳播路徑的頻率特性密切相關(guān)[1]。1981年,SANIIE J[2]提出了利用高斯回波模型模擬超聲信號(hào)的方法; 2001年,DEMIRLI R[3]利用高斯牛頓算法迭代出了高斯回波模型的參數(shù),具有十分高的精度。但是高斯牛頓法過(guò)度依賴(lài)于初值點(diǎn)的選取,若初值點(diǎn)選取偏差較大將會(huì)影響算法整體的收斂性,而且該算法從理論上來(lái)說(shuō)獲取的是局部最優(yōu)解,并非全局最優(yōu)解,所以這種方法在實(shí)際應(yīng)用中具有很大局限性。
針對(duì)以上問(wèn)題,本文基于超聲回波信號(hào)的高斯模型,提出改進(jìn)的差分進(jìn)化算法對(duì)超聲信號(hào)的參數(shù)進(jìn)行估計(jì)。
差分進(jìn)化算法是基于群體智能理論的優(yōu)化算法[4- 6]。按照算法分析指標(biāo)對(duì)DE算法分析如下:(1)在編碼標(biāo)準(zhǔn)方面,DE算法與粒子群優(yōu)化算法類(lèi)似,都是采用實(shí)數(shù)編碼;(2)在參數(shù)設(shè)置問(wèn)題上,研究表明DE算法的參數(shù)設(shè)置較少,對(duì)結(jié)果影響不是很明顯,優(yōu)于遺傳算法和粒子群優(yōu)化算法;(3)在高維問(wèn)題中,DE算法和粒子群優(yōu)化算法能夠很好地解決問(wèn)題,而且DE算法收斂快且精確;(4)收斂性能上,粒子群優(yōu)化算法容易陷入局部最優(yōu)解并且不穩(wěn)定,DE算法不存在這樣的問(wèn)題。差分進(jìn)化算法的基本原理是:采用對(duì)個(gè)體進(jìn)行方向擾動(dòng),以達(dá)到個(gè)體的函數(shù)值下降的目的。同其他進(jìn)化算法一樣,差分進(jìn)化算法不利用函數(shù)的梯度信息,因此對(duì)函數(shù)的可導(dǎo)性甚至連續(xù)性沒(méi)有要求,適用性很強(qiáng)。DE 算法的搜索性能取決于算法全局探索和局部開(kāi)發(fā)能力的平衡 ,而這在很大程度上依賴(lài)于算法的控制參數(shù)的選取 ,包括種群規(guī)模、縮放比例因子和交叉概率等[7]。但是實(shí)驗(yàn)表明,在噪聲環(huán)境下差分進(jìn)化算法的優(yōu)化功能遠(yuǎn)遠(yuǎn)低于進(jìn)化算法。本文所應(yīng)用的改進(jìn)差分進(jìn)化算法能有效處理噪聲環(huán)境下函數(shù)優(yōu)化。
超聲換能器的脈沖響應(yīng)可以模擬成高斯信號(hào)[8],響應(yīng)的脈沖幅度可以表示為:
x(t)=βe-αt2cos(2πfct)
(1)
高斯回波的一般模型為:
s(θ;t)=βe-α(t-τ)2cos[2πfc(t-τ)+φ]
(2)
參數(shù)向量θ=[βατfcφ]為超聲回波的參數(shù)向量。其中,α為帶寬,τ為到達(dá)時(shí)間,fc為中心頻率,φ為相位,β為幅度系數(shù)。
回波模型做加性白噪聲處理得到:
y(t)=s(θ;t)+v(t)
(3)
其中,v(t)為加性高斯白噪聲。
為了準(zhǔn)確估計(jì)超聲信號(hào)模型參數(shù),首先應(yīng)用最小二乘法構(gòu)造目標(biāo)待優(yōu)化函數(shù)[9]如下:
(4)
問(wèn)題就轉(zhuǎn)化為求J(θ)為最小值時(shí)的參數(shù)向量θ,J(θ)為目標(biāo)待優(yōu)化函數(shù),y為實(shí)際信號(hào),s(θ)為信號(hào)模型,θ為目標(biāo)待估測(cè)參數(shù)向量。
差分進(jìn)化算法的參數(shù)定義[10-11]如下:
NP:群體大小,即種群中個(gè)體數(shù)目,NP越大,種群多樣性越強(qiáng),活動(dòng)最優(yōu)解的概率越大,但是計(jì)算量也越大,一般選擇在5D~10D之間;
F:變異因子,經(jīng)典DE算法中它是一個(gè)實(shí)常數(shù)因數(shù),它決定偏差向量的放大比例,本文方法中將F定義如下:
圖1 超聲估計(jì)信號(hào)
F=0.5×(1+rand(0,1))
(5)
比列因子的均值保持在0.75,因子的隨機(jī)變化有利于搜索的多樣性,在真正的最優(yōu)解搜索到之前不易陷入停滯,在噪聲背景下,尤其適用。
CR:交叉因子,經(jīng)典DE算法中它是一個(gè)范圍在[0,1]的實(shí)數(shù),控制一個(gè)試驗(yàn)向量來(lái)自隨機(jī)選擇的變異向量而不是原來(lái)向量的概率的參數(shù)。
G:最大迭代次數(shù),G越大,獲得的最優(yōu)解越精確,計(jì)算量增加;
D:?jiǎn)栴}維數(shù);
Xi(t):群體中每個(gè)目標(biāo)個(gè)體;
Vi(t):變異的目標(biāo)個(gè)體;
Ui(t):交叉操作后的試驗(yàn)個(gè)體。
算法具體步驟如下:
(1)初始化種群,同時(shí)初始化算法中的各參數(shù);
(2)計(jì)算每個(gè)個(gè)體的適應(yīng)值,適應(yīng)度函數(shù)根據(jù)具體問(wèn)題決定;
(3)判斷是否滿(mǎn)足終止條件(終止條件是指找到最優(yōu)解或者達(dá)到最大迭代次數(shù)),若是則進(jìn)化終止得到最優(yōu)解,若不滿(mǎn)足終止條件則跳到步驟(4);
(4)執(zhí)行變異操作,將式(5)代入式(6)中:
Vi(t)=Xr1(t)+F×(Xr2(t)-Xr3(t))
(6)
產(chǎn)生變異個(gè)體Vi(t),其中i、r1、r2、r3是取值在1~NP之間的4個(gè)不同整數(shù),執(zhí)行步驟(5);
(5)執(zhí)行交叉操作,按照式(7)生成試驗(yàn)個(gè)體Ui(t);
(7)
(6)執(zhí)行選擇操作,采用的是貪婪選擇的策略,即選擇較優(yōu)的個(gè)體作為新的個(gè)體:
(8)
f(*)即目標(biāo)待優(yōu)化函數(shù)J(θ)。
在仿真試驗(yàn)中,超聲信號(hào)的模型參數(shù)向量取值為θ=[3.2 25 3.5 5 1.8]。超聲信號(hào)的采樣頻率為100 MHz,采樣信號(hào)保存為s(θ),分別在該信號(hào)中添加30 dB、20 dB、10 dB的高斯白噪聲作為實(shí)測(cè)信號(hào)。估計(jì)出的信號(hào)如圖1所示。
在差分進(jìn)化算法中設(shè)置種群NP=50,交叉因子CR=0.3,最大迭代次數(shù)為500,仿真結(jié)果如表1所示。
從表1可以看出,無(wú)噪聲情況下高斯法可以100%的估計(jì)得到原始信號(hào),本文方法有誤差,但是與真實(shí)值很接近,如圖1(a)所示。在有噪聲信號(hào)的影響下,高斯法根本無(wú)法迭代估計(jì)參數(shù),本文方法可以得到優(yōu)化解。但是隨著噪聲影響增大,估測(cè)得到參數(shù)的誤差也在加大,但是超聲信號(hào)特征參數(shù)的3個(gè)重要參數(shù)估測(cè)的精確度依然很高,誤差分析如表2所示。
表1 參數(shù)估計(jì)結(jié)果
表2 誤差分析(%)
本文提出利用最小二乘法構(gòu)造目標(biāo)待優(yōu)化函數(shù),并應(yīng)用DE算法的變異因子的特點(diǎn)提出改進(jìn)的DE算法估計(jì)噪聲情況下的參數(shù),該方法充分利用DE算法對(duì)初值依賴(lài)性不大、不需計(jì)算梯度、善于在全局范圍內(nèi)搜索等優(yōu)點(diǎn)。實(shí)驗(yàn)結(jié)果表明,與高斯牛頓法相比,該算法收斂速度很快而且結(jié)果很精確,具有很高的工程應(yīng)用性。
[1] DEMIRLI R. Model based estimation of ultrasonic echoes: analysis, algorithms, and applications[M]. Illinois Institute of Technology, 2001.
[2] SANIIE J. Ultrasonic signal processing: system identification and parameter estimation of reverberant and inhomogeneous targets[J]. Ultrasonic Signal Processing: System Identification
and Parameter EST, 1981.
[3] DEMIRLI R,SANIIE J.Model-based estimation of ultrasonic echoes. Part I: analysis and algorithms[J]. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 2001,48(3):787-802.
[4] Fan Huiyuan, LAMPINEN J. A trigonometric mutation operation to differential evolution[J]. Journal of Global Optimization, 2003, 27(1): 105-129.
[5] 胡中波.差分演化算法及其在函數(shù)優(yōu)化中的應(yīng)用研究[D]. 武漢:武漢理工大學(xué),2006.
[6] 李凱斌.智能進(jìn)化優(yōu)化算法的研究與應(yīng)用[D]. 杭州:浙江大學(xué),2008.
[7] KRINK T, FILIPIC B, FOGEL G B, et al. Noisy optimization problems-a particular challenge for differential evolution[C].Proceedings of the 2004 Congress on Evolutionary Computation,2004:332-339.
[8] 周方,張小鳳,張光斌.超聲回波參數(shù)的蟻群算法估計(jì)[J]. 陜西師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012,40(2):35-40.
[9] 周西峰, 朱文文, 郭前崗. 基于遺傳算法和高斯牛頓法的超聲回波信號(hào)參數(shù)估計(jì)[J]. 解放軍理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 13(3): 247-251.
[10] 汪艷, 張小鳳, 張光斌, 等. 基于高斯模型的多重超聲回波信號(hào)重?cái)?shù)估計(jì)[J]. 陜西師范大學(xué)學(xué)報(bào) (自然科學(xué)版), 2015,43(1):7.
[11] 潘欣,高曉智.基于差分演化的果蠅優(yōu)化算法[J].微型機(jī)與應(yīng)用,2015,34(1):23-25,28.
Parameters estimation of ultrasonic echo signal based on improved DE algorithm method
Zhang Qiang,Zhou Xifeng,Guo Qian′gang
(School of Automation, Nanjing University of Posts andTelecommuninations,Nanjing 210046,China)
This paper uses the least square method to make the ultrasonic echo parameter estimation problem converted into the optimization problem, and the improved differential evolution algorithm (DE algorithm) is used to optimize the parameters. Simulation results show that the proposed method is effective in the noise environment, does not depend on the initial value selection, without the need to calculate the gradient, and can be in the global range search.
DE algorithm; Gauss-Newton method; ultrasonic echo signal; signal processing
TN911.6
A
1674-7720(2016)04-0088-03
張檣,周西峰,郭前崗.基于改進(jìn)的DE算法的超聲回波參數(shù)估計(jì)[J].微型機(jī)與應(yīng)用,2016,35(4):88-90.
2015-10-31)
張檣(1992-),通信作者,女,碩士研究生,主要研究方向:智能檢測(cè)與控制技術(shù)。E-mail: 63876496@qq.com。
周西峰(1960-),男,本科,碩士生導(dǎo)師,主要研究方向:計(jì)算機(jī)智能監(jiān)測(cè)與控制技術(shù)等。