基于協(xié)方差矩陣重構(gòu)和導(dǎo)向矢量估計的穩(wěn)健自適應(yīng)波束形成

謝菊蘭　李昕亞　李會勇　王旭

(電子科技大學(xué),成都 611731)

?

基于協(xié)方差矩陣重構(gòu)和導(dǎo)向矢量估計的穩(wěn)健自適應(yīng)波束形成

謝菊蘭李昕亞李會勇王旭

(電子科技大學(xué),成都 611731)

摘要實(shí)際應(yīng)用中,當(dāng)假定的與真實(shí)的期望信號導(dǎo)向矢量之間存在一定誤差時,波束形成器的性能會急劇下降,特別是當(dāng)期望信號功率很強(qiáng)的時候.為解決這個問題,提出了一種新的算法.當(dāng)信源數(shù)小于陣元數(shù)時,干擾加噪聲協(xié)方差矩陣具有稀疏性.新方法首先利用該特性重構(gòu)干擾加噪聲協(xié)方差矩陣并由此得到與干擾導(dǎo)向矢量正交的子空間,使接收的數(shù)據(jù)通過該子空間得到只含有期望信號和噪聲的混合信號,然后，對該混合信號基于最大化輸出功率原理估計期望信號導(dǎo)向矢量,最后，把得到的導(dǎo)向矢量和正交子空間來構(gòu)造陣列加權(quán)值.仿真結(jié)果表明：該算法分別在假定的期望信號導(dǎo)向矢量存在誤差、期望信號很強(qiáng)和低快拍數(shù)時仍然具有良好的性能.

關(guān)鍵詞穩(wěn)健自適應(yīng)波束形成;導(dǎo)向矢量誤差;期望信號功率；矩陣重構(gòu)

引言

自適應(yīng)波束形成技術(shù)是陣列信號處理中的關(guān)鍵技術(shù),具有廣泛的應(yīng)用前景．傳統(tǒng)的波束形成算法往往是建立在我們已經(jīng)知道準(zhǔn)確的期望信號導(dǎo)向矢量這個條件下來抑制干擾的．然而在強(qiáng)期望信號背景下,由于采樣快拍中含有期望信號成分,使得這些算法的性能對系統(tǒng)中出現(xiàn)的失配量非常敏感[1-3],這些失配量可能來自信號傳播環(huán)境、陣元分布、快拍數(shù)有限等．任何一種失配都將對系統(tǒng)性能造成嚴(yán)重的影響．因此,近年來,穩(wěn)健的自適應(yīng)波束形成算法吸引了許多研究者的關(guān)注,一些穩(wěn)健的自適應(yīng)波束形成算法也被陸續(xù)提出[4-5]．

一類常見穩(wěn)健的自適應(yīng)波束形成算法是特征分解波束形成法(Eigen-Subspace Beamforming, ESB)[6],這種方法基于特征值分解利用了信號子空間的特性．該類算法的不足是當(dāng)信噪比較低的時候性能會受影響并且需已知信源個數(shù)．還有一些穩(wěn)健的自適應(yīng)波束形成算法通過給樣本協(xié)方差矩陣加載一個對角矩陣改善系統(tǒng)穩(wěn)健性,這類算法稱為對角加載算法[7-8],它們的缺點(diǎn)是在不同的背景下無法確定最優(yōu)的加載量．第三類穩(wěn)健波束形成算法是基于最差情況的最優(yōu)化問題準(zhǔn)則的算法[9-10],這種算法實(shí)際上與第二種算法是等價的[8]．此外大多數(shù)此類穩(wěn)健波束形成算法是基于二階錐規(guī)劃問題(Second Order Cone Programming, SOCP)的,求解需要使用一些特殊的最優(yōu)化工具箱[11],這類算法容易陷入局部尋優(yōu)的問題且不利于工程實(shí)現(xiàn),因此其實(shí)際應(yīng)用價值有限．最近,有一種算法吸引了很多研究者的關(guān)注[12-18],這種算法在已知期望信號來波方向的空域范圍內(nèi)積分,通過最小化輸出功率估計真實(shí)的期望信號導(dǎo)向矢量．然而大部分此類算法是基于二次約束二次規(guī)劃(Quadratically Constrained Quadratic Programming, QCQP)的,通過凸優(yōu)化工具箱比如CVX[19]求解．此類算法同樣面臨局部尋優(yōu)的問題,且在弱期望信號時往往期望信號導(dǎo)向矢量估計不準(zhǔn)確．以上提到的大部分穩(wěn)健的自適應(yīng)波束形成算法在強(qiáng)期望信號背景下性能都會嚴(yán)重衰退,不能對干擾進(jìn)行有效抑制．

本文提出了一種穩(wěn)健的自適應(yīng)波束形成算法,這種算法首先利用矩陣的稀疏特性重構(gòu)干擾加噪聲協(xié)方差矩陣并由此得到與干擾信號導(dǎo)向矢量正交的矩陣,利用該矩陣的正交性去除干擾信號之后,再利用旋轉(zhuǎn)子空間的原理通過最大化輸出功率估計導(dǎo)向矢量,將得到的導(dǎo)向矢量作為陣列加權(quán)值,最終使波束形成器在強(qiáng)期望信號背景下具有對抗期望信號導(dǎo)向矢量失配的穩(wěn)健性．該方法基本不受期望信號輸入信噪比(Signal-to-Noise Ratio, SNR)的影響,并且在期望信號來波方向失配量很大時也可以保持很好的性能,同時還具有快速的收斂速率．

1信號模型

假設(shè)有N+1個相互獨(dú)立窄帶信號入射到M元均勻陣列,入射信號與噪聲相互獨(dú)立,其中期望信號s0(t)的導(dǎo)向矢量記為a0,干擾信號的導(dǎo)向矢量為ad(d=1,2,…,N),接收信號記為

x(t)=As(t)+n(t).

(1)

最小方差無失真響應(yīng)(Minimum Variance Dis-

tortionless Response, MVDR)波束形成算法基于以下最優(yōu)問題[20]為

(2)

得出如下的最優(yōu)解為

(3)

(4)

(5)

2基于干擾加噪聲協(xié)方差矩陣重構(gòu)和導(dǎo)向矢量估計的穩(wěn)健波束形成

2.1干擾加噪聲協(xié)方差矩陣稀疏重構(gòu)

(6)

將式(6)中的壓縮感知問題分解為兩個獨(dú)立的子問題:通過分析可用的訓(xùn)練數(shù)據(jù)計算來波方向(Direction of Arrival, DOA);再利用得到的來波方向和不等式條件估計信源的功率．這兩步合起來便得到了式(6)的解．

得到DOA之后,公式(6)的壓縮感知問題就可以分解為一個不等式約束的最小均方問題

(7)

忽略限制條件,公式(7)的解可以表示為

(8)

式中:

G≡ vec(a(θp,1)aH(θp,1)),…,[

(9)

(10)

估計的空間譜為

(11)

(12)

式中:

(13)

,

(14)

(15)

修改后的功率譜為

(16)

用得到的稀疏譜重構(gòu)干擾加噪聲協(xié)方差矩陣:

(17)

y(t)=Px(t)=Pa0s0(t)+Pn(t).

(18)

從式(18)可以看出輸入信號與矩陣P相乘后只剩下信號成分和部分噪聲．

2.2利用旋轉(zhuǎn)子空間求解導(dǎo)向矢量

期望信號來波方向在一個已知的空域范圍Θ內(nèi),并且干擾信號不在這個范圍內(nèi)．因此,我們在這個范圍內(nèi)建立正交子空間,通過旋轉(zhuǎn)子空間來估計真實(shí)的導(dǎo)向矢量[17]．首先通過積分構(gòu)造正定矩陣:

C=∫Θa(θ)aH(θ)dθ.

(19)

a0=Ur.

(20)

式中,矢量r∈CK×1是所需的旋轉(zhuǎn)矢量．

由于矩陣P與干擾信號導(dǎo)向矢量ad正交,令權(quán)向量為

w=Pa0=PUr.

(21)

此時,權(quán)向量滿足條件:

wHad=0.

(22)

陣列輸出為

z(t)=wHy(t).

(23)

s.t. rHr=M.

(24)

式中:

(25)

(26)

利用Lagrange乘子法,最小化下面的函數(shù)為

(27)

(28)

2.3本算法主要步驟

為便于更好地理解本文算法,現(xiàn)就所提算法的步驟作如下歸納:

3仿真結(jié)果

為了驗(yàn)證上述分析,進(jìn)行了仿真驗(yàn)證,所有仿真實(shí)驗(yàn)都基于500次獨(dú)立實(shí)驗(yàn)．仿真是基于陣元間距為半波長的16元均勻線陣．入射信號包括一個期望信號和兩個功率相等的干擾信號．對兩個干擾信號,干噪比(Interference-to-NoiseRario,INR)為30dB,其來波方向分別為θ1=-30°和θ2=50°．假定期望信號的來波方向?yàn)?5°,期望信號可能的來波方向范圍為Θ=[θs-5°,θs+5°]．

最優(yōu)波束形成(Optimal)[24-25]是基于最大輸出信干噪比(SignaltoInterferenceplusNoiseRario,SINR)準(zhǔn)則得到的波束形成,是衡量波束形成算法性能的重要指標(biāo)．將本文提出的算法與下述算法進(jìn)行對比:傳統(tǒng)的MVDR算法[20],稀疏重構(gòu)(SparseReconstruction,SR)法[21],旋轉(zhuǎn)子空間(RotatingSubspace,RS)估計導(dǎo)向矢量法[22],基于特征子空間的投影法(ProjectionMethod,PM)[23]．其中RS估計導(dǎo)向矢量法中構(gòu)造信號子空間時主特征值個數(shù)取4個[22]．本文提出的算法中構(gòu)造矩陣U時主特征值的個數(shù)取4個．

仿真實(shí)驗(yàn)1: 在這個實(shí)驗(yàn)中,我們對比幾種算法的歸一化波束圖．這里快拍數(shù)取5 000．假設(shè)觀測角度誤差為-5°,也就是說真實(shí)的期望信號來波方向?yàn)?0°．期望信號對應(yīng)的信噪比為20 dB．圖1為幾種算法歸一化波束圖的對比．

由圖1可以看出,本文提出的算法(記為Proposed)不僅在干擾信號的來波方向上形成了較深的零陷,而且主瓣也對準(zhǔn)了真實(shí)的期望信號來波方向．作為對比的其他幾種算法,其中MVDR算法在強(qiáng)期望信號背景下,期望信號導(dǎo)向矢量存在誤差時,不僅主瓣處出現(xiàn)了下陷,也沒有在干擾信號的來波方向上形成零陷;SR法在干擾信號的來波方向上形成了較深的零陷,但是其主瓣并沒有對準(zhǔn)期望信號的來波方向;RS估計導(dǎo)向矢量法盡管主瓣對準(zhǔn)了期望信號的來波方向,同時也在干擾信號的來波方向上形成了零陷,但是其旁瓣過高,在強(qiáng)期望信號背景下主瓣甚至低于旁瓣;基于特征子空間的PM法的主瓣對準(zhǔn)了期望信號的來波方向,但是在干擾方向上并沒有形成零陷．綜上所述,本文提出的算法在強(qiáng)期望信號背景下,期望信號導(dǎo)向矢量存在誤差時,性能優(yōu)于其他幾種算法,在強(qiáng)期望信號背景下和期望信號導(dǎo)向矢量存在誤差時具有一定的穩(wěn)健性．這是由于與之對比的其他幾種算法只考慮到了期望信號過強(qiáng)時直接使用樣本協(xié)方差矩陣作為干擾加噪聲協(xié)方差矩陣的估計會因?yàn)榛烊肫谕盘柍煞謱?dǎo)致主瓣形成下陷;或是只考慮到期望信號導(dǎo)向矢量存在誤差時的情況．因此當(dāng)期望信號過強(qiáng)并且期望信號導(dǎo)向矢量存在誤差時,這些算法的性能都不如本文提出算法的性能．

圖1　輸入信噪比SNR為20 dB時的波束圖對比

仿真實(shí)驗(yàn)2: 在這個實(shí)驗(yàn)中,我們觀察幾種算法輸入SNR對輸出SINR的影響．這里,期望信號來波方向的隨機(jī)誤差滿足[-5°,5°]的均勻分布,期望信號的來波方向滿足[15°,25°]的均勻分布．期望信號的信噪比從-30 dB到50 dB變化．快拍數(shù)為5 000．

從圖2可以看出,本文提出的算法在輸入SNR較低的時候輸出SINR低于最優(yōu)的SINR,同時隨輸入SNR的增大而增大,隨著輸入SNR的逐漸增大,輸出SINR基本接近最優(yōu)的輸出SINR．作為對比的幾種算法,其中MVDR算法的輸出SINR一開始隨輸入SNR的增大而增大,當(dāng)輸入SNR增加到-10 dB左右的時候,輸出SINR隨輸入SNR的增大而減小;SR法的輸出SINR基本隨輸入SNR增大而線性增大,一開始比本文提出的算法的輸出SINR要高,但是隨著輸入SNR的增加最終要略低于本文提出的算法的輸出SINR; RS估計導(dǎo)向矢量法輸出SINR一開始隨輸入SNR的增大而增大,當(dāng)輸入SNR增大到15 dB左右,輸出SINR隨輸入SNR的增大而減?。鳳M法輸出SINR一開始隨輸入SNR的增大而增大,當(dāng)輸入SNR增大到20 dB左右,輸出SINR隨輸入SNR的增大而減小．綜上所述,當(dāng)輸入SNR不是特別弱的時候,本文所提出的算法性能優(yōu)于其他幾種算法．

圖2　輸出SINR隨期望信號輸入SNR變化曲線圖

仿真實(shí)驗(yàn)3: 在這個實(shí)驗(yàn)中,我們將期望信號來波方向的角度失配量固定,觀察不同算法的輸出SINR．這里,期望信號輸入SNR取20 dB,期望信號來波方向的角度失配量從-5°到5°變化．快拍數(shù)為5 000．

從圖3可以看出,當(dāng)期望信號來波方向失配量極大時,本文所提的算法依然與最優(yōu)的輸出SINR接近．而SR法只在失配量較小的時候接近最優(yōu)的輸出SINR．PM法的輸出SINR隨失配量絕對值的增大而減?。甊S估計導(dǎo)向矢量法SINR盡管不隨期望信號來波方向角度失配量的變化而變化,但是其輸出SINR與最優(yōu)的SINR有一定距離．傳統(tǒng)的MVDR算法對期望信號來波方向角度失配最敏感,從圖中可以看出極小的失配也會造成輸出SINR的極大下降．綜上可知,本文所提算法性能優(yōu)于其他幾種算法．

圖3　輸出SINR隨期望信號來波方向失配角度變化曲線圖

仿真實(shí)驗(yàn)4: 在第四個實(shí)驗(yàn)中我們研究快拍數(shù)對輸出SINR的影響．期望信號來波方向的隨機(jī)誤差滿足[-5°,5°]的均勻分布,期望信號的來波方向滿足[15°,25°]的均勻分布．期望信號對應(yīng)的信噪比為20 dB．這里隨機(jī)的期望信號來波方向每次蒙特卡洛實(shí)驗(yàn)都是隨機(jī)的,但對應(yīng)的快拍數(shù)固定．圖4為幾種算法輸出SINR隨快拍數(shù)變化的對比圖．從圖4可以看出,本文所提出的算法輸出SINR是最接近最優(yōu)的SINR的,并且除SR法外其輸出SINR的收斂速度是最快的．盡管SR法輸出SINR的收斂速度與本文所提算法相同,但是其輸出SINR總是低于本文所提算法的輸出SINR．PM與RS估計導(dǎo)向矢量法的輸出SINR的收斂趨勢比較接近,只是PM的輸出SINR要高于RS估計導(dǎo)向矢量法．而傳統(tǒng)的MVDR算法的性能是最差的．綜上可知,本文所提的算法不僅輸出SINR的收斂速度優(yōu)于其他算法,同時輸出SINR也是最接近最優(yōu)的SINR的,其性能最佳．

圖4　輸出SINR隨快拍數(shù)變化曲線圖

4結(jié)論

本文針對強(qiáng)期望信號背景下,期望信號導(dǎo)向矢量存在誤差時波束形成器性能下降的問題,提出了一種結(jié)合了稀疏重構(gòu)干擾加噪聲協(xié)方差矩陣和旋轉(zhuǎn)子空間估計期望信號導(dǎo)向矢量兩種算法的新算法．它不僅解決了期望信號過強(qiáng)時波束圖的主瓣處出現(xiàn)下陷的問題,同時還解決了期望信號導(dǎo)向矢量出現(xiàn)誤差時的穩(wěn)健性問題,此外它的收斂速度也很好．通過仿真對比分析,該算法比其他幾種算法性能更好,具有一定的理論意義和實(shí)用價值．

參考文獻(xiàn)

[1] WAX M, ANU Y. Performance analysis of the minimum variance beamformer[J]. IEEE transactions on signal processing, 1996, 44 (4): 928-937.

[2] WAX M, ANU Y. Performance analysis of the minimum variance beamformer in the presence of steering vector errors[J]. IEEE transactions on signal processing, 1996, 44(4):938-947.

[3] VAN TREES H L. Optimum array processing[M]. New York: John Wiley & Sons, 2002.

[4]李洪濤, 賀亞鵬, 朱曉華, 等. 基于譜分析的穩(wěn)健自適應(yīng)波束形成算法[J]. 電波科學(xué)學(xué)報, 2012, 27(1): 147-151.

LI H T, HE Y P, ZHU X H, et al. Spectral analysis based robust adaptive beamforming algorithm[J]. Chinese journal of radio science, 2012, 27(1): 147-151.(in Chinese)

[5] 李洪濤, 陳誠, 曾文浩, 等. 基于導(dǎo)向矢量估計的魯棒自適應(yīng)波束形成算法[J]. 電波科學(xué)學(xué)報, 2015, 30(1): 188-193.

LI H T, CHEN C, ZENG W H, et al. Steering vector estimation based robust adaptive beamforming algorithm[J]. Chinese journal of radio science, 2015, 30(1): 188-193. (in Chinese)

[6] FELDMAN D D, GRIFFITHS L J. A projection approach for robust adaptive beamforming[J]. IEEE transactions on signal processing, 1994, 42(4): 867-876.

[7] CARLSON B D. Covariance matrix estimation errors and diagonal loading in adaptive arrays[J]. IEEE transactions on aerospace and electronic system, 1988, 24(4): 397-401.

[8] LI J, STOICA P, WANG Z. On robust Capon beamforming and diagonal loading[J]. IEEE transactions on signal processing, 2003, 51(7): 1702-1715.

[9]VOROBYOV S A, GERSHMAN A B, LUO Z. Robust adaptive beamforming using worst-case performance optimization: a solution to the signal mismatch problem[J]. IEEE transactions on signal processing, 2003, 51(2): 313-324.

[10]LORENZ R G, BOYD S P. Robust minimum variance beamforming[J]. IEEE transactions on signal processing, 2005, 53(5):1684-1696.

[11]STURM J F. Using sedumi 1.02, a MATLAB toolbox for optimization over symmetric cones[J]. Optim methods software, 1999,11(1/4): 625-653.

[12]HASSANIEN A, VOROBYOV S A, WONG K M. Robust adaptive beamforming using sequential quadratic programming: an iterative solution to the mismatch problem[J]. IEEE signal processing letters, 2008, 15: 733-736.

[13]KHABBAZIBASMENJ A, VOROBYOV S A, HASSANIEN A. Robust adaptive beamforming via estimating steering vector based on semidefinite relaxation[C]//Conference on Signals, Systems & Computers. Pacific Grove: IEEE, 2010:1102-1106.

[14]MALLIPEDDI R, LIE J P, SUGANTHAN P N, et al. Robust adaptive beamforming based on covariance Matrix reconstruction for look direction mismatch[J]. Progress in electromagnetics research letters, 2011, 25: 37-46.

[15]KHABBAZIBASMENJ A, VOROBYOV S A, HASSANIEN A. Robust adaptive beamforming based on steering vector estimation with as little as possible prior information[J]. IEEE transactions on signal processing, 2012, 60(6): 2974-2987.

[16]GU Y, LESHEM A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation[J]. IEEE transactions on signal processing, 2012, 60(7): 3881-3885.

[17]ZHANG W, WANG J, WU S. Robust capon beamforming against large DOA mismatch[J]. Signal processing, 2013, 93(4):804-810.

[18]ZHUANG J, MANIKAS A. Interference cancellation beamforming robust to pointing errors[J]. IET signal processing, 2013, 7(2): 120-127.

[19]GRANT M, BOYD S, YE Y Y. CVX: MATLAB software for disciplined convex programming[EB/OL]. 2007[2015-05-26]. http://www.stanford.edu/boyd/cvx/V.1.0RC3.

[20]CAPON J. High-resolution frequency-wavenumber spectrum analysis[J]. Proceedings of the IEEE, 1969, 57(8):1408-1418.

[21]GU Y, GOODMAN N A, HONG S, et al. Robust adaptive beamforming based on interference covariance matrix sparse reconstruction[J]. Signal processing, 2014, 96(5): 375-381.

[22]ZHANG W, WANG J, WU S. Robust Capon beamforming against large DOA mismatch[J]. Signal processing, 2013, 93(4): 804-810.

[23]劉聰峰. 穩(wěn)健自適應(yīng)波束形成算法[M]. 西安:西安電子科技大學(xué)出版社, 2012: 36-99.

[24]龔耀寰. 自適應(yīng)濾波[M]. 北京: 電子工業(yè)出版社, 2003: 27-29, 295-297.

[25]曾浩, 楊士中, 馮文江. 基于SVD的特征空間波束合成器[J]. 電波科學(xué)學(xué)報, 2008, 23(6): 1194-1198.

ZENG H, YANG S Z, FENG W J. Eigenspace beamforming based on SVD[J]. Chinese journal of radio science, 2008, 23(6): 1194-1198. (in Chinese)

Robust adaptive beamforming based on reconstruction of the covariance matrix and the estimation of steering vector

XIE JulanLI XinyaLI HuiyongWANG Xu

(UniversityofElectronicScienceandTechnology,Chengdu611731,China)

AbstractIn the practical application, the error between the assumed steering vector of the desired signal and the real one causes the sharp degrading of the beamforming performance, especially when the desired signal power is strong. To solve this problem, a novel robust algorithm is proposed. The interference-plus-noise covariance matrix is sparse when the number of the signal source is less than the array element number. At first, by using this sparse property, the interference-plus-noise covariance matrix can be reconstructed, from which a subspace orthogonal to the steering vectors of interference signals can be obtained. Then, a mixed signal only containing the desired signal and the noise can be gained by mapping the received data through to this orthogonal subspace. The real steering vector of the desired signal can be estimated based on maximizing the out power of the mixed signal. The final beamforming weight can be constructed by using the orthogonal subspace and the real steering vector of the desired signal. Simulation results demonstrate the efficiency of the proposed algorithm in the case of the mismatch of the steering vector of the desired signal, the strong desired signal and the low snapshot number.

Keywordsrobust adaptive beamforming; the error of steering vector; the desired signal power; matrix reconstruction

收稿日期：2015-05-30

中圖分類號TN9117

文獻(xiàn)標(biāo)志碼A

文章編號1005-0388(2016)02-0297-07

DOI10.13443/j.cjors.2015053001

作者簡介

謝菊蘭(1981-),女,湖南人,電子科技大學(xué)電子工程學(xué)院副教授,研究方向?yàn)閿?shù)字波束形成與波束到達(dá)角估計等．

李昕亞(1990-),男,四川人,碩士研究生,研究方向?yàn)閿?shù)字波束形成與抗干擾算法等．

李會勇(1975-),男,湖北人,電子科技大學(xué)電子工程學(xué)院教授,研究方向?yàn)樽赃m應(yīng)陣列處理和MIMO系統(tǒng)．

王旭(1987-),男,四川人,博士研究生,研究方向?yàn)殛嚵锌諘r處理中的干擾抑制和波束形成．

謝菊蘭, 李昕亞, 李會勇, 等. 基于協(xié)方差矩陣重構(gòu)和導(dǎo)向矢量估計的穩(wěn)健自適應(yīng)波束形成[J]. 電波科學(xué)學(xué)報,2016,31(2):297-303. DOI: 10.13443/j.cjors.2015053001

XIE J L, LI X Y, LI H Y, et al. Robust adaptive beamforming based on reconstruction of the covariance matrix and the estimation of steering vector[J]. Chinese journal of radio science,2016,31(2):297-303. (in Chinese). DOI: 10.13443/j.cjors.2015053001

資助項(xiàng)目: 國家自然科學(xué)基金(61301262,61371184); 中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)基金(ZYGX2013J022)

聯(lián)系人： 李昕亞 E-mail：lixinya24@126.com