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

    基于希爾伯特變換的兩棲車輛非線性橫搖運(yùn)動(dòng)模型辨識(shí)研究

    2016-08-05 03:36:05馬新謀郝一凝常列珍潘玉田
    關(guān)鍵詞:包絡(luò)線

    馬新謀, 郝一凝, 常列珍, 潘玉田

    (1. 中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051; 2. 中北大學(xué) 理學(xué)院, 山西 太原 030051)

    ?

    基于希爾伯特變換的兩棲車輛非線性橫搖運(yùn)動(dòng)模型辨識(shí)研究

    馬新謀1, 郝一凝1, 常列珍2, 潘玉田1

    (1. 中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051; 2. 中北大學(xué) 理學(xué)院, 山西 太原 030051)

    摘要:為了識(shí)別兩棲車輛非線性橫搖運(yùn)動(dòng)方程, 首先采用希爾伯特變換法把自由橫搖衰減信號(hào)變換成解析信號(hào), 該解析信號(hào)包含了自由橫搖運(yùn)動(dòng)的瞬時(shí)幅值信息和瞬時(shí)頻率信息, 并給出了兩棲車輛非線性橫搖運(yùn)動(dòng)方程阻尼函數(shù)和恢復(fù)函數(shù)與瞬時(shí)幅值和頻率的關(guān)系式, 最后由最小二乘法擬合擬合即可獲得兩棲車輛非線性橫搖運(yùn)動(dòng)方程的非線性阻尼項(xiàng)和恢復(fù)項(xiàng)系數(shù). 應(yīng)用所提出的方法針對(duì)數(shù)值仿真實(shí)例和某兩棲車輛1∶4模型自由橫搖衰減實(shí)驗(yàn)數(shù)據(jù), 成功地估計(jì)出了非線性參數(shù), 估計(jì)值與真值及實(shí)驗(yàn)值吻合很好. 估計(jì)結(jié)果表明所提出的方法可以用于估計(jì)兩棲車輛非線性橫搖運(yùn)動(dòng)方程, 而且精度滿足工程需求.

    關(guān)鍵詞:希爾伯特變換法; 兩棲車輛; 非線性橫搖運(yùn)動(dòng); 包絡(luò)線; 瞬時(shí)頻率

    兩棲車輛在水上的運(yùn)動(dòng)為六自由度的運(yùn)動(dòng), 其中橫搖運(yùn)動(dòng)與船舶橫搖運(yùn)動(dòng)相似, 是一個(gè)非線性振動(dòng)學(xué)問題, 因此兩棲車輛非線性橫搖運(yùn)動(dòng)也成為影響兩棲車輛水上安全性、 成員舒適性的重要因素[1-3]. 眾所周知, 船舶非線性橫搖是影響船舶性能與安全性的重要因素, 特別是在突風(fēng)載荷、 波浪載荷和流聯(lián)合激勵(lì)下, 船舶會(huì)發(fā)生奇異傾覆[4-9]. 所以有必要針對(duì)兩棲車輛非線性橫搖運(yùn)動(dòng)進(jìn)行深入研究. 要開展兩棲車輛非線性橫搖運(yùn)動(dòng)特性的首要任務(wù)是建立非線性橫搖運(yùn)動(dòng)方程. 到目前為止, 研究船舶橫搖運(yùn)動(dòng)的工作很多, 常用來建立船舶非線性橫搖運(yùn)動(dòng)方程的方法有經(jīng)驗(yàn)公式法[5]和船模自由橫搖衰減試驗(yàn)法[4,6-7,9]兩種. 對(duì)于經(jīng)驗(yàn)公式法僅對(duì)特定的船型適用; 而船模自由橫搖衰減試驗(yàn)法適用性強(qiáng), 應(yīng)用范圍廣, 由此演化而來的建立船舶非線性橫搖運(yùn)動(dòng)方程的方法有最小二乘法[2]、 傅汝德能量法[1]、 改進(jìn)的能量法[1,4,6]、 遺傳算法[6]、 反問題法[7]、 消滅線法、 支持向量機(jī)法[9]等, 這些方法都要從自由橫搖衰減曲線來識(shí)別非線性橫搖運(yùn)動(dòng)方程的參數(shù).

    另一方面, 希爾伯特變換(Hilberttransform,HT)在信號(hào)處理、 故障診斷和參數(shù)識(shí)別領(lǐng)域都有了廣泛的應(yīng)用, 其中關(guān)于參數(shù)辨識(shí)領(lǐng)域的文獻(xiàn)也較多.Feldman[10-11]采用希爾伯特變換對(duì)非線性系統(tǒng)的各項(xiàng)參數(shù)進(jìn)行了識(shí)別研究, 并針對(duì)一般形式的運(yùn)動(dòng)方程提出了一種自由振動(dòng)分析方法, 是近期關(guān)于希爾伯特變換的經(jīng)典之作.Shin等[12]把希爾伯特變換在信號(hào)處理領(lǐng)域進(jìn)行了系統(tǒng)的闡述, 給出了不同特征信號(hào)的處理技巧.Ignacio等[13]成功用希爾伯特變換估計(jì)了鐵路橋振動(dòng)系統(tǒng)的模型參數(shù).Shi等[14]用希爾伯特變換和經(jīng)驗(yàn)?zāi)B(tài)公式法從強(qiáng)迫振動(dòng)信號(hào)里估計(jì)出了多自由度模型的參數(shù), 表明希爾伯特變換法不僅對(duì)單自由度系統(tǒng)適用, 對(duì)多自由度系統(tǒng)也適用.Kim[15]用希爾伯特變換法對(duì)浮式生產(chǎn)儲(chǔ)/卸油裝置(FPSO)的自由橫搖衰減曲線進(jìn)行了處理, 成功辨識(shí)了FPSO橫搖運(yùn)動(dòng)方程的非線性阻尼項(xiàng)和非線性恢復(fù)項(xiàng)系數(shù), 精度滿足工程需求. 從這些成功的范例來看, 希爾伯特變換法也可以用來估計(jì)兩棲車輛的非線性橫搖運(yùn)動(dòng)模型參數(shù).

    基于上述原因, 本文借鑒船舶橫搖運(yùn)動(dòng)的研究成果和希爾伯特變換理論, 擬采用希爾伯特變換法來處理兩棲車輛自由橫搖衰減曲線以估計(jì)兩棲車輛非線性橫搖運(yùn)動(dòng)方程的參數(shù).

    1非線性橫搖方程參數(shù)估計(jì)理論

    借鑒船舶非線性橫搖運(yùn)動(dòng)研究理論, 不考慮與其他方向運(yùn)動(dòng)的耦合效應(yīng), 則靜水里的兩棲車輛單自由度的非線性橫搖運(yùn)動(dòng)方程可以用二階常微分方程[1]表示為

    (1)

    為了方便描述, 式(1)可以改寫為

    (2)

    (3)

    (4)

    采用希爾伯特變換法估計(jì)非線性橫搖參數(shù)的基本理論是建立在如下的解析信號(hào)基礎(chǔ)上.

    (5)

    (6)

    式中:pv為柯西主值.

    將式(5)改寫成幅值和相位表示的復(fù)指數(shù)形式

    (7)

    式中:A(t)是瞬時(shí)包絡(luò)線; ψ(t)是瞬時(shí)相位角.

    由歐拉公式可知, 存在如下關(guān)系

    (8)

    (9)

    而且

    (10)

    (11)

    (12)

    (13)

    式中:ω(t)是瞬時(shí)頻率, 是時(shí)間的實(shí)函數(shù).

    考慮到式(12)和(13), 計(jì)算式(5)的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)

    (14)

    (15)

    對(duì)式(2)進(jìn)行希爾伯特變換, 可得

    (16)

    將式(16) 乘以虛數(shù)單位i再加上式(2)可得用解析信號(hào)Y(t)表達(dá)的微分方程

    (17)

    將式(14)和(15)代入式(17), 合并實(shí)部和虛部可得

    (18)

    要使式(18)恒成立, 需要實(shí)部和虛部的系數(shù)恒為零, 即

    (19)

    (20)

    2數(shù)值仿真實(shí)驗(yàn)

    在獲得上述信息的基礎(chǔ)上, 如圖 1 所示, 存在邊界效應(yīng), 為了使計(jì)算方便和可靠, 可以截取信號(hào)質(zhì)量較好的中間部分, 再根據(jù)式(10)和(13)可以計(jì)算瞬時(shí)包絡(luò)A(t)和瞬時(shí)頻率ω(t), 如圖 2 和圖 3 所示. 進(jìn)而可以計(jì)算A(t)和ω(t)的一階導(dǎo)數(shù)和二階導(dǎo)數(shù).

    圖 1 橫搖位移信號(hào)曲線Fig.1 Signal curve of roll motion

    圖 2 自由橫搖曲線和包絡(luò)線Fig.2 Free decay curve and envelope

    圖 3 瞬時(shí)頻率曲線Fig.3 Instantaneous frequency curve

    圖 4 阻尼函數(shù)和擬合曲線Fig.4 Damping coefficients function and fitted curve

    圖 5 頻率函數(shù)和擬合曲線Fig.5 Frequency function and fitted curve

    估計(jì)出的阻尼項(xiàng)和恢復(fù)項(xiàng)系數(shù)d1, d2, k1和k3如表 1 所示. 由表 1 可知, 采用希爾伯特變換法和最小二乘擬合相結(jié)合可以快速有效地估計(jì)兩棲車輛非線性橫搖運(yùn)動(dòng)方程的系數(shù), 最大估計(jì)參數(shù)的相對(duì)誤差為0.313%.

    表 1 參數(shù)估計(jì)結(jié)果表

    將估計(jì)的系數(shù)代入式(1), 用四階龍格庫塔法求解, 與原系數(shù)的結(jié)果對(duì)比如圖 6 所示.

    圖 6 仿真和估計(jì)的衰減曲線Fig.6 Simulated decay curve and estimated decay curve

    圖 6 表明估計(jì)的角位移曲線與原數(shù)據(jù)的角位移曲線高度吻合. 因此可以證明用所提出估計(jì)非線性橫搖運(yùn)動(dòng)方程的方法可行性, 精度滿足工程需求.

    需要特別指出, 上面仿真實(shí)例所生成的自由橫搖衰減曲線是光滑的, 不含噪聲, 因此估計(jì)的結(jié)果精度很高. 如果自由橫搖衰減曲線包含噪聲, 則精度會(huì)有所下降. 因此在使用該方法時(shí), 必須先對(duì)測試所得的自由橫搖信號(hào)進(jìn)行濾波和平滑, 再應(yīng)用本文提出的方法估計(jì).

    3希爾伯特變換法估計(jì)橫搖運(yùn)動(dòng)方程參數(shù)

    3.1兩棲車輛自由橫搖衰減實(shí)驗(yàn)

    為了研究兩棲車輛的非線性橫搖運(yùn)動(dòng), 根據(jù)某兩棲車輛的技術(shù)參數(shù), 設(shè)計(jì)了一個(gè)1∶4的縮小模型, 使其滿足幾何相似、 運(yùn)動(dòng)相似和動(dòng)力相似, 如圖 7 所示. 在試驗(yàn)水池里進(jìn)行自由橫搖衰減實(shí)驗(yàn), 采樣頻率100Hz, 測得的自由橫搖衰減曲線如圖 8 所示.

    圖 7 兩棲車輛拖模(1∶4)Fig.7 Scaled amphibious vehicle model(1∶4)

    為了減小估計(jì)誤差, 必須先對(duì)測試獲取的自由橫搖衰減曲線上的噪聲進(jìn)行濾波和平滑處理. 由圖8可知, 前面和后面各有一段不規(guī)則的數(shù)據(jù), 為了估計(jì)方便可以預(yù)先把自由橫搖衰減曲線掐頭去尾, 并認(rèn)為起始時(shí)間為0, 以避免出現(xiàn)邊界效應(yīng).

    圖 8 兩棲車輛拖模自由橫搖衰減曲線Fig.8 Test free rolling decay curve of amphibious vehicle model

    在此基礎(chǔ)上進(jìn)行濾波, 濾波時(shí)分為兩步. 第一步先用MATLAB里的butter函數(shù)設(shè)計(jì)一個(gè)Butterworth低通濾波器以去除高頻噪聲; 第二步進(jìn)行小波分解重構(gòu)法濾波, 小波基選擇“rbio6.8”, 分解尺度為2; 最后選用三點(diǎn)平滑公式

    (21)

    進(jìn)行數(shù)據(jù)平滑, 平滑后的自由橫搖衰減曲線如圖 9 所示.

    圖 9 濾波后的自由橫搖衰減曲線Fig.9 Free rolling decay curve after filter processing

    3.2橫搖運(yùn)動(dòng)參數(shù)估計(jì)

    設(shè)兩棲車輛非線性橫搖的控制方程如方程(1)所示, 根據(jù)上述自由橫搖衰減曲線估計(jì)d1, d2, k1和k3. 首先運(yùn)用希爾伯特變換法把平滑后的自由橫搖衰減曲線進(jìn)行變換, 構(gòu)造解析信號(hào), 如圖 10 所示.

    圖 10 濾波后的橫搖位移信號(hào)曲線Fig.10 Signal curve of roll motion after filter processing

    (22)

    4估計(jì)結(jié)果與討論

    由圖 11 可知, 估計(jì)橫搖角度和測試的橫搖角度吻合得較好. 與圖 6 的對(duì)比結(jié)果可知, 該估計(jì)方法可以估計(jì)非線性橫搖運(yùn)動(dòng)方程的參數(shù), 但是對(duì)信號(hào)的噪聲相應(yīng)比較敏感, 因此要提高使用該方法估計(jì)非線性橫搖參數(shù)的精度有兩個(gè)途徑:

    圖 11 測試和估計(jì)的衰減曲線Fig.11 Measured decay curve and estimated decay curve

    1) 提高測試系統(tǒng)的測試精度和信噪比, 使測試結(jié)果更準(zhǔn)確.

    2) 運(yùn)用硬件和軟件進(jìn)行信號(hào)濾波, 并對(duì)測試信號(hào)進(jìn)行平滑處理, 恢復(fù)理想的信號(hào).

    文獻(xiàn)[1]中, 馬新謀等采用能量法估計(jì)了兩棲車輛的非線性阻尼系數(shù), 用流體靜力學(xué)理論計(jì)算了兩棲車輛的GZ曲線, 并根據(jù)GZ曲線給了恢復(fù)項(xiàng)系數(shù), 最終給出的兩棲車輛非線性橫搖運(yùn)動(dòng)方程為

    (23)

    為了對(duì)比文獻(xiàn)[1]和本文估計(jì)的非線性橫搖方程的系數(shù), 把它們列入表中, 如表 2 所示.

    表 2 兩種方法的參數(shù)估計(jì)結(jié)果

    由表 2 可知, 兩種方法獲取的阻尼系數(shù)差距較大, 用能量法估計(jì)的阻尼值比希爾伯特法估計(jì)的小.GZ曲線擬合的線性恢復(fù)項(xiàng)系數(shù)和立方項(xiàng)恢復(fù)項(xiàng)系數(shù)的絕對(duì)值都比希爾伯特法估計(jì)的大, 但兩個(gè)方法估計(jì)的非線性橫搖運(yùn)動(dòng)方程都與實(shí)驗(yàn)數(shù)據(jù)吻合較好. 把文獻(xiàn)[1]估計(jì)的非線性方程預(yù)測和本文預(yù)測的自由橫搖衰減曲線和實(shí)驗(yàn)數(shù)據(jù)放在一起對(duì)比, 如圖 12 所示.

    由圖 12 可知, 本文估計(jì)的固有頻率相對(duì)于用GZ曲線法估計(jì)的固有頻率而言更接近于試驗(yàn)數(shù)據(jù). 分析原因是由于GZ曲線是根據(jù)設(shè)計(jì)圖紙通過三維建模, 通過流體靜力學(xué)計(jì)算出來的, 是圖定數(shù)據(jù), 為兩棲車輛水上橫搖運(yùn)動(dòng)的理論固有頻率; 而本文估計(jì)的固有頻率則是通過試驗(yàn)得到, 雖然該實(shí)驗(yàn)的拖模是滿足幾何相似、 運(yùn)動(dòng)相似和動(dòng)力相似的, 但是不可避免地存在誤差, 加之測試也存在噪聲誤差. 因此, 用希爾伯特變換法估計(jì)的非線性橫搖運(yùn)動(dòng)從固有頻率角度看更接近于實(shí)驗(yàn)值. 本文和文獻(xiàn)[1]中僅用實(shí)驗(yàn)數(shù)據(jù)估計(jì)了阻尼項(xiàng), 由圖12可知, 本文估計(jì)的結(jié)果在幅值上更接近測試數(shù)據(jù), 因而阻尼系數(shù)比文獻(xiàn)[1]給出的結(jié)果大更符合測試結(jié)果. 從固有頻率角度看誤差比本文的大是正常的. 但是, 兩種方法獲得的固有頻率非常接近, 無論用哪種方法來估計(jì)兩棲車輛非線性橫搖的恢復(fù)項(xiàng)都是可信的.

    圖 12 兩種方法所得測試和估計(jì)的衰減曲線Fig.12 Measured decay curve and estimated decay curve by two methods

    5結(jié)論

    通過上述分析, 可以得出以下結(jié)論:

    1) 希爾伯特變換法與最小二乘法相結(jié)合, 可以估計(jì)兩棲車輛水上非線性橫搖運(yùn)動(dòng)方程的系數(shù).

    2) 該方法不僅可以估計(jì)出阻尼項(xiàng), 也可以估計(jì)出恢復(fù)項(xiàng).

    3) 該方法對(duì)測試信號(hào)的噪聲敏感, 因此在使用該方法時(shí)必須先對(duì)自由橫搖衰減試驗(yàn)信號(hào)進(jìn)行濾波和平滑.

    4) 該方法也可以用來估計(jì)船舶非線性橫搖運(yùn)動(dòng)的模型, 也可以用來估計(jì)其他非線性動(dòng)力學(xué)系統(tǒng)的模型.

    參考文獻(xiàn):

    [1]馬新謀, 潘玉田, 常列珍. 兩棲車輛非線性橫搖阻尼系數(shù)估計(jì)[J]. 船舶力學(xué), 2013, 17(5): 488-493.

    MaXinmou,PanYutian,ChangLiezhen.Estimatednonlineardampingcoefficientsofamphibiousvehiclerollmotion[J].JournalofShipMechanics, 2013, 17(5): 488-493. (inChinese)

    [2]馬新謀, 常列珍, 侯宏花. 兩棲車輛非線性橫搖動(dòng)力學(xué)參數(shù)辨識(shí)研究[J]. 中北大學(xué)學(xué)報(bào)(自然版), 2015, 36(2): 134-139.

    MaXinmou,ChangLiezhen,HouHonghua.Aninvestigationondynamicparametersidentificationforamphibiousvehiclenonlinearrollmotion[J].JournalofNorthUniversityofChina(NaturalScienceEdition), 2015, 36(2): 134-139. (inChinese)

    [3]徐國英, 薛勁櫓. 兩棲車輛在波浪中航行時(shí)的搖蕩分析及解決方法[J]. 兵工學(xué)報(bào), 2010, 31(5): 541-546.

    XuGuoying,XueJinlu.Analysisandsolutionofamphibiousvehicle’stossmotioninwave[J].ActaArmamentarII, 2010, 31(5): 541-546. (inChinese)

    [4]RobertsJB.Estimationofnonlinearshiprolldampingfromfree-decaydata[J].JournalofShipResearch, 1985, 29(2):127-138.

    [5]ChakrabartiS.Empiricalcalculationofrolldampingforshipsandbarges[J].OceanEngineering, 2001, 28(7): 915-932.

    [6]KimihikoU,NobuoK,KiyoshiA.Estimationofcoefficientsoftheequationofnonlinearrollmotionforfishingboatsbyimprovedenergymethodandgeneticalgorithm[J].JanpanJ.Indust.Appl.Math., 2003, 20: 155-192.

    [7]JangTS,ChoiHS,HanSL.Anewmethodfordetectingnon-lineardampingandrestoringforcesinnon-linearoscillationsystemsfromtransientdata[J].InternationalJournalofNon-LinearMechanics, 2009, 44(7): 801-808.

    [8]TaylanM.Theeffectofnonlineardampingandrestoringinshiproll[J].OceanEngineering, 2000, 27(9): 921-932.

    [9]HouXR,ZouZJ.SVR-basedidentificationofnonlinearrollmotionequationforFPSOsinregularwaves[J].OceanEngineering, 2015, 109: 531-538.

    [10]FeldmanM.Hilberttransforminvibrationanalysis[J].MechanicalSystemsandSignalProcessing, 2011, 25(3): 735- 802.

    [11]FeldmanM.Hilberttransformmethodsfornonparametricidentificationofnonlineartimevaryingvibrationsystems[J].MechanicalSystemsandSignalProcessing, 2014, 47(3): 66 -77.

    [12]ShinK,HammondJ.Fundamentalsofsignalprocessingforsoundandvibrationengineers[M].NewJersey:JohnWiley&Sons,Ltd., 2008.

    [13]IgnacioG,RaidK.AnalysisoftheannualvariationsinthedynamicbehaviorofaballastedrailwaybridgeusingHilberttransform[J].EngineeringStructures, 2014, 60: 126-132.

    [14]ShiZY,LawSS,XuX.Identificationoflineartime-varyingMDOFdynamicsystemsfromforcedexcitationusingHilberttransformandEMDmethod[J].JournalofSoundandVibration. 2009, 321(3-5): 572-589.

    [15]KimY,ParkMJ.IdentificationofthenonlinearrolldampingandrestoringmomentofaFPSOusingHilberttransform[J].OceanEngineering. 2015, 109: 381-388.

    文章編號(hào):1673-3193(2016)04-0361-07

    收稿日期:2016-03-31

    基金項(xiàng)目:中北大學(xué)科學(xué)基金; 國防預(yù)研課題資助

    作者簡介:馬新謀(1978-), 男, 博士, 講師, 主要從事火炮自動(dòng)武器與彈藥工程的研究.

    中圖分類號(hào):TJ811+.6

    文獻(xiàn)標(biāo)識(shí)碼:A

    doi:10.3969/j.issn.1673-3193.2016.04.008

    Identification for the Nonlinear Roll Motion Model ofAmphibiousVehicleUsingHilbertTransform

    MA Xin-mou1, HAO Yi-ning1, CHANG Lie-zhen2, PAN Yu-tian1

    (1. School of Mechatronics Engineering, North University of China, Taiyuan 030051, China;2.SchoolofScience,NorthUniversityofChina,Taiyuan030051,China)

    Abstract:To identify the coefficients of nonlinear roll motion modelling of amphibious vehicle, the Hilbert Transform (HT) method is applied for a free roll decay signal in time domain to obtain an analytical signal. The analytical signal provides information on the instantaneous amplitude and instantaneous frequency. Then, the relationship between the instantaneous amplitude/frequency and the nonlinear parameters can be derived. Least square method is used to estimate unknown parameters of nonlinear roll motion of amphibious vehicle. The present method is applied to a numerical experiment data and then to real free roll decay data obtained from a 1∶4 scale model of an amphibious vehicle to identify nonlinear parameters. The obtained nonlinear damping coefficients and nonlinear restoring coefficients are good agreement with the true value and tested data. So, the present method can be used to identify nonlinear roll motion equation of amphibious vehicle with enough accuracy.

    Key words:Hilbert transform method; amphibious vehicle; nonlinear roll motion; envelope; instantaneous frequency

    猜你喜歡
    包絡(luò)線
    平面曲線運(yùn)動(dòng)中的包絡(luò)線問題探討
    平方反比有心力場中的3個(gè)拋體包絡(luò)線及其證明
    基于ISO 14692 標(biāo)準(zhǔn)的玻璃鋼管道應(yīng)力分析
    由橢圓張角為直角的弦所在直線形成的“包絡(luò)”
    拋體的包絡(luò)線方程的推導(dǎo)
    一種用于故障隔離的參數(shù)區(qū)間包絡(luò)線計(jì)算方法
    包絡(luò)線在經(jīng)濟(jì)學(xué)中的應(yīng)用
    中國商論(2016年22期)2016-11-30 01:29:35
    斜拉索振動(dòng)信號(hào)的包絡(luò)線稀疏復(fù)原算法
    EMD 包絡(luò)線擬合算法改進(jìn)及在泡沫尺寸趨勢提取中的應(yīng)用
    函數(shù)最值中的包絡(luò)線
    久久久精品欧美日韩精品| 国产综合懂色| 成人av在线播放网站| 国产黄a三级三级三级人| 手机成人av网站| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频| 综合色av麻豆| 在线观看66精品国产| 成人鲁丝片一二三区免费| 一个人免费在线观看电影 | 免费在线观看亚洲国产| 午夜免费观看网址| 亚洲欧美日韩东京热| 麻豆国产av国片精品| 成人特级av手机在线观看| 身体一侧抽搐| 成人亚洲精品av一区二区| 嫩草影院精品99| 悠悠久久av| 九色国产91popny在线| 免费人成视频x8x8入口观看| 少妇丰满av| 黑人欧美特级aaaaaa片| 亚洲精品久久国产高清桃花| www日本黄色视频网| www.熟女人妻精品国产| 免费人成视频x8x8入口观看| 神马国产精品三级电影在线观看| 欧美最黄视频在线播放免费| 日本免费a在线| 国内精品美女久久久久久| 精品久久久久久久久久免费视频| 18禁黄网站禁片午夜丰满| 国产精品久久久av美女十八| 中文字幕高清在线视频| 桃色一区二区三区在线观看| 99在线视频只有这里精品首页| 国产成人福利小说| av天堂在线播放| 国产aⅴ精品一区二区三区波| 精品久久久久久久毛片微露脸| 国产精品综合久久久久久久免费| 在线观看舔阴道视频| 久久性视频一级片| 欧美一区二区国产精品久久精品| 国产极品精品免费视频能看的| 国产高清三级在线| 精品国产乱码久久久久久男人| av天堂中文字幕网| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品粉嫩美女一区| 国产精品亚洲一级av第二区| 亚洲欧美激情综合另类| 欧美一区二区国产精品久久精品| 丁香欧美五月| 一级黄色大片毛片| 国产人伦9x9x在线观看| 俄罗斯特黄特色一大片| 一区福利在线观看| 老熟妇乱子伦视频在线观看| 欧美日本视频| 国产真实乱freesex| 欧美成人一区二区免费高清观看 | 成年女人永久免费观看视频| 18禁观看日本| 色老头精品视频在线观看| 欧美黄色淫秽网站| 中国美女看黄片| 免费av毛片视频| 精品一区二区三区视频在线观看免费| 少妇的逼水好多| 男插女下体视频免费在线播放| 午夜免费激情av| 免费大片18禁| 在线视频色国产色| 好看av亚洲va欧美ⅴa在| 亚洲第一欧美日韩一区二区三区| 黄频高清免费视频| 日韩三级视频一区二区三区| 午夜视频精品福利| 在线国产一区二区在线| 精品一区二区三区四区五区乱码| 国产精品女同一区二区软件 | 国产精品,欧美在线| 又粗又爽又猛毛片免费看| 欧美最黄视频在线播放免费| 黄色丝袜av网址大全| 女人被狂操c到高潮| av女优亚洲男人天堂 | 久久精品国产综合久久久| 国产av一区在线观看免费| 桃色一区二区三区在线观看| 18禁黄网站禁片午夜丰满| 国产精华一区二区三区| 欧美成狂野欧美在线观看| 国产午夜精品久久久久久| 国产精品久久久人人做人人爽| 免费电影在线观看免费观看| 久久精品国产清高在天天线| 日本免费一区二区三区高清不卡| 不卡一级毛片| 一个人看的www免费观看视频| 身体一侧抽搐| 黄色女人牲交| 天堂√8在线中文| 日本 欧美在线| 性色avwww在线观看| 国产精品久久视频播放| 久久中文字幕一级| 国产伦精品一区二区三区视频9 | 最新美女视频免费是黄的| 亚洲一区高清亚洲精品| 亚洲人成伊人成综合网2020| 黄片大片在线免费观看| 国产成人精品久久二区二区91| 哪里可以看免费的av片| 久久天堂一区二区三区四区| 亚洲精品一区av在线观看| 亚洲电影在线观看av| 欧美日韩综合久久久久久 | 亚洲男人的天堂狠狠| 午夜激情福利司机影院| 午夜免费激情av| a在线观看视频网站| 免费看美女性在线毛片视频| 狠狠狠狠99中文字幕| 国产一区二区三区在线臀色熟女| 一个人免费在线观看的高清视频| 国内久久婷婷六月综合欲色啪| 欧美黑人巨大hd| cao死你这个sao货| 国产一区二区激情短视频| АⅤ资源中文在线天堂| 成在线人永久免费视频| 精品久久久久久久毛片微露脸| 久久精品国产亚洲av香蕉五月| 亚洲av免费在线观看| 女人高潮潮喷娇喘18禁视频| 久久婷婷人人爽人人干人人爱| 天天添夜夜摸| 国产精品九九99| 国产三级黄色录像| 18禁国产床啪视频网站| 久久中文看片网| 免费电影在线观看免费观看| 久久久久久久久久黄片| 床上黄色一级片| 波多野结衣高清作品| 日韩欧美三级三区| 久久婷婷人人爽人人干人人爱| 美女高潮喷水抽搐中文字幕| 国产精品,欧美在线| 丰满的人妻完整版| 婷婷六月久久综合丁香| 啦啦啦免费观看视频1| 一区福利在线观看| 欧美在线黄色| 中文字幕人成人乱码亚洲影| 午夜影院日韩av| 国产精品一及| 亚洲国产欧洲综合997久久,| 岛国视频午夜一区免费看| 久久久精品欧美日韩精品| 女生性感内裤真人,穿戴方法视频| 美女黄网站色视频| 日韩欧美三级三区| 麻豆国产97在线/欧美| 哪里可以看免费的av片| 免费看美女性在线毛片视频| 18禁美女被吸乳视频| 亚洲精品久久国产高清桃花| 无人区码免费观看不卡| 国产麻豆成人av免费视频| 色av中文字幕| 成人18禁在线播放| 99久久综合精品五月天人人| 欧美极品一区二区三区四区| 婷婷丁香在线五月| 日本在线视频免费播放| 国产aⅴ精品一区二区三区波| 中文在线观看免费www的网站| 女人被狂操c到高潮| 蜜桃久久精品国产亚洲av| 免费在线观看视频国产中文字幕亚洲| 亚洲专区字幕在线| 久久国产精品影院| 日韩 欧美 亚洲 中文字幕| a级毛片在线看网站| 在线播放国产精品三级| 又黄又爽又免费观看的视频| 国产一区在线观看成人免费| 欧美性猛交╳xxx乱大交人| 亚洲国产精品成人综合色| 国产成人av教育| 黑人欧美特级aaaaaa片| 亚洲精品一卡2卡三卡4卡5卡| 少妇人妻一区二区三区视频| 最近视频中文字幕2019在线8| 一本一本综合久久| 久久午夜综合久久蜜桃| 黄色日韩在线| 99热精品在线国产| av福利片在线观看| av欧美777| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av熟女| xxxwww97欧美| 看片在线看免费视频| 不卡一级毛片| 免费高清视频大片| 亚洲五月婷婷丁香| 级片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 岛国在线免费视频观看| 女警被强在线播放| 亚洲精品乱码久久久v下载方式 | 国产在线精品亚洲第一网站| 亚洲av成人一区二区三| 一区二区三区高清视频在线| 午夜免费观看网址| 亚洲精品中文字幕一二三四区| 免费在线观看亚洲国产| 免费在线观看日本一区| 午夜两性在线视频| 精华霜和精华液先用哪个| 嫩草影院精品99| 夜夜爽天天搞| av片东京热男人的天堂| 国产又黄又爽又无遮挡在线| 成人永久免费在线观看视频| 国内精品美女久久久久久| 日本精品一区二区三区蜜桃| 日日干狠狠操夜夜爽| 亚洲成人久久爱视频| 日本一本二区三区精品| 99视频精品全部免费 在线 | 桃色一区二区三区在线观看| 可以在线观看毛片的网站| 亚洲熟女毛片儿| 国产免费av片在线观看野外av| 18美女黄网站色大片免费观看| 国产激情久久老熟女| 韩国av一区二区三区四区| 午夜视频精品福利| 亚洲熟妇中文字幕五十中出| 亚洲国产精品成人综合色| 国产v大片淫在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 身体一侧抽搐| 婷婷六月久久综合丁香| 99国产精品一区二区蜜桃av| 狂野欧美激情性xxxx| 黑人操中国人逼视频| 18禁观看日本| 亚洲精品美女久久久久99蜜臀| 18禁裸乳无遮挡免费网站照片| 精品99又大又爽又粗少妇毛片 | 亚洲专区国产一区二区| www日本黄色视频网| 久久精品综合一区二区三区| 国产美女午夜福利| 成人精品一区二区免费| 精品国产超薄肉色丝袜足j| 亚洲av成人av| 黄片大片在线免费观看| 久久热在线av| 不卡一级毛片| 天堂动漫精品| 18禁裸乳无遮挡免费网站照片| 91麻豆av在线| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| av在线天堂中文字幕| а√天堂www在线а√下载| 真人一进一出gif抽搐免费| 色综合欧美亚洲国产小说| 日日夜夜操网爽| 动漫黄色视频在线观看| 国产欧美日韩一区二区三| 成人午夜高清在线视频| 精品午夜福利视频在线观看一区| 91在线观看av| 波多野结衣高清无吗| 欧美性猛交黑人性爽| 一级作爱视频免费观看| 夜夜躁狠狠躁天天躁| 成年女人看的毛片在线观看| 亚洲av免费在线观看| a级毛片在线看网站| www.精华液| 99视频精品全部免费 在线 | 色综合欧美亚洲国产小说| 757午夜福利合集在线观看| 97碰自拍视频| 久久香蕉国产精品| 午夜激情欧美在线| 国产精品一区二区三区四区免费观看 | 综合色av麻豆| 真人一进一出gif抽搐免费| 99热这里只有是精品50| 国产极品精品免费视频能看的| 男女床上黄色一级片免费看| 男人的好看免费观看在线视频| 桃色一区二区三区在线观看| 亚洲中文日韩欧美视频| 99国产极品粉嫩在线观看| 中文字幕熟女人妻在线| 非洲黑人性xxxx精品又粗又长| 亚洲国产日韩欧美精品在线观看 | 成人精品一区二区免费| 97人妻精品一区二区三区麻豆| 啪啪无遮挡十八禁网站| 看黄色毛片网站| 一区福利在线观看| 久99久视频精品免费| 成人国产综合亚洲| 日本 av在线| www.www免费av| 中文在线观看免费www的网站| 日韩欧美国产一区二区入口| 嫩草影院入口| 怎么达到女性高潮| 亚洲精品在线观看二区| 九九热线精品视视频播放| 欧美丝袜亚洲另类 | 啦啦啦韩国在线观看视频| 伊人久久大香线蕉亚洲五| 国产精品av久久久久免费| 久久久久性生活片| 久久久久免费精品人妻一区二区| 精品久久久久久,| 观看美女的网站| 欧美中文日本在线观看视频| 在线观看舔阴道视频| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 成年女人毛片免费观看观看9| 欧美乱码精品一区二区三区| 免费看a级黄色片| 夜夜爽天天搞| 91老司机精品| 久久香蕉精品热| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 亚洲av中文字字幕乱码综合| 男人舔女人的私密视频| 亚洲性夜色夜夜综合| e午夜精品久久久久久久| 欧美绝顶高潮抽搐喷水| 97人妻精品一区二区三区麻豆| 51午夜福利影视在线观看| 欧美日韩亚洲国产一区二区在线观看| 小说图片视频综合网站| 国产黄a三级三级三级人| 非洲黑人性xxxx精品又粗又长| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 久久久久国内视频| 一区二区三区国产精品乱码| 中文字幕人成人乱码亚洲影| 亚洲欧美日韩卡通动漫| 国产人伦9x9x在线观看| 动漫黄色视频在线观看| 国产三级黄色录像| 精品99又大又爽又粗少妇毛片 | 国语自产精品视频在线第100页| 国产精品久久久人人做人人爽| 给我免费播放毛片高清在线观看| 成人一区二区视频在线观看| 人人妻,人人澡人人爽秒播| 国产伦人伦偷精品视频| 成人三级做爰电影| 夜夜看夜夜爽夜夜摸| 亚洲中文av在线| 久久草成人影院| 成年版毛片免费区| 久久中文看片网| 国产成人av教育| 天天一区二区日本电影三级| 精品国产美女av久久久久小说| 1024手机看黄色片| 亚洲专区字幕在线| 国产免费av片在线观看野外av| 美女 人体艺术 gogo| 日本熟妇午夜| 欧美黄色片欧美黄色片| 亚洲成人久久爱视频| 少妇熟女aⅴ在线视频| 宅男免费午夜| 两人在一起打扑克的视频| 白带黄色成豆腐渣| 中文字幕人妻丝袜一区二区| 91麻豆av在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲九九香蕉| svipshipincom国产片| 久久久久久九九精品二区国产| 日韩欧美免费精品| 一本精品99久久精品77| 久久精品人妻少妇| 两个人视频免费观看高清| 免费在线观看日本一区| 国产免费av片在线观看野外av| 99久久综合精品五月天人人| 夜夜看夜夜爽夜夜摸| 叶爱在线成人免费视频播放| 色尼玛亚洲综合影院| 国产爱豆传媒在线观看| 母亲3免费完整高清在线观看| 啦啦啦韩国在线观看视频| 在线观看一区二区三区| av片东京热男人的天堂| 三级国产精品欧美在线观看 | 一级毛片高清免费大全| 宅男免费午夜| 亚洲精品456在线播放app | 此物有八面人人有两片| 母亲3免费完整高清在线观看| 国产精品电影一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲五月天丁香| 亚洲av成人一区二区三| 99久久精品热视频| 黄色片一级片一级黄色片| av天堂中文字幕网| 一个人免费在线观看的高清视频| 首页视频小说图片口味搜索| 小说图片视频综合网站| 久久久久久久午夜电影| 国产高清videossex| 久久精品人妻少妇| 男插女下体视频免费在线播放| 香蕉久久夜色| 中文字幕人妻丝袜一区二区| 久久久久国产一级毛片高清牌| 精品国产亚洲在线| 免费观看人在逋| 久久久精品大字幕| 亚洲成人精品中文字幕电影| 成人无遮挡网站| 一区福利在线观看| 成年版毛片免费区| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 好男人电影高清在线观看| 91久久精品国产一区二区成人 | 给我免费播放毛片高清在线观看| 亚洲国产日韩欧美精品在线观看 | 中文字幕人妻丝袜一区二区| 首页视频小说图片口味搜索| 午夜视频精品福利| 一本久久中文字幕| 国产精品国产高清国产av| 少妇裸体淫交视频免费看高清| 日韩欧美在线乱码| 亚洲18禁久久av| 国产精品一区二区免费欧美| 草草在线视频免费看| 精品久久蜜臀av无| av中文乱码字幕在线| 国产成人av教育| 国产极品精品免费视频能看的| 欧美黑人巨大hd| 两性午夜刺激爽爽歪歪视频在线观看| 哪里可以看免费的av片| 日本三级黄在线观看| 国产三级中文精品| 久久久久国内视频| 国产午夜精品论理片| 哪里可以看免费的av片| 中文在线观看免费www的网站| 欧美在线黄色| 国产成人av教育| 国产一级毛片七仙女欲春2| cao死你这个sao货| 在线观看舔阴道视频| 国产高清视频在线观看网站| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 男人和女人高潮做爰伦理| 国内精品一区二区在线观看| 久久久久久久午夜电影| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 色视频www国产| 精品一区二区三区视频在线观看免费| 亚洲精品中文字幕一二三四区| 成人欧美大片| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 日本免费一区二区三区高清不卡| 俺也久久电影网| 亚洲va日本ⅴa欧美va伊人久久| АⅤ资源中文在线天堂| 高清在线国产一区| 国产精品久久久av美女十八| 在线观看免费午夜福利视频| 曰老女人黄片| 观看美女的网站| 一区二区三区国产精品乱码| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| ponron亚洲| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 黄频高清免费视频| 国产成人精品久久二区二区免费| 国产成年人精品一区二区| 99热这里只有精品一区 | 免费高清视频大片| 亚洲第一电影网av| 国产精品香港三级国产av潘金莲| 亚洲一区高清亚洲精品| 老熟妇仑乱视频hdxx| 18美女黄网站色大片免费观看| 国产免费男女视频| 久久久国产精品麻豆| 久久人妻av系列| 一区福利在线观看| 又爽又黄无遮挡网站| 亚洲国产精品999在线| 亚洲第一欧美日韩一区二区三区| 老熟妇仑乱视频hdxx| 国产成人影院久久av| 欧美日韩综合久久久久久 | 午夜a级毛片| 美女高潮的动态| 午夜福利在线在线| 欧美日韩一级在线毛片| 99视频精品全部免费 在线 | 国产毛片a区久久久久| or卡值多少钱| 亚洲av片天天在线观看| 99在线人妻在线中文字幕| 国产人伦9x9x在线观看| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 一本综合久久免费| 91麻豆精品激情在线观看国产| 国产视频内射| 日韩人妻高清精品专区| 亚洲色图 男人天堂 中文字幕| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区| 亚洲自拍偷在线| 免费看日本二区| 免费观看精品视频网站| 欧美又色又爽又黄视频| 后天国语完整版免费观看| 日韩中文字幕欧美一区二区| 两性夫妻黄色片| 久久精品91蜜桃| 丝袜人妻中文字幕| av福利片在线观看| 久久精品91无色码中文字幕| 久久久久九九精品影院| 精品国产亚洲在线| 久99久视频精品免费| av天堂在线播放| 99久久综合精品五月天人人| 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 黄色视频,在线免费观看| 亚洲av片天天在线观看| 综合色av麻豆| 日韩欧美在线乱码| 欧美性猛交黑人性爽| 9191精品国产免费久久| 国产精品自产拍在线观看55亚洲| 亚洲中文日韩欧美视频| 国产精品av久久久久免费| 搡老熟女国产l中国老女人| 99国产综合亚洲精品| 一个人看视频在线观看www免费 | 18禁黄网站禁片免费观看直播| 又爽又黄无遮挡网站| av欧美777| 色尼玛亚洲综合影院| 黄片小视频在线播放| 久久久久久久久免费视频了| 国产欧美日韩精品亚洲av| 丁香六月欧美| 久久久久久大精品| 熟妇人妻久久中文字幕3abv| 亚洲七黄色美女视频| 国内精品一区二区在线观看| 免费看十八禁软件| 国内精品美女久久久久久| 亚洲成人久久爱视频| 999久久久国产精品视频| 久99久视频精品免费| 午夜福利高清视频| 在线观看免费午夜福利视频| 久久久久久人人人人人| 在线免费观看不下载黄p国产 | 日韩成人在线观看一区二区三区| 三级毛片av免费| 美女大奶头视频| 亚洲国产精品久久男人天堂| 精品乱码久久久久久99久播| 国产在线精品亚洲第一网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲男人的天堂狠狠| 他把我摸到了高潮在线观看| 男女那种视频在线观看| 十八禁网站免费在线| 波多野结衣巨乳人妻| 嫩草影院精品99| 日本熟妇午夜| 国产精品 国内视频|