• <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| 极品少妇高潮喷水抽搐| 亚洲情色 制服丝袜| 老汉色av国产亚洲站长工具| 免费日韩欧美在线观看| 国产一区亚洲一区在线观看| 精品免费久久久久久久清纯 | tube8黄色片| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩一级在线毛片| 伦理电影大哥的女人| 亚洲专区中文字幕在线 | 精品一区二区三卡| 王馨瑶露胸无遮挡在线观看| 又大又黄又爽视频免费| 国产毛片在线视频| 女人被躁到高潮嗷嗷叫费观| 中文字幕最新亚洲高清| 超色免费av| 久久精品国产亚洲av涩爱| 久久久精品94久久精品| 你懂的网址亚洲精品在线观看| 9色porny在线观看| 国产97色在线日韩免费| 久久亚洲国产成人精品v| 1024香蕉在线观看| 最近手机中文字幕大全| 自拍欧美九色日韩亚洲蝌蚪91| 男女下面插进去视频免费观看| 18禁国产床啪视频网站| 水蜜桃什么品种好| 国产精品免费大片| 一区二区三区精品91| 亚洲中文av在线| 久久精品熟女亚洲av麻豆精品| 韩国av在线不卡| 免费观看性生交大片5| 十分钟在线观看高清视频www| 51午夜福利影视在线观看| 蜜桃在线观看..| 一区二区三区激情视频| 国精品久久久久久国模美| 桃花免费在线播放| 日韩大码丰满熟妇| 久久久国产精品麻豆| 成年av动漫网址| 国产精品99久久99久久久不卡 | 天天影视国产精品| 国产精品免费视频内射| 黄色怎么调成土黄色| 99热网站在线观看| 亚洲精品视频女| 国产高清不卡午夜福利| av在线app专区| 9热在线视频观看99| 亚洲精品av麻豆狂野| 成人国语在线视频| 日韩一本色道免费dvd| 一本一本久久a久久精品综合妖精| 久久久国产欧美日韩av| 久久久久国产精品人妻一区二区| 国产日韩一区二区三区精品不卡| 国产国语露脸激情在线看| 中文字幕人妻丝袜制服| 欧美人与善性xxx| 久久久国产欧美日韩av| 777久久人妻少妇嫩草av网站| xxxhd国产人妻xxx| 成人国语在线视频| 制服诱惑二区| 男女免费视频国产| 久久女婷五月综合色啪小说| 久久精品aⅴ一区二区三区四区| 老司机靠b影院| www日本在线高清视频| 天美传媒精品一区二区| 国产成人欧美在线观看 | 亚洲精品国产一区二区精华液| 母亲3免费完整高清在线观看| 19禁男女啪啪无遮挡网站| 亚洲国产精品国产精品| 久热爱精品视频在线9| 精品少妇内射三级| 多毛熟女@视频| 国产麻豆69| 1024视频免费在线观看| 波野结衣二区三区在线| 午夜福利一区二区在线看| 啦啦啦啦在线视频资源| 日韩免费高清中文字幕av| 色视频在线一区二区三区| 国产精品人妻久久久影院| 亚洲av国产av综合av卡| 天堂中文最新版在线下载| 久久天堂一区二区三区四区| 国产一区有黄有色的免费视频| 黄网站色视频无遮挡免费观看| 亚洲欧美一区二区三区久久| www.熟女人妻精品国产| 一区二区三区精品91| 亚洲精品aⅴ在线观看| 亚洲国产毛片av蜜桃av| 亚洲精品在线美女| av国产久精品久网站免费入址| 一级毛片电影观看| 精品一区二区三区四区五区乱码 | 秋霞在线观看毛片| 我要看黄色一级片免费的| 久久久久久久久免费视频了| 18禁裸乳无遮挡动漫免费视频| 国产黄频视频在线观看| a级毛片在线看网站| 成年美女黄网站色视频大全免费| 免费高清在线观看日韩| 人妻一区二区av| 精品视频人人做人人爽| 九九爱精品视频在线观看| 性色av一级| 久热爱精品视频在线9| 热99国产精品久久久久久7| 99久国产av精品国产电影| 国语对白做爰xxxⅹ性视频网站| 午夜免费鲁丝| 9色porny在线观看| 精品一区二区三区av网在线观看 | 免费观看av网站的网址| 欧美日本中文国产一区发布| 国产毛片在线视频| 亚洲精品一二三| 国产成人午夜福利电影在线观看| h视频一区二区三区| 午夜免费观看性视频| 国产一区有黄有色的免费视频| 亚洲av在线观看美女高潮| 男人舔女人的私密视频| 女的被弄到高潮叫床怎么办| 999精品在线视频| 亚洲精品国产区一区二| 韩国高清视频一区二区三区| 免费黄色在线免费观看| 久久久久精品性色| 曰老女人黄片| 高清欧美精品videossex| 中文字幕精品免费在线观看视频| 久久人人爽av亚洲精品天堂| 成人漫画全彩无遮挡| 国产男女内射视频| 成人国产av品久久久| 女性被躁到高潮视频| 亚洲婷婷狠狠爱综合网| 国精品久久久久久国模美| 久久久久人妻精品一区果冻| 免费久久久久久久精品成人欧美视频| 日韩一卡2卡3卡4卡2021年| 亚洲熟女精品中文字幕| 十八禁人妻一区二区| 国产精品 国内视频| 亚洲国产精品一区二区三区在线| 人体艺术视频欧美日本| 国产在线免费精品| 国产精品二区激情视频| 国产午夜精品一二区理论片| 国产精品麻豆人妻色哟哟久久| 国产亚洲精品第一综合不卡| 国产成人啪精品午夜网站| 中文字幕色久视频| 日韩一卡2卡3卡4卡2021年| 国产成人91sexporn| 菩萨蛮人人尽说江南好唐韦庄| 69精品国产乱码久久久| 国产精品一国产av| 国产精品成人在线| 如日韩欧美国产精品一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲美女黄色视频免费看| 亚洲美女黄色视频免费看| 波多野结衣一区麻豆| 一区在线观看完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 国产极品粉嫩免费观看在线| 久久这里只有精品19| 久久这里只有精品19| 日日撸夜夜添| 一本色道久久久久久精品综合| 色婷婷av一区二区三区视频| 亚洲免费av在线视频| e午夜精品久久久久久久| 色婷婷av一区二区三区视频| 可以免费在线观看a视频的电影网站 | 国产欧美亚洲国产| 亚洲av在线观看美女高潮| av国产久精品久网站免费入址| 久久精品国产综合久久久| 老司机深夜福利视频在线观看 | 亚洲欧美成人精品一区二区| a级片在线免费高清观看视频| 一本一本久久a久久精品综合妖精| 蜜桃在线观看..| 亚洲国产日韩一区二区| 国产欧美亚洲国产| 国产一区二区三区综合在线观看| 国产毛片在线视频| 深夜精品福利| 十八禁网站网址无遮挡| 亚洲中文av在线| 成人亚洲精品一区在线观看| 欧美日韩综合久久久久久| 国产有黄有色有爽视频| 色94色欧美一区二区| 国产成人精品在线电影| 99国产综合亚洲精品| 国产有黄有色有爽视频| 成人毛片60女人毛片免费| 免费高清在线观看视频在线观看| 亚洲天堂av无毛| 国产一区二区在线观看av| 女人精品久久久久毛片| 精品午夜福利在线看| 多毛熟女@视频| 国产一区二区 视频在线| 日日摸夜夜添夜夜爱| av国产久精品久网站免费入址| 国产精品久久久久成人av| 亚洲精品日本国产第一区| 久久久久精品性色| 欧美日韩亚洲国产一区二区在线观看 | 国产精品女同一区二区软件| 免费观看a级毛片全部| 日本av免费视频播放| 99热国产这里只有精品6| 成人国语在线视频| 亚洲精品国产一区二区精华液| 久久午夜综合久久蜜桃| 国产免费视频播放在线视频| 国产精品熟女久久久久浪| 亚洲欧美日韩另类电影网站| 亚洲精品成人av观看孕妇| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠躁躁| svipshipincom国产片| 七月丁香在线播放| 777米奇影视久久| 大片电影免费在线观看免费| 搡老乐熟女国产| 久久韩国三级中文字幕| 老司机深夜福利视频在线观看 | 宅男免费午夜| 亚洲av日韩精品久久久久久密 | 国产亚洲一区二区精品| 免费人妻精品一区二区三区视频| 制服诱惑二区| 色综合欧美亚洲国产小说| av天堂久久9| 丝袜美足系列| av女优亚洲男人天堂| 一级黄片播放器| 黑人巨大精品欧美一区二区蜜桃| 日本91视频免费播放| 国产日韩欧美亚洲二区| 另类亚洲欧美激情| 国产欧美亚洲国产| 国产成人精品久久久久久| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 一边亲一边摸免费视频| 久久久久人妻精品一区果冻| 亚洲精品av麻豆狂野| 亚洲自偷自拍图片 自拍| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久精品古装| 高清在线视频一区二区三区| 国产精品香港三级国产av潘金莲 | 少妇精品久久久久久久| 久久亚洲国产成人精品v| 午夜福利免费观看在线| 国产97色在线日韩免费| 91精品伊人久久大香线蕉| 又黄又粗又硬又大视频| 亚洲精品在线美女| 国产乱人偷精品视频| 国产成人系列免费观看| 亚洲视频免费观看视频| 国产探花极品一区二区| 亚洲美女黄色视频免费看| 中文精品一卡2卡3卡4更新| 色视频在线一区二区三区| 天天躁日日躁夜夜躁夜夜| 亚洲国产精品成人久久小说| 悠悠久久av| 精品亚洲成a人片在线观看| 免费看不卡的av| 女性生殖器流出的白浆| 午夜福利网站1000一区二区三区| 国产精品无大码| 精品午夜福利在线看| 午夜免费鲁丝| 久久久久网色| 亚洲国产成人一精品久久久| 精品亚洲成国产av| 熟女av电影| 午夜福利在线免费观看网站| 超色免费av| 一本—道久久a久久精品蜜桃钙片| svipshipincom国产片| 中文字幕高清在线视频| 少妇人妻精品综合一区二区| 丝瓜视频免费看黄片| 亚洲精华国产精华液的使用体验| 亚洲精品美女久久av网站| 中文字幕人妻丝袜制服| 别揉我奶头~嗯~啊~动态视频 | 热re99久久精品国产66热6| 亚洲少妇的诱惑av| 久久av网站| 久久女婷五月综合色啪小说| 亚洲美女黄色视频免费看| 美女扒开内裤让男人捅视频| 天天躁日日躁夜夜躁夜夜| 亚洲精品中文字幕在线视频| 久久精品aⅴ一区二区三区四区| 波野结衣二区三区在线| 亚洲精品,欧美精品| 精品午夜福利在线看| 少妇猛男粗大的猛烈进出视频| 国产精品一二三区在线看| 搡老岳熟女国产| 人人妻人人澡人人爽人人夜夜| 日韩av不卡免费在线播放| 女人精品久久久久毛片| 国精品久久久久久国模美| 韩国av在线不卡| 在线观看一区二区三区激情| 国产欧美亚洲国产| 亚洲精品视频女| 中文精品一卡2卡3卡4更新| 王馨瑶露胸无遮挡在线观看| 国产精品99久久99久久久不卡 | 91国产中文字幕| 欧美国产精品一级二级三级| 成年人午夜在线观看视频| 乱人伦中国视频| 黄色视频不卡| 国产黄频视频在线观看| 爱豆传媒免费全集在线观看| 黄色毛片三级朝国网站| 韩国高清视频一区二区三区| 午夜免费观看性视频| 人人妻人人澡人人爽人人夜夜| 日本av免费视频播放| 久久午夜综合久久蜜桃| 下体分泌物呈黄色| av线在线观看网站| 十八禁网站网址无遮挡| 婷婷色麻豆天堂久久| 18禁裸乳无遮挡动漫免费视频| 日韩成人av中文字幕在线观看| 亚洲综合精品二区| 女人高潮潮喷娇喘18禁视频| 91精品三级在线观看| 国产又色又爽无遮挡免| 午夜免费男女啪啪视频观看| 99国产综合亚洲精品| 亚洲美女黄色视频免费看| av网站在线播放免费| 亚洲精品国产av成人精品| 欧美成人精品欧美一级黄| 午夜免费观看性视频| 精品国产一区二区三区久久久樱花| 哪个播放器可以免费观看大片| 中文字幕色久视频| 亚洲婷婷狠狠爱综合网| 精品国产超薄肉色丝袜足j| 亚洲视频免费观看视频| 精品福利永久在线观看| 十分钟在线观看高清视频www| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 欧美av亚洲av综合av国产av | 国产成人欧美在线观看 | 久久精品熟女亚洲av麻豆精品| 亚洲中文av在线| 日本色播在线视频| 国产成人av激情在线播放| 午夜激情av网站| 久久国产精品男人的天堂亚洲| 免费观看人在逋| 日韩av免费高清视频| 日韩 亚洲 欧美在线| av视频免费观看在线观看| 51午夜福利影视在线观看| 国产精品亚洲av一区麻豆 | 久久国产精品男人的天堂亚洲| 亚洲成人av在线免费| 欧美日韩成人在线一区二区| 国产精品偷伦视频观看了| 国产精品秋霞免费鲁丝片| 777米奇影视久久| 日韩成人av中文字幕在线观看| 好男人视频免费观看在线| 国产精品久久久久成人av| 老司机在亚洲福利影院| 亚洲熟女毛片儿| 久热这里只有精品99| 久久久国产精品麻豆| 男女下面插进去视频免费观看| 天天躁夜夜躁狠狠久久av| 午夜av观看不卡| 操美女的视频在线观看| 免费av中文字幕在线| 我的亚洲天堂| 亚洲 欧美一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区精品视频观看| 精品国产超薄肉色丝袜足j| 欧美97在线视频| 十八禁网站网址无遮挡| 欧美激情高清一区二区三区 | 免费久久久久久久精品成人欧美视频| 校园人妻丝袜中文字幕| 欧美日韩亚洲高清精品| 精品人妻熟女毛片av久久网站| 亚洲一卡2卡3卡4卡5卡精品中文| av在线播放精品| 丝瓜视频免费看黄片| 丝袜美腿诱惑在线| 国产精品免费大片| 亚洲熟女毛片儿| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 国产免费视频播放在线视频| 久久久久久人妻| 在线看a的网站| 亚洲av福利一区| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 成人亚洲欧美一区二区av| 一级片'在线观看视频| 丰满饥渴人妻一区二区三| 伦理电影大哥的女人| 亚洲成国产人片在线观看| 最近的中文字幕免费完整| 99re6热这里在线精品视频| 少妇人妻 视频| 成年人午夜在线观看视频| 在线精品无人区一区二区三| 日韩不卡一区二区三区视频在线| 一本大道久久a久久精品| 亚洲精品一二三| 日韩制服骚丝袜av| 97人妻天天添夜夜摸| 国产欧美亚洲国产| 七月丁香在线播放| 亚洲欧美精品综合一区二区三区| 国产一区二区三区综合在线观看| 狂野欧美激情性xxxx| 久久久精品区二区三区| 女人久久www免费人成看片| 欧美日韩亚洲国产一区二区在线观看 | 少妇人妻久久综合中文| 久久久久网色| 成年人免费黄色播放视频| 各种免费的搞黄视频| 久久久久久免费高清国产稀缺| 欧美最新免费一区二区三区| 一本—道久久a久久精品蜜桃钙片| 女性被躁到高潮视频| 1024香蕉在线观看| 极品人妻少妇av视频| 午夜福利,免费看| 国产黄色免费在线视频| 久久精品久久精品一区二区三区| 亚洲美女搞黄在线观看| 国产女主播在线喷水免费视频网站| av国产精品久久久久影院| 韩国精品一区二区三区| 亚洲av日韩在线播放| 一边亲一边摸免费视频| 亚洲成人一二三区av| 色网站视频免费| 久久影院123| 久久人人97超碰香蕉20202| 国产av码专区亚洲av| 国产成人精品在线电影| 久久av网站| 蜜桃国产av成人99| 嫩草影院入口| 日本午夜av视频| 欧美激情高清一区二区三区 | av电影中文网址| 在线精品无人区一区二区三| 国产成人免费观看mmmm| 黄网站色视频无遮挡免费观看| 亚洲成色77777| 日韩中文字幕视频在线看片| 婷婷成人精品国产| 免费女性裸体啪啪无遮挡网站| 观看美女的网站| 国产又爽黄色视频| 制服诱惑二区| 国产老妇伦熟女老妇高清| 19禁男女啪啪无遮挡网站| 亚洲国产中文字幕在线视频| av线在线观看网站| 日韩熟女老妇一区二区性免费视频| 一级毛片 在线播放| 毛片一级片免费看久久久久| 999久久久国产精品视频| 人成视频在线观看免费观看| 亚洲精品久久午夜乱码| 国产日韩欧美在线精品| 欧美精品一区二区免费开放| 人妻人人澡人人爽人人| 不卡视频在线观看欧美| 国产精品免费大片| av线在线观看网站| 日韩大码丰满熟妇| 午夜日韩欧美国产| 一本一本久久a久久精品综合妖精| 日韩中文字幕视频在线看片| 韩国av在线不卡| 亚洲av电影在线进入| 久久综合国产亚洲精品| 欧美精品高潮呻吟av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 人妻人人澡人人爽人人| 亚洲男人天堂网一区| 久久久久国产一级毛片高清牌| 免费看av在线观看网站| 狠狠婷婷综合久久久久久88av| 国产99久久九九免费精品| 国产淫语在线视频| 精品国产一区二区久久| 九色亚洲精品在线播放| 亚洲中文av在线| 亚洲av电影在线观看一区二区三区| 亚洲第一av免费看| 国产麻豆69| 久久久久久久国产电影| 国产精品偷伦视频观看了| 美国免费a级毛片| 国产在线视频一区二区| 亚洲国产精品一区三区| 久久久精品免费免费高清| 一本—道久久a久久精品蜜桃钙片| 午夜精品国产一区二区电影| 亚洲精品视频女| 日韩欧美精品免费久久| 美女高潮到喷水免费观看| 亚洲欧洲精品一区二区精品久久久 | 美女福利国产在线| 亚洲,欧美,日韩| 欧美国产精品一级二级三级| 巨乳人妻的诱惑在线观看| 99久久精品国产亚洲精品| 天天躁夜夜躁狠狠躁躁| 麻豆av在线久日| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频 | 人人妻人人添人人爽欧美一区卜| 日本91视频免费播放| 精品一区二区三区四区五区乱码 | 国产亚洲最大av| 精品视频人人做人人爽| 国产成人免费观看mmmm| 伊人久久国产一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲三区欧美一区| 18在线观看网站| 不卡av一区二区三区| 中文字幕人妻丝袜制服| 精品久久久精品久久久| 成人亚洲精品一区在线观看| 亚洲精品av麻豆狂野| 黄片小视频在线播放| 国产人伦9x9x在线观看| 另类亚洲欧美激情| 男女高潮啪啪啪动态图| 国产精品久久久人人做人人爽| 啦啦啦在线观看免费高清www| 国产精品麻豆人妻色哟哟久久| 成人国产av品久久久| 美女国产高潮福利片在线看| 男女之事视频高清在线观看 | 中文欧美无线码| 亚洲欧美激情在线| 亚洲一区二区三区欧美精品| 亚洲欧美中文字幕日韩二区| 人成视频在线观看免费观看| 男人爽女人下面视频在线观看| 狠狠婷婷综合久久久久久88av| 一级毛片我不卡| 久久国产精品大桥未久av| 丝袜美足系列| 久久女婷五月综合色啪小说| 卡戴珊不雅视频在线播放| 伦理电影大哥的女人| 91精品伊人久久大香线蕉| videosex国产| 精品国产国语对白av| 最近中文字幕2019免费版| 亚洲精品国产av成人精品| 在线天堂最新版资源| 亚洲欧美清纯卡通| 免费女性裸体啪啪无遮挡网站| 老司机深夜福利视频在线观看 | 观看美女的网站| 欧美人与善性xxx| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| 国产成人免费无遮挡视频| 日韩熟女老妇一区二区性免费视频| 黄色毛片三级朝国网站| 人人妻人人澡人人爽人人夜夜| 99久久99久久久精品蜜桃| 99精国产麻豆久久婷婷|