專利名稱：無(wú)線信號(hào)接收系統(tǒng)的陣列綜合測(cè)向方法
技術(shù)領(lǐng)域：
本發(fā)明屬無(wú)線接收范疇，涉及無(wú)線信號(hào)接收系統(tǒng)的陣列綜合測(cè)向方法，尤其涉及 由一部分相位同步的子陣列和另一部分相位不同步的子陣列組成的陣列的高分辨率 測(cè)向方法。
背景技術(shù)：
在無(wú)線信號(hào)接收系統(tǒng)中，當(dāng)由多個(gè)陣元組成的陣列接收電》茲波、聲波等無(wú)線信號(hào) 時(shí)，各個(gè)陣元接收的來(lái)波信號(hào)與某個(gè)指定的參考陣元之間存在一定的相位差。由于信 號(hào)的波達(dá)方向與該相位差之間存在一一對(duì)應(yīng)的關(guān)系，因此可以利用該關(guān)系設(shè)計(jì)陣列測(cè) 向技術(shù)，通過(guò)對(duì)陣列接收的信號(hào)進(jìn)行處理，實(shí)現(xiàn)陣列測(cè)向。在基于這種原理的現(xiàn)有陣 列測(cè)向方法中，典型的是具有高分辨率的多信號(hào)分類(lèi)(MUSIC, Multiple Signals Classification)方法。
這種方法的缺點(diǎn)是，當(dāng)各個(gè)陣元的信號(hào)接收機(jī)本身的相位未保持同步時(shí)，將在由 信號(hào)波達(dá)方向決定的相位差上面附加一個(gè)額外的、未知的、隨機(jī)的相位差，此時(shí)如果 仍利用相位差與波達(dá)方向之間的關(guān)系設(shè)計(jì)陣列測(cè)向方法，將導(dǎo)致測(cè)向性能出現(xiàn)顯著的 惡化。同時(shí)，對(duì)多個(gè)目標(biāo)同時(shí)進(jìn)行高分辨率測(cè)向時(shí)，需要使用較大孔徑和陣元數(shù)較多 的陣列。由于各個(gè)陣元的信號(hào)接收機(jī)的相位保持同步需要同一本振，當(dāng)陣列較大、陣 元數(shù)較多時(shí)，陣列測(cè)向系統(tǒng)的成本將顯著增加。

發(fā)明內(nèi)容
本發(fā)明的目的是為無(wú)線接收系統(tǒng)中相位不完全同步的陣列系統(tǒng)(由一部分相位同 步的子陣列和另一部分相位不同步的子陣列組成的陣列)提供一種較低成本的高分辨 率測(cè)向方法。本方法綜合利用參考信號(hào)的方向和其他信號(hào)與參考信號(hào)的方向差，估計(jì) 其他信號(hào)的波達(dá)方向。利用這種方法能在陣列較大、陣元數(shù)較多的情況下，盡可能降 低對(duì)各個(gè)陣元信號(hào)接收機(jī)之間的相位同步的要求，將在保持高分辨率測(cè)向性能的同 時(shí)，有效的降低陣列測(cè)向系統(tǒng)的成本。 本發(fā)明的目的是這樣達(dá)到的
首先利用不受相位不完全同步影響的陣列接收信號(hào)的模向量，得到其它信號(hào)與參考信號(hào)的方向差估計(jì)，再利用相位完全同步子陣列接收信號(hào)的復(fù)向量和完整的子空間 正交特性，得到參考信號(hào)的方向，從而完成所有信號(hào)的波達(dá)方向估計(jì)。 具體步驟是
第一步，確定相位不完全同步陣列接收信號(hào)的模向量及其自相關(guān)矩陣，對(duì)模自相 關(guān)矩陣進(jìn)行特征分解，確定相位不完全同步陣列的噪聲子空間；
第二步，確定相位不完全同步陣列的搜索向量，利用噪聲子空間與相位不完全同 步陣列的搜索向量之間的正交特性，確定相位不完全同步陣列的方向差空間謙及其峰 值位置，"^值位置確定其他信號(hào)與參考信號(hào)的方向差；
第三步，確定相位完全同步子陣列接收信號(hào)的自相關(guān)矩陣，對(duì)自相關(guān)矩陣進(jìn)行特 征分解，確定相位完全同步子陣列的噪聲子空間；
第四步，確定相位完全同步子陣列的搜索向量，利用噪聲子空間與相位完全同步 子陣列的搜索向量之間的完整正交特性確定參考信號(hào)空間鐠及其峰值位置，峰值位置 確定參考信號(hào)的方向；
最后，利用參考信號(hào)的方向和其它信號(hào)與參考信號(hào)的方向差估計(jì)其它信號(hào)的波達(dá) 方向。
參考信號(hào)指方向最小的信號(hào)。 所勤目位不完全同步陣列接收信號(hào)的模向量為
為f時(shí)刻相位不完全同步陣列的接收信號(hào)向量，
M為相位不完全同步陣列的陣元數(shù)，I |2表示取復(fù)數(shù)的模； 所述相位不完全同步的陣列接收樣本的模自相關(guān)矩陣為ifv =丄1>(/)/(0
其中，
其中，[f表示向量轉(zhuǎn)置，iV為陣列接收數(shù)據(jù)的快攝數(shù)。式中丄可以為任何非零值
計(jì)算相位不完全同步的陣列模自相關(guān)矩陣的特征分解為:
其中，矩陣力是以特征值為對(duì)角元素的對(duì)角矩陣，
<formula>formula see original document page 6</formula>
4矩陣2是以對(duì)應(yīng)的特征向量為列向量的矩陣；
所i^目位不完全同步的陣列噪聲子空間為
2 =[&+1 …《M]其中，矩陣^是相位不完全同步的陣列模自相關(guān)矩
陣^的M-Z個(gè)最小特征值對(duì)應(yīng)的特征向量構(gòu)成的矩陣。
所勤目位不完全同步的陣列搜索向量為 l
<formula>formula see original document page 7</formula>
其中，d為相臨陣元之間的間距，義為接收信號(hào)的波長(zhǎng)，上兩式的右邊均 可乘以任何非零常數(shù)；
所勤目位不完全同步的陣列的方向差空間譜為
<formula>formula see original document page 7</formula>
所勤目位同步子陣列接收樣本的自相關(guān)矩陣為及2 =丄5>(/)^( )
式中的右邊可乘以任何非負(fù)值。
其中，柳=
柳
為相位同步子陣列的接收信號(hào)向量，[f表示向量共軛轉(zhuǎn)置；
相位同步子陣列的自相關(guān)矩陣的特征分解為
及2=尸12尸" 其中，矩陣i2是以特征值為對(duì)角元素的對(duì)角矩陣，<formula>formula see original document page 7</formula>
，^ ^…2^ ,矩陣"是以對(duì)應(yīng)的特征向量為列向量的矩陣；
<formula>formula see original document page 7</formula>
所勤目位同步子陣列的噪聲子空間為/ =[仏+1 pA+2…pm]
其中，矩陣戶 是相位同步子陣列的自相關(guān)矩陣i z的w-個(gè)最小特征值對(duì)應(yīng)的特征向量構(gòu)成的矩陣； 所述相位同步子陣列的搜索向量為:
<formula>formula see original document page 8</formula>式中的右邊可以乘以任何非零常數(shù)。 所述參考信號(hào)方向空間譜為gSP^)=
<formula>formula see original document page 8</formula>
式中的右邊可
<formula>formula see original document page 8</formula>
以乘以任何非零值。
相位不完全相同陣列的方向差空間譜搜索式gISP =
的右邊乘以任何非負(fù)數(shù)時(shí)得到相位不完全相同陣列的方向差空間語(yǔ)的峰值位置，峰值位
置確定其他信號(hào)與參考信號(hào)的方向差A(yù) (z、l,2,…，D-1)，搜索式的右邊乘以任何負(fù)數(shù) 時(shí)是搜索空間譜的谷值位置得到的對(duì)應(yīng)方向。式中，下標(biāo)ISP表示的是相位不完全同步
的陣列，
參考信號(hào)方向空間譜搜索式gSP^)=^:
<formula>formula see original document page 8</formula>的右邊乘以任何非
負(fù)參考信號(hào)方向差空間譜的峰值位置，峰值位置確定參考信號(hào)的方向，搜索式的右邊
乘以任何負(fù)數(shù)時(shí)是搜索空間譜的谷值位置得到的對(duì)應(yīng)方向。式中，下標(biāo)SP表示相位不完 全同步陣列的同步子陣列。
所述相位不完全同步陣列是均勻線陣列或非均勻線陣列或圓陣列或隨機(jī)分布的陣 列，不受形態(tài)限制。
本發(fā)明的突出優(yōu)點(diǎn)是
1 、綜合確定相位不完全相同陣列的方向差空間語(yǔ)及其峰值位置和相位同步子陣 列的方向最小的參考信號(hào)空間譜及其峰值位置，利用相位不完全同步陣列的方向差空 間譜的峰值位置確定其他信號(hào)與參考信號(hào)的方向差；利用相位同步子陣列的參考方向 空間謙的峰值位置確定參考信號(hào)的方向。進(jìn)而利用參考信號(hào)的方向和其他信號(hào)與參考 信號(hào)的方向差估計(jì)其他信號(hào)的波達(dá)方向，確保達(dá)到高分辨率的效果 2、 能在陣列較大、陣元數(shù)較多的情況下，盡可能降低對(duì)各個(gè)陣元信號(hào)接收機(jī)之 間的相位同步的要求，做到較低成本。
3、 本方法使用范圍廣闊、不受接收系統(tǒng)陣列形狀的限制?？蓮V泛應(yīng)用于增強(qiáng)現(xiàn) 有陣列測(cè)向系統(tǒng)的分辨率，從而滿足智能天線陣、聲納陣列、無(wú)線電成像陣列等陣列 系統(tǒng)對(duì)高分辨率測(cè)向性能的要求。
4、根據(jù)本發(fā)明的綜合方法進(jìn)行實(shí)驗(yàn)表明，將本方法應(yīng)用于相位同步的的子陣列陣 元數(shù)為6個(gè)、相位不同步的子陣列陣元數(shù)為26、相鄰陣元間隔都為半波長(zhǎng)的相位不 完全同步的均勻線陣，在快攝數(shù)等于256的情況下，可以清晰的分辨信噪比都為 6. OdB、方向分別等于7. 0、 18.0、 22. 0度的三個(gè)信號(hào)，滿足測(cè)向系統(tǒng)對(duì)高分辨率性 能的要求。而在相同條件下常用的MUSIC方法已經(jīng)無(wú)法分辨這三個(gè)信號(hào)。


圖1示出本發(fā)明的方框流程圖。
圖2示出本發(fā)明的相位不完全同步陣列的結(jié)構(gòu)圖。
圖3示出3個(gè)信號(hào)情況下相位不完全同步陣列系統(tǒng)采用已有的陣列測(cè)向方法得到 的歸一化空間i普估計(jì)結(jié)果圖。
圖4示出3個(gè)信號(hào)情況下相位不完全同步陣列系統(tǒng)采用本發(fā)明的綜合測(cè)向方法得 到的歸一化方向差空間譜-估計(jì)結(jié)果圖。
圖5示出3個(gè)信號(hào)情況下相位不完全同步陣列系統(tǒng)采用本發(fā)明的綜合測(cè)向方法得 到的歸一化參考信號(hào)空間譜估計(jì)結(jié)果圖。
具體實(shí)施例方式
附圖也即實(shí)施例。下面結(jié)合附圖對(duì)本發(fā)明進(jìn)行詳細(xì)說(shuō)明。 以圖2所示的相位不完全同步的均勻線陣列為例子，陣列的陣元素為M，相臨陣元之 間的間距為"其中w個(gè)相臨的陣元構(gòu)成相位同步的子陣列，另外M-w個(gè)相臨的陣 元構(gòu)成相位不同步的子陣列。假設(shè)空間有D個(gè)非相干的窄帶信號(hào)A(O, K，A ，
到達(dá)此線陣，與線陣法線方向的夾角為《。以相位不完全同步的陣列的第一個(gè)陣元為 參考陣元，f時(shí)刻陣列的接收信號(hào)向量為
<formula>formula see original document page 9</formula>1
々(m-l)rfsin(e)
)與(m+l)"n((9)+処
2;r
(2)
(3)
(4)
仍(/ = 1,2，...，^-附)表示相位不同步子陣列各陣元接》|^的隨^目位，[f表示向量轉(zhuǎn)置，
;i為信號(hào)的波長(zhǎng)，v(o為加性白噪聲向量，與信號(hào)*0)相互獨(dú)立。
對(duì)應(yīng)相位不完全同步陣列的相位同步子陣列和相位不同步子陣列，^時(shí)刻相位不完全 同步陣列的接收信號(hào)向量又可以表示為
其中，z(0 =
竭=
柳 一)
—國(guó)
為相位不同
為相位同步子陣列的接收信號(hào)向量，M^): 步子陣列的接收信號(hào)向量。
參見(jiàn)附圖1。流程開(kāi)始于步驟201。在步驟2021,先確定相位不完全同步的陣列 接收樣本的模向量為
k(《 |W)|2
(5)
其中，I f表示取復(fù)數(shù)的模，W為陣列接收數(shù)據(jù)的快攝數(shù)。
在步驟2022，確定相位不完全同步的陣列接收樣本的模自相關(guān)矩陣為:<formula>formula see original document page 11</formula>
在步驟2023,計(jì)算相位不完全同步的陣列模自相關(guān)矩陣的特征分解為
<formula>formula see original document page 11</formula>
其中，矩陣J是以特征值為對(duì)角元素的對(duì)角矩陣， 、0 … 0 0 ^ … 0
<formula>formula see original document page 11</formula>
矩陣2是以對(duì)應(yīng)的特征向量為列向量的矩陣。
在步驟2024,確定相位不完全同步的陣列噪聲子空間為
<formula>formula see original document page 11</formula>
其中，矩陣2 是相位不完全同步的陣列模自相關(guān)矩陣^的M-￡個(gè)最小特征值對(duì)應(yīng)
的特征向量構(gòu)成的矩陣。
在步驟2031,確定相位不完全同步的陣列搜索向量為<formula>formula see original document page 11</formula>
和
湖=
其中<formula>formula see original document page 11</formula>
在步驟2032,先確定本發(fā)明設(shè)計(jì)的綜合測(cè)向技術(shù)的方向差空間譜:<formula>formula see original document page 12</formula>(11)
其中下標(biāo)'ISP，表示本發(fā)明的陣列測(cè)向方法是針對(duì)相位不完全同步(Incompletely Synchronized Phase )的陣列。然后通過(guò)搜索式(11)表示的方向差空間譜的峰值位置確 定其它信號(hào)與參考信號(hào)(指方向最小的信號(hào))的方向差-,(/ = 1，2，...，￡>-1)。
在步驟2041，確定相位同步子陣列接收樣本的自相關(guān)矩陣為:
<formula>formula see original document page 12</formula>
其中，[f表示向量共軛轉(zhuǎn)置。
在步驟2042，計(jì)算相位同步子陣列的自相關(guān)矩陣的特征分解為 i z=i^, (13)
其中，矩陣O是以特征值為對(duì)角元素的對(duì)角矩陣，
、0 ... 0
0 72 … 0 . . .
.00…"m
矩陣/2是以對(duì)應(yīng)的特征向量為列向量的矩陣。
在步驟2043，確定相位同步子陣列的噪聲子空間為
…A] (")
其中，矩陣/> 是相位同步子陣列的自相關(guān)矩陣Uz的w-個(gè)最小特征值對(duì)應(yīng)的特征向
量構(gòu)成的矩陣。
在步驟2051，確定相位同步子陣列的搜索向量為
1
v與—
(15)
其中，-2^0^三，/ = 1，2廣-，￡>-1 2 2在步驟2052,先確定本發(fā)明設(shè)計(jì)的綜合測(cè)向技術(shù)的參考信號(hào)方向空間譜
gSPW)=^l-^- (16)
'■=i
其中下標(biāo)'SP，表示相位不完全同步陣列的同步子陣列(Synchronous Phase )。然后通 過(guò)搜索式(16)表示的參考信號(hào)方向空間譜的峰值位置確定參考信號(hào)的方向，
在步驟206,利用參考信號(hào)的方向和其它信號(hào)與參考信號(hào)的方向差估計(jì)其它信號(hào) 的波達(dá)方向。
根據(jù)本方法的流程結(jié)束于步驟207。 下面對(duì)本方法的上述實(shí)施方式進(jìn)行更詳細(xì)的比較和舉例說(shuō)明。 已有方法輸出的空間鐠等于
gMUSIC(^=Z^^^ (17)
其中下標(biāo)'MUSIC，分別表示多信號(hào)分類(lèi)方法(MUSIC), R為樣本自相關(guān)矩陣及的 M-Z)個(gè)最小特征值對(duì)應(yīng)的特征向量構(gòu)成的矩陣。
參見(jiàn)圖3、圖4。在相位同步陣元數(shù)6，相位異步陣元數(shù)26，信號(hào)數(shù)3，波達(dá)方向(度)7.0,18.0,22.0j言""桑比(db)6.0,6.0,6.0的牙目同'清;兄下，與已有的MUSIC才目比，本發(fā) 明設(shè)計(jì)的綜合測(cè)向技術(shù)具有比MUSIC更高的方向分辨率。從圖3可以看出，已有的MUSIC 存在無(wú)法分辨的信號(hào)，而利用本發(fā)明能夠清晰的分辨所有信號(hào)。
對(duì)本發(fā)明的綜合測(cè)向方法雖然已經(jīng)參考附圖以舉例方式進(jìn)行了描述，但是本發(fā)明不 限于上述這些細(xì)節(jié)，本申請(qǐng)含蓋權(quán)利要求范圍之內(nèi)的各種變型或改變。
權(quán)利要求
1、一種無(wú)線信號(hào)接收系統(tǒng)的陣列綜合測(cè)向方法，其特征在于首先利用不受相位不完全同步影響的陣列接收信號(hào)的模向量，得到其它信號(hào)與參考信號(hào)的方向差估計(jì)，再利用相位完全同步子陣列接收信號(hào)的復(fù)向量和完整的子空間正交特性，得到參考信號(hào)的方向，完成所有信號(hào)的波達(dá)方向估計(jì)，具體步驟是第一步，確定相位不完全同步陣列接收信號(hào)的模向量及其自相關(guān)矩陣，對(duì)模自相關(guān)矩陣進(jìn)行特征分解，確定相位不完全同步陣列的噪聲子空間；第二步，確定相位不完全同步陣列的搜索向量，利用噪聲子空間與相位不完全同步陣列的搜索向量之間的正交特性，確定相位不完全同步陣列的方向差空間譜及其峰值位置，峰值位置確定其他信號(hào)與參考信號(hào)的方向差；第三步，確定相位完全同步子陣列接收信號(hào)的自相關(guān)矩陣，對(duì)自相關(guān)矩陣進(jìn)行特征分解，確定相位完全同步子陣列的噪聲子空間；第四步，確定相位完全同步子陣列的搜索向量，利用噪聲子空間與相位完全同步子陣列的搜索向量之間的完整正交特性，確定參考信號(hào)空間譜及其峰值位置，峰值位置確定參考信號(hào)的方向；最后，利用參考信號(hào)方向和其它信號(hào)與參考信號(hào)的方向差估計(jì)其它信號(hào)的波達(dá)方向。
2. 根據(jù)權(quán)利要求l所述的陣列綜合測(cè)向方法，其特征在于所述參考信號(hào)是指方向 最小的信號(hào)。
3. 根據(jù)權(quán)利要求1或2或3所述的陣列綜合測(cè)向方法，其特征在于所述相位不完 全同步陣列接收信號(hào)的模向量為'為/時(shí)刻相位不完全同步陣列的接收信號(hào)向量，M為相位不完全同步陣列的陣元數(shù)，I |2表示取復(fù)數(shù)的模； 所勤目位不完全同步的陣列接收樣本的模自相關(guān)矩陣為J v =丄|^(,)/(0<formula>formula see original document page 2</formula>其中，[f表示向量轉(zhuǎn)置，iV為陣列接收數(shù)據(jù)的快攝數(shù)，式中丄可以為任何非零值； 計(jì)算相位不完全同步的陣列模自相關(guān)矩陣的特征分解為其中，矩陣J是以特征值為對(duì)角元素的對(duì)角矩陣，<formula>formula see original document page 3</formula>矩陣g是以對(duì)應(yīng)的特征向量為列向量的矩陣；所勤目位不完全同步的陣列噪聲子空間為O^kw …其中，矩陣込是相位不完全同步的陣列模自相關(guān)矩陣A的M-Z個(gè)最小特征值對(duì)應(yīng)的特征向量構(gòu)成的矩陣。
4.根據(jù)權(quán)利要求1或2所述的陣列綜合測(cè)向方法，其特征在于所述相位不完全同 步的陣列搜索向量為<formula>formula see original document page 3</formula>和<formula>formula see original document page 3</formula>其中，d為相臨陣元之間的間距，義為接收信號(hào)的波長(zhǎng)，os^si，上兩式的右邊均可乘以任何非零常數(shù)；所勤目位不完全同步的陣列的方向差空間譜為<formula>formula see original document page 3</formula>式中的右邊可乘以任何非負(fù)值.
5.根據(jù)權(quán)利要求1或2所述的陣列綜合測(cè)向方法，其特征在于所勤目位同步子陣列接收樣本的自相關(guān)矩陣為及2=丄|> ^(0 其中，洲=列的接收信號(hào)向量，[f表示向量共軛轉(zhuǎn)置；為相位同步子陣相位同步子陣列的自相關(guān)矩陣的特征分解為<formula>formula see original document page 0</formula>其中，矩陣D是以特征值為對(duì)角元素的對(duì)角矩陣，<formula>formula see original document page 4</formula>所勤目位同步子陣列的噪聲子空間為尸 =[^+1 A+2…Pm] 其中，矩陣戶 是相位同步子陣列的自相關(guān)矩陣&的m- A個(gè)最小特征值對(duì)應(yīng)的特征向量 構(gòu)成的矩陣。
6.根據(jù)權(quán)利要求5所述的陣列綜合測(cè)向方法，其特征在于所述相位同步子陣列的 搜索向量為<formula>formula see original document page 4</formula>式中的右邊可以乘以任何非零常數(shù)。
7.根據(jù)權(quán)利要求1或2所述的陣列綜合測(cè)向方法，其特征在于所勤目位不完全相同的陣列的方向差空間i普搜索式<formula>formula see original document page 4</formula>的右邊乘以任何非負(fù)數(shù)時(shí)得到相位不完全相同信號(hào)的峰值位置，峰值位置確定對(duì)應(yīng)信號(hào)與參考信號(hào)的方向差是^ (/ = 1，2,-.，"-1);搜索式的右邊乘以任何負(fù)數(shù)時(shí)，是搜索空間諳的谷值位置得到的對(duì)應(yīng)方向。
8.根據(jù)權(quán)利要求1或2所述的陣列綜合測(cè)向方法，其特征在于所述參考信號(hào)方1向空間譜搜索式<formula>formula see original document page 4</formula>的右邊乘以任何非負(fù)數(shù)時(shí)得到參考信號(hào)方向空間譜的峰值位置，峰值位置確定參考信號(hào)的方向；搜索式的右邊乘以任 何負(fù)數(shù)時(shí)，是搜索空間譜的谷值位置得到的對(duì)應(yīng)方向。
9.根據(jù)權(quán)利要求4所述的陣列綜合測(cè)向方法，其特征在于所勤目位不完全同步陣 列是均勻線陣列或非均勻線陣列或圓陣列或隨機(jī)分布的陣列。
全文摘要
無(wú)線信號(hào)接收系統(tǒng)的陣列綜合測(cè)向方法屬信號(hào)接收范疇。步驟確定相位不完全同步陣列接收信號(hào)的模向量及其自相關(guān)矩陣，通過(guò)模自相關(guān)矩陣特征分解確定噪聲子空間；利用噪聲子空間與相位不完全同步陣列搜索向量間的正交特性確定方向差空間譜，由其峰值位置確定其它信號(hào)與參考信號(hào)的方向差；確定相位完全同步子陣列接收信號(hào)自相關(guān)矩陣，通過(guò)自相關(guān)矩陣特征分解確定噪聲子空間；利用噪聲子空間與相位完全同步子陣列搜索向量之間完整的正交特性確定參考信號(hào)方向空間譜，由其峰值位置確定參考信號(hào)方向；利用參考信號(hào)方向和其它信號(hào)與參考信號(hào)的方向差估計(jì)其它信號(hào)的方向。本方法能降低各信號(hào)接收機(jī)間的相位同步要求，保持高分辨率，有效降低系統(tǒng)成本。
文檔編號(hào)H04B7/04GK101309101SQ200710049079
公開(kāi)日2008年11月19日 申請(qǐng)日期2007年5月14日 優(yōu)先權(quán)日2007年5月14日
發(fā)明者群 萬(wàn), 楊萬(wàn)麟, 楊瑞明, 沈曉峰, 衡 竇 申請(qǐng)人:電子科技大學(xué)