• 
    

    
    

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

      多發(fā)送端拓撲結構中網(wǎng)絡鏈路延遲推斷方法

      2010-05-31 06:10:20許鑫何涇沙石恒華
      中南大學學報(自然科學版) 2010年3期
      關鍵詞:多播估計值數(shù)據(jù)流

      許鑫,何涇沙,石恒華

      (1. 北京工業(yè)大學 計算機學院,北京,100124;2. 北京工業(yè)大學 軟件學院,北京,100124)

      網(wǎng)絡透視(Network tomography)[1]作為新興的網(wǎng)絡測量方法已受到很多研究機構[2-7]的關注。其基本思想是利用統(tǒng)計學方法,分析端到端的網(wǎng)絡路徑級屬性信息并推斷網(wǎng)絡鏈路級屬性信息,以更好地檢測網(wǎng)絡行為及診斷網(wǎng)絡的性能變化。Chen等[8]認為這種方法將成為復雜網(wǎng)絡性能診斷和評估中最受關注的方法之一。在對基于網(wǎng)絡透視的鏈路延遲分布推斷研究中,Lo等[9]提出使用端到端多播測量法建立基于邏輯多播樹的離散延遲推斷模型。Shih等[10]提出使用端到端單播測量法建立有限的混合模型,采用最大似然估計(Maximum likelihood estimation, MLE)和期望最大化(Expectation maximization, EM)算法來推斷網(wǎng)絡鏈路延遲情況。Xia等[11]對固定和隨機這2種網(wǎng)絡鏈路延遲情況分別進行了研究,并對隨機情況采用混合指數(shù)模型進行建模。雖然這些方法可以較精確地推斷出網(wǎng)絡鏈路延遲情況,但是,這些都是針對單一發(fā)送端的情況進行研究的,對于復雜拓撲結構的網(wǎng)絡鏈路性能研究還較少。Bu等[12]提出采用最小方差權值平均(Minimum variance weighted average, MVWA)方法對多發(fā)送端拓撲結構的網(wǎng)絡鏈路丟包率和延遲分布進行研究。通過對每一個邏輯多播樹進行單獨推斷,根據(jù)所獲得的探測包數(shù)量進行加權得到融合結果。然而,這種方法沒有充分利用來自各邏輯多播樹的路徑數(shù)據(jù),忽略了其中的關聯(lián)關系,造成最終推斷結果不精確。Rabbat等[13]通過研究網(wǎng)絡拓撲結構變化,探測包到達序列及時間同步等來對復雜拓撲結構的網(wǎng)絡鏈路丟包率進行推斷,該方法的分析過程較復雜且實現(xiàn)較困難。由于多播測量能有效節(jié)省帶寬且易于統(tǒng)計推斷[9],本文作者采用端到端多播測量法進行推斷研究。在盡可能充分利用路徑數(shù)據(jù)的基礎上,提出一種簡單易行的針對多發(fā)送端拓撲結構的網(wǎng)絡鏈路延遲推斷方法。在滿足網(wǎng)絡平穩(wěn)性、網(wǎng)絡鏈路延遲的時間獨立性和空間獨立性的假設下,將復雜的多發(fā)送端拓撲結構的網(wǎng)絡分解成多個分解單元,并根據(jù)各分解單元所含鏈路數(shù)量的升序依次推斷各分解單元中的網(wǎng)絡鏈路延遲分布。

      1 多發(fā)送端拓撲結構分解

      在滿足網(wǎng)絡平穩(wěn)性、網(wǎng)絡鏈路延遲的時間獨立性和空間獨立性的假設下,本文作者提出分解排序法,即提出將多發(fā)送端拓撲結構的網(wǎng)絡分解為多個單發(fā)送端拓撲結構的分解單元,并按照各分解單元所含鏈路數(shù)量的升序進行排序。

      1.1 假設條件

      基于如下假設[2,9-10]對網(wǎng)絡延遲分布進行推斷研究:

      (1) 網(wǎng)絡平穩(wěn)性。每次測量時網(wǎng)絡保持平穩(wěn)。

      (2) 時間獨立性。不同數(shù)據(jù)包在同一鏈路上的延遲統(tǒng)計獨立,并且同分布。

      (3) 空間獨立性。同一數(shù)據(jù)包在不同鏈路上的延遲統(tǒng)計獨立。

      網(wǎng)絡流量突發(fā)或長時間發(fā)生延遲會造成不同數(shù)據(jù)包之間發(fā)生關聯(lián),此時,時間獨立性不成立;經過同一路徑的不同數(shù)據(jù)流會發(fā)生相互影響,此時,空間獨立性不成立。然而,Lo等[9,14]證明了當上述假設條件不能滿足時,對網(wǎng)絡性能推斷的影響程度很小。同時,這些假設又有助于分析網(wǎng)絡屬性和結構信息,因此,本文作者基于上面的假設條件進行推斷研究。

      1.2 分解單元

      選擇邏輯多播樹[9]作為研究多發(fā)送端拓撲結構中的基本分解單元。定義邏輯多播樹為T = (V, L)(其中:V為節(jié)點集;L為鏈路集)。發(fā)送端節(jié)點作為邏輯多播樹的唯一根節(jié)點。節(jié)點對(k, j)表示從節(jié)點 k到節(jié)點 j的鏈路j。這里,鏈路j對應的節(jié)點就是節(jié)點j,節(jié)點k是節(jié)點j的父節(jié)點。

      分解單元如圖1所示,其中:節(jié)點0是發(fā)送端節(jié)點,節(jié)點1是分支節(jié)點,并且是節(jié)點2和節(jié)點3的父節(jié)點;節(jié)點2和節(jié)點3是互為兄弟的葉節(jié)點。發(fā)送端鏈路1與鏈路2和鏈路3分別組成路徑1和路徑2。因此,基于網(wǎng)絡透視的延遲分布推斷模型 Y=AX[2],根據(jù)可測的路徑i上的延遲Yi(i=1, 2)和已知的拓撲結構A,來推斷鏈路j上的延遲Xj(j=1, 2, 3)。

      圖1 分解單元Fig.1 Decomposition unit

      1.3 分解排序法

      設多發(fā)送端拓撲結構網(wǎng)絡中的發(fā)送端節(jié)點個數(shù)為M,接收端節(jié)點個數(shù)為R,節(jié)點集為V,鏈路集為L。該網(wǎng)絡將分解出M個分解單元,每個分解單元都包括唯一發(fā)送端鏈路和多個數(shù)據(jù)流共享鏈路。多發(fā)送端拓撲結構分解過程如圖2所示,具有4個發(fā)送端拓撲結構的網(wǎng)絡將會被分解為4個分解單元,其中:節(jié)點對(S1, 1),(S2, 1),(S3, 1)和(S4, 3)分別為發(fā)送端鏈路 S1,S2,S3和S4。由于發(fā)送端節(jié)點S1,S2和S3均向接收端節(jié)點2,4,5,7和8發(fā)送數(shù)據(jù)流,發(fā)送端節(jié)點S4向接收端節(jié)點 4,5,7和 8發(fā)送數(shù)據(jù)流,因此,鏈路j(j=2, …, 8)為數(shù)據(jù)流共享鏈路。節(jié)點1和節(jié)點3為網(wǎng)絡的數(shù)據(jù)流匯聚節(jié)點。

      設Sm(m=1, …, M)表示發(fā)送端節(jié)點,Rm表示發(fā)送端節(jié)點Sm發(fā)送的數(shù)據(jù)包所到達的接收端節(jié)點個數(shù),且Rm≤R,Tm表示以 Sm為發(fā)送端節(jié)點的分解單元,Vm表示分解單元Tm的節(jié)點集,Lm表示分解單元Tm的鏈路集,并且 Vm?V,Lm?L,Cm表示 Lm中鏈路個數(shù)。分解排序法描述如下:

      步驟1 標記發(fā)送端節(jié)點并確定分解單元個數(shù)M。

      步驟2 劃分出M個分解單元。對于每個發(fā)送端節(jié)點Sm進行如下操作:從Sm出發(fā),尋找由Sm發(fā)送的數(shù)據(jù)包達到對應的所有 Rm個接收端節(jié)點所經過的節(jié)點集Vm和鏈路集Lm,即得到分解單元Tm= (Vm, Lm)。

      步驟3 依照每個分解單元Tm所包含的鏈路個數(shù)Cm的升序方式排列各分解單元,完成分解排序過程。

      由于C4<C1=C2=C3,按照上面的分解排序法,圖2中的分解單元排序為T4T1T2T3。每個數(shù)據(jù)包都標記其發(fā)送端節(jié)點,同時,保證其接收端節(jié)點能確認該數(shù)據(jù)包的發(fā)送端地址;因此,在多發(fā)送端拓撲結構的網(wǎng)絡中,每個分解單元Tm都可以采集到對應Rm條路徑延遲數(shù)據(jù)。通過分析這些路徑延遲數(shù)據(jù)來計算出分解單元Tm的網(wǎng)絡鏈路延遲分布。除對第1個分解單元計算外,其余分解單元的計算都是基于前面計算結果依次進行的。在使用相同測量次數(shù)的延遲數(shù)據(jù)時,分解單元Tm含有鏈路個數(shù)Cm越少,計算結果越精確;因此,采用鏈路個數(shù)較少的分解單元作為最先推斷對象,使得后續(xù)計算能在較好的計算基礎上進行修正。

      2 分解單元的鏈路延遲分布計算方法

      采用離散延遲模型[10]對網(wǎng)絡中的各分解單元進行鏈路延遲分布的計算。在對分解單元中的所有鏈路進行拓撲簡化后,使用最大似然估計法[15-16]對該分解單元中的網(wǎng)絡鏈路延遲分布參數(shù)進行估計。

      2.1 離散延遲模型

      設q為離散延遲模型中的量化單元,主要用來量化所有鏈路 j上的延遲 Xj(j=1, …, J)。當延遲 Xj在內時,則該延遲被量化為dq。其中:量化索引d=0, …, D;D為1個正整數(shù)。當d=0時,量化區(qū)間為 ( 0 , 0.5q]。

      2.2 分解單元的計算方法

      2.2.1 拓撲簡化

      圖2 多發(fā)送端拓撲結構分解Fig.2 Decomposition of multiple-source topology

      為便于對分解單元中的各條鏈路延遲分布進行參數(shù)估計,提出1種拓撲簡化法來分析和處理各鏈路在所屬分解單元中的拓撲結構。即對于某條鏈路,在其所屬分解單元中尋找通過該鏈路的所有路徑,通過進行拓撲簡化使得每條鏈路只屬于該分解單元中的某一條路徑。若僅有1條路徑通過,則不需要進行拓撲簡化處理。

      對于某條鏈路j的拓撲簡化描述如下。

      步驟1 尋找子孫節(jié)點。尋找并標記鏈路j對應節(jié)點的所有子孫節(jié)點。

      步驟 2 自底向上,逐步合并。在這些子孫節(jié)點中,從深度最大(距離鏈路 j所對應節(jié)點的跳數(shù)最多)的葉節(jié)點開始,先對同父的兄弟節(jié)點進行合并,然后,與其父節(jié)點合并。依此類推,直至下一個將合并的節(jié)點為鏈路j對應的節(jié)點為止,此時得到1條路徑。

      以圖2分解單元T1中的鏈路3為例,拓撲簡化過程如圖 3所示。在步驟 1中尋找鏈路 3的子孫節(jié)點,鏈路3對應節(jié)點3,節(jié)點3的所有子孫節(jié)點為節(jié)點4、5,6,7和8。在步驟2中,先合并兄弟葉節(jié)點7和8為節(jié)點(7, 8),然后,與其父節(jié)點6合并為節(jié)點(6, (7, 8)),再與節(jié)點6的兄弟節(jié)點4,5合并為節(jié)點(4, 5, (6, (7, 8))),此時,它們的父節(jié)點為節(jié)點3,是所求鏈路3的對應節(jié)點??梢姡兄雇負浜喕?,得到路徑S1-1-3-(4, 5, (6, (7, 8)))。

      圖3中節(jié)點(4, 5, (6, (7, 8)))對應鏈路的延遲概率可以表達為:

      其中:dj為鏈路j上延遲的量化索引;y1,y2,y3和y4分別為到葉節(jié)點 4,5,7和 8的路徑延遲值,函數(shù)index(x)表示延遲值為x時所對應的量化索引。

      2.2.2 參數(shù)估計

      采用最大似然估計法[15-16]對該分解單元中的網(wǎng)絡鏈路延遲分布進行參數(shù)估計。

      令Y為可測的端到端路徑延遲數(shù)據(jù)。第j條鏈路上延遲Xj的概率為:

      因此,在給定完整延遲數(shù)據(jù)X1, …, XT的情況下,X的對數(shù)似然函數(shù)表示為:

      通過EM算法求解網(wǎng)絡鏈路延遲分布的最大似然估計值。

      令 )(k

      α 表示第k次估計,最大化第k+1次目標方程為:

      計算出第 k+1次估計值 α(k+1),循環(huán)計算,直至得到滿足一定精度的網(wǎng)絡鏈路延遲分布的估計值為止。

      3 多發(fā)送端拓撲結構中網(wǎng)絡鏈路延遲推斷方法

      對于整體的多發(fā)送端拓撲結構,提出的推斷過程描述如下。

      步驟 1 按照分解排序法分解多發(fā)送端拓撲結構的網(wǎng)絡,得到分解單元的排序;

      步驟2 依次選擇分解單元Tm,按照分解單元計算方法對Tm中所有鏈路進行拓撲簡化,使用其對應的Rm條路徑延遲數(shù)據(jù)進行最大似然估計,計算發(fā)送端鏈路和數(shù)據(jù)流共享鏈路的延遲分布,此時,將數(shù)據(jù)流共享鏈路的延遲分布作為計算下一個分解單元延遲分布的初始值,即下一次計算過程是修正前面的計算結果。

      圖3 拓撲簡化Fig.3 Topology simplification

      步驟3 重復步驟2,直到完成對所有分解單元的計算為止,從而得到多發(fā)送端拓撲結構中所有鏈路的延遲分布情況。

      其中,發(fā)送端鏈路只包含于1個分解單元中,因此,只使用所在分解單元的路徑延遲數(shù)據(jù)進行1次計算,而數(shù)據(jù)流共享鏈路可以使用包含它的所有分解單元的路徑延遲數(shù)據(jù)進行多次計算。

      4 仿真實驗

      4.1 實驗平臺

      使用 OPNET仿真軟件進行多發(fā)送端拓撲結構的網(wǎng)絡延遲實驗。網(wǎng)絡拓撲如圖4所示。節(jié)點對(S1, 1),(S2, 1),(S3, 1)和(S4, 3)分別為發(fā)送端鏈路 S1,S2,S3和S4,其余鏈路均為數(shù)據(jù)流共享鏈路。

      4.2 實驗設置

      圖4中各鏈路的數(shù)據(jù)傳送速率如表1所示,各發(fā)送端節(jié)點發(fā)送的數(shù)據(jù)包屬性如表2所示。實驗仿真時間為1 000 s,選擇量化單元q=0.1,共有40個量化區(qū)間。在每條鏈路上采集1 000個點作為該鏈路延遲的真實值,并采用同一量化單元q=0.1進行量化,統(tǒng)計得到各鏈路延遲的真實分布情況,用于與本文作者提出的方法的推斷結果進行對比。

      根據(jù)分解排序法,得到分解單元的推斷順序為T4T1T2T3,即:首先計算分解單元T4得到鏈路S4及數(shù)據(jù)流共享鏈路4,5,6,7和8上的延遲分布估計值;將這些數(shù)據(jù)流共享鏈路上的延遲分布估計值作為計算分解單元T1中對應鏈路上的延遲分布初始值,以得到鏈路S1及數(shù)據(jù)流共享鏈路2,3,4,5,6,7和8上的延遲分布估計值;將這些數(shù)據(jù)流共享鏈路上的延遲分布估計值作為下一個分解單元 T2中對應鏈路上的延遲分布初始值。依此類推,計算得出所有鏈路上的延遲分布估計值。

      圖4 多發(fā)送端的網(wǎng)絡拓撲Fig.4 Multiple-source network topology

      4.3 實驗結果

      4.3.1 有效性驗證

      選擇非數(shù)據(jù)流共享鏈路 S1以及數(shù)據(jù)流共享鏈路2,4和8作為代表鏈路進行分析。圖5所示為這4條鏈路上的延遲分布真實值和估計值的對比結果。

      如圖5(a)所示,對于非數(shù)據(jù)流共享鏈路S1,只使用分解單元T1進行網(wǎng)絡鏈路延遲分布估計,雖然沒有多次估計修正過程,但是網(wǎng)絡鏈路延遲分布估計值已經能較好模擬真實網(wǎng)絡的變化情況。

      設Pd和Qd分別表示當量化索引為d時數(shù)據(jù)流共享鏈路 j上的延遲分布真實值和估計值,根據(jù)計算鏈路j的Pd和Qd之間的差異,其中,數(shù)據(jù)流共享鏈路4和8共進行了4次估計,結果如表3所示。

      表1 各鏈路上數(shù)據(jù)傳送速率Table 1 Data rate of each link

      表2 各發(fā)送端節(jié)點發(fā)送的數(shù)據(jù)包屬性Table 2 Packet properties of each source

      圖5 鏈路延遲分布真實值和估計值的對比Fig.5 Comparisons between true and estimated values of links delay distribution

      對于數(shù)據(jù)流共享鏈路,對真實值與多次修正的估計值進行對比,結果如圖 5(b)~(d)所示。可見:每次修正結果總體上都更趨近于真實的延遲分布情況。鏈路4和鏈路8上的延遲真實值和估計值之間的差異如表3所示,可見:鏈路4和8的每次估計都使得真實值和估計值的差異減小。圖5也顯示出每次修正后的變化趨勢。

      表3 鏈路4和鏈路8上的延遲真實值和估計值之間的差異Table 3 Difference between true and estimated values of delay for Link 4 and Link 8

      在對多發(fā)送端拓撲結構的網(wǎng)絡鏈路延遲分布進行推斷研究時,要充分考慮所有分解單元中的數(shù)據(jù)流對網(wǎng)絡性能的影響。從圖4可見:該網(wǎng)絡中的4個發(fā)送端節(jié)點所產生的數(shù)據(jù)流共同影響其數(shù)據(jù)流共享鏈路,因此,只使用其中某個分解單元的路徑數(shù)據(jù)進行分析,則推斷結果只是受該分解單元中數(shù)據(jù)流影響的粗略估計值。但是,若在滿足網(wǎng)絡平穩(wěn)性、網(wǎng)絡鏈路延遲時間獨立性和空間獨立性的假設條件下,使用其余分解單元的路徑延遲數(shù)據(jù)對該分解單元的粗略估計值進行多次迭代修正,則會使推斷結果更符合網(wǎng)絡內部的真實情況。

      4.3.2 與MVWA方法的比較

      采用本文作者提出的方法與 MVWA方法[12]所得的鏈路延遲分布比較結果如圖6所示。鏈路延遲分布概率越大,精度越高。從圖6可見:作者提出的方法優(yōu)于MVWA方法。由于MVWA方法只是通過對每個邏輯多播樹進行單獨計算,雖然考慮了所有數(shù)據(jù)流的影響,但是,這種較簡單的通過加權得到估計值的方法卻降低了精度。而作者提出的方法在沒有增加計算復雜度的情況下,充分利用了路徑數(shù)據(jù),通過多次估計修正達到較高的精度。

      圖6 本文方法與MVWA方法鏈路延遲分布概率的比較Fig.6 Comparisons of probability of links delay distribution between method in this paper and MVWA method

      5 結論

      (1) 提出一種網(wǎng)絡鏈路延遲分布推斷方法來解決復雜的多發(fā)送端拓撲結構中網(wǎng)絡內部延遲推斷問題。該方法將復雜的多發(fā)送端拓撲結構的網(wǎng)絡分解成多個簡單的分解單元,并按照各分解單元所含鏈路個數(shù)的升序計算各分解單元網(wǎng)絡鏈路延遲分布。

      (2) 每個分解單元的計算是在完成拓撲簡化處理后,使用最大似然估計法進行網(wǎng)絡鏈路延遲分布的估計。由于對每個分解單元的數(shù)據(jù)流共享鏈路延遲分布的計算,都是基于前面估計結果進行修正,因此,數(shù)據(jù)流共享鏈路延遲分布的真實值和估計值之間的差異逐漸減小,即推斷結果逐漸趨近于網(wǎng)絡內部的真實情況。

      (3) 與MVWA方法相比,本文作者所提出的方法能更精確地推斷出多發(fā)送端拓撲結構中網(wǎng)絡內部鏈路延遲情況,具有較強的實用性。

      [1] Vardi Y. Network tomography: estimating source-destination traffic intensities from link data[J]. American Statistical Association, 1996, 91: 365-377.

      [2] Castro R, Coates M J, LIANG Gang, et al. Network tomography:recent developments[J]. Statistical Science, 2004, 52(3):499-517.

      [3] Duffield N G. Network tomography of binary network performance characteristics[J]. IEEE Transactions on Information Theory, 2006, 52(12): 5373-5388.

      [4] Shih M, Hero A O. Hierarchical inference of unicast network topologies based on end-to-end measurements[J]. IEEE Transactions on Signal Processing, 2007, 55(5): 1708-1718.

      [5] JIN Xing, TU Wang-qing, Gary Chan S H. Scalable and efficient end-to-end network topology inference[J]. IEEE Transaction on Parallel and Distributed System, 2008, 19(6): 837-850.

      [6] Duffield N G, Lo Presti F, Paxson V, et al. Network loss tomography using striped unicast probes[J]. IEEE/ACM Transaction on Networking, 2006, 14(4): 691-710.

      [7] Sharma G, Jaggi S, Dey B K. Network tomography via network coding[C]//Masimo F. IEEE Information Theory and Applications Workshop. Piscataway: IEEE, 2008: 151-157.

      [8] CHEN Ai-you, CAO Jin, BU Tian. Network tomography:identifiability and Fourier domain estimation[C]//IEEE INFOCOM. Anchorage: IEEE, 2007: 1875-1883.

      [9] Lo P F, Duffield N G, Horowitz J, et al. Multicast-based inference of network-internal delay distributions[J]. IEEE/ACM Transactions on Networking, 2002, 10(6): 761-775.

      [10] Shih M, Hero A O. Unicast-based inference of network link delay distributions with finite mixture models[J]. IEEE Transactions on Signal Processing, 2003, 51(8): 2219-2228.

      [11] Xia Y, Tse D. Inference of link delay in communication networks[J]. IEEE Journal on Selected Areas in Communications,2006, 24(12): 2235-2248.

      [12] Bu T, Duffield N G, Lo Presti F, et al. Network tomography on general topologies[C]//Richard M. Proceeding of the ACM SIGMETRICS. Marina Del Rey: ACM Press, 2002: 21-30.

      [13] Rabbat M G, Coates M J, Nowak R D. Multiple-source Internet tomography[J]. IEEE Journal on Selected Areas in Communications, 2006, 24(12): 2221-2234.

      [14] Caceres R, Duffield N G, Horowitz J, et al. Multicast-based inference of network-internal loss characteristics[J]. IEEE Transactions on Information Theory, 1999, 45(7): 2462-2480.

      [15] LIANG Gang, YU Bin. Maximum pseudo likelihood estimation in network tomography[J]. IEEE Transactions on Signal Processing, 2003, 51(8): 2043-2053.

      [16] Tsang Y, Nowak R. Network tomography in theory and practice[D]. Houston: Rice University. Department of Electrical and Computer Engineering, 2005.

      猜你喜歡
      多播估計值數(shù)據(jù)流
      胖樹拓撲中高效實用的定制多播路由算法
      用于超大Infiniband網(wǎng)絡的負載均衡多播路由
      InfiniBand中面向有限多播表條目數(shù)的多播路由算法
      汽車維修數(shù)據(jù)流基礎(下)
      一道樣本的數(shù)字特征與頻率分布直方圖的交匯問題
      統(tǒng)計信息
      2018年4月世界粗鋼產量表(續(xù))萬噸
      一種提高TCP與UDP數(shù)據(jù)流公平性的擁塞控制機制
      基于數(shù)據(jù)流聚類的多目標跟蹤算法
      北醫(yī)三院 數(shù)據(jù)流疏通就診量
      洪洞县| 屏山县| 泗洪县| 平阳县| 玉林市| 阜新市| 万宁市| 高尔夫| 汾阳市| 镶黄旗| 福清市| 临湘市| 宽城| 米易县| 犍为县| 阳东县| 囊谦县| 布尔津县| 精河县| 吴忠市| 汉阴县| 色达县| 浏阳市| 新龙县| 辉县市| 桐柏县| 确山县| 枝江市| 项城市| 清丰县| 溧阳市| 聂拉木县| 民勤县| 五莲县| 广饶县| 皮山县| 调兵山市| 略阳县| 锡林浩特市| 都江堰市| 库伦旗|