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

    卡爾曼濾波在海洋浮標(biāo)數(shù)據(jù)預(yù)處理中的應(yīng)用

    2024-01-26 08:19:26張新文林冠英劉同木周保成
    關(guān)鍵詞:原始數(shù)據(jù)浮標(biāo)卡爾曼濾波

    張新文,林冠英,劉同木,周保成

    (1.自然資源部南海調(diào)查中心,廣東 廣州 510300;2.自然資源部海洋環(huán)境探測(cè)技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室,廣東 廣州 510300)

    海洋浮標(biāo)是利用無(wú)動(dòng)力漂浮載體獲取海洋環(huán)境信息的自動(dòng)化監(jiān)測(cè)系統(tǒng),具有“海洋上的地球同步衛(wèi)星”和“海上天氣偵察兵”等稱(chēng)號(hào),是海洋環(huán)境監(jiān)測(cè)、海洋資源調(diào)查評(píng)估和海洋科學(xué)研究的重要平臺(tái)[1-2]。受海洋環(huán)境噪聲、海洋生物附著和傳感器元件老化等因素的影響,浮標(biāo)搭載的傳感器測(cè)量值不可避免地存在環(huán)境噪聲,直接影響數(shù)據(jù)應(yīng)用效果[3]。當(dāng)前,海洋浮標(biāo)受到衛(wèi)星通信帶寬限制,一般在浮標(biāo)端需要對(duì)原始數(shù)據(jù)按照固定時(shí)間將觀測(cè)數(shù)據(jù)融合后回傳。因此,在浮標(biāo)端實(shí)時(shí)開(kāi)展數(shù)據(jù)異常值檢測(cè)和測(cè)量噪聲預(yù)處理,是提高數(shù)據(jù)融合精度的可靠手段之一。

    國(guó)內(nèi)外學(xué)者在環(huán)境監(jiān)測(cè)數(shù)據(jù)預(yù)處理方面進(jìn)行了相關(guān)研究,主要分為以下三類(lèi)預(yù)處理方法:(1)基于數(shù)值統(tǒng)計(jì)規(guī)律的數(shù)據(jù)預(yù)處理方法,該類(lèi)方法通過(guò)計(jì)算數(shù)據(jù)組的分布規(guī)律識(shí)別離散異常值[4-5],但通常閾值設(shè)定需要人工干預(yù),在延時(shí)數(shù)據(jù)預(yù)處理方面應(yīng)用較好[6-7]。(2)基于時(shí)空特征的數(shù)據(jù)預(yù)處理方法,該類(lèi)方法將觀測(cè)要素的基于時(shí)間/空間的變化特征進(jìn)行相關(guān)性分析,使得預(yù)處理結(jié)果符合時(shí)間/空間變化規(guī)律,主要包括卡爾曼濾波算法[8-10]、自回歸移動(dòng)平均模型(ARMA)[11]、基于時(shí)序相關(guān)性分析方法[12]等。(3)基于人工智能的數(shù)據(jù)預(yù)處理方法,該類(lèi)方法以大數(shù)據(jù)量的歷史數(shù)據(jù)挖掘分析為主,具有數(shù)據(jù)量大、算法復(fù)雜度高特點(diǎn),適用于監(jiān)測(cè)數(shù)據(jù)的延時(shí)預(yù)處理過(guò)程[13-14]。綜上,基于數(shù)值統(tǒng)計(jì)規(guī)律的數(shù)據(jù)預(yù)處理方法在數(shù)據(jù)異常值剔除方面具有良好的表現(xiàn),但不能實(shí)現(xiàn)浮標(biāo)原始數(shù)據(jù)的實(shí)時(shí)去噪聲處理,而浮標(biāo)數(shù)據(jù)采集器的計(jì)算能力無(wú)法滿足基于人工智能方法的高復(fù)雜度、大計(jì)算量的需求?;跁r(shí)空特征的數(shù)據(jù)預(yù)處理方法計(jì)算開(kāi)銷(xiāo)較低且符合浮標(biāo)測(cè)量要素隨時(shí)間序列變化特征,但在工程應(yīng)用上,浮標(biāo)端尚未實(shí)現(xiàn)原始數(shù)據(jù)的自動(dòng)化預(yù)處理,國(guó)內(nèi)外的相關(guān)研究也以數(shù)據(jù)異常值檢測(cè)或延時(shí)數(shù)據(jù)的批量質(zhì)控處理為主,對(duì)于在浮標(biāo)端開(kāi)展原始測(cè)量數(shù)據(jù)實(shí)時(shí)自動(dòng)化預(yù)處理方面的研究較少。因此,本研究結(jié)合海洋浮標(biāo)的數(shù)據(jù)特征,提出一種基于改進(jìn)的自適應(yīng)卡爾曼濾波自動(dòng)化數(shù)據(jù)預(yù)處理方法,該方法采用箱形圖法對(duì)原始數(shù)據(jù)中的異常值進(jìn)行檢測(cè)和修正,以期通過(guò)改進(jìn)的卡爾曼濾波方法實(shí)現(xiàn)浮標(biāo)觀測(cè)數(shù)據(jù)的噪聲濾波處理,并在廣東某海域進(jìn)行實(shí)際的海洋工程應(yīng)用實(shí)驗(yàn),以驗(yàn)證該方法的有效性和可行性。

    1 數(shù)據(jù)預(yù)處理方法

    1.1 浮標(biāo)數(shù)據(jù)特征

    在確定數(shù)據(jù)預(yù)處理方法前,對(duì)海洋浮標(biāo)搭載的傳感器采樣原理、數(shù)據(jù)特征和采樣/融合頻率等進(jìn)行分析,以保證數(shù)據(jù)處理方法的合理應(yīng)用,提高數(shù)據(jù)預(yù)處理效果。

    依據(jù)GB/T 14914.3—2021《海洋觀測(cè)規(guī)范 第3部分:浮標(biāo)潛標(biāo)觀測(cè)》[15]和HY/T 147.6—2013《海洋監(jiān)測(cè)技術(shù)規(guī)程第6 部分:海洋水文、氣象與海冰》[16]等海洋環(huán)境監(jiān)測(cè)相關(guān)標(biāo)準(zhǔn),各監(jiān)測(cè)要素在到達(dá)融合和回傳時(shí)間后,浮標(biāo)端將監(jiān)測(cè)期間的原始數(shù)據(jù)進(jìn)行提取并融合成該時(shí)段的測(cè)量值,浮標(biāo)數(shù)據(jù)具有采樣與融合頻率固定、單參數(shù)融合、隨時(shí)間序列演變等特征。因此,必須在融合前對(duì)該時(shí)段內(nèi)所有監(jiān)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)一的去噪聲處理。同時(shí),浮標(biāo)監(jiān)測(cè)參數(shù)在時(shí)間序列上符合環(huán)境變化規(guī)律,在設(shè)計(jì)自動(dòng)化數(shù)據(jù)處理程序時(shí)可以根據(jù)參數(shù)的融合頻率,采用基于預(yù)測(cè)與修正的實(shí)時(shí)濾波去噪方法。

    1.2 卡爾曼濾波算法

    卡爾曼濾波算法是一種具有較強(qiáng)容錯(cuò)能力,以最小均方誤差為準(zhǔn)則的遞推無(wú)偏估計(jì)算法,可有效去除系統(tǒng)中的噪聲干擾,在無(wú)人導(dǎo)航、電氣控制、傳感器數(shù)據(jù)融合等多個(gè)領(lǐng)域得到廣泛應(yīng)用[17-18]。傳統(tǒng)的卡爾曼濾波是在標(biāo)準(zhǔn)條件下獲得的,但海洋浮標(biāo)是搭載多傳感器的實(shí)時(shí)環(huán)境監(jiān)測(cè)平臺(tái),在工程應(yīng)用中由于噪聲統(tǒng)計(jì)特性不準(zhǔn)確而產(chǎn)生濾波精度降低甚至產(chǎn)生濾波發(fā)散問(wèn)題。為解決卡爾曼濾波發(fā)散問(wèn)題,有學(xué)者提出Sage-Husa 自適應(yīng)卡爾曼濾波方法[19],該方法通過(guò)時(shí)變?cè)肼暪烙?jì)器,在實(shí)時(shí)處理數(shù)據(jù)過(guò)程中,更新系統(tǒng)噪聲和測(cè)量噪聲,從而抑制濾波發(fā)散。其狀態(tài)空間模型如下:

    式(1)中,xk為系統(tǒng)k時(shí)刻的n維狀態(tài)量;Φk,k-1為系統(tǒng)從k-1 時(shí)刻到k時(shí)刻的n×n維狀態(tài)轉(zhuǎn)移矩陣;zk為系統(tǒng)k時(shí)刻的m維觀測(cè)量;Hk為系統(tǒng)k時(shí)刻的m×n維測(cè)量矩陣;wk-1為系統(tǒng)k-1 時(shí)刻的過(guò)程噪聲;vk為m維測(cè)量噪聲,滿足如下統(tǒng)計(jì)特性:

    式(2)中,j,k=1,2,3,...,當(dāng)k=j時(shí),δkj=1,否 則δkj=0;qk和qj為wk和wj的數(shù)學(xué)期望;rk為vk的數(shù)學(xué)期望;Qk和Rk為系統(tǒng)噪聲和測(cè)量噪聲方差矩陣。

    Sage-Husa自適應(yīng)卡爾曼濾波算法可描述為

    1.3 改進(jìn)的自適應(yīng)卡爾曼濾波

    受海洋環(huán)境噪聲、海洋生物附著和傳感器元件老化等因素的影響,浮標(biāo)數(shù)據(jù)不可避免的存在異常測(cè)量值,這些異常值將導(dǎo)致算法在式(3)求狀態(tài)值時(shí)被放大并引入到狀態(tài)估計(jì)中,特別是當(dāng)異常值偏離較大時(shí),殘差無(wú)法正確修正,會(huì)嚴(yán)重影響濾波精度,因而在使用卡爾曼濾波器進(jìn)行濾波時(shí),必須對(duì)測(cè)量異常數(shù)據(jù)進(jìn)行剔除并修正。同時(shí),測(cè)量噪聲方差在計(jì)算過(guò)程中,當(dāng)初始噪聲設(shè)置過(guò)大或?qū)嶋H海洋環(huán)境噪聲與理論建模值相比較小時(shí),卡爾曼濾波算法將出現(xiàn)濾波發(fā)散問(wèn)題。為提高算法的穩(wěn)定性,需要解決異常值對(duì)濾波的干擾和抑制濾波發(fā)散問(wèn)題。

    1.3.1 異常值修正 基于統(tǒng)計(jì)規(guī)律的方法在異常值檢測(cè)方面應(yīng)用最為廣泛,如拉依達(dá)準(zhǔn)則和狄克遜準(zhǔn)則等基于測(cè)量值正態(tài)分布規(guī)律的方法。但基于正態(tài)分布的疏失誤差檢測(cè)方法是以假定數(shù)據(jù)服從正態(tài)分布為前提的,其判斷異常值的標(biāo)準(zhǔn)是以計(jì)算數(shù)據(jù)組的均值和標(biāo)準(zhǔn)差為基礎(chǔ)的,異常值本身會(huì)對(duì)它們產(chǎn)生較大影響,進(jìn)而導(dǎo)致方法的耐抗性較小。同時(shí),海洋浮標(biāo)數(shù)據(jù)往往難以服從正態(tài)分布規(guī)律,在使用基于正態(tài)分布的異常值檢測(cè)方法時(shí)存在漏剔除和誤剔除現(xiàn)象?;诮y(tǒng)計(jì)規(guī)律的還有箱形圖法,如圖1所示。

    圖1 箱形圖結(jié)構(gòu)Fig.1 Box plot Structure

    該方法利用數(shù)據(jù)序列中的5個(gè)統(tǒng)計(jì)量:下限、下四分位數(shù)(Q1)、中位數(shù)(Q2)、上四分位數(shù)(Q3)與上限來(lái)描述數(shù)據(jù)離散程度的方法。假設(shè)預(yù)處理數(shù)據(jù)集含n個(gè)測(cè)量值,則箱形圖Q1取1×(n+1)/4 位置的數(shù);若n為奇數(shù),則Q2取(n+1)/2 位置的數(shù),否則取n/2 和n/ 2+1 位置數(shù)的均值;Q3取3 × (n+1)/ 4 位置的數(shù)。因此,四分位距IQR 為Q3—Q1,上限Q3+1.5 × IQR 和下限Q1—1.5 × IQR 作為異常值截?cái)帱c(diǎn),測(cè)量值分布在該區(qū)間內(nèi)的判定為有效測(cè)量值,否則判定為異常測(cè)量值。箱形圖判斷異常值的標(biāo)準(zhǔn)以四分位數(shù)和四分位距為基礎(chǔ),具有一定的耐抗性,能夠客觀地識(shí)別數(shù)據(jù)中的異常值。在識(shí)別到異常值后,為使得數(shù)據(jù)維數(shù)保持一致,本研究使用臨近數(shù)據(jù)求取的平均值進(jìn)行插補(bǔ)。

    1.3.2 方差受限抑制濾波發(fā)散 為避免初始噪聲設(shè)置過(guò)大或?qū)嶋H系統(tǒng)噪聲與理論建模值相比較小,將造成測(cè)量噪聲方差失去正定性而導(dǎo)致濾波發(fā)散問(wèn)題。本研究設(shè)計(jì)采用方差受限的方法解決測(cè)量噪聲導(dǎo)致的濾波發(fā)散問(wèn)題,這種改變雖然降低了一定的濾波精度,但提高了濾波的穩(wěn)定性。

    式(4)中,βk作為方差加權(quán)系數(shù),其取值為βk=βk-1/(βk-1+b),其中β0=1,β∞=1 -b,b為漸消因子,取值區(qū)間為[0.90,0.99];Rmin為設(shè)置的方差下限,保證為正,Rmax為設(shè)置的方差上限,能夠快速降低測(cè)量值的可信度并判定是否異常,若異常則放棄本次更新。

    2 數(shù)據(jù)處理流程

    結(jié)合海洋浮標(biāo)數(shù)據(jù)特征和處理程序,提出基于改進(jìn)擴(kuò)展卡爾曼濾波的自動(dòng)化數(shù)據(jù)處理方法,如圖2所示。

    圖2 數(shù)據(jù)預(yù)處理流程Fig.2 Data preprocessing process

    在到達(dá)參數(shù)數(shù)據(jù)融合/回傳時(shí)間后,浮標(biāo)端將監(jiān)測(cè)時(shí)段內(nèi)的原始數(shù)據(jù)集提取作為數(shù)據(jù)預(yù)處理方法的輸入,對(duì)數(shù)據(jù)集依次進(jìn)行常規(guī)檢驗(yàn)、箱型圖異常值檢測(cè)以及卡爾曼濾波去噪處理。主要步驟如下:

    (1)判斷是否到達(dá)參數(shù)數(shù)據(jù)融合與回傳時(shí)間。若到達(dá)則提取該監(jiān)測(cè)時(shí)段內(nèi)參數(shù)的所有采樣數(shù)據(jù),作為預(yù)處理數(shù)據(jù)集D;若未到達(dá)則繼續(xù)進(jìn)行采樣監(jiān)測(cè)。

    (2)數(shù)據(jù)常規(guī)校驗(yàn)。對(duì)數(shù)據(jù)集D內(nèi)的各測(cè)量值進(jìn)行遍歷處理,判斷數(shù)據(jù)是否符合格式要求,對(duì)不符合要求的數(shù)據(jù)進(jìn)行修正或剔除,同時(shí)判斷測(cè)量值是否在理論區(qū)間,對(duì)于超過(guò)傳感器理論量程的數(shù)據(jù)予以剔除,輸出數(shù)據(jù)集D1。

    (3)異常值檢測(cè)。完成常規(guī)校驗(yàn)后,以數(shù)據(jù)集D1作為箱形圖輸入,統(tǒng)計(jì)該時(shí)段內(nèi)的數(shù)據(jù)分布規(guī)律,對(duì)于異常值進(jìn)行剔除并使用臨近數(shù)據(jù)的平均值進(jìn)行插補(bǔ),保證數(shù)據(jù)序列隨時(shí)間的變化特征,輸出數(shù)據(jù)集D2。

    (4)初始化卡爾曼濾波算法的參數(shù)??紤]到浮標(biāo)監(jiān)測(cè)系統(tǒng)噪聲由傳感器元件自身的特性所決定的,如風(fēng)速采樣頻率為3 s,浮標(biāo)觀測(cè)系統(tǒng)相對(duì)穩(wěn)定,相鄰間隔的測(cè)量值變化不大,因此根據(jù)k—1 時(shí)刻的最優(yōu)估計(jì)值來(lái)預(yù)測(cè)k時(shí)刻的預(yù)測(cè)值時(shí),我們由經(jīng)驗(yàn)認(rèn)為,測(cè)量噪聲和系統(tǒng)噪聲近似為高斯噪聲,變量Φ和H初始設(shè)置為1,方差初始設(shè)置為q=0.05和r=0.50。

    (5)去噪聲處理。采用改進(jìn)的自適應(yīng)卡爾曼濾波方法進(jìn)行去噪聲處理,以數(shù)據(jù)集D2各測(cè)量值按時(shí)間序列作為算法輸入,輸出預(yù)處理后的數(shù)據(jù)集D3。

    3 應(yīng)用與分析

    3.1 海洋工程應(yīng)用

    2021年6月28日至2021年8月16日,本研究數(shù)據(jù)預(yù)處理方法的海洋浮標(biāo)在廣東汕尾海域進(jìn)行實(shí)際工程應(yīng)用,該浮標(biāo)搭載2 臺(tái)風(fēng)速風(fēng)向儀、1 臺(tái)聲學(xué)多普勒流速剖面儀(ADCP)、1 臺(tái)波浪儀進(jìn)行海洋風(fēng)、水溫、海流、波浪的監(jiān)測(cè),實(shí)時(shí)回傳數(shù)據(jù)為廣東某地的系統(tǒng)工程建設(shè)提供數(shù)據(jù)支撐。站點(diǎn)位置與傳感器配置情況詳見(jiàn)表1。

    表1 浮標(biāo)系統(tǒng)部署與配置情況Table 1 Deployment and configuration of buoy system

    3.2 應(yīng)用結(jié)果與分析

    3.2.1 系統(tǒng)運(yùn)行穩(wěn)定性 傳感器完成數(shù)據(jù)采集后實(shí)時(shí)發(fā)送至浮標(biāo)端的數(shù)據(jù)采集與控制系統(tǒng),由該系統(tǒng)按照數(shù)據(jù)融合/回傳時(shí)間,提取原始數(shù)據(jù)進(jìn)行異常值修正和濾波降噪處理,并同步存儲(chǔ)在浮標(biāo)本地。海上實(shí)驗(yàn)結(jié)果表明,在整個(gè)監(jiān)測(cè)周期內(nèi),部署本研究自動(dòng)化數(shù)據(jù)預(yù)處理程序的浮標(biāo)系統(tǒng),數(shù)據(jù)采集率達(dá)100%,云數(shù)據(jù)中心數(shù)據(jù)接收率達(dá)97%。因此,本研究提出的自動(dòng)化數(shù)據(jù)預(yù)處理程序計(jì)算開(kāi)銷(xiāo)較低且不影響浮標(biāo)系統(tǒng)的正常采集和數(shù)據(jù)融合工作,系統(tǒng)整體運(yùn)行穩(wěn)定。

    3.2.2 常規(guī)檢驗(yàn)及異常值檢測(cè)情況 2021 年8 月16日回收浮標(biāo)后,將原始數(shù)據(jù)和處理后的數(shù)據(jù)下載后進(jìn)行統(tǒng)計(jì)分析。2021 年7 月18 日臺(tái)風(fēng)“查帕卡”期間浮標(biāo)監(jiān)測(cè)的風(fēng)、海流流速數(shù)據(jù)異常值檢測(cè)前后對(duì)比情況如圖3所示。

    圖3 數(shù)據(jù)異常值檢測(cè)Fig.3 Data outlier detection

    圖3 顯示了數(shù)據(jù)異常值的檢測(cè)過(guò)程,由圖3(a、c、e)可知,風(fēng)、海流的原始數(shù)據(jù)中存在明顯的異常測(cè)量值,這是由于傳感器受海洋浮標(biāo)搖晃、電信號(hào)不穩(wěn)定及環(huán)境噪聲干擾產(chǎn)生的。由于海洋浮標(biāo)測(cè)量要素具有隨時(shí)間有序變化的特征,如圖3(b、d、f)所示,在到達(dá)數(shù)據(jù)融合時(shí)間后,通過(guò)常規(guī)檢驗(yàn)和箱形圖法能夠很好地統(tǒng)計(jì)區(qū)間范圍內(nèi)的數(shù)據(jù)分布情況和識(shí)別離群異常值。2021 年6 月28 日至8 月16 日期間,浮標(biāo)數(shù)據(jù)采集器未檢測(cè)到傳感器測(cè)量值超出理論閾值情況,但識(shí)別出部分?jǐn)?shù)據(jù)存在格式錯(cuò)誤和異常測(cè)量值。如表2 所示,浮標(biāo)數(shù)據(jù)采集器檢測(cè)到各傳感器錯(cuò)誤率均較低,最高的ADCP 錯(cuò)誤率為8.90‰,主要原因是串口通信存在一定的通信誤碼率,而剖面海流的測(cè)量參數(shù)較多且單組數(shù)據(jù)量較大,單個(gè)參數(shù)格式有誤則直接判定為格式錯(cuò)誤。浮標(biāo)數(shù)據(jù)采集器對(duì)格式錯(cuò)誤進(jìn)行初步修正后存儲(chǔ)至本地,并對(duì)箱型圖檢測(cè)出的異常測(cè)量值使用臨近數(shù)據(jù)求取平均值進(jìn)行插補(bǔ)。

    表2 常規(guī)檢驗(yàn)及異常值檢測(cè)情況Table 2 Routine inspection and detection of outliers

    3.2.3 卡爾曼濾波數(shù)據(jù)預(yù)處理與應(yīng)用分析 在完成常規(guī)校驗(yàn)和異常值檢測(cè)后,浮標(biāo)端將在每個(gè)融合/數(shù)據(jù)回傳時(shí)刻,提取該時(shí)間段內(nèi)的原始數(shù)據(jù)并采用改進(jìn)的自適應(yīng)卡爾曼濾波方法進(jìn)行去噪聲處理。如圖4所示,浮標(biāo)運(yùn)行期間水溫、海流流速和風(fēng)速的原始數(shù)據(jù)經(jīng)過(guò)濾波處理前后的對(duì)比。

    圖4 原始數(shù)據(jù)與卡爾曼濾波處理結(jié)果對(duì)比Fig.4 Comparison of raw data and Kalman filter processing results

    根據(jù)圖4,原始數(shù)據(jù)經(jīng)卡爾曼濾波方法處理后,有效降低了水溫、海流流速和風(fēng)速的數(shù)據(jù)離散度,數(shù)據(jù)曲線更為平滑且符合隨時(shí)間序列變化的規(guī)律,濾波效果良好。由于卡爾曼濾波算法“預(yù)測(cè)+修正”的特征,原始數(shù)據(jù)經(jīng)濾波后,不僅能夠消除測(cè)量噪聲,也過(guò)濾了原始測(cè)量數(shù)據(jù)的高頻信號(hào)。但部分海上施工作業(yè)時(shí)甚至只關(guān)注極值、最大值和海流流速等高頻信號(hào)的變化,使得該方法在部分海洋工程領(lǐng)域存在應(yīng)用局限性。因此,各融合時(shí)段的平均風(fēng)速、平均風(fēng)向、水溫、氣壓、氣溫、平均波高等需要數(shù)據(jù)融合的參數(shù),可直接采用濾波后的數(shù)據(jù)進(jìn)行算數(shù)或矢量平均;極大風(fēng)、最大風(fēng)、最大波高等極值統(tǒng)計(jì)參數(shù),則需要結(jié)合原始測(cè)量數(shù)據(jù)進(jìn)行計(jì)算。針對(duì)大剖面、長(zhǎng)序列的海流數(shù)據(jù),研究人員一般更關(guān)注高頻信號(hào)的變化以開(kāi)展中大尺度的海洋動(dòng)力過(guò)程研究,如內(nèi)波實(shí)時(shí)監(jiān)測(cè)與預(yù)警,通過(guò)計(jì)算原始數(shù)據(jù)與濾波后數(shù)據(jù)的殘差,獲取剖面海流流速的高頻信號(hào),進(jìn)而通過(guò)海流高頻變化特征實(shí)現(xiàn)浮標(biāo)端內(nèi)波的自動(dòng)化識(shí)別[20-21]。因此,本研究提出的改進(jìn)的自適應(yīng)卡爾曼方法在有效消除測(cè)量噪聲的同時(shí),也能夠間接應(yīng)用于海流等參數(shù)高頻信號(hào)的提取,驗(yàn)證了本研究方法的可行性。

    3.3 與其他方法的仿真對(duì)比

    為進(jìn)一步驗(yàn)證算法的準(zhǔn)確性、可靠性和抑制濾波發(fā)散能力,將本研究算法與標(biāo)準(zhǔn)卡爾曼濾波方法[8]、自適應(yīng)卡爾曼濾波方法[9]和改進(jìn)的Sage-Husa自適應(yīng)卡爾曼濾波方法[22]的數(shù)據(jù)處理效果進(jìn)行對(duì)比分析,運(yùn)行硬件為具有一定計(jì)算能力的低功耗數(shù)據(jù)采集設(shè)備(Raspberry Pi,運(yùn)行內(nèi)存2 GB)。如圖5(a、b)所示,該實(shí)驗(yàn)數(shù)據(jù)采用已知正弦信號(hào)作為模擬的實(shí)測(cè)值,添加高斯噪聲(sigma=0.5)和突變數(shù)據(jù)的正弦信號(hào)作為傳感器的測(cè)量值作為仿真實(shí)驗(yàn)的原始數(shù)據(jù),因?yàn)檎鎸?shí)值已知,通過(guò)將預(yù)處理后的數(shù)據(jù)與真實(shí)測(cè)量值進(jìn)行誤差比對(duì)分析,能夠很好地評(píng)價(jià)數(shù)據(jù)預(yù)處理效果和驗(yàn)證算法的穩(wěn)定性。3 種算法的數(shù)據(jù)預(yù)處理效果對(duì)比如圖5(c)所示。

    圖5 不同預(yù)處理方法仿真對(duì)比Fig.5 Simulation comparison of different preprocessing methods

    通過(guò)計(jì)算處理后數(shù)據(jù)與正弦信號(hào)真實(shí)值之間的均方根誤差(RMSE)、平均絕對(duì)誤差(MAE)、相關(guān)系數(shù)(R)及算法運(yùn)行耗時(shí)等指標(biāo),對(duì)比分析不同方法的處理效果。如表3 所示,文獻(xiàn)[8]方法受到數(shù)據(jù)噪聲及突變數(shù)據(jù)的影響,即噪聲過(guò)大及噪聲與理論建模值相比較小時(shí),隨著誤差累積出現(xiàn)了濾波發(fā)散現(xiàn)象。文獻(xiàn)[9]、文獻(xiàn)[22]及本研究方法能夠有效的降低異常值對(duì)卡爾曼濾波發(fā)散的影響,濾波處理后的數(shù)據(jù)與真實(shí)數(shù)據(jù)曲線擬合效果較好,且本文方法處理后的數(shù)據(jù)均方根誤差、平均絕對(duì)誤差最低,與原始數(shù)據(jù)集的相關(guān)系數(shù)達(dá)0.98,數(shù)據(jù)預(yù)處理結(jié)果的準(zhǔn)確性最優(yōu)。另外,本研究方法利用箱型圖法處理突變數(shù)據(jù)和限定方差上限,在有效降低異常值對(duì)濾波精度影響的同時(shí),能夠快速降低測(cè)量值的可信度,降低了算法的計(jì)算開(kāi)銷(xiāo),其計(jì)算耗時(shí)相比文獻(xiàn)[9]和文獻(xiàn)[22]降低了50%,計(jì)算耗時(shí)最少。

    表3 對(duì)添加噪聲及突變數(shù)據(jù)的預(yù)處理效果評(píng)價(jià)對(duì)比Table 3 Comparison of preprocessing performance evaluation on adding noise and mutation data

    4 結(jié)論

    結(jié)合海洋浮標(biāo)數(shù)據(jù)采集特征,本文提出了基于改進(jìn)自適應(yīng)卡爾曼濾波的自動(dòng)化數(shù)據(jù)預(yù)處理方法,主要結(jié)論如下:

    (1)該方法通過(guò)箱形圖進(jìn)行異常測(cè)量值的檢測(cè),并采用方差受限的方法解決測(cè)量噪聲導(dǎo)致的濾波發(fā)散問(wèn)題。

    (2)海上工程應(yīng)用結(jié)果表明,浮標(biāo)端常規(guī)數(shù)據(jù)采集器能夠滿足該方法的計(jì)算開(kāi)銷(xiāo),數(shù)據(jù)采集及自動(dòng)化數(shù)據(jù)預(yù)處理程序運(yùn)行穩(wěn)定,浮標(biāo)系統(tǒng)有效數(shù)據(jù)接收率達(dá)97%。根據(jù)統(tǒng)計(jì),常規(guī)檢驗(yàn)和箱形圖法能夠有效地識(shí)別格式錯(cuò)誤和離群異常值,且各時(shí)段測(cè)量數(shù)據(jù)經(jīng)過(guò)異常值修正和濾波降噪處理后,數(shù)據(jù)曲線更為平滑且符合隨時(shí)間序列變化的規(guī)律,其濾波處理后的數(shù)據(jù)可應(yīng)用于各測(cè)量要素的數(shù)據(jù)融合與高頻信號(hào)的間接提取。

    (3)與同類(lèi)方法的仿真對(duì)比發(fā)現(xiàn),本研究方法能夠有效的降低異常值對(duì)卡爾曼濾波發(fā)散的影響,濾波處理后的數(shù)據(jù)與真實(shí)數(shù)據(jù)曲線擬合效果較好,且計(jì)算耗時(shí)最小。

    本研究成果豐富了我國(guó)海洋浮標(biāo)數(shù)據(jù)預(yù)處理與質(zhì)控手段,可為其他領(lǐng)域的基于時(shí)間序列變化的數(shù)據(jù)處理提供參考和借鑒。

    猜你喜歡
    原始數(shù)據(jù)浮標(biāo)卡爾曼濾波
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    家教世界(2023年7期)2023-03-22 12:11:24
    受了委屈的浮標(biāo)君
    家教世界(2023年4期)2023-03-04 07:31:28
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    受特定變化趨勢(shì)限制的傳感器數(shù)據(jù)處理方法研究
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    全新Mentor DRS360 平臺(tái)借助集中式原始數(shù)據(jù)融合及直接實(shí)時(shí)傳感技術(shù)實(shí)現(xiàn)5 級(jí)自動(dòng)駕駛
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    基于擴(kuò)展卡爾曼濾波的PMSM無(wú)位置傳感器控制
    videos熟女内射| 午夜福利免费观看在线| 最新的欧美精品一区二区| 国产精品亚洲av一区麻豆| 男男h啪啪无遮挡| 韩国高清视频一区二区三区| 1024香蕉在线观看| 亚洲一区二区三区欧美精品| 一级片'在线观看视频| 99久久综合免费| 不卡一级毛片| 19禁男女啪啪无遮挡网站| 亚洲五月婷婷丁香| 美女扒开内裤让男人捅视频| 在线 av 中文字幕| 午夜福利乱码中文字幕| 1024视频免费在线观看| av一本久久久久| 欧美日韩av久久| 国产亚洲精品第一综合不卡| 黄色视频不卡| 精品国产一区二区久久| 看免费av毛片| 久久人人97超碰香蕉20202| 99re6热这里在线精品视频| 亚洲一区二区三区欧美精品| 黑人欧美特级aaaaaa片| 老司机靠b影院| 人妻一区二区av| 嫁个100分男人电影在线观看| √禁漫天堂资源中文www| 成年人午夜在线观看视频| 无遮挡黄片免费观看| 成人亚洲精品一区在线观看| 国产97色在线日韩免费| 国产老妇伦熟女老妇高清| 丝袜喷水一区| 欧美精品人与动牲交sv欧美| 黑人猛操日本美女一级片| 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区激情| 啪啪无遮挡十八禁网站| 中文字幕最新亚洲高清| 另类精品久久| 成人免费观看视频高清| 欧美黄色淫秽网站| 9热在线视频观看99| 国产真人三级小视频在线观看| 大陆偷拍与自拍| 精品一区二区三区四区五区乱码| 精品少妇内射三级| tube8黄色片| 精品国产一区二区三区久久久樱花| 久久这里只有精品19| 在线观看免费高清a一片| 18禁观看日本| 女性生殖器流出的白浆| 成年女人毛片免费观看观看9 | 免费高清在线观看视频在线观看| 久久99热这里只频精品6学生| 欧美亚洲日本最大视频资源| 亚洲国产日韩一区二区| 精品熟女少妇八av免费久了| 国产免费视频播放在线视频| 国产野战对白在线观看| cao死你这个sao货| 多毛熟女@视频| 国产精品一区二区精品视频观看| 国产精品二区激情视频| 国产精品成人在线| 少妇猛男粗大的猛烈进出视频| 亚洲专区中文字幕在线| 考比视频在线观看| 久9热在线精品视频| 亚洲色图综合在线观看| 新久久久久国产一级毛片| 19禁男女啪啪无遮挡网站| 在线永久观看黄色视频| 天堂8中文在线网| 精品少妇久久久久久888优播| 亚洲成人手机| 咕卡用的链子| 欧美 亚洲 国产 日韩一| 黄网站色视频无遮挡免费观看| 成年av动漫网址| 国产主播在线观看一区二区| 宅男免费午夜| 一本一本久久a久久精品综合妖精| 精品久久久久久久毛片微露脸 | 国产欧美日韩综合在线一区二区| 亚洲国产欧美在线一区| 天堂8中文在线网| 国产精品久久久久成人av| 少妇 在线观看| 伊人久久大香线蕉亚洲五| 国产av一区二区精品久久| 90打野战视频偷拍视频| 久久精品亚洲av国产电影网| 精品欧美一区二区三区在线| 91精品国产国语对白视频| www.自偷自拍.com| 69av精品久久久久久 | 免费日韩欧美在线观看| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 国产免费福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 99热全是精品| 国产免费现黄频在线看| 欧美日韩亚洲国产一区二区在线观看 | 日本五十路高清| 少妇裸体淫交视频免费看高清 | 12—13女人毛片做爰片一| 色94色欧美一区二区| 亚洲免费av在线视频| 亚洲自偷自拍图片 自拍| 精品亚洲乱码少妇综合久久| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 自线自在国产av| 欧美激情高清一区二区三区| 男女高潮啪啪啪动态图| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 大香蕉久久网| 老司机福利观看| av免费在线观看网站| a级片在线免费高清观看视频| 久久久国产成人免费| 精品卡一卡二卡四卡免费| 欧美精品高潮呻吟av久久| 啪啪无遮挡十八禁网站| 国产免费视频播放在线视频| 精品一品国产午夜福利视频| 狠狠婷婷综合久久久久久88av| 18禁黄网站禁片午夜丰满| 国产三级黄色录像| 老汉色av国产亚洲站长工具| 91成年电影在线观看| 久久久久久久久久久久大奶| 国产区一区二久久| 久久 成人 亚洲| 欧美亚洲日本最大视频资源| 久久av网站| 18在线观看网站| 高清av免费在线| 高潮久久久久久久久久久不卡| 老熟妇乱子伦视频在线观看 | 欧美在线一区亚洲| 亚洲综合色网址| 国产成人a∨麻豆精品| 999精品在线视频| 国产又爽黄色视频| videosex国产| 中文欧美无线码| 淫妇啪啪啪对白视频 | 久久精品人人爽人人爽视色| 久久九九热精品免费| 麻豆国产av国片精品| 成人18禁高潮啪啪吃奶动态图| 黄网站色视频无遮挡免费观看| 亚洲欧美色中文字幕在线| 国产精品久久久av美女十八| 肉色欧美久久久久久久蜜桃| 亚洲黑人精品在线| 久久久久久亚洲精品国产蜜桃av| 午夜91福利影院| 妹子高潮喷水视频| 国产一卡二卡三卡精品| 男女边摸边吃奶| 丰满人妻熟妇乱又伦精品不卡| 午夜激情av网站| av欧美777| 满18在线观看网站| 狂野欧美激情性xxxx| 亚洲欧美清纯卡通| 精品亚洲成a人片在线观看| 亚洲av日韩在线播放| 亚洲国产欧美在线一区| 免费日韩欧美在线观看| 国产真人三级小视频在线观看| 性色av乱码一区二区三区2| 777米奇影视久久| 侵犯人妻中文字幕一二三四区| 啦啦啦免费观看视频1| 午夜福利免费观看在线| 亚洲国产欧美一区二区综合| 老司机午夜十八禁免费视频| 午夜成年电影在线免费观看| 美国免费a级毛片| 欧美 日韩 精品 国产| 欧美午夜高清在线| 免费观看人在逋| 日本vs欧美在线观看视频| 久久青草综合色| 欧美成人午夜精品| avwww免费| 亚洲性夜色夜夜综合| 日本一区二区免费在线视频| 精品第一国产精品| 老司机在亚洲福利影院| 啦啦啦免费观看视频1| 亚洲国产精品一区二区三区在线| 极品人妻少妇av视频| 美国免费a级毛片| 少妇被粗大的猛进出69影院| 亚洲第一欧美日韩一区二区三区 | 亚洲精品第二区| 国产欧美亚洲国产| av一本久久久久| 精品亚洲乱码少妇综合久久| 亚洲五月婷婷丁香| 亚洲一区中文字幕在线| 亚洲性夜色夜夜综合| 久久午夜综合久久蜜桃| 巨乳人妻的诱惑在线观看| 制服诱惑二区| 久久国产精品影院| 狠狠婷婷综合久久久久久88av| 人人澡人人妻人| 亚洲成人免费av在线播放| 91国产中文字幕| av不卡在线播放| 一本一本久久a久久精品综合妖精| 欧美黄色淫秽网站| 日韩制服丝袜自拍偷拍| 亚洲人成电影观看| 老熟妇乱子伦视频在线观看 | 制服人妻中文乱码| 美女午夜性视频免费| 午夜精品久久久久久毛片777| 成年人午夜在线观看视频| 国产精品九九99| 欧美中文综合在线视频| 丝瓜视频免费看黄片| 中国美女看黄片| 91九色精品人成在线观看| 美女高潮喷水抽搐中文字幕| 男女高潮啪啪啪动态图| 美女中出高潮动态图| 视频区图区小说| 国产精品免费视频内射| 成年美女黄网站色视频大全免费| 中文字幕高清在线视频| 欧美午夜高清在线| 可以免费在线观看a视频的电影网站| 亚洲精品国产区一区二| 18在线观看网站| 亚洲自偷自拍图片 自拍| 亚洲成av片中文字幕在线观看| 少妇的丰满在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 亚洲自偷自拍图片 自拍| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 久久久欧美国产精品| 国产97色在线日韩免费| av线在线观看网站| 欧美日韩福利视频一区二区| 秋霞在线观看毛片| 国产成人精品久久二区二区91| 国产1区2区3区精品| 国产精品久久久人人做人人爽| 十分钟在线观看高清视频www| 欧美日韩黄片免| 电影成人av| 丰满少妇做爰视频| 十八禁网站免费在线| 亚洲综合色网址| 亚洲伊人色综图| 如日韩欧美国产精品一区二区三区| 曰老女人黄片| 亚洲久久久国产精品| 国产精品 国内视频| 国产精品熟女久久久久浪| 欧美精品高潮呻吟av久久| 一区二区三区精品91| 丝袜美足系列| 日韩欧美一区视频在线观看| tocl精华| 国产成人av激情在线播放| 亚洲 国产 在线| 老司机影院毛片| 18禁裸乳无遮挡动漫免费视频| www日本在线高清视频| 久久精品aⅴ一区二区三区四区| 最新的欧美精品一区二区| 脱女人内裤的视频| 亚洲国产成人一精品久久久| 国产一区二区三区在线臀色熟女 | 国产精品一区二区在线不卡| 欧美亚洲日本最大视频资源| 欧美日韩中文字幕国产精品一区二区三区 | 一边摸一边抽搐一进一出视频| 亚洲av电影在线进入| 精品一区二区三区av网在线观看 | 黄色视频不卡| 国产av国产精品国产| 91麻豆精品激情在线观看国产 | 亚洲精品美女久久久久99蜜臀| 一级片'在线观看视频| 男人舔女人的私密视频| 最近中文字幕2019免费版| 成人影院久久| 午夜福利视频在线观看免费| 又大又爽又粗| 亚洲欧美一区二区三区久久| 男人添女人高潮全过程视频| 亚洲一区二区三区欧美精品| 黄片播放在线免费| 老司机靠b影院| 亚洲av电影在线进入| av网站在线播放免费| 狂野欧美激情性bbbbbb| 99国产精品免费福利视频| 日本wwww免费看| 欧美国产精品va在线观看不卡| 久久毛片免费看一区二区三区| 一区在线观看完整版| 咕卡用的链子| 久热这里只有精品99| 麻豆国产av国片精品| 老司机亚洲免费影院| 欧美老熟妇乱子伦牲交| 99久久精品国产亚洲精品| 天天躁夜夜躁狠狠躁躁| 久久久国产成人免费| 妹子高潮喷水视频| 免费观看a级毛片全部| 激情视频va一区二区三区| 交换朋友夫妻互换小说| 欧美国产精品va在线观看不卡| 久久中文字幕一级| 欧美 亚洲 国产 日韩一| 99精国产麻豆久久婷婷| 午夜日韩欧美国产| 国产成人精品在线电影| 咕卡用的链子| 天堂中文最新版在线下载| 亚洲成人手机| 91麻豆精品激情在线观看国产 | 一区二区三区四区激情视频| 精品人妻1区二区| 久久人妻熟女aⅴ| 久久影院123| 久久人人97超碰香蕉20202| 亚洲av日韩在线播放| 精品一区二区三区av网在线观看 | 午夜两性在线视频| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片| 午夜激情av网站| 日韩,欧美,国产一区二区三区| 最新在线观看一区二区三区| 国产精品1区2区在线观看. | 成年女人毛片免费观看观看9 | 男女边摸边吃奶| 老熟妇乱子伦视频在线观看 | 男男h啪啪无遮挡| 欧美黄色片欧美黄色片| 黑人欧美特级aaaaaa片| 国产精品二区激情视频| 日日摸夜夜添夜夜添小说| 在线精品无人区一区二区三| 久久青草综合色| 亚洲精品美女久久av网站| 久久青草综合色| 91精品三级在线观看| xxxhd国产人妻xxx| 久久久水蜜桃国产精品网| 午夜免费鲁丝| 女人久久www免费人成看片| 青春草亚洲视频在线观看| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 十八禁网站网址无遮挡| 一级毛片女人18水好多| 国产成人精品在线电影| 高清欧美精品videossex| 久久中文字幕一级| 免费日韩欧美在线观看| 交换朋友夫妻互换小说| 秋霞在线观看毛片| 波多野结衣一区麻豆| 成人av一区二区三区在线看 | 国产一区有黄有色的免费视频| 可以免费在线观看a视频的电影网站| 国产精品久久久av美女十八| 国产亚洲午夜精品一区二区久久| 国产91精品成人一区二区三区 | 免费在线观看影片大全网站| 久久精品久久久久久噜噜老黄| 成年av动漫网址| 国产一区二区在线观看av| 日本五十路高清| 一本—道久久a久久精品蜜桃钙片| 波多野结衣av一区二区av| 老司机午夜福利在线观看视频 | 久久久久久久国产电影| 久久久久久久久久久久大奶| 免费一级毛片在线播放高清视频 | 水蜜桃什么品种好| 日韩,欧美,国产一区二区三区| 国产亚洲欧美精品永久| 另类精品久久| 99国产精品一区二区蜜桃av | 免费高清在线观看视频在线观看| 一本久久精品| 亚洲av美国av| 国产成人欧美在线观看 | 欧美黑人精品巨大| 性少妇av在线| 在线观看www视频免费| 亚洲av日韩精品久久久久久密| svipshipincom国产片| 国产真人三级小视频在线观看| 成人三级做爰电影| 另类亚洲欧美激情| 我的亚洲天堂| 久久国产精品大桥未久av| 又黄又粗又硬又大视频| 69av精品久久久久久 | 国产伦理片在线播放av一区| 精品国内亚洲2022精品成人 | 黑人巨大精品欧美一区二区蜜桃| 可以免费在线观看a视频的电影网站| 国产成人系列免费观看| 大片电影免费在线观看免费| 国产精品久久久人人做人人爽| 亚洲情色 制服丝袜| 国产精品九九99| 国产成人a∨麻豆精品| 最近中文字幕2019免费版| 老司机亚洲免费影院| av网站免费在线观看视频| 欧美成人午夜精品| 视频区图区小说| 97在线人人人人妻| 2018国产大陆天天弄谢| 日本精品一区二区三区蜜桃| 亚洲欧美日韩高清在线视频 | 美女主播在线视频| 欧美日韩亚洲高清精品| 性少妇av在线| 久久狼人影院| 青青草视频在线视频观看| 999久久久精品免费观看国产| 欧美精品一区二区免费开放| 99国产极品粉嫩在线观看| 大码成人一级视频| 黄片小视频在线播放| 久久 成人 亚洲| 肉色欧美久久久久久久蜜桃| 亚洲欧美日韩另类电影网站| 久久香蕉激情| 欧美av亚洲av综合av国产av| 国产精品一区二区在线不卡| videosex国产| 五月天丁香电影| 中国美女看黄片| 欧美日韩亚洲高清精品| 狂野欧美激情性xxxx| 久久毛片免费看一区二区三区| 久久青草综合色| 国内毛片毛片毛片毛片毛片| 黄色视频在线播放观看不卡| 久久久国产成人免费| 精品亚洲成a人片在线观看| www.精华液| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 人成视频在线观看免费观看| 国产成人欧美| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 亚洲 欧美一区二区三区| 香蕉丝袜av| 少妇 在线观看| 免费在线观看黄色视频的| 天天躁夜夜躁狠狠躁躁| bbb黄色大片| 亚洲精品国产av蜜桃| 国产精品.久久久| tocl精华| 久久人人爽人人片av| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 亚洲精品美女久久av网站| 国产精品一区二区在线不卡| 宅男免费午夜| 久久人人爽av亚洲精品天堂| 一二三四社区在线视频社区8| av在线播放精品| 久久女婷五月综合色啪小说| 国产1区2区3区精品| 日韩 欧美 亚洲 中文字幕| 操出白浆在线播放| 久久久久久免费高清国产稀缺| 亚洲av欧美aⅴ国产| 岛国在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 在线观看舔阴道视频| 成在线人永久免费视频| 高清欧美精品videossex| 成在线人永久免费视频| 国产片内射在线| 久久毛片免费看一区二区三区| 成人三级做爰电影| 精品国产一区二区三区四区第35| 欧美一级毛片孕妇| 搡老岳熟女国产| 欧美中文综合在线视频| 日本91视频免费播放| 法律面前人人平等表现在哪些方面 | 极品少妇高潮喷水抽搐| 中文字幕色久视频| 日韩 欧美 亚洲 中文字幕| www.精华液| 亚洲视频免费观看视频| 久久综合国产亚洲精品| h视频一区二区三区| 国产成人啪精品午夜网站| 无遮挡黄片免费观看| 免费不卡黄色视频| 日韩制服丝袜自拍偷拍| 国产成人一区二区三区免费视频网站| 国产高清视频在线播放一区 | 国产无遮挡羞羞视频在线观看| 久久精品久久久久久噜噜老黄| 久久久久网色| 国产黄色免费在线视频| 日韩中文字幕欧美一区二区| 亚洲美女黄色视频免费看| 黑人巨大精品欧美一区二区mp4| 亚洲精品中文字幕一二三四区 | 国产99久久九九免费精品| 欧美老熟妇乱子伦牲交| 久久久久精品国产欧美久久久 | 嫁个100分男人电影在线观看| 免费在线观看影片大全网站| 后天国语完整版免费观看| 国产视频一区二区在线看| 一区二区三区精品91| 国产视频一区二区在线看| 亚洲人成电影观看| 黄网站色视频无遮挡免费观看| 欧美日韩一级在线毛片| 欧美+亚洲+日韩+国产| 一本—道久久a久久精品蜜桃钙片| 成年动漫av网址| 久久久久久久久免费视频了| 欧美黑人欧美精品刺激| 久久久久久久久免费视频了| 午夜福利免费观看在线| 亚洲va日本ⅴa欧美va伊人久久 | 欧美国产精品一级二级三级| 久久中文看片网| xxxhd国产人妻xxx| 午夜福利免费观看在线| 一区二区三区乱码不卡18| 亚洲一码二码三码区别大吗| 777久久人妻少妇嫩草av网站| 桃红色精品国产亚洲av| 蜜桃国产av成人99| 热99re8久久精品国产| 久久国产精品男人的天堂亚洲| 久久久国产一区二区| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 日韩熟女老妇一区二区性免费视频| 亚洲va日本ⅴa欧美va伊人久久 | 搡老乐熟女国产| 国产亚洲欧美在线一区二区| 亚洲精品国产av成人精品| 狠狠婷婷综合久久久久久88av| 老司机在亚洲福利影院| 亚洲五月婷婷丁香| 黑人巨大精品欧美一区二区蜜桃| 大片电影免费在线观看免费| 80岁老熟妇乱子伦牲交| av在线播放精品| 久久国产精品大桥未久av| 午夜激情av网站| 中文字幕人妻熟女乱码| 欧美在线黄色| 免费在线观看视频国产中文字幕亚洲 | bbb黄色大片| 国产xxxxx性猛交| 熟女少妇亚洲综合色aaa.| 欧美黄色片欧美黄色片| 亚洲色图综合在线观看| 美女高潮到喷水免费观看| 一区二区三区乱码不卡18| 中国美女看黄片| 99久久人妻综合| 法律面前人人平等表现在哪些方面 | 青草久久国产| 午夜精品国产一区二区电影| 操美女的视频在线观看| 80岁老熟妇乱子伦牲交| 青青草视频在线视频观看| 国内毛片毛片毛片毛片毛片| 黑人巨大精品欧美一区二区mp4| 亚洲人成电影免费在线| 嫩草影视91久久| 老司机午夜福利在线观看视频 | 每晚都被弄得嗷嗷叫到高潮| 日本撒尿小便嘘嘘汇集6| 国产xxxxx性猛交| 国产日韩一区二区三区精品不卡| 成人亚洲精品一区在线观看| www日本在线高清视频|