郭翠娟,龔楚云,榮 鋒,宋雅琪
(天津工業(yè)大學(xué)電子與信息工程學(xué)院,天津300387)
基于熵權(quán)法的實(shí)時(shí)振動(dòng)信號(hào)存儲(chǔ)管理技術(shù)
郭翠娟,龔楚云,榮 鋒,宋雅琪
(天津工業(yè)大學(xué)電子與信息工程學(xué)院,天津300387)
針對(duì)實(shí)時(shí)振動(dòng)信號(hào)的存儲(chǔ)問題,提出一種基于熵權(quán)法的實(shí)時(shí)振動(dòng)信號(hào)存儲(chǔ)管理模式.根據(jù)振動(dòng)信號(hào)的波動(dòng)特點(diǎn),將采集到的信號(hào)波形從時(shí)域和頻域2個(gè)方面分別計(jì)算出相關(guān)的指標(biāo)參數(shù),并采用熵權(quán)法對(duì)各參數(shù)進(jìn)行加權(quán)計(jì)算,按照加權(quán)結(jié)果對(duì)各段波形進(jìn)行排序,從理論角度對(duì)比甄選出含有有用信息較多的振動(dòng)信號(hào)波形,舍棄較為無用的波形信號(hào),從而節(jié)省數(shù)據(jù)庫(kù)存儲(chǔ)空間,方便后續(xù)信號(hào)的存儲(chǔ),而且使得有用波形可以完整攜帶并保存與機(jī)械運(yùn)行狀態(tài)相關(guān)的信號(hào).
振動(dòng)信號(hào);大數(shù)據(jù);熵權(quán)法;數(shù)據(jù)存儲(chǔ)
對(duì)旋轉(zhuǎn)機(jī)械的振動(dòng)信號(hào)進(jìn)行監(jiān)測(cè)與處理有助于實(shí)時(shí)了解設(shè)備的運(yùn)行狀態(tài)[1].采集的振動(dòng)信號(hào)通常被實(shí)時(shí)存儲(chǔ)到數(shù)據(jù)庫(kù)系統(tǒng)中,便于后續(xù)的分析處理.當(dāng)前,國(guó)內(nèi)外主流的工業(yè)實(shí)時(shí)數(shù)據(jù)庫(kù)產(chǎn)品一般都提供了歷史數(shù)據(jù)存儲(chǔ)功能,能較長(zhǎng)時(shí)間地保存?zhèn)浞莠F(xiàn)場(chǎng)所采集到的數(shù)據(jù)[2].由于機(jī)械需要長(zhǎng)期運(yùn)行,隨著采集時(shí)間的增加,采集點(diǎn)增加,采集時(shí)間周期又較短,因此要保存的振動(dòng)信號(hào)的數(shù)據(jù)量也隨之不斷增大,而數(shù)據(jù)庫(kù)容量有限,如果將這些數(shù)據(jù)直接存儲(chǔ),不僅會(huì)占用大量的系統(tǒng)存儲(chǔ)空間,無法滿足長(zhǎng)時(shí)間存儲(chǔ)的需求,降低數(shù)據(jù)庫(kù)的實(shí)時(shí)性,而且數(shù)據(jù)庫(kù)的安全性和穩(wěn)定性也會(huì)受到影響,數(shù)據(jù)的傳輸、查詢將變得十分困難[3].
針對(duì)海量信息存儲(chǔ)問題,傳統(tǒng)的解決方案大多采用網(wǎng)絡(luò)存儲(chǔ),然而網(wǎng)絡(luò)存儲(chǔ)需要專用服務(wù)器和專用磁盤陣列,成本昂貴,并且磁盤運(yùn)行壽命有限,故障率高[4].而對(duì)數(shù)據(jù)庫(kù)進(jìn)行邏輯劃分,即將數(shù)據(jù)分布到不同的服務(wù)器中進(jìn)行存儲(chǔ),需要在數(shù)據(jù)增加時(shí)重新進(jìn)行劃分,進(jìn)而要求修改程序并丟棄模型的獨(dú)立性[5].對(duì)關(guān)系數(shù)據(jù)進(jìn)行非規(guī)范化處理增加了數(shù)據(jù)的冗余,雖然有利于數(shù)據(jù)分布到多個(gè)節(jié)點(diǎn)上進(jìn)行并行操作,卻沒能節(jié)省數(shù)據(jù)庫(kù)容量.目前常用的數(shù)據(jù)壓縮技術(shù)即對(duì)數(shù)據(jù)進(jìn)行壓縮再存儲(chǔ),可以提供較高的壓縮率,因此歷史數(shù)據(jù)占用的硬盤空間少,在存儲(chǔ)及調(diào)用歷史數(shù)據(jù)時(shí)對(duì)硬盤的讀寫減少,能夠在保持?jǐn)?shù)據(jù)庫(kù)高實(shí)時(shí)性等特點(diǎn)的同時(shí),盡可能地提高數(shù)據(jù)庫(kù)的容量,從而節(jié)省了系統(tǒng)資源[6].數(shù)據(jù)壓縮技術(shù)主要有有損壓縮算法和無損壓縮算法兩種.有損算法依據(jù)數(shù)據(jù)的特征,通常采用了特殊舍點(diǎn)算法,該種算法會(huì)帶來不可恢復(fù)的損失與誤差,變換后數(shù)據(jù)變得很平滑,高值和低值被濾掉,也就不能準(zhǔn)確地由波形信息推測(cè)出設(shè)備運(yùn)行狀態(tài).而無損算法更不能從根本上解決數(shù)據(jù)庫(kù)容量有限的問題.由此觀之,在大數(shù)據(jù)背景下,信號(hào)處理的瓶頸已由“信號(hào)采集量太少”轉(zhuǎn)移到“如何對(duì)海量信號(hào)數(shù)據(jù)進(jìn)行高速有效處理”[7].
本文提出一種振動(dòng)信號(hào)存儲(chǔ)管理模式.根據(jù)振動(dòng)信號(hào)的波動(dòng)特點(diǎn),計(jì)算所采集的信號(hào)在時(shí)域與頻域中的相關(guān)指標(biāo)參數(shù),對(duì)這些參數(shù)按照熵權(quán)法[8]賦予權(quán)值進(jìn)行加權(quán)計(jì)算,再將所采集的信號(hào)按照加權(quán)結(jié)果進(jìn)行排序,對(duì)比甄選出含有有用信息較多的振動(dòng)信號(hào),節(jié)省數(shù)據(jù)庫(kù)存儲(chǔ)空間.
信號(hào)是信息的載體,選擇適當(dāng)?shù)男盘?hào)處理方法對(duì)信號(hào)進(jìn)行處理,有利于后續(xù)的狀態(tài)監(jiān)測(cè)以及故障診斷工作的開展.通常,振動(dòng)信號(hào)時(shí)域特征參數(shù)有振幅、周期、相位等,頻域特征則主要包含在頻率、能量信息中.采集板卡采集到振動(dòng)信號(hào)之后,時(shí)域連續(xù)信號(hào)經(jīng)過采樣成為離散信號(hào),上位機(jī)通過以太網(wǎng)接收離散振動(dòng)信號(hào),保存在數(shù)據(jù)庫(kù)中.通過對(duì)離散信號(hào)進(jìn)行分析處理,結(jié)合機(jī)械工作原理,可以對(duì)設(shè)備運(yùn)行狀態(tài)做出粗略的判斷,為狀態(tài)維修提供維修策略.
1.1 信號(hào)的時(shí)域指標(biāo)參數(shù)計(jì)算
中心頻率在信號(hào)的時(shí)域中描述幅值隨時(shí)間的變化關(guān)系稱作幅值域分析.幅值域分析是信號(hào)處理中最常用的信號(hào)分析方法.信號(hào)的幅值域參數(shù)主要包括峰峰值、峰值、均值、有效值、方差、峭度指標(biāo)等:
(1)峰峰值,主要用來描述振動(dòng)的幅值,即機(jī)械振動(dòng)的位移量,由波形上最高點(diǎn)與最低點(diǎn)之差表示,計(jì)算公式為:
(2)峰值,是指波形上與零線最大的偏移量,振動(dòng)信號(hào)的加速度值常用峰值描述.根據(jù)經(jīng)驗(yàn),一般在檢測(cè)剝落、裂紋、壓痕、凹坑等原因造成的沖擊性振動(dòng)時(shí),峰值會(huì)急劇增大,計(jì)算公式為:
(3)均值,用來描述信號(hào)的平均水平,也稱數(shù)學(xué)期望或一次矩,反映了信號(hào)變化的中心趨勢(shì),計(jì)算公式為:
(4)有效值,即信號(hào)的均方根值,可以描述振動(dòng)的強(qiáng)烈程度,計(jì)算公式為:
(5)方差,反映了信號(hào)繞均值的波動(dòng)程度,也是反映信號(hào)的離散程度,計(jì)算公式為:
(6)峭度指標(biāo),可以敏感捕捉信號(hào)中的沖擊成分,是描述波形尖峰度的指標(biāo),計(jì)算公式為:
峭度指標(biāo)的計(jì)算公式中分子為x的4次方,而分母是一個(gè)平均量,這就必然導(dǎo)致分子的增加快于分母,若峭度指標(biāo)的值上升很快,說明故障已出現(xiàn),從而大大提高故障診斷準(zhǔn)確度.
利用時(shí)域指標(biāo)來分析振動(dòng)信號(hào),優(yōu)點(diǎn)是直觀、便于理解和直接計(jì)算,但只能做定性分析.從理論上講,對(duì)信號(hào)來說,時(shí)域分析對(duì)時(shí)間的分辨率是無窮的,但對(duì)頻率的分辨率為零;而頻域分析對(duì)頻率的分辨率是無窮的,但對(duì)時(shí)間的分辨率為零.因此,還需要結(jié)合信號(hào)的頻域分析,才能更準(zhǔn)確地評(píng)估出信號(hào)所包含有用信息量的多少,才能更有利于正確評(píng)估信號(hào)的有用價(jià)值.
1.2 信號(hào)的頻域指標(biāo)參數(shù)計(jì)算
頻域和時(shí)域表明了動(dòng)態(tài)信號(hào)的2個(gè)觀察面,即這2種觀察信號(hào)方法以不同的角度揭示了信號(hào)的物理特征,而傅里葉變換建立起它們之間的聯(lián)系.
傅里葉變換可以看作是時(shí)間函數(shù)在頻率域上的表示.由傅里葉變換給出的頻率域包含的信息和原函數(shù)時(shí)間域內(nèi)所包含的完全相同,不同的僅是信息的表示形式.根據(jù)振動(dòng)信號(hào)采集已采樣成為離散點(diǎn)存儲(chǔ)在數(shù)據(jù)庫(kù)的實(shí)際情況,本文選擇更適合的離散傅里葉變換(DFT)來做離散信號(hào)的頻域分析.DFT的表達(dá)方式如下:
式中:DFS[*]表示離散傅里葉正變換;IDFS[*]表示離散傅里葉逆變換.利用傅里葉變換,可以得到幅度譜.在時(shí)間計(jì)算時(shí),可以使用快速傅里葉變換(FFT),以減少運(yùn)算量.
另外,功率譜也是一個(gè)常用的頻域指標(biāo).功率譜反映了隨機(jī)信號(hào)各頻率成份功率能量的分布情況.當(dāng)機(jī)械設(shè)備發(fā)生故障時(shí),會(huì)在某些特征頻率附近出現(xiàn)明顯的峰值,功率譜反映了各頻率成分能量的分布關(guān)系,更容易得到峰值所對(duì)應(yīng)的頻率,與設(shè)備特征頻率進(jìn)行對(duì)比分析,有助于對(duì)設(shè)備進(jìn)行故障診斷.對(duì)于連續(xù)信號(hào)來說,當(dāng)采樣長(zhǎng)度為T0時(shí),其功率譜密度為:
當(dāng)采樣間隔為T、采樣點(diǎn)為N時(shí),T0=NT,由連續(xù)傅里葉變換和離散傅里葉變換關(guān)系得:
便可以得到離散功率譜密度表達(dá)式:
通過以上公式,可以將信號(hào)的使用價(jià)值大小定量化,方便作為后續(xù)排序的依據(jù).
對(duì)振動(dòng)信號(hào)分析中各指標(biāo)參數(shù)的關(guān)心程度通常不盡相同,對(duì)于在監(jiān)測(cè)機(jī)械運(yùn)轉(zhuǎn)與檢測(cè)機(jī)械故障方面更能直觀、方便和準(zhǔn)確地提供信息的參數(shù)則應(yīng)該多給予關(guān)心.本文方法的基本思想是:對(duì)于已經(jīng)計(jì)算出來的各參數(shù),對(duì)應(yīng)地給予一組非負(fù)數(shù)權(quán)值ωi(i=1~n),然后將各指標(biāo)參數(shù)計(jì)算值乘以相應(yīng)的加權(quán)系數(shù),求和,得到振動(dòng)信號(hào)各項(xiàng)指標(biāo)的加權(quán)結(jié)果.基于多指標(biāo)與客觀性考慮,本文選用熵權(quán)法進(jìn)行權(quán)值的確定[9].
熵是信息論中的一個(gè)概念,熵權(quán)法的基本思路是根據(jù)指標(biāo)變異性的大小來確定客觀權(quán)重,在多指標(biāo)綜合評(píng)定方面幾乎不受主觀因素干擾,很有優(yōu)勢(shì)[10].一般來說,若某個(gè)指標(biāo)的信息熵越小,表明指標(biāo)值的變異程度越大,提供的信息量越多,在綜合評(píng)價(jià)中所能起到的作用也越大,其權(quán)重也就越大.相反,某個(gè)指標(biāo)的信息熵越大,表明指標(biāo)值的變異程度越小,提供的信息量越少,在綜合評(píng)價(jià)中所起到的作用越小,其權(quán)重也就越小[11-12].
利用熵權(quán)法,結(jié)合實(shí)際情況,本文所提出的振動(dòng)信號(hào)指標(biāo)參數(shù)加權(quán)計(jì)算過程如下:
(1)根據(jù)實(shí)際情況需要,將每一天的波形數(shù)據(jù)分為n組進(jìn)行檢測(cè),即確定n個(gè)對(duì)象;
(2)計(jì)算出n組波形的時(shí)域以及頻域各指標(biāo)參數(shù)值,共有m項(xiàng)指標(biāo).從而得到矩陣:
(3)對(duì)指標(biāo)參數(shù)矩陣進(jìn)行標(biāo)準(zhǔn)化處理得到新的矩陣:
式中:rij為第i個(gè)對(duì)象在第j個(gè)評(píng)價(jià)指標(biāo)上的標(biāo)準(zhǔn)值.其標(biāo)準(zhǔn)化計(jì)算公式為:
公式(14)適用于指標(biāo)數(shù)值越大越好的情況下,公式(15)則適用于指標(biāo)數(shù)值越小越好的情況.
(4)計(jì)算第j項(xiàng)指標(biāo)的熵:
(5)計(jì)算各指標(biāo)權(quán)重:
(6)根據(jù)各指標(biāo)的熵與權(quán)值,便可以得到最終的加權(quán)結(jié)果Vi:
以上便是利用熵權(quán)法對(duì)振動(dòng)信號(hào)的各項(xiàng)指標(biāo)進(jìn)行加權(quán)計(jì)算的全部過程.根據(jù)得到的加權(quán)結(jié)果Vi,便可以對(duì)抽取出來的波形信號(hào)信息量進(jìn)行排序,從而刪除包含信息量最少的信號(hào),保留富含信息量的信號(hào),方便后續(xù)的信號(hào)分析處理.
振動(dòng)信號(hào)存儲(chǔ)在數(shù)據(jù)庫(kù)中之后,根據(jù)實(shí)驗(yàn)需要,可以先抽取5組數(shù)據(jù),按照上述方法進(jìn)行計(jì)算.
本次實(shí)驗(yàn)利用振動(dòng)信號(hào)采集系統(tǒng)分別采集5組波形,其中前4組分別為正常運(yùn)行的波形產(chǎn)生的參數(shù),第5組為故障波形產(chǎn)生的參數(shù).正常振動(dòng)時(shí)產(chǎn)生的波形以及故障時(shí)產(chǎn)生的波形分別如圖1和圖2所示.
圖1 正常振動(dòng)波形圖Fig.1 Normal vibration signal
圖2 振動(dòng)故障波形圖Fig2Vibration signal with bugs
根據(jù)前邊所述過程,計(jì)算結(jié)果如表1所示.
表1 5組振動(dòng)波形各項(xiàng)指標(biāo)計(jì)算結(jié)果Tab.1 Indexes of five vibration signals
對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化后,得到結(jié)果如表2所示.
表2 5組振動(dòng)波形各項(xiàng)指標(biāo)計(jì)算標(biāo)準(zhǔn)化后結(jié)果Tab.2 Indexes of five vibration signals after standardization
根據(jù)信息熵的計(jì)算公式,得到各指標(biāo)信息熵如表3所示.
表3 各項(xiàng)指標(biāo)信息熵表Tab.3 Entropy weight of each index
根據(jù)得到的各指標(biāo)信息熵,計(jì)算出各指標(biāo)權(quán)重如表4所示.
表4 各指標(biāo)權(quán)重Tab.4 Weight of each index
根據(jù)計(jì)算出的指標(biāo)權(quán)重,5組波形加權(quán)結(jié)果如表5所示.
表5 5組波形評(píng)定結(jié)果Tab.5 Result of five wave caculate
根據(jù)最終的加權(quán)結(jié)果,可知第2組波形變化較為平緩,相應(yīng)時(shí)段機(jī)械運(yùn)轉(zhuǎn)并無太大異常,所以并不需要對(duì)這一時(shí)段產(chǎn)生的波形做進(jìn)一步地分析來得知機(jī)械是否故障,因此可以舍棄,節(jié)省數(shù)據(jù)庫(kù)的存儲(chǔ)空間.
針對(duì)實(shí)時(shí)振動(dòng)信號(hào)采集系統(tǒng)的數(shù)據(jù)庫(kù)存儲(chǔ)容量有限,數(shù)據(jù)不能長(zhǎng)期穩(wěn)定且安全地存儲(chǔ)在數(shù)據(jù)庫(kù)中,而壓縮數(shù)據(jù)又不能完整還原波形信號(hào)的問題,提出一種新的振動(dòng)信號(hào)存儲(chǔ)管理模式,將采集到的振動(dòng)信號(hào)波形從時(shí)域與頻域2個(gè)方面分別計(jì)算出相關(guān)的波形指標(biāo)參數(shù),將計(jì)算得到的參數(shù)按照關(guān)心程度賦予權(quán)值,進(jìn)行加權(quán)計(jì)算,便可以對(duì)各段波形按照權(quán)值進(jìn)行排序,從理論角度對(duì)比甄選出含有有用信息較多的振動(dòng)信號(hào)波形,繼續(xù)保存在數(shù)據(jù)庫(kù)中,舍棄較為無用的波形信號(hào).這樣不但節(jié)省了數(shù)據(jù)庫(kù)存儲(chǔ)空間,而且使得有用波形可以完整攜帶并保存與機(jī)械運(yùn)行狀態(tài)相關(guān)的信號(hào).
本文提出的基于熵權(quán)法的實(shí)時(shí)振動(dòng)信號(hào)存儲(chǔ)管理模式,是基于實(shí)驗(yàn)室已有的硬件振動(dòng)設(shè)備平臺(tái),研究平臺(tái)上位機(jī)系統(tǒng)如何能應(yīng)對(duì)硬件系統(tǒng)采集的大量振動(dòng)信號(hào)數(shù)據(jù)而創(chuàng)新得出的一種大數(shù)據(jù)管理模式.由于該數(shù)據(jù)管理模式是基于信號(hào)特征參數(shù)計(jì)算得出的,而模擬和數(shù)字信號(hào)都有著諸如峰峰值、均值、有效值、方差和功率譜等相同的特征參量,因此理論上看似可以推廣應(yīng)用到其他信號(hào)如語(yǔ)音、視頻、圖像等信號(hào)的存儲(chǔ).但實(shí)際上,并不是所有信號(hào)都適用這一方法.因?yàn)樵谠O(shè)備監(jiān)測(cè)領(lǐng)域,一旦設(shè)備發(fā)生故障,設(shè)備故障部位產(chǎn)生的振動(dòng)信號(hào)會(huì)發(fā)生突變,其各方面特征參量也有別于正常運(yùn)作時(shí)產(chǎn)生的振動(dòng)信號(hào),因此通過算法計(jì)算,很容易找到突變信號(hào)而了解到故障所在.但語(yǔ)音、視頻和圖像信息的獲取目的卻與設(shè)備監(jiān)測(cè)不同,或許并不是為了監(jiān)測(cè)故障,而是為了進(jìn)行圖像或者語(yǔ)音信號(hào)的處理,那么按照本文方法并不能得出加權(quán)值大的信號(hào)就有利于信號(hào)處理的結(jié)論.非同一目的,不能混淆,這也是本文所提出的基于熵權(quán)法的實(shí)時(shí)振動(dòng)信號(hào)存儲(chǔ)管理模式的局限之處,還需作進(jìn)一步研究.
[1]彭暢.旋轉(zhuǎn)機(jī)械軸承振動(dòng)信號(hào)分析方法研究[D].重慶:重慶大學(xué),2014.
[2]賴達(dá)波.某齒輪箱故障振動(dòng)信號(hào)特征提取及分析技術(shù)研究[D].成都:電子科技大學(xué),2013.
[3]徐國(guó)風(fēng).實(shí)時(shí)數(shù)據(jù)庫(kù)關(guān)鍵技術(shù)研究[D].西安:西安建筑科技大學(xué),2006.
[4]吳劍強(qiáng).流程工業(yè)實(shí)時(shí)數(shù)據(jù)庫(kù)系統(tǒng)研究與開發(fā)[D].杭州:浙江大學(xué),2004.
[5]欒振華.工業(yè)實(shí)時(shí)數(shù)據(jù)庫(kù)的應(yīng)用設(shè)計(jì)及數(shù)據(jù)處理研究[D].杭州:浙江大學(xué),2008.
[6]HOURANI B,SALIM R.A performance analysis framework for the design of DSP systems.North Carolina State University,2008(4):100-104.
[7]朱林.基于特征加權(quán)與特征選擇的數(shù)據(jù)挖掘算法研究[D].上海:上海交通大學(xué),2013.
[8]劉志,李軍華,胡克澤,等.熵權(quán)法在企業(yè)績(jī)效綜合評(píng)價(jià)中的應(yīng)用[J].石油化工管理干部學(xué)院學(xué)報(bào),2008,4(10):63-69.
[9]王昆,宋海洲.三種客觀權(quán)重賦權(quán)法的比較分析[J].經(jīng)濟(jì)技術(shù)與管理研究,2003(6):48-50.
[10]章穗,張梅,遲國(guó)泰.基于熵權(quán)法的科學(xué)技術(shù)評(píng)價(jià)模型及其實(shí)證研究[J].管理學(xué)報(bào),2010,7(1):34-41.
[11]高孝偉.熵權(quán)法在教學(xué)評(píng)優(yōu)中的應(yīng)用研究[J].中國(guó)地質(zhì)教育,2008(4):100-104.
[12]俞立平,武夷山.學(xué)術(shù)期刊客觀賦權(quán)評(píng)價(jià)新方法:指標(biāo)難度賦權(quán)法[J].現(xiàn)代圖書情報(bào)技術(shù),2011,27(4):64-70.
[13]MENDAL Jery M.Tutorial on higher-order statistics spectra in signal processing and system theory:theoretial results and some applications[J].Proe IEEE,1991,49(3):278-305.
[14]PORTNOFF M R.Time-frequency representation of digital signals and system based on short time Fourier analysis[J]. IEEE Trans on ASSP,1980,28:55-69.
[15]MALLAT S.Singularity deteetion and processing with wavelet [J].IEEE Trans Information Theory,1992,38(2):617-693.
[16]MALLAT S G.A theory for multisolution signal decomposition:The wavelet representation[J].IEEE Transaetion on Pattern Analysis and Maehine Intelligence,1989,11(3):674-693.
Storage management method of real-time vibration signal based on entropy weight technology
GUO Cui-juan,GONG Chu-yun,RONG Feng,SONG Ya-qi
(School of Electronics and Information Engineering,Tianjin Polytechnic University,Tianjin 300387,China)
Aiming at the stored problem of the real-time vibration signal,a storage management method of real-time vibration signal was presented based on entropy weight technology.The parameters of signal were calculated from timedomain and frequency-domain respectively,and each part of signal was evaluated with different weight valued parameters,then these signals were ordered,the useful signals with higher weighed result were kept and less useful signals with lower weighed results were deleted.This method can not only save the storage space of database for subsequent signals,but also store useful signals completely with useful information of machinery's states.
vibration signal;big data;entropy weight;data storage
TB535;TP311
A
1671-024X(2015)06-0067-05
10.3969/j.issn.1671-024x.2015.06.014
2015-07-07
國(guó)家自然科學(xué)基金資助項(xiàng)目(61405144)
郭翠娟(1975—),女,博士,副教授,主要研究方向?yàn)楝F(xiàn)代通信網(wǎng)絡(luò).E-mail:guocuijuan@tjpu.edu.cn