徐先峰 張華竹 段晨東

（長安大學(xué)電子與控制工程學(xué)院 西安 710064）

1 引言

實(shí)現(xiàn)大型土木工程結(jié)構(gòu)損傷特征的精確提取，是解決土木工程應(yīng)用問題的關(guān)鍵點(diǎn)之一。同時，特征提取也是信號處理研究課題中的重要步驟［1～5］。獨(dú) 立 量 分 析［6］（Independent Component Analysis，ICA）是一類非常重要的特征提取方法，它假設(shè)各個源信號之間具有統(tǒng)計獨(dú)立性，這為后來的許多算法奠定了基礎(chǔ)［7～10］。但在土木工程應(yīng)用問題中，由于實(shí)際工程結(jié)構(gòu)具有較高的阻尼，導(dǎo)致ICA算法的穩(wěn)健性降低，很大程度上限制了算法的分離精度，使得模態(tài)參數(shù)不能有效地識別，為結(jié)構(gòu)損傷特征提取帶來了較大困難［11～13］。

針對傳統(tǒng)ICA算法在結(jié)構(gòu)損傷特征提取中的不足，本文引入RobustICA算法［14］，并將該算法應(yīng)用于實(shí)測IASC-ASCE結(jié)構(gòu)損傷數(shù)據(jù)的特征頻率提取。實(shí)驗(yàn)結(jié)果顯示RobustICA算法分離出的各分量之間獨(dú)立性較傳統(tǒng)FastICA算法更高，并在對各分量的固有頻率進(jìn)行識別后發(fā)現(xiàn)，RobustICA算法兼有較高的特征提取精度。

2 RobustICA算法

2.1 信號模型

考慮有n個未知的源信號，經(jīng)過線性混合之后可以得出如下信號模型：

其中，向量 x（t）=[x1（t），x2（t），…，xm（t）]T表示 m 個觀測信號組成的觀測向量；s（t）=[s1（t），s2（t），…， sn（t）]T是由n個未知的源信號所組成的向量；矩陣H表示未知的混合矩陣（m×n），若假設(shè)矩陣H行滿秩，那么其中每一個元素hij（hij=[H]ij）表示第 j個源信號在第i個觀測信號上的混合系數(shù)。符號（?）T表示轉(zhuǎn)置運(yùn)算符。在下文中，為了方便表示，省略時間t，將式（1）記為

盲源分離所要解決的問題就是從已知的觀測信號x中估計出源信號s。但是，由于混合矩陣H未知，所以需要引入一個分離矩陣W（n×m），使得式（2）轉(zhuǎn)換為

通過估計分離矩陣W ，使得輸出的y=Wx包含對源信號的估計。則矩陣W中的每一行都相當(dāng)于用來提取單個源的一個空間濾波器。

2.2 算法理論

在求解分離矩陣W時，需要選擇一個合適的目標(biāo)函數(shù)作為非高斯性的度量，使得W作用于觀測信號x后能將分量的非高斯性最大化。在ICA算法中，通常使用到的一類目標(biāo)函數(shù)為峭度表達(dá)式。零均值隨機(jī)變量y的峭度一般定義式為

在基于峭度的FastICA算法中，通常采用式（4）的簡化形式：

相對于一般定義式（4），簡化定義式（5）雖然計算簡便，但會引入計算誤差；不僅如此，F(xiàn)astICA算法在尋求最優(yōu)解時采用的是基于固定點(diǎn)的迭代法，容易陷入局部收斂，導(dǎo)致求得的解非全局最優(yōu)解。

而Zarzoso等提出的RobustICA算法［14］則是將峭度的一般定義式（4）的絕對值形式作為目標(biāo)函數(shù)，使用精確線性搜索方法來獲得最優(yōu)步長μ。假設(shè)含噪觀測矩陣為x，由于RobustICA算法不需要對其進(jìn)行白化處理，僅要求其均值為零，則分離后的一個輸出 y=wi?x（i=1，2，…，m）的峭度表達(dá)式為

式（6）同樣為峭度函數(shù)的一般定義式，是將y=wi?x（i=1，2，…，m）代入式（4）后的結(jié)果，不影響其一般性。如果對源信號非高斯性的假設(shè)不變，那么則可以將式（6）的絕對值形式作為目標(biāo)函數(shù)，得到最優(yōu)步長的精確線性搜索方法：

其中g(shù)為搜索方向，通常取梯度。至此，峭度準(zhǔn)則便可以表示成如下關(guān)于步長μ的四次多項(xiàng)式

再對 wi?（k）進(jìn)行歸一化處理：

最后求得對源信號的估計：

如上步驟為提取第 i個（i=1，2，…，m）分離矢量wi?。反復(fù)上述步驟即可獲得完整的分離矩陣W。

本文所做的主要工作如圖1所示。

3 特征分離方法

為了評價算法性能，將平均相似度rˉ作為相關(guān)性的評價指標(biāo)，先計算各分量頻域之間的相關(guān)系數(shù)：對于隨機(jī)向量bi，bj，其相關(guān)系數(shù)可以表示為ri，j：

其中cov（?X，Y）表示變量 X和變量Y的協(xié)方差。ri，j評價相關(guān)性可分為以下三種相關(guān)等級：

圖1 本文所做的主要工作

1）當(dāng) ri，j∈［0.8，1.0］時，表明分量 i與 j具有強(qiáng)相關(guān)性；

2）當(dāng) ri，j∈［0.4，0.8］時，表明分量 i與 j具有一般相關(guān)性；

3）當(dāng) ri，j∈［0，0.4］時，表明分量 i與 j具有弱相關(guān)性或無相關(guān)性。

計算出ri，j后，再利用平均相似度 rˉ對分量的整體相關(guān)性予以衡量：

rˉ越小，說明各分量總體呈現(xiàn)較好的獨(dú)立性；反之，則說明各分量之間相關(guān)性較大，沒有實(shí)現(xiàn)完全分離。

3.1 數(shù)值實(shí)驗(yàn)

實(shí)驗(yàn)采用ASCE的SHM Benchmark模型，如圖2（b）所示共設(shè)置16個節(jié)點(diǎn)，對應(yīng)設(shè)置16個傳感器分別安裝在每一層的2，6，8，4的柱子上，模擬6種損傷模式［15］。

該模型具有模擬對稱性，且激勵模擬為作用于每層y方向上的環(huán)境風(fēng)載，因此，在x，y方向上各選取4個作為檢測節(jié)點(diǎn)，分別為11、13、20、22、29、31、38、40。選取12自由度對稱模型下第一層所有斜支撐斷裂工況的數(shù)據(jù)進(jìn)行分析，數(shù)據(jù)長度為16384，采樣頻率為1000Hz。

由于源信號與混合矩陣均未知，且僅有傳感器接收到的信號。所以，該實(shí)際問題便可以看作是盲源分離模型 x（t）=Hs（t）。8個傳感器接采集到的信號 便 可 看 作 8 路 觀 測 信 號 x（t）=[x1（t），x2（t），…，x8（t）]T，記為（m為觀測信號個數(shù)）。

圖2 結(jié)構(gòu)模型的部件標(biāo)注與檢測節(jié)點(diǎn)分布

圖3和圖4分別是利用FastICA算法和Robust-ICA算法對8路觀測信號的盲源分離結(jié)果的頻譜。

圖3 FastICA分離后的各分量頻譜

圖4 RobustICA分離后的各分量頻譜

由圖3、圖4可以看出，觀測信號經(jīng)FastICA算法和RobustICA算法分離后均得到四階固有頻率f1、f2、f3、f4。表1為兩種算法分離出的12自由度對稱模型下第一層所有斜支撐斷裂工況下的模態(tài)固有頻率與文獻(xiàn)［15］計算出的該工況下模態(tài)固有頻率的對比。

表1 FastICA與RobustICA分離出的模態(tài)固有頻率與理論值的對比

表1顯示，F(xiàn)astICA與RobustICA提取出的特征頻率與理論值的均方根誤差分別為0.2230和0.0415，可見RobustICA算法在提取精度上優(yōu)于FastICA算法；并且，對比圖3與圖4可以看出，F(xiàn)astICA算法分離出的特征頻率出現(xiàn)混疊現(xiàn)象，而RobustICA分離出的各分量具有更高的獨(dú)立性。為了更有效地對這些分量進(jìn)行分析，現(xiàn)利用式（12）對分量兩兩之間的相關(guān)系數(shù)進(jìn)行計算，結(jié)果如表2。

表2 FastICA和RobustICA分離出的各分量之間相關(guān)系數(shù)

表2顯示，混合信號經(jīng)FastICA算法分離后各分量之間相關(guān)系數(shù) ri，j全部集中在［0.5，0.97］范圍內(nèi)，且通過式（13）可以計算出分量的平均相似度rˉ1=0.776∈［0.4，0.8］，說明FastICA算法分離出的各分量總體上具有較強(qiáng)的相關(guān)性，出現(xiàn)了嚴(yán)重的頻率混疊，沒有實(shí)現(xiàn)對源信號的通縮分離。相反，混合信號經(jīng)RobustICA算法分離后，各分量相關(guān)系數(shù)的最大值為0.96＞0.8，認(rèn)為分量b2和b3具有強(qiáng)相關(guān)性，屬于同一分量；次大值為0.50∈［0.4，0.8］，則分量b4和b6具有一般相關(guān)性，近似為同一分量；其余各分量之間的相關(guān)系數(shù)均小于0.3，可以認(rèn)為不具有相關(guān)性，b5、b7、b8為噪聲信號；且平均相似度 rˉ2=0.125＜ rˉ1=0.776。由此可以看出 RobustICA算法可以大大降低各分量的相關(guān)性，將主要的特征頻率單獨(dú)分離出來，近似實(shí)現(xiàn)通縮分離，分離精度顯著提高。

3.2 工程實(shí)驗(yàn)

利用RobustICA算法對Phase II IASC-ASCE實(shí)測結(jié)構(gòu)數(shù)據(jù)進(jìn)行損傷特征頻率提取。Phase II實(shí)測結(jié)構(gòu)共有16個數(shù)據(jù)采集通道，通道16為損傷結(jié)構(gòu)的激勵（隨機(jī)環(huán)境、力錘、振蕩器等）信號，通道1～15為加速度傳感器采集的損傷結(jié)構(gòu)響應(yīng)信號。傳感器布置及檢測節(jié)點(diǎn)如圖5所示。

圖5 Phase II鋼結(jié)構(gòu)實(shí)物模型及檢測節(jié)點(diǎn)分布

將環(huán)境噪聲作為激勵源，選取無損傷和去掉東側(cè)所有斜支撐兩種工況。采集每層中心節(jié)點(diǎn)2、5、8、11、14和激勵源節(jié)點(diǎn)16的響應(yīng)信號，數(shù)據(jù)長度為4096，采樣頻率為200Hz，利用RobustICA算法對這6個通道的信號進(jìn)行盲源分離。

1）工況一：無損傷

圖6為環(huán)境噪聲激勵下工況一的頻譜。

圖6 環(huán)境噪聲激勵下工況一的頻譜

將環(huán)境激勵下的工況一記為config1.ambient。表3為環(huán)境激勵下RobustICA分離出工況一的模態(tài)固有頻率與文獻(xiàn)［16］計算出的模態(tài)固有頻率的對比。

表3 環(huán)境激勵下工況一的模態(tài)固有頻率

表3顯示，無損傷工況下，RobustICA算法分離出四階模態(tài)頻率，與文獻(xiàn)［16］采用的貝葉斯兩步計算法計算出的四階模態(tài)頻率一致，具有較好的分離效果。

2）工況二：去掉東側(cè)所有斜支撐

圖7為環(huán)境噪聲激勵下工況二的頻譜。

圖7 環(huán)境噪聲激勵下工況二的頻譜

將環(huán)境激勵下的工況二記為config2.ambient。表4為環(huán)境激勵下RobustICA分離出工況二的模態(tài)固有頻率與文獻(xiàn)［16］計算出的模態(tài)固有頻率的對比。

表4 環(huán)境激勵下工況二的模態(tài)固有頻率

表4顯示，去掉東側(cè)所有斜支撐工況下，RobustICA算法分離出三階模態(tài)頻率，與文獻(xiàn)［16］采用的貝葉斯兩步計算法計算出的三階模態(tài)頻率一致，具有較好的分離效果。

4 結(jié)語

為了精確提取土木工程中結(jié)構(gòu)損傷特征頻率，將一種穩(wěn)健的盲源分離算法RobustICA應(yīng)用于實(shí)際工程中，對采集到的ASCE實(shí)測結(jié)構(gòu)的混合信號進(jìn)行分離，提取損傷模態(tài)的固有頻率，具有較好的分離精度，并能實(shí)現(xiàn)對模態(tài)頻率的通縮分離，解決了頻率混疊的問題。