• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于空間插值算法的2014—2018年河南省PM2.5時(shí)空分布及演變特征研究

    2019-09-10 04:25:34張馨文
    河南科技 2019年8期
    關(guān)鍵詞:克里插值河南省

    張馨文

    摘 要:本文選取河南省2014—2018年監(jiān)測(cè)網(wǎng)絡(luò)中269個(gè)監(jiān)測(cè)站點(diǎn)的空氣質(zhì)量資料,在地理統(tǒng)計(jì)學(xué)和GIS的支持下,分別用四種不同的空間插值方法對(duì)各年P(guān)M2.5年均濃度進(jìn)行插值,通過(guò)模擬精度的對(duì)比以及制圖分析發(fā)現(xiàn):泛克里金和普通克里金插值精度較高;而可視化方面普通克里金具有最佳效果。將普通克里金插值結(jié)果與前期探索性數(shù)據(jù)進(jìn)行整合,得出結(jié)論:PM2.5濃度大致從由中北部向四周遞減變?yōu)橛杀辈肯蚰喜窟f減,在東西方向上呈現(xiàn)中部高兩側(cè)低的特點(diǎn),2014年以來(lái)河南省大氣質(zhì)量整體在改善。

    關(guān)鍵詞:PM2.5;空間插值;時(shí)空分布

    Abstract: In this paper, selected the air quality data of 269 monitoring stations in the monitoring network of Henan Province from 2014 to 2018, and the annual average concentration of PM2.5 was interpolated by four different spatial interpolation methods with the support of geostatistics and GIS. Through the comparison of simulation precision and cartographic analysis, it was found that Ordinary Kriging and Universal Kriging have the highest accuracy, while the Ordinary Kriging has the best effect in the visual. Combining the Ordinary Kriging interpolation with the ESDA, the result showed that the concentration of PM2.5 is changing from the pattern of decline from north-central to the sidesinto the north to the south, whilein the east-west direction the middle is higher than the sides, and the overall air quality has improved in Henan Province since 2014.

    Keywords: PM2.5;spatial interpolation;spatial and temporal distribution

    空氣質(zhì)量與每個(gè)人的生活息息相關(guān)。但是,隨著我國(guó)城市化和工業(yè)化進(jìn)程不斷加快,大范圍霧霾天氣頻頻發(fā)生,嚴(yán)重影響了居民的健康水平、經(jīng)濟(jì)與環(huán)境的可持續(xù)發(fā)展[1]。作為重點(diǎn)大氣污染源[2],PM2.5一方面粒徑細(xì)小,能較長(zhǎng)時(shí)間懸浮于空氣中,顯著減弱大氣能見(jiàn)度;另一方面,其化學(xué)組成復(fù)雜,富含大量有害物質(zhì),極易損害人體健康。黨的十九大報(bào)告中指出“污染防治”是我國(guó)全面建成小康社會(huì)實(shí)現(xiàn)第一個(gè)一百年奮斗目標(biāo)必須打好的“三大攻堅(jiān)戰(zhàn)”之一,以防治PM2.5為重點(diǎn)的大氣污染防治工作在我國(guó)必然被著力推進(jìn)。然而,近年來(lái),河南省空氣質(zhì)量形勢(shì)不容樂(lè)觀,尤其到18年末,多地持續(xù)出現(xiàn)霧霾紅色預(yù)警,17個(gè)地級(jí)市中就有12個(gè)頻頻登上全國(guó)城市空氣質(zhì)量排行的末位。由此,結(jié)合區(qū)域特征,利用更適宜的方法探討河南省PM2.5的時(shí)空分布規(guī)律,成為優(yōu)化大氣防治措施的關(guān)鍵。

    國(guó)內(nèi)有關(guān)PM2.5時(shí)空分布的研究區(qū)域主體往往以全國(guó)為重點(diǎn)[3-4],而區(qū)域尺度上,多數(shù)研究基本關(guān)注京津冀、長(zhǎng)三角等污染核心地帶[5-6],尚未發(fā)現(xiàn)以河南省為區(qū)域主體的研究。同時(shí),研究數(shù)據(jù)主要來(lái)自遙感反演和地面監(jiān)測(cè)[7]。由于反演所用的光學(xué)氣溶膠厚度值(AOD)時(shí)間序列上不穩(wěn)定,只能通過(guò)分辨率較粗的氣象數(shù)據(jù)修正[8],反演模型精度難以提升。針對(duì)時(shí)間尺度不斷增長(zhǎng)的需求,越來(lái)越多的學(xué)者選擇將研究數(shù)據(jù)主體轉(zhuǎn)向具有穩(wěn)定時(shí)序的監(jiān)測(cè)數(shù)據(jù)。

    本文選擇河南省具有完整監(jiān)測(cè)數(shù)據(jù)的2014—2018年為研究時(shí)段,基于392個(gè)監(jiān)測(cè)站實(shí)時(shí)監(jiān)測(cè)的空氣質(zhì)量數(shù)據(jù),利用4種不同空間插值方法探索出PM2.5模擬分布的最佳結(jié)果,進(jìn)而獲得2014—2018年來(lái)河南省PM2.5質(zhì)量濃度時(shí)空分布格局和演變規(guī)律,為相關(guān)研究和河南省大氣污染防治政策優(yōu)化提供科學(xué)的參考依據(jù)。

    1 研究數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    河南省位于中國(guó)中部,界于北緯31°23'—36°22',東經(jīng)10°21'—116°39',東接安徽、山東,北接河北、山西,西連陜西,南臨湖北,總面積16.7萬(wàn)km2,大部分地區(qū)屬北亞熱帶向暖溫帶過(guò)渡的大陸性季風(fēng)氣候,南北過(guò)渡性明顯。經(jīng)濟(jì)方面,河南省地處沿海開(kāi)放地區(qū)與中西部地區(qū)的結(jié)合部,是中國(guó)經(jīng)濟(jì)活動(dòng)承東啟西、通南達(dá)北的中間地帶。截至2018年,河南省區(qū)域總?cè)丝诤湍闓DP總量分別位于全國(guó)省份排行的第三和第四名。

    1.2 數(shù)據(jù)來(lái)源與處理

    本實(shí)驗(yàn)的數(shù)據(jù)主要包括:①?gòu)恼鏆饩W(wǎng)(https://www.zq12369.com/)上爬取的河南省17個(gè)地級(jí)市共包括國(guó)控點(diǎn)和省控點(diǎn)在內(nèi)的392個(gè)監(jiān)測(cè)站2014年1月1日至2018年12月31日監(jiān)測(cè)的空氣質(zhì)量數(shù)據(jù);②來(lái)自2014—2018年《河南統(tǒng)計(jì)年鑒》的各縣(市)年末常住人口、GDP數(shù)據(jù);③由Esri公司提供的ChinaData矢量數(shù)據(jù)集(包括省級(jí)行政區(qū)、地級(jí)市、縣城等)。

    針對(duì)各站點(diǎn)監(jiān)測(cè)的PM2.5逐時(shí)濃度數(shù)據(jù),剔除了數(shù)據(jù)不完整的123個(gè)站點(diǎn),余下269個(gè)監(jiān)測(cè)站作為后續(xù)插值的采樣點(diǎn)(如圖1所示),逐步求得各采樣點(diǎn)24h日平均質(zhì)量濃度和年平均質(zhì)量濃度。

    1.3 探索性空間數(shù)據(jù)分析

    探索性空間數(shù)據(jù)分析(Exploratory Spatial Data Analysis,ESDA)不僅能顯示空間數(shù)據(jù)的內(nèi)在規(guī)律,而且是后續(xù)插值方法相關(guān)參數(shù)和模型選擇的依據(jù)。

    1.3.1 空間自相關(guān)分析。空間插值建立在空間相關(guān)性分析的基礎(chǔ)上,若分析檢測(cè)的是同一屬性變量,稱(chēng)為空間自相關(guān),反映的是一個(gè)區(qū)域單元上的某種地理現(xiàn)象或某一屬性值與鄰近區(qū)域單元上同一現(xiàn)象或?qū)傩灾档南嚓P(guān)程度[9]。常用Moran’s I系數(shù)來(lái)描述空間鄰近區(qū)域單元的相似程度,計(jì)算結(jié)果在-1~1,大于0呈正相關(guān)關(guān)系,屬性值高或低的地區(qū)都有聚集現(xiàn)象;若結(jié)果小于0,呈負(fù)相關(guān)關(guān)系,數(shù)據(jù)空間分布呈現(xiàn)高低間隔分布的狀態(tài);絕對(duì)值越接近1說(shuō)明相關(guān)性越強(qiáng),若結(jié)果趨近0,則相鄰區(qū)域單元相關(guān)性低,在研究區(qū)呈隨機(jī)分布狀態(tài)。

    1.3.2 同步統(tǒng)計(jì)檢驗(yàn)。同步統(tǒng)計(jì)檢驗(yàn)可以完成檢驗(yàn)數(shù)據(jù)分布、尋找離群值以及進(jìn)行全局趨勢(shì)分析等一系列任務(wù)。直方圖可以對(duì)大量采樣點(diǎn)進(jìn)行整理加工,找出其數(shù)據(jù)分布的規(guī)律和形態(tài),以便對(duì)其總體分布特征進(jìn)行推斷。QQ圖是一種用于用于檢驗(yàn)數(shù)據(jù)分布的相似性的統(tǒng)計(jì)圖,被檢驗(yàn)的采樣點(diǎn)數(shù)據(jù)越符合所指定的分布,越聚集在一條直線上。趨勢(shì)分析圖用X、Y平面表示采樣點(diǎn)位置,用Z軸表示數(shù)據(jù)屬性值,通過(guò)觀察Z軸數(shù)據(jù)值投影到X、Y平面散點(diǎn)曲線的趨勢(shì),判斷數(shù)據(jù)是否存在全局趨勢(shì)以及其趨勢(shì)形態(tài),若存在,可使用某種確定性?xún)?nèi)插方法生成表面,或通過(guò)范克里金或協(xié)同克里金模型移除全局趨勢(shì)。

    1.4 插值方法

    采樣點(diǎn)的數(shù)量有限,空間分布也不夠均衡,要獲得區(qū)域內(nèi)整體數(shù)據(jù)的分布認(rèn)知,必然要進(jìn)行空間插值??臻g插值的方法多樣,每種方法都有其特點(diǎn)和適用性。用多種方法對(duì)PM2.5濃度插值,在保證預(yù)測(cè)結(jié)果適用性的同時(shí)還可對(duì)比出其誤差規(guī)律[10],并探索與PM2.5有較強(qiáng)相關(guān)性的因素,為河南省PM2.5濃度預(yù)測(cè)和調(diào)控工作提供參考。

    空間插值的公式為:

    1.4.1 確定性空間數(shù)據(jù)插值。反距離加權(quán)法(Inverse Distanceweight,IDW)認(rèn)為,預(yù)測(cè)點(diǎn)的數(shù)據(jù)值是其鄰近采樣點(diǎn)綜合貢獻(xiàn)的結(jié)果,貢獻(xiàn)程度與距離成反比。確定權(quán)重的公式為:

    1.4.2 克里金(Kriging)插值??死锝鸩逯捣ú粌H考慮了預(yù)測(cè)點(diǎn)與鄰近采樣點(diǎn)的空間距離關(guān)系,還考慮了其位置關(guān)系。該方法建立在區(qū)域化變量理論分析的基礎(chǔ)上,以半變異函數(shù)作為工具,對(duì)有限區(qū)域內(nèi)的區(qū)域化變量取值進(jìn)行線性無(wú)偏最優(yōu)估計(jì)。

    半變異函數(shù)是一個(gè)關(guān)于數(shù)據(jù)點(diǎn)的半變異值與數(shù)據(jù)點(diǎn)間距離的函數(shù),計(jì)

    1.4.2.2 泛克里金(Universal Kriging,UK)。對(duì)不滿足平穩(wěn)性假設(shè)的數(shù)據(jù),UK可以對(duì)空間坐標(biāo)表示的變量按照自相關(guān)進(jìn)行回歸分析,即線性無(wú)偏最優(yōu)估計(jì)時(shí)考慮其趨勢(shì),其模型可表示為:

    1.4.2.3 協(xié)同克里金法(Collaborative Kriging,CK)。CK是OK區(qū)域化變量單個(gè)向多個(gè)的拓展,通過(guò)計(jì)算屬性間的半方差函數(shù)和交叉半方差函數(shù)利用變量間的空間連續(xù)性、相關(guān)關(guān)系[11],較為復(fù)雜。

    1.5 插值檢驗(yàn)

    1.5.1 交叉驗(yàn)證。本實(shí)驗(yàn)選擇ArcMap10.4軟件Geostatistical Analyst模塊中的交叉驗(yàn)證工具調(diào)整插值方法的參數(shù),根據(jù)平均誤差最接近于零,均方根預(yù)測(cè)誤差最小,平均標(biāo)準(zhǔn)誤差最接近于均方根預(yù)測(cè)誤差,標(biāo)準(zhǔn)均方根預(yù)測(cè)誤差最接近1的標(biāo)準(zhǔn)判斷出克里金插值最合適的半變異函數(shù)模型。

    1.5.2 子集驗(yàn)證。實(shí)驗(yàn)通過(guò)子集驗(yàn)證的方式比較四種不同插值方法的精度。先從研究區(qū)的269個(gè)采樣點(diǎn)中隨機(jī)選擇20個(gè)作為檢驗(yàn)站點(diǎn),剩余249個(gè)采樣點(diǎn)作為分析插值數(shù)據(jù)重新進(jìn)行插值,各點(diǎn)分布如圖1所示。對(duì)不同插值結(jié)果的精度選擇絕對(duì)平均絕對(duì)誤差(Mean Absolute Error,MAE)和均方根誤差(Root-Mean-Square Error,RMSE)進(jìn)行比較。二者越小插值精度越高;MAE相等時(shí),RMSE越小越好。

    2 插值應(yīng)用過(guò)程

    首先觀察ESDA結(jié)果(見(jiàn)圖2),直方圖顯示年均PM2.5濃度的中位數(shù)基本都與平均值較為接近,說(shuō)明數(shù)據(jù)接近正態(tài)分布;而正態(tài)QQ圖則可以顯示出數(shù)據(jù)總體是傾向于正態(tài)分布的,插值前無(wú)需進(jìn)行相應(yīng)的數(shù)據(jù)轉(zhuǎn)換。通過(guò)趨勢(shì)分析可以發(fā)現(xiàn),河南省各市的年均PM2.5濃度具有一定的空間相關(guān)性,其基本趨勢(shì)為北高南低,東西向?yàn)橹胁扛邇蓚?cè)低。

    OK不考慮移除趨勢(shì),其精度主要受半變異函數(shù)的影響。本文在保證變異函數(shù)的“步長(zhǎng)個(gè)數(shù)[×]步長(zhǎng)大小=研究區(qū)域長(zhǎng)度的一半”的情況下,分別選擇球面、高斯、指數(shù)、三角函數(shù)四種半變異函數(shù)模型進(jìn)行插值,對(duì)比交叉驗(yàn)證結(jié)果后,除2015年選擇三角函數(shù)外,其余各年都選擇高斯函數(shù)做OK的半變異函數(shù)模型。UK綜合了探索性數(shù)據(jù)分析結(jié)果去除數(shù)據(jù)的趨勢(shì)特征,選擇趨勢(shì)的移除階數(shù)為1,確定高斯函數(shù)作為適合表面的核函數(shù),其半變異函數(shù)建模過(guò)程與普通克里金保持一致。CK考慮霧霾與經(jīng)濟(jì)發(fā)展的緊密關(guān)聯(lián),分析二者的正交協(xié)方差函數(shù)云圖后,選擇各市縣的年末GDP作為增加的區(qū)域化變量,半變異函數(shù)建模過(guò)程依舊與OK保持一致。IDW需要確定的參數(shù)極少,這里經(jīng)調(diào)試后統(tǒng)一設(shè)置指數(shù)為2,權(quán)重取自15個(gè)相鄰要素。

    3 實(shí)驗(yàn)結(jié)果分析

    3.1 2014—2018年河南PM2.5濃度的空間集聚趨勢(shì)

    實(shí)驗(yàn)利用GeoDa軟件中的Arc Distance距離算法生成空間權(quán)重矩陣計(jì)算各年的Moran’s I,結(jié)果如圖3所示。

    從圖3可以看出,PM2.5的分布在2014和2017年具有一定的空間自相關(guān)性,存在集聚現(xiàn)象;在2016和2018年存在空間相關(guān)性,但相關(guān)性不高;而2015年幾乎沒(méi)有自相關(guān)性。觀測(cè)Moran散點(diǎn)圖可以看出濃度值的聚集類(lèi)型。相關(guān)性最強(qiáng)的是2014年,其散點(diǎn)圖中顯示,大部分測(cè)點(diǎn)在一、三象限,說(shuō)明在2014年河南省內(nèi)PM2.5多呈冷、熱點(diǎn)分明的分布模式。

    3.2 插值預(yù)測(cè)結(jié)果分析

    3.2.1 不同方法預(yù)測(cè)結(jié)果比較。對(duì)插值結(jié)果制圖得到不同方法模擬的PM2.5濃度空間分布圖像(圖略)。IDW預(yù)測(cè)結(jié)果非常平滑,但能明顯看出PM2.5濃度的極值只出現(xiàn)在采樣點(diǎn)處,受采樣點(diǎn)分布的影響較大,尤其是在河南省中北部地區(qū)產(chǎn)生很多孤立異常的圖斑,與實(shí)際分布情況有差距。OK所得濃度結(jié)果過(guò)渡性好,比較符合實(shí)際,可以體現(xiàn)出空間變化趨勢(shì),制圖的視覺(jué)效果也很理想,但其存在預(yù)測(cè)的最高值偏低,而最低值偏高的問(wèn)題。UK也可以反映出PM2.5的空間變化趨勢(shì),但成圖存在塊狀效果。CK明顯受到GDP變量的影響,使得模擬的PM2.5濃度在部分地區(qū)空間變異過(guò)大。GDP極值處常有孤立異常的圖斑存在,與實(shí)際情況有區(qū)別,這也說(shuō)明實(shí)際的PM2.5濃度分布受眾多因素影響,僅考慮GDP作為協(xié)同變量還不夠??梢暬螼K最優(yōu),可以準(zhǔn)確地表現(xiàn)PM2.5的空間分布特征與變化規(guī)律。

    3.2.2 2014—2018年河南PM2.5的空間分布及動(dòng)態(tài)變化。4種插值方法均大致顯示出河南省的PM2.5濃度分布格局與變化規(guī)律,綜合趨勢(shì)面分析結(jié)果:近年來(lái),PM2.5濃度大致從由中北部向四周遞減變?yōu)橛杀辈肯蚰喜窟f減,在東西方向上依然呈現(xiàn)中部高兩側(cè)低的特點(diǎn)。此外,插值預(yù)測(cè)到的PM2.5濃度的年均值的域值逐年降低,說(shuō)明2014年以來(lái)河南省大氣質(zhì)量整體上是逐年改善的。

    3.3 插值精度分析

    由于IDW沒(méi)有對(duì)預(yù)測(cè)誤差產(chǎn)生估計(jì),并不適用于交叉驗(yàn)證,實(shí)驗(yàn)轉(zhuǎn)而利用子集驗(yàn)證的方式進(jìn)行誤差比對(duì),對(duì)比結(jié)果如表1所示。首先確定實(shí)測(cè)值與各種插值方法的估算值無(wú)顯著差異。

    綜合兩指標(biāo)的計(jì)算值及其排序情況得出結(jié)論:河南省2014和2016年的年均PM2.5濃度用UK插值所得結(jié)果精度最高,2015、2017年和2018年用OK插值所得結(jié)果精度最高。雖然不同插值方法插值精度沒(méi)有特別顯著的差別,但克里金法要普遍優(yōu)于確定性空間插值方法。

    4 討論與總結(jié)

    本文選取了2014—2018年河南省17個(gè)地級(jí)市269個(gè)監(jiān)測(cè)站點(diǎn)的空氣質(zhì)量資料,以此為采樣點(diǎn)檢驗(yàn)數(shù)據(jù)分布特征,利用IDW、OK、UK、CK共4種方法分別進(jìn)行PM2.5年均濃度的插值,制作了各方法模擬各年P(guān)M2.5空間分布的圖像,通過(guò)交叉驗(yàn)證確定高斯函數(shù)為OK擬合區(qū)域化變量的最優(yōu)函數(shù),接著進(jìn)行子集驗(yàn)證,結(jié)合ESDA與OK的結(jié)果,最后得出結(jié)論:①2014和2016年UK的精度最高,2015、2017年和2018年OK精度最高,視覺(jué)上OK最能反映空間變化規(guī)律;②PM2.5年均濃度大致從由中北部向四周遞減變?yōu)橛杀辈肯蚰喜窟f減,在東西方向上依然呈現(xiàn)出中部高兩側(cè)低的特點(diǎn);③插值預(yù)測(cè)到的PM2.5濃度的年均值的域值在逐年降低,說(shuō)明2014年以來(lái)河南省大氣質(zhì)量整體上是逐年改善的。

    此外,CK的應(yīng)用過(guò)程表明,近年來(lái),經(jīng)濟(jì)的發(fā)展與空氣質(zhì)量的沖突在逐漸減弱,體現(xiàn)了河南省經(jīng)濟(jì)結(jié)構(gòu)正逐步向與環(huán)境協(xié)調(diào)化發(fā)展轉(zhuǎn)型的特征。深入研究PM2.5的影響因素,如風(fēng)力風(fēng)向、濕度、汽車(chē)尾氣排放,從而確定最適宜研究區(qū)的多協(xié)同變量的克里金模型也是今后研究的主要方向。

    參考文獻(xiàn):

    [1]Lelieveld J, Evans J S, Fnais M, et al. The contribution of outdoor air pollution sources to premature mortality on a global scale[J]. Nature,2015(7569):367-371.

    [2]盛黎,繆育聰,徐權(quán)周,等.對(duì)話黃順祥:走近中國(guó)大氣污染防治[J].科學(xué)通報(bào),2018(16):1515-1518.

    [3]王振波,方創(chuàng)琳,許光,等.2014年中國(guó)城市PM2.5濃度的時(shí)空變化規(guī)律[J].地理學(xué)報(bào),2015(11):1720-1734.

    [4]李名升,任曉霞,于洋,等.中國(guó)大陸城市PM污染時(shí)空分布規(guī)律[J].中國(guó)環(huán)境科學(xué),2016(3):641-650.

    [5]楊興川,趙文吉,熊秋林,等.2016年京津冀地區(qū)PM2.5時(shí)空分布特征及其與氣象因素的關(guān)系[J].生態(tài)環(huán)境學(xué)報(bào),2017(10):1747-1754.

    [6]戴昭鑫,張?jiān)浦?,胡云鋒,等.基于地面監(jiān)測(cè)數(shù)據(jù)的2013—2015年長(zhǎng)三角地區(qū)PM2.5時(shí)空特征[J].長(zhǎng)江流域資源與環(huán)境,2016(5):1004-8227.

    [7]熊歡歡,梁龍武,曾贈(zèng),等.中國(guó)城市PM2.5時(shí)空分布的動(dòng)態(tài)比較分析[J].資源科學(xué),2017(1):136-146.

    [8]張?jiān)浦?基于地面監(jiān)測(cè)和遙感反演的京津冀地區(qū)PM2.5時(shí)空分布研究[D].呼和浩特:內(nèi)蒙古師范大學(xué),2016.

    [9]劉湘南,王平,關(guān)麗,盧浩,張春曉.GIS空間分析[M].3版.北京:科學(xué)出版社,2008.

    [10]朱求安,張萬(wàn)昌,余鈞輝.基于GIS的空間插值方法研究[J].江西師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2004(2):183-188.

    [11]章清,張海濤,郭龍,等.基于主成分分析的協(xié)同克里格插值模型對(duì)土壤銅含量的空間分布預(yù)測(cè)[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào),2016(1):60-68.

    猜你喜歡
    克里插值河南省
    今晚不能去你家玩啦!
    我可以咬一口嗎?
    知識(shí)窗(2023年2期)2023-03-05 11:28:27
    河南省樹(shù)人教育交流中心
    河南省樹(shù)人教育交流中心
    河南省樹(shù)人教育交流中心
    河南省樹(shù)人教育交流中心
    你今天真好看
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    你今天真好看
    讀者(2018年24期)2018-12-04 03:01:34
    一種改進(jìn)FFT多譜線插值諧波分析方法
    国内精品宾馆在线| 国产亚洲最大av| 久久久午夜欧美精品| 日本午夜av视频| 国产精品人妻久久久影院| 欧美一区二区亚洲| 综合色av麻豆| 18禁在线无遮挡免费观看视频| 18禁在线无遮挡免费观看视频| 99久国产av精品国产电影| 中文乱码字字幕精品一区二区三区 | 亚洲精品成人久久久久久| 日韩精品青青久久久久久| 国产精品综合久久久久久久免费| 国产精品1区2区在线观看.| 国产精品人妻久久久久久| 天堂影院成人在线观看| av女优亚洲男人天堂| 久久6这里有精品| 欧美性感艳星| 免费观看无遮挡的男女| 男插女下体视频免费在线播放| 午夜福利在线观看吧| av国产久精品久网站免费入址| 欧美日韩视频高清一区二区三区二| 在线观看美女被高潮喷水网站| 欧美一区二区亚洲| 亚洲最大成人av| 乱人视频在线观看| 淫秽高清视频在线观看| 高清av免费在线| 狂野欧美激情性xxxx在线观看| 精品一区在线观看国产| 狂野欧美激情性xxxx在线观看| 一级毛片久久久久久久久女| 免费看美女性在线毛片视频| 天美传媒精品一区二区| 色播亚洲综合网| 亚洲欧美日韩卡通动漫| 2021少妇久久久久久久久久久| 精品一区二区三区视频在线| 午夜爱爱视频在线播放| 日韩一区二区三区影片| 老女人水多毛片| 久久久久免费精品人妻一区二区| 精品国产露脸久久av麻豆 | 欧美日韩视频高清一区二区三区二| 日本一二三区视频观看| 欧美3d第一页| 日本一本二区三区精品| 十八禁网站网址无遮挡 | 搡老乐熟女国产| 一区二区三区乱码不卡18| 亚洲成人一二三区av| 国产免费又黄又爽又色| 亚洲在久久综合| 国产免费福利视频在线观看| 国产 一区精品| 最近最新中文字幕免费大全7| 国产欧美另类精品又又久久亚洲欧美| 99九九线精品视频在线观看视频| 成人无遮挡网站| 人妻少妇偷人精品九色| 韩国高清视频一区二区三区| 亚洲美女搞黄在线观看| 人妻一区二区av| 亚洲成人av在线免费| 亚洲av日韩在线播放| 哪个播放器可以免费观看大片| 国产精品一区二区三区四区久久| 国精品久久久久久国模美| 能在线免费观看的黄片| 亚洲不卡免费看| 午夜免费男女啪啪视频观看| 国产av码专区亚洲av| 亚洲熟女精品中文字幕| 国产亚洲91精品色在线| 嫩草影院入口| 建设人人有责人人尽责人人享有的 | 久久鲁丝午夜福利片| av网站免费在线观看视频 | 禁无遮挡网站| 日日摸夜夜添夜夜添av毛片| 少妇的逼水好多| 日本免费在线观看一区| 日韩av在线免费看完整版不卡| av国产免费在线观看| 亚洲内射少妇av| 18禁在线无遮挡免费观看视频| 韩国高清视频一区二区三区| 免费观看精品视频网站| 三级男女做爰猛烈吃奶摸视频| 自拍偷自拍亚洲精品老妇| 国产亚洲91精品色在线| 大片免费播放器 马上看| 国产精品美女特级片免费视频播放器| av线在线观看网站| 麻豆成人av视频| 日韩av在线免费看完整版不卡| 日韩大片免费观看网站| 亚洲av电影在线观看一区二区三区 | 久久这里只有精品中国| 国产黄a三级三级三级人| 日韩强制内射视频| 日本av手机在线免费观看| 亚洲最大成人手机在线| 美女高潮的动态| 中文字幕人妻熟人妻熟丝袜美| 99久久精品一区二区三区| 国产毛片a区久久久久| 麻豆久久精品国产亚洲av| 亚洲精品日韩av片在线观看| 国产成人91sexporn| 国产亚洲午夜精品一区二区久久 | 欧美日韩国产mv在线观看视频 | 亚洲自偷自拍三级| 欧美日韩视频高清一区二区三区二| 干丝袜人妻中文字幕| 亚洲人成网站高清观看| 国产av码专区亚洲av| 99热网站在线观看| 超碰av人人做人人爽久久| 晚上一个人看的免费电影| 亚洲av二区三区四区| 我的女老师完整版在线观看| 久久久久精品性色| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 男插女下体视频免费在线播放| 日韩一区二区视频免费看| 2021天堂中文幕一二区在线观| 91在线精品国自产拍蜜月| 久久精品国产鲁丝片午夜精品| 日本-黄色视频高清免费观看| 色哟哟·www| 国产亚洲一区二区精品| 日韩国内少妇激情av| 免费av毛片视频| 日韩亚洲欧美综合| 欧美高清成人免费视频www| 日本色播在线视频| 18禁裸乳无遮挡免费网站照片| 夜夜爽夜夜爽视频| 2021天堂中文幕一二区在线观| 亚洲国产精品成人综合色| 又黄又爽又刺激的免费视频.| 亚洲精品aⅴ在线观看| 乱人视频在线观看| 国产乱人偷精品视频| 一区二区三区高清视频在线| 最近中文字幕高清免费大全6| 国产久久久一区二区三区| 国产精品一区二区性色av| 嫩草影院精品99| 午夜福利成人在线免费观看| 欧美三级亚洲精品| 亚洲国产成人一精品久久久| 欧美不卡视频在线免费观看| 亚洲自偷自拍三级| 一级av片app| 美女脱内裤让男人舔精品视频| 国产精品三级大全| 亚洲人成网站在线播| 国产一级毛片七仙女欲春2| 日韩强制内射视频| 国产久久久一区二区三区| 日韩欧美 国产精品| 青青草视频在线视频观看| 亚洲美女搞黄在线观看| ponron亚洲| 日本免费a在线| 色尼玛亚洲综合影院| 纵有疾风起免费观看全集完整版 | 天堂av国产一区二区熟女人妻| 亚洲精品日韩av片在线观看| 日韩亚洲欧美综合| 看免费成人av毛片| 一边亲一边摸免费视频| 极品教师在线视频| 麻豆久久精品国产亚洲av| 国产综合精华液| 国产大屁股一区二区在线视频| 亚洲精品自拍成人| 丰满少妇做爰视频| 亚洲精品一二三| 亚洲av成人精品一二三区| 五月天丁香电影| 人妻一区二区av| 99久久中文字幕三级久久日本| 亚洲欧美一区二区三区国产| 看十八女毛片水多多多| 成人特级av手机在线观看| 中文字幕免费在线视频6| 精品一区二区三区人妻视频| 久久精品久久精品一区二区三区| 中文字幕久久专区| 精品久久久久久久久亚洲| 在线 av 中文字幕| 亚洲最大成人中文| 国产亚洲午夜精品一区二区久久 | 一个人看的www免费观看视频| 精品久久久久久久久av| 欧美极品一区二区三区四区| 一区二区三区高清视频在线| 91精品伊人久久大香线蕉| 免费观看精品视频网站| 免费播放大片免费观看视频在线观看| 欧美xxxx性猛交bbbb| h日本视频在线播放| 国产精品三级大全| av黄色大香蕉| 精品人妻视频免费看| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 国产男人的电影天堂91| 国产精品三级大全| 99久国产av精品| 日本黄大片高清| 一夜夜www| 联通29元200g的流量卡| 日韩不卡一区二区三区视频在线| 97精品久久久久久久久久精品| 成人亚洲精品一区在线观看 | 午夜日本视频在线| 丰满少妇做爰视频| 日韩国内少妇激情av| 久久人人爽人人爽人人片va| 高清毛片免费看| 三级国产精品欧美在线观看| 亚洲国产欧美人成| 一二三四中文在线观看免费高清| 两个人视频免费观看高清| 欧美丝袜亚洲另类| 国产综合懂色| 精品酒店卫生间| 99热网站在线观看| 精品人妻熟女av久视频| 精品人妻一区二区三区麻豆| 在线天堂最新版资源| 少妇的逼好多水| a级毛色黄片| 特级一级黄色大片| 国产一区二区三区av在线| 亚洲av福利一区| 亚洲精品国产av成人精品| 一级爰片在线观看| 日本-黄色视频高清免费观看| 国产精品一及| 亚洲婷婷狠狠爱综合网| 亚洲图色成人| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美成人综合另类久久久| 日日摸夜夜添夜夜爱| 色尼玛亚洲综合影院| 色综合色国产| 日韩制服骚丝袜av| 欧美zozozo另类| av线在线观看网站| 亚洲成人av在线免费| 国产不卡一卡二| 色网站视频免费| 免费电影在线观看免费观看| 69人妻影院| 国产精品嫩草影院av在线观看| 亚洲天堂国产精品一区在线| 国产精品99久久久久久久久| 少妇丰满av| 在线观看免费高清a一片| 国产在线一区二区三区精| 高清午夜精品一区二区三区| 免费av毛片视频| 春色校园在线视频观看| 国产乱人偷精品视频| 看十八女毛片水多多多| 淫秽高清视频在线观看| 国产免费又黄又爽又色| 最近最新中文字幕免费大全7| 国产精品久久久久久久久免| 免费高清在线观看视频在线观看| 久久久久精品性色| 永久网站在线| 日韩欧美精品免费久久| 国产精品.久久久| 一级毛片电影观看| 久久精品熟女亚洲av麻豆精品 | 最后的刺客免费高清国语| 久久鲁丝午夜福利片| 嫩草影院精品99| 精品国内亚洲2022精品成人| 97超视频在线观看视频| 伊人久久国产一区二区| 26uuu在线亚洲综合色| 国产精品一二三区在线看| 我的老师免费观看完整版| 国语对白做爰xxxⅹ性视频网站| 伊人久久精品亚洲午夜| 听说在线观看完整版免费高清| 一级毛片aaaaaa免费看小| 一级毛片 在线播放| 免费黄网站久久成人精品| av播播在线观看一区| 狂野欧美白嫩少妇大欣赏| 亚洲欧美一区二区三区黑人 | freevideosex欧美| 51国产日韩欧美| 永久网站在线| 日本一二三区视频观看| 最近中文字幕2019免费版| 嫩草影院新地址| 亚洲精品一二三| 男女边吃奶边做爰视频| 欧美潮喷喷水| 一边亲一边摸免费视频| 啦啦啦韩国在线观看视频| 国产 亚洲一区二区三区 | 亚洲精品456在线播放app| 午夜日本视频在线| 午夜激情福利司机影院| 亚洲精品视频女| 简卡轻食公司| 中文天堂在线官网| 午夜爱爱视频在线播放| 午夜福利高清视频| 男人爽女人下面视频在线观看| 亚洲国产高清在线一区二区三| 欧美成人a在线观看| 十八禁网站网址无遮挡 | 精品国产露脸久久av麻豆 | 国产一区二区三区综合在线观看 | 男女国产视频网站| 国产午夜精品一二区理论片| 亚洲欧洲日产国产| 久久这里有精品视频免费| 午夜激情久久久久久久| 国产免费又黄又爽又色| 高清av免费在线| 国产日韩欧美在线精品| 久久人人爽人人爽人人片va| 免费人成在线观看视频色| 亚洲av免费在线观看| 又爽又黄a免费视频| 国产伦精品一区二区三区视频9| 18禁动态无遮挡网站| 97人妻精品一区二区三区麻豆| 女人被狂操c到高潮| av播播在线观看一区| a级毛色黄片| 日日啪夜夜撸| a级毛色黄片| 2022亚洲国产成人精品| 青春草国产在线视频| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 亚洲熟女精品中文字幕| 亚洲av不卡在线观看| 亚洲国产精品sss在线观看| 建设人人有责人人尽责人人享有的 | 国产单亲对白刺激| 精华霜和精华液先用哪个| 成人鲁丝片一二三区免费| 人人妻人人澡欧美一区二区| 国产精品伦人一区二区| 少妇被粗大猛烈的视频| av线在线观看网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 九色成人免费人妻av| 日韩欧美国产在线观看| 亚洲18禁久久av| 99热这里只有是精品50| 精品午夜福利在线看| 亚洲在线自拍视频| 国产成人精品婷婷| 赤兔流量卡办理| 亚洲国产欧美人成| 中国美白少妇内射xxxbb| 成人漫画全彩无遮挡| 国产精品不卡视频一区二区| 综合色av麻豆| 欧美+日韩+精品| 一夜夜www| 高清视频免费观看一区二区 | 国产有黄有色有爽视频| videossex国产| 在线a可以看的网站| 亚洲综合精品二区| 亚洲精品aⅴ在线观看| 人妻少妇偷人精品九色| 久久99精品国语久久久| 国产v大片淫在线免费观看| 3wmmmm亚洲av在线观看| 日韩一区二区三区影片| 成年女人看的毛片在线观看| av一本久久久久| 老女人水多毛片| 最新中文字幕久久久久| 99久国产av精品| 在线观看一区二区三区| 伊人久久国产一区二区| 中文精品一卡2卡3卡4更新| 国产综合懂色| 日韩av免费高清视频| 精品久久久久久久人妻蜜臀av| a级毛片免费高清观看在线播放| 久久久欧美国产精品| 亚洲最大成人av| 国产成人一区二区在线| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 卡戴珊不雅视频在线播放| 老司机影院毛片| 亚洲综合精品二区| 国产爱豆传媒在线观看| 一级黄片播放器| 亚洲精品视频女| 国产69精品久久久久777片| 蜜臀久久99精品久久宅男| 麻豆精品久久久久久蜜桃| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 久久人人爽人人爽人人片va| 在线播放无遮挡| 69人妻影院| 亚洲经典国产精华液单| 日本欧美国产在线视频| 男女啪啪激烈高潮av片| 如何舔出高潮| 国产成人免费观看mmmm| 亚洲在线观看片| 亚洲怡红院男人天堂| 亚洲精品456在线播放app| 又大又黄又爽视频免费| 国产成人a区在线观看| 亚洲欧洲日产国产| 一个人看视频在线观看www免费| 丰满乱子伦码专区| 久久久久免费精品人妻一区二区| 国产成人a∨麻豆精品| 一级毛片黄色毛片免费观看视频| 国产黄色视频一区二区在线观看| 少妇人妻一区二区三区视频| 精品久久久久久久久久久久久| 日韩三级伦理在线观看| 国产精品人妻久久久久久| 日本爱情动作片www.在线观看| 亚洲经典国产精华液单| 一级毛片我不卡| 色综合色国产| 欧美zozozo另类| 精品酒店卫生间| 国产在视频线精品| 麻豆乱淫一区二区| 久久99精品国语久久久| 欧美成人午夜免费资源| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中国美白少妇内射xxxbb| 一个人免费在线观看电影| 亚洲精品日本国产第一区| 一个人看视频在线观看www免费| 久久久国产一区二区| 别揉我奶头 嗯啊视频| 中文欧美无线码| 亚洲精品日本国产第一区| 肉色欧美久久久久久久蜜桃 | av在线天堂中文字幕| 天堂av国产一区二区熟女人妻| 国产欧美另类精品又又久久亚洲欧美| 高清在线视频一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久久久免| 国产黄频视频在线观看| 少妇人妻精品综合一区二区| 两个人的视频大全免费| 国产亚洲精品久久久com| 大陆偷拍与自拍| 美女内射精品一级片tv| 大香蕉久久网| 日日干狠狠操夜夜爽| 不卡视频在线观看欧美| 乱人视频在线观看| 久久久欧美国产精品| 淫秽高清视频在线观看| 欧美不卡视频在线免费观看| 国产老妇伦熟女老妇高清| 国产黄频视频在线观看| 午夜久久久久精精品| 久久久久久久午夜电影| 少妇高潮的动态图| 在线免费十八禁| 美女大奶头视频| 日韩欧美国产在线观看| 噜噜噜噜噜久久久久久91| 看十八女毛片水多多多| 乱系列少妇在线播放| av国产免费在线观看| 日本一本二区三区精品| 青青草视频在线视频观看| 三级男女做爰猛烈吃奶摸视频| 精品人妻熟女av久视频| 日韩精品有码人妻一区| 国产欧美日韩精品一区二区| av在线播放精品| 国产成人精品久久久久久| 久久精品久久久久久久性| 亚洲怡红院男人天堂| 久久这里有精品视频免费| 久久97久久精品| 国产精品一区二区性色av| 又爽又黄无遮挡网站| 最近最新中文字幕免费大全7| 国产又色又爽无遮挡免| 18禁在线无遮挡免费观看视频| ponron亚洲| 成人亚洲精品av一区二区| 日韩三级伦理在线观看| 特级一级黄色大片| 亚洲av免费高清在线观看| 免费看a级黄色片| 天堂影院成人在线观看| 亚洲av电影在线观看一区二区三区 | 一级毛片黄色毛片免费观看视频| 青春草国产在线视频| 亚洲av.av天堂| 国产又色又爽无遮挡免| 成人无遮挡网站| 日日啪夜夜撸| 中文字幕人妻熟人妻熟丝袜美| 噜噜噜噜噜久久久久久91| 国产高清不卡午夜福利| 麻豆精品久久久久久蜜桃| 久久久久久伊人网av| 久久精品熟女亚洲av麻豆精品 | 日日摸夜夜添夜夜添av毛片| 插阴视频在线观看视频| 免费看a级黄色片| 老师上课跳d突然被开到最大视频| 中国国产av一级| 精品国产一区二区三区久久久樱花 | 国产不卡一卡二| 成人一区二区视频在线观看| 卡戴珊不雅视频在线播放| 午夜福利高清视频| av线在线观看网站| 青春草亚洲视频在线观看| 国产有黄有色有爽视频| 女人被狂操c到高潮| 五月玫瑰六月丁香| 看免费成人av毛片| 国产精品一区www在线观看| 男女下面进入的视频免费午夜| 国产伦精品一区二区三区四那| av卡一久久| 非洲黑人性xxxx精品又粗又长| 亚洲伊人久久精品综合| 久久亚洲国产成人精品v| 男插女下体视频免费在线播放| 国产午夜精品久久久久久一区二区三区| 十八禁国产超污无遮挡网站| 一级毛片久久久久久久久女| 国产成人精品久久久久久| 欧美成人精品欧美一级黄| 国产黄a三级三级三级人| 三级经典国产精品| 久久久久性生活片| 美女黄网站色视频| 国产一区二区三区综合在线观看 | 精品欧美国产一区二区三| eeuss影院久久| 成人午夜精彩视频在线观看| 欧美成人一区二区免费高清观看| av网站免费在线观看视频 | 精华霜和精华液先用哪个| 欧美潮喷喷水| videossex国产| 国产午夜精品久久久久久一区二区三区| 国产成人精品一,二区| 丝袜美腿在线中文| freevideosex欧美| 成人美女网站在线观看视频| 一区二区三区四区激情视频| 国内少妇人妻偷人精品xxx网站| 九九爱精品视频在线观看| 久久久精品欧美日韩精品| 国产熟女欧美一区二区| 极品教师在线视频| 久久久精品94久久精品| 看十八女毛片水多多多| 青春草国产在线视频| 国产男人的电影天堂91| 成人亚洲精品一区在线观看 | 日日干狠狠操夜夜爽| 男女啪啪激烈高潮av片| 国产综合懂色| 久久6这里有精品| 神马国产精品三级电影在线观看| 日本一本二区三区精品| 日韩av在线大香蕉| 免费看a级黄色片| 精品国内亚洲2022精品成人| 18禁在线无遮挡免费观看视频| 黄片无遮挡物在线观看| 成人欧美大片| 一本久久精品| 欧美bdsm另类| 中国美白少妇内射xxxbb| 成年人午夜在线观看视频 | 国产v大片淫在线免费观看| 97在线视频观看| 免费观看在线日韩| 欧美高清性xxxxhd video| 亚洲av电影在线观看一区二区三区 | 久久久午夜欧美精品| 天堂网av新在线| 午夜福利在线观看吧| 哪个播放器可以免费观看大片| 少妇的逼水好多|