• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      面向智能電網(wǎng)的無線傳感器網(wǎng)絡(luò)目標探測算法

      2015-03-16 06:03:34王曉玲王昕彬宋志明王紅梅辛俊峰
      湖南電力 2015年2期
      關(guān)鍵詞:向量無線矩陣

      王曉玲,王昕彬,宋志明,王紅梅,辛俊峰

      (1.國網(wǎng)技術(shù)學(xué)院,山東 泰安 271000;2.江蘇電力公司南京供電公司,江蘇 南京 210000;3.山東電力公司濟寧供電公司,山東 濟寧 272000)

      面向智能電網(wǎng)的無線傳感器網(wǎng)絡(luò)目標探測算法

      王曉玲1,王昕彬2,宋志明1,王紅梅3,辛俊峰1

      (1.國網(wǎng)技術(shù)學(xué)院,山東 泰安 271000;2.江蘇電力公司南京供電公司,江蘇 南京 210000;3.山東電力公司濟寧供電公司,山東 濟寧 272000)

      文中提出一種結(jié)合壓縮感知 (CS)路由算法進行目標源定位的方法。該路由算法通過對網(wǎng)絡(luò)節(jié)點進行分簇,將簇內(nèi)節(jié)點的信息集中在簇頭上。然后對無線傳感器網(wǎng)絡(luò)進行拓撲優(yōu)化,并采用多跳路由策略完成簇頭節(jié)點與匯聚節(jié)點間的通信。最后利用次梯度投影算法,從來自簇頭的少量數(shù)據(jù)中恢復(fù)匯聚節(jié)點信號源,用于目標源信號的定位。采用一種自適應(yīng)處理機制優(yōu)化壓縮感知算法的信號重構(gòu)性能。仿真結(jié)果表明,算法能準確探測目標,減少計算量,并具有良好的重構(gòu)性能。

      智能電網(wǎng);壓縮感知;次梯度投影;目標探測;無線傳感器網(wǎng)絡(luò)

      當前,電力行業(yè)的發(fā)展面臨投資效率低、供電可靠性不高、靈活應(yīng)變能力不強、環(huán)境破壞嚴重等諸多問題的挑戰(zhàn)〔1〕?;谌录夹g(shù)和架構(gòu)的智能電網(wǎng) (smart grid)〔2-3〕是當前各國公認的最佳解決途徑。智能電網(wǎng)是將先進的傳感量測技術(shù)、信息通訊技術(shù)、分析策略技術(shù)、自動控制技術(shù)、能源電力技術(shù)以及電網(wǎng)技術(shù)設(shè)施高度集中而形成的新型現(xiàn)代化電網(wǎng)。先進的無線傳感器網(wǎng)絡(luò) (Wireless Sensor Network,WSN)和通信技術(shù)是其中的基礎(chǔ)和關(guān)鍵技術(shù)。因此,要實現(xiàn)智能電網(wǎng)開放、互動的目標,必須具備一個可靠、高效的通信網(wǎng)絡(luò)。智能電網(wǎng)涵蓋了電力傳輸?shù)母鱾€環(huán)節(jié),但是發(fā)電、輸電、變電、配電和電力用戶等不同主體之間的信息傳輸、各種應(yīng)用系統(tǒng)間的數(shù)據(jù)傳輸,都必須依賴通信載體。無線傳感器網(wǎng)絡(luò)作為通信技術(shù),具有大規(guī)模性、自組織性、動態(tài)性、可靠性、多跳路由等特點〔4〕。無線傳感器網(wǎng)絡(luò)通過分布式處理大量的采集信息能夠提高監(jiān)測的精確度;大量節(jié)點能夠增大覆蓋的監(jiān)測區(qū)域,便于遙控遙測;大量冗余節(jié)點的存在,使得系統(tǒng)具有很強的容錯性能。無線傳感器網(wǎng)絡(luò)在多目標、短距離通信領(lǐng)域被廣泛應(yīng)用。因此,該項技術(shù)對于建設(shè)智能電網(wǎng)有著非常重要的作用。

      但是,隨著無線通信技術(shù)的發(fā)展和網(wǎng)絡(luò)規(guī)模的擴張,如何實現(xiàn)對網(wǎng)絡(luò)中大規(guī)模異質(zhì)數(shù)據(jù)的有效采集,降低信息存儲、處理和傳輸?shù)某杀?,成為亟待解決的焦點問題。壓縮感知理論〔5-7〕在無線傳感器網(wǎng)絡(luò)中的應(yīng)用研究能促進相關(guān)技術(shù)的發(fā)展和壓縮感知理論自身的完善。需要解決的問題主要有2個方面:首先,伴隨著網(wǎng)絡(luò)發(fā)展,急速擴充的網(wǎng)絡(luò)規(guī)模給數(shù)據(jù)的處理、傳輸和存儲帶來了巨大的壓力。這點和壓縮感知減低數(shù)據(jù)量的設(shè)計初衷相吻合。其次,壓縮感知理論還有很多值得研究的地方,尤其是無線傳感器網(wǎng)絡(luò)方面,包括算法效率、魯棒性、可靠性等方面都有待提高。

      目前壓縮感知理論己廣泛應(yīng)用于無線通信網(wǎng)絡(luò)的物理層、網(wǎng)絡(luò)層、應(yīng)用層、數(shù)據(jù)鏈路層等各個層面〔8〕, 為解決信號檢測〔9-12〕、 信道估計〔13〕、 成簇算法設(shè)計、 網(wǎng)絡(luò)數(shù)據(jù)處理〔14-15〕、 網(wǎng)絡(luò)監(jiān)管〔16-17〕等方面的問題提供了有效的途徑。在大規(guī)模分布式的無線傳感器網(wǎng)絡(luò)中,所需處理的數(shù)據(jù)量非常大,但是大部分數(shù)據(jù)源產(chǎn)生的信息是一樣的,故該信息具有稀疏性。然而,壓縮感知理論的有效應(yīng)用也面臨2個方面的挑戰(zhàn):如何得到網(wǎng)絡(luò)數(shù)據(jù)的稀疏化表示;如何得到有效的觀測值。一般解決的思路是將壓縮感知理論和網(wǎng)絡(luò)拓撲及路由相結(jié)合,通過匯聚各個隨機節(jié)點的數(shù)據(jù)信息,在匯聚節(jié)點處得到觀測值集合。而壓縮感知理論中的觀測矩陣可用加權(quán)信息表示。所以,經(jīng)過路由構(gòu)建的觀測矩陣本身是稀疏的,因此利用壓縮感知理論重構(gòu)序列可有效提高信號的傳輸效率。

      文中使用路由算法首先通過對網(wǎng)絡(luò)節(jié)點進行分簇,簇的信息集中在一個簇頭節(jié)點。然后,通過選擇簇頭并將收集到的信息進行傳輸,可大大減少傳輸?shù)絽R聚節(jié)點的數(shù)據(jù)量。最后,基于次梯度投影算法,匯聚節(jié)點可將來自信號源的少量簇頭數(shù)據(jù)進行目標源恢復(fù)。由于采用了一種自適應(yīng)機制的算法,優(yōu)化了壓縮感知重構(gòu)算法的效率。仿真結(jié)果表明,該算法能夠準確地檢測目標,體現(xiàn)了該算法在定位精度和抗干擾性上的優(yōu)勢。

      1 壓縮重構(gòu)理論

      在CS理論中,如果一個信號的投影域和變換域的基是不相干的,而且該信號經(jīng)過域變換后是稀疏的,則該信號可以通過其在投影域的投影近似無損的獲得〔1〕。

      假設(shè)有離散信號X=[x(1),x(2),…,x(N)]T,其采樣點數(shù)為 N。變換矩陣 Ψ=[φ1,φ2,…,φN],其列向量φi為N×1的向量,且相互正交。此時,信號X為:

      式中 θi為列向量φi的對應(yīng)系數(shù)。當‖Θ‖0=K時,則稱信號X為K-稀疏的。由于Ψ為正交矩陣,則有Θ=ΨTX。此時,存在一個平穩(wěn)的、與變換基Ψ不相關(guān)的M×N維觀測矩陣Φ,使得觀測集合Y=ΦΘ=ΦΨTX成立。針對此等式,如果想求解X,那么可利用0-范數(shù)意義下的優(yōu)化問題求解其精確值或近似逼近值X^。

      考慮到0-范數(shù)的優(yōu)化問題是個NP問題,因此常將式 (2)轉(zhuǎn)換為1-范數(shù)的優(yōu)化問題〔2〕,即:

      令s=ΨTX,則式 (3)改為

      上述問題是一個典型的非線性凸優(yōu)化問題,通過求解該最優(yōu)化問題,可得到變換域上的s,再經(jīng)反變換就可得到原始信號估計X^。

      2 基于壓縮感知的數(shù)據(jù)采集方法

      在包含目標源的傳感器網(wǎng)絡(luò)中,網(wǎng)絡(luò)節(jié)點通過與全體傳感器節(jié)點的信息交互,可以獲得對目標源的探測。但是,受傳感器節(jié)點的尺寸、功耗等限制,節(jié)點的處理能力不足以獲取全部節(jié)點的數(shù)據(jù)信息。為了通過少量節(jié)點獲得整個網(wǎng)絡(luò)數(shù)據(jù),首先進行節(jié)點融合。假設(shè)無線信號傳輸過程中受到干擾,干擾源為高斯白噪聲,因此匯聚節(jié)點的信息收集模型的數(shù)學(xué)形式為〔18〕:

      式 (5)中 y為匯聚節(jié)點的接收向量為M維匯聚節(jié)點信息向量;n為高斯白噪聲向量,n∈RM;v為目標源向量;Φ為觀測矩陣;Ψ為稀疏基矩陣。

      目標源數(shù)目K相比子區(qū)域的數(shù)量N幾乎可以忽略。目標源向量v中非零元素的數(shù)量同維數(shù)相比幾乎為0,即目標源向量v是稀疏的。故由式 (5)可得,在無線傳感器網(wǎng)絡(luò)中,目標源探測問題實際上就是從信息向量 y中恢復(fù)出目標源向量 v的問題。

      按照壓縮感知理論可知,式 (5)符合壓縮感知理論的采樣觀測模型,即y相當于觀測值向量,矩陣Φ為觀測矩陣,矩陣Ψ為稀疏基矩陣,目標源向量v是稀疏系數(shù)向量,高斯白噪聲n為觀測噪聲。由于,y的維數(shù)M遠小于v的維數(shù)N,而且v具有稀疏性,根據(jù)壓縮感知理論模型,v的求解問題轉(zhuǎn)變?yōu)槿绾螛?gòu)建觀測矩陣Φ,即如何對傳感器節(jié)點進行信息融合操作的問題。

      通過采集各個簇點的信息,匯聚節(jié)點對信息進行匯總,從而判斷目標源參數(shù)信息。此時,信息向量y可表示為:

      式 (6)中 yi為匯聚節(jié)點所接受到的來自第i個簇的簇頭節(jié)點的信息;ki為第i個簇所包含的節(jié)點數(shù)量,包括簇頭節(jié)點和普通非簇頭節(jié)點;wi,k為第i個簇中的第k個節(jié)點的信息權(quán)值系數(shù);xi,k為第i個簇中的第k個節(jié)點上的本地信息;ni為匯聚節(jié)點所接收到的來自第i個簇的信息所受到的噪聲干擾,假設(shè)ni滿足零均值高斯分布N(0,β-1)。

      構(gòu)建權(quán)值矩陣W,令wi,k為權(quán)值矩陣W的第i行第k列元素。當且僅當?shù)趇行第k列元素所對應(yīng)的傳感器節(jié)點屬于第i個簇的范圍時,wi,k的值才為1;否則,取值為0。此時,整個無線傳感器網(wǎng)絡(luò)的所有簇集合表示為:

      此處,矩陣W大小為M×N,M為無線傳感器網(wǎng)絡(luò)中簇數(shù),n為網(wǎng)絡(luò)中傳感器節(jié)點數(shù)目。

      對比可知,權(quán)值矩陣W可以看作是觀測矩陣Φ,此時觀測值數(shù)量也就是簇的數(shù)量。在傳感器網(wǎng)絡(luò)中,簇頭節(jié)點的作用大,負責采集信息給匯聚節(jié)點。簇頭節(jié)點中包含整個網(wǎng)絡(luò)的數(shù)據(jù),從簇頭節(jié)點可以恢復(fù)整個網(wǎng)絡(luò)。但是,在傳感器網(wǎng)絡(luò)中,簇頭節(jié)點的數(shù)量同網(wǎng)絡(luò)中的節(jié)點總數(shù)相比可以忽略不計。因此,網(wǎng)絡(luò)傳感器節(jié)點信息的恢復(fù)算法通過普通算法是無法求解的,這是一個欠定問題。求解這個問題的算法比較多,其中壓縮感知算法就是一個較好的方法。該算法可利用原始信號的稀疏性,在欠定情況下,恢復(fù)原始信號。

      3 基于自適應(yīng)次梯度投影算法的數(shù)據(jù)重構(gòu)方法

      假設(shè)不同目標信號源的發(fā)射信號強度相同,則目標源向量中元素的值可以設(shè)置為1或者0。為了利用上述目標信號源向量v的相關(guān)先驗信息,文中采用自適應(yīng)次梯度投影重構(gòu)算法〔19〕。其具體步驟為:

      第1步,令Φ=W,s=Ψv,Y=y(tǒng),則式 (7)變?yōu)?/p>

      第2步,求解式 (8),可轉(zhuǎn)換為

      上式是一個最優(yōu)化求解的問題,求解這個問題可以得到變量s,從而通過反變換得到原始信號x。

      第3步,引入向量u和v代替s,將式 (9)做正則化處理這里, ui=(si)+, vi=(-si)+, i=1, 2, …,n,操作符(·)+表示取向量中的正數(shù)部分,負數(shù)部分置零。 所以, ‖s‖1=1Tnu+1Tnv, 這里 1Tn=[1,1,…,1]T是長度為n的1向量。式 (9)可改寫為

      第4步,將式 (11)轉(zhuǎn)化為標準的受限二次規(guī)劃形式

      這里

      第5步,根據(jù)式 (12),定義凸集 {Ck,k=1,...,m}滿足約束條件 (14)

      此時,z的求解問題就等同于向凸集Ck的投影問題,當滿足z∈Ιk∈ICk≠φ時,該求解問題能保證收斂,即

      第6步,定義凸函數(shù)

      以及梯度算子

      由上可知,膨脹系數(shù)ρ的取值直接影響系統(tǒng)性能,其取值和F(z)、和信號的稀疏度存在直接的關(guān)聯(lián)。這些參數(shù)的誤差決定了系統(tǒng)的性能,因此,要想提高系統(tǒng)性能,必須將ρ的取值與這些參數(shù)進行關(guān)聯(lián),自適應(yīng)的取值。由式 (10)可知,ρ的取值和凸集Ck(ρ)面積的取值成正比。而凸集Ck(ρ)面積的取值與迭代的次數(shù)成反比。但是迭代次數(shù)小并不意味系統(tǒng)的性能就好,因此當凸集Ck(ρ)面積大時,∈Ck(ρ)的概率增大,從而使收斂容易,但是也會引入較大的誤差,從而影響系統(tǒng)的性能。因此,如何綜合系統(tǒng)性能和收斂速度間的關(guān)系是一個非常關(guān)鍵的因素。

      為此,文中采用了一種自適應(yīng)的調(diào)節(jié)算法。算法的思路是:初始時,ρstart取一個相對大的值,并且?Ck(ρ),然后通過設(shè)置迭代使k∈Ck(ρ)。然后以步長Δ逐步減小ρ,使得Ck(ρ)面積不斷減小,每次迭代的誤差逐步減少。最后,通過設(shè)置最終膨脹系數(shù)值 ρstop, 求得最優(yōu)解 zopt。此時,迭代結(jié)束。否則,當k=M(M為終止迭代次數(shù))時,迭代過程也結(jié)束。

      4 仿真與驗證

      仿真實驗的設(shè)置如下:傳感器節(jié)點個數(shù)為100,目標源的個數(shù)為K,分布區(qū)域的大小為100× 100。在分布的區(qū)域里,所有節(jié)點都是隨機分布,彼此通信,匯聚節(jié)點位于整個區(qū)域的中心。目標源的信號強度均為1。系統(tǒng)干擾為均值為0,方差為0.005的高斯白噪聲。為了評價目標源探測的性能,引入定位誤差ε,定義為:

      式中 vi和 v′i分別代表原始目標和估計目標。D(g)代表歐式距離。

      不同算法下目標源探測性能對比結(jié)果如圖1所示。比較的算法包括自適應(yīng)次梯度投影算法、基本追蹤算法和傳統(tǒng)算法〔20〕。由圖1可見,對于傳統(tǒng)算法來說,要想實現(xiàn)重構(gòu)恢復(fù),所需的節(jié)點簇數(shù)幾乎與節(jié)點總數(shù)相等。而CS算法只需少量簇就能精確探測目標源,因此CS算法的應(yīng)用具有一定的必要性。在兩種CS算法比較中,自適應(yīng)次梯度算法的性能更佳。

      圖1 CS算法和傳統(tǒng)算法的性能對比

      信號源對文中所提出的算法性能的影響見圖2,圖2中M代表信號源個數(shù)。在同樣的簇數(shù)下,算法的重構(gòu)性能隨著目標源數(shù)的增加而降低。因此,目標源越多,就需要越多的簇來實現(xiàn)目標源的精確探測。

      圖2 不同信號源數(shù)下的算法性能對比

      5 結(jié)束語

      針對智能電網(wǎng)中無線傳感器網(wǎng)絡(luò)應(yīng)用存在的問題,文中提出了一種基于壓縮感知的無線傳感器網(wǎng)絡(luò)進行目標源探測的算法。算法利用自適應(yīng)次梯度算法從來自簇頭的少量數(shù)據(jù)中恢復(fù)出目標源。該算法能有效地減少網(wǎng)絡(luò)中待存儲和傳輸?shù)臄?shù)據(jù)量,從而給無線通信網(wǎng)絡(luò)的發(fā)展提供了有效地保障。

      仿真實驗結(jié)果表明,文中算法同傳統(tǒng)算法相比,能更精確恢復(fù)目標源。但是,算法性能仍存在不足,后繼研究應(yīng)主要面對壓縮感知算法進行改進。研究內(nèi)容包括觀測矩陣的稀疏化和重構(gòu)算法,并推動壓縮感知理論在無線通信網(wǎng)絡(luò)中的應(yīng)用研究,提高算法的運算效率。

      〔1〕Smart Grid Working Group.Challenge and opportunity:charting a new energy future,Appendix A:working group reports〔R〕. Washington,DC,USA:Energy Future Coalition,2003.

      〔2〕Commission E.European SmartGrids technology platform:vision and strategy for europe’ s electricity networks of the future〔R〕. European Commission, Brussels, Belgium, Tech. Rep. EUR22040,2006.

      〔3〕Gungor V C,Lambert F C.A survey on communication networks for electric system automation〔J〕.Elsevier Computer Networks,2006,50(7):877-897.

      〔4〕任豐原,黃海寧,林闖.無線傳感器網(wǎng)絡(luò) 〔J〕.軟件學(xué)報,2003,14(7):1 282-1 291.

      〔5〕Candes E J,Romberg J.Quantitative robust uncertainty principles and optimally sparse decompositions〔J〕.Foundations of Computational Mathematics,2006,6(2):227-254.

      〔6〕Candes E J,Tao T.Near-Optimal Signal Recovery From Random Projections:Universal Encoding Strategies?〔J〕.IEEE Transactions on Information Theory,2006,52(12):5 406-5 425.

      〔7〕 Donoho D L.Compressed sensing〔J〕.IEEE Transactions on Information Theory,2006,52(4):1 289-1 306.

      〔8〕王曉.壓縮感知在無線通信網(wǎng)絡(luò)數(shù)據(jù)采集中的應(yīng)用 〔D〕. 2011,浙江大學(xué).

      〔9〕Jia Y, Pengwu T, Hongyi Y.The identification of frequency hopping signal using compressive sensing〔J〕.Communications and Network,2009,1(1):52-56.

      〔10〕Fletcher A K,Rangan S,Goyal V K.On-off random access channels:A compressed sensing framework〔J〕.Arxiv preprint ArXiv:0903.1022,2009.

      〔11〕Al-Safadi E B,Al-Naffouri T Y.On reducing the complexity of tone-reservation based PAPR reduction schemes by compressive sensing〔C〕.//IEEE Conference on Global Telecommunications,2009:1-6.

      〔12〕Duarte M F,Sarvotham S,Baron D,et al.Distributed compressed sensing of jointly sparse signals〔C〕.Asilomar Conf.Signals,Sys.,Comput.2005:1 537-1 541.

      〔13〕Bajwa W U,Sayeed A,Nowak R.Compressed sensing of wireless channels in time,frequency,and space〔C〕.IEEE Conference on Signals,Systems and Computers,2008:2 048-2 052.

      〔14〕Berger C R,Zhou S,Preisig J C,Willett P.Sparse channel estimation for multicarrier underwater acoustic communication: From subspace methods tocompressed sensing〔J〕.IEEE Transactions on Signal Processing,2010,58(3):1 708-1721.

      〔15〕Haupt J,Nowak R,Castro R.Adaptive sensing for sparse signal recovery〔C〕.IEEE 13th Digital Signal Processing Workshop and 5th IEEE SignalProcessing Education Workshop, 2009: 702-707.

      〔16〕Coates M,Pointurier Y,Rabbat M.Compressed network monitoring〔C〕.IEEE/SP 14th Workshop on Statistical Signal Processing,2007,SSP'07:418-422.

      〔17〕CoatesM, PointurierY, RabbatM.Compressed network monitoring for IP and all-optical networks.ACM Proceedings of the 7th ACM SIGCOMM conference on Internet measurement,2007:241-252.

      〔18〕Heinzelman W R,Chandrakasan A,Balakrishnan H.Energyefficient communication protocol for wireless microsensor networks〔C〕.IEEE Conference on Proceedings of the 33rd Annual Hawaii International,2000:10.

      〔19〕梁瑞宇,鄒采榮,王青云,等.基于自適應(yīng)次梯度投影算法的壓縮感知信號重構(gòu) 〔J〕.信號處理,2010,26(12):1 883-1 889.

      〔20〕李志剛,屈玉貴,劉桂英.用無線傳感器網(wǎng)絡(luò)探測跟蹤目標〔J〕.信息安全與通信保密,2006(12):72-74.

      Target detection algorithm of wireless sensor network in smart grid

      WANG Xiao-ling1,WANG Xin-bin2,SONG Zhi-ming1,WANG Hong-mei3,XIN Jun-feng1
      (1.State Grid of China Technology College,Taian 271000,China;2.Jiangsu Electric Power Company Nanjing Power Supply Company,Nanjing 210000,China;3.Shandong Electric Power Company Jining Power Supply Company,Nanjing 272000,China)

      In order to effectively detect the signal target in the monitoring area of the wireless sensor network,a routine algorithm combining with the compressive sensing(CS)is proposed in this paper,in which network clustering routing algorithm and CS are integrated.In this paper,the network nodes are clustered by the routine algorithm,and information within the cluster nodes is concentrated in cluster-head.After the optimization of wireless sensor network topology,communication between the cluster head node and the sink node is realized by multi-hop routing.Then subgradient projection method is used to restore source signal from few cluster head data.To improve the signal reconstruction performance of CS,a new adaptive processing mechanism is also proposed in this paper.Simulations demonstrate that the proposed algorithm which has excellent good reconstruction performance can accurately detect the target and the computational complexity is reduced effectively.

      smart grid;compression sensing;subgradient projection;target detection;wireless sensor network

      10.3969/j.issn.1008-0198.2015.02.002

      TM727;TN929.5;TP212.6

      A

      1008-0198(2015)02-0005-05

      王曉玲(1969),女,山東魚臺人,國網(wǎng)技術(shù)學(xué)院高級講師,學(xué)士,主要研究方向為智能電網(wǎng)、變電運維等。

      2014-07-07 改回日期:2015-01-15

      猜你喜歡
      向量無線矩陣
      向量的分解
      《無線互聯(lián)科技》征稿詞(2021)
      聚焦“向量與三角”創(chuàng)新題
      無線追蹤3
      基于ARM的無線WiFi插排的設(shè)計
      電子制作(2018年23期)2018-12-26 01:01:08
      ADF7021-N在無線尋呼發(fā)射系統(tǒng)中的應(yīng)用
      電子制作(2016年15期)2017-01-15 13:39:03
      初等行變換與初等列變換并用求逆矩陣
      向量垂直在解析幾何中的應(yīng)用
      向量五種“變身” 玩轉(zhuǎn)圓錐曲線
      矩陣
      南都周刊(2015年4期)2015-09-10 07:22:44
      赣榆县| 安仁县| 寿宁县| 澜沧| 綦江县| 保德县| 乌拉特中旗| 长沙县| 搜索| 如皋市| 桐柏县| 滕州市| 岳阳县| 中阳县| 内丘县| 临沂市| 剑川县| 鞍山市| 盱眙县| 长寿区| 东宁县| 怀安县| 乐清市| 皋兰县| 阜平县| 石棉县| 昌平区| 三门县| 南涧| 洛宁县| 五家渠市| 定远县| 会昌县| 安龙县| 淄博市| 舞阳县| 昌图县| 谢通门县| 北宁市| 五华县| 四会市|