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

    QR分解與GLS改進算法在壓力傳感器動態(tài)性能改善中的應用

    2016-12-15 12:31:57楊文杰張志杰李巖峰穆欣榮
    傳感技術學報 2016年11期
    關鍵詞:階次階躍補償

    楊文杰,張志杰,李巖峰,穆欣榮

    (中北大學儀器科學與動態(tài)測試教育部重點試驗室,太原030051)

    QR分解與GLS改進算法在壓力傳感器動態(tài)性能改善中的應用

    楊文杰,張志杰*,李巖峰,穆欣榮

    (中北大學儀器科學與動態(tài)測試教育部重點試驗室,太原030051)

    將矩形矩陣的正交(QR)分解和改進的最小二乘(GLS)算法相結合,實現(xiàn)傳感器的動態(tài)建模,運用零極點配置理論改善傳感器動態(tài)性能。利用傳感器動態(tài)校準實驗數(shù)據(jù),首先構造矩形數(shù)據(jù)矩陣并對矩陣進行正交分解確定傳感器模型階次,繼而利用改進的最小二乘算法建立傳感器的數(shù)學模型,而后運用零極點理論對模型零極點進行重新配置改善傳感器的動態(tài)性能。實驗結果表明:該方法簡潔,建模精度高,能夠較準確得到反映傳感器動態(tài)特性的數(shù)學模型,且基于該方法的傳感器動態(tài)性能改善效果明顯。

    壓力傳感器動態(tài)建模正交分解最小二乘零極點

    隨著科學技術的發(fā)展,在科學研究、技術開發(fā)、工程實踐等領域,越來越多的要求進行動態(tài)測試,要求定量、深入地了解瞬態(tài)過程參數(shù)的變化規(guī)律[1-2]。動態(tài)測試產(chǎn)生的誤差主要是由于傳感器動態(tài)性能不佳引起的,因此改善傳感器的動態(tài)性能是減小測試誤差的一條有效途徑。

    改善傳感器動態(tài)性能的方法根據(jù)是否依賴傳感器動態(tài)模型可分兩類:依賴傳感器動態(tài)模型類型和不依賴傳感器動態(tài)模型類型。零極點配置法[3],反濾波法[4]等屬于典型的依賴傳感器動態(tài)模型類型的改善方法。利用PSO算法直接對硅加速度傳感器進行動態(tài)補償?shù)姆椒ǎ?],基于神經(jīng)網(wǎng)絡逆系統(tǒng)的傳感器動態(tài)補償方法[6-7]等都屬于不依賴傳感器動態(tài)模型方法。不依賴傳感器動態(tài)模型的改善方法較依賴傳感器動態(tài)模型的方法少了一個求解傳感器系統(tǒng)特性的環(huán)節(jié),顯得簡便,但是卻得不到傳感器系統(tǒng)特性,并且傳感器動態(tài)性能改善程度也很難定量表示。兩種類型的改善方法共同面臨著階次確定的問題,即依賴傳感器動態(tài)模型類型面臨著模型階次判定問題,不依賴傳感器動態(tài)模型類型面臨著補償環(huán)節(jié)階次確定問題。本文用矩陣正交分解來判定傳感器動態(tài)模型階次,用改進的最小二乘算法建立傳感器動態(tài)模型,用零極點配置理論對傳感器動態(tài)模型的零極點進行重新配置以改善傳感器的動態(tài)性能。

    1 壓力傳感器動態(tài)校準

    動態(tài)校準的首要問題是要有頻帶能夠充分覆蓋被校傳感器模態(tài)的動態(tài)激勵信號發(fā)生器,這樣才能將被校準傳感器的主要模態(tài)激發(fā)出來。本文選用激波管作為壓力傳感器動態(tài)校準實驗所施加激勵信號發(fā)生器。選用激波管作為激勵源是因為激波管可以產(chǎn)生上升沿很短(可以為納秒級),持續(xù)時間為數(shù)個毫秒的理想階躍信號[8]。校準時,將壓力傳感器齊平安裝在激波管末端端面上,如圖1所示。

    整個校準系統(tǒng)由高壓氣源、激波管、測速傳感器、信號調(diào)理電路、多通道數(shù)據(jù)采集卡等組成。標定激波管的主體為高壓區(qū)和低壓區(qū)兩部分,中間使用膜片隔開,校準時打開高壓氣瓶,當兩室壓差達到破膜壓力時膜片破裂,高壓室氣體沖向低壓室形成入射激波[9]。

    圖1 激波管動態(tài)校準系統(tǒng)原理圖

    2 QR分解判定動態(tài)模型階次

    矩形矩陣的正交分解又稱QR分解,是把一個m×n的矩陣A分解為一個正交矩陣Q和一個上三角矩陣R的乘積,即A=QR。

    利用實驗數(shù)據(jù){(u(k),y(k))|k=1,2,…,N}構造如下矩陣:

    矩陣D中的v為待考察的階次,N為數(shù)據(jù)長度,對D陣進行矩陣QR分解,得到如下上三角陣。

    式中,M=2v+2,×為任意元素。R′是一個M×M維的上三角方陣。則R′陣對角線偶元素平方和是各階差分方程模型對應于最小二乘估計的殘差平方和[10],因此可以通過對比各階模型對應的殘差平方和的大小來判斷模型的階次。

    3 改進GLS算法進行線性動態(tài)建模的原理與實現(xiàn)

    單輸入、單輸出線性定常系統(tǒng)可以用下面差分方程描述

    如果考慮噪聲,則式(1)可寫為

    式(2)簡寫為

    式中:u(k)為系統(tǒng)輸入數(shù)據(jù),y(k)為輸出數(shù)據(jù),d-1為延遲算子,ε(k)為擬合殘差,有

    建模就是根據(jù)標定實驗得到的輸入、輸出序列{(u(k),y(k))|k=1,2,…,N}確定模型階次n,并對模型參數(shù)ai,bi,i=0,1,…,n(a0=1)做出估計。模型的最小二乘格式

    式中,數(shù)據(jù)向量和參數(shù)向量定義為

    定義最小二乘優(yōu)化算法的指標函數(shù)

    由于建模時輸入是假定的理想階躍信號序列,此時輸入端噪聲可不予以考慮。如果假定輸出端是零均值的白噪聲,式(6)可變?yōu)?/p>

    則式(3)可寫為

    式(15)可以用最小二乘求得參數(shù) ai,bi,i=0,1,…,n(a0=1)的無偏估計。但由于在ai,bi尚未估計出之前,是未知的,因而需利用迭代法加以估計,其步驟如下:

    Step 1:由輸入、輸出數(shù)據(jù)(u(k),y(k)),用最小二乘法求出參數(shù)的最小二乘估計值,作為迭代運算的初值。

    Step 2:設迭代第l次時,求得

    計算

    Step 4:令l=l+1,回到第二步,直到迭代收斂或者達到迭代次數(shù)為止。當?shù)諗繒r,可以得到下式:

    4 傳感器動態(tài)補償方法

    傳感器的動態(tài)特性與其傳遞函數(shù)的零極點尤其是極點位置密切相關[11],對傳遞函數(shù)的零極點進行分析,觀察它們對傳感器動態(tài)性能的影響,繼而采用零極點相消補償原理,即串接一個補償環(huán)節(jié),而將原來不符合要求的零極點消去,配以新的零極點,使系統(tǒng)動態(tài)特性得以改善,這便是零極點相消補償法。本文應用參考文獻[12]補償傳感器動態(tài)特性的基本思路對傳感器動態(tài)特性進行改善。

    5 仿真實驗

    5.1 模型辨識方法驗證

    設壓力傳感器為二階系統(tǒng),模型為

    式(18)中假定,k=1,ξ=0.01,ω0=3.6×105π。傳感器的動態(tài)特性一般用傳遞函數(shù)描述,為了便于數(shù)據(jù)采集、處理和存儲,通常用離散傳遞函數(shù)描述。離散傳遞函數(shù)的一般表示形式為

    采用雙線性變換法,可以將連續(xù)傳遞函數(shù)離散化,變換公式為

    式中:T為采樣時間,本文假定為0.5 μs。則式(18)在Z域的表達式為

    圖2 壓力傳感器模型階躍響應曲線

    從圖2可以比較直觀看出,該傳感器模型對輸入階躍信號的響應存在超調(diào)量大,響應時間長等不足,在實際應用中會帶來較大的動態(tài)誤差,必須改善傳感器的動態(tài)性能。

    改善傳感器的動態(tài)性能之前往往需要獲取反映傳感器動態(tài)特性的動態(tài)指標,欲獲取動態(tài)指標需要建立全面描述傳感器動態(tài)特性的動態(tài)數(shù)學模型,而建立模型之前需要確定傳感器的階次。運用矩陣QR分解求得模型階次與殘差平方和的關系如圖3所示。

    圖3 殘差平方和隨模型階次變化曲線

    分析圖3數(shù)據(jù)知,前1、2階模型對應的殘差平方和分別為19.042 0、0.193 7,之后的變化都比較緩慢,因此可以認定該壓力傳感器模型為2階。

    當模型階次確定后,利用改進的最小二乘求得系統(tǒng)離散傳遞函數(shù)為

    圖4反映的是辨識系統(tǒng)與原始系統(tǒng)對階躍信號的響應,從圖4可以看出,兩者吻合較好,幾乎重疊在一起。

    圖4 原始系統(tǒng)和辨識系統(tǒng)階躍響應曲線

    5.2 改善模型特性方法驗證

    為了使系統(tǒng)滿足測試要求,必須要改善其動態(tài)特性,具體是拓寬傳感器的工作頻帶,因為對于傳感器而言,工作頻帶是其最重要的動態(tài)指標。參照零極點相消補償方法設計步驟[3,12],設計補償濾波器

    該濾波器串聯(lián)在原系統(tǒng)后面便能夠拓寬工作頻帶,圖5是傳感器系統(tǒng)頻帶拓寬效果。

    圖5 補償前后傳感器系統(tǒng)幅頻特性曲線

    圖5從頻率域表明,補償后傳感器的工作頻帶被大大拓寬,且補償后系統(tǒng)抗干擾性能提高。進一步分析圖5中的數(shù)據(jù)知,辨識系統(tǒng)在176 kHz附近發(fā)生諧振,接近之前設定的系統(tǒng)諧振頻率180 kHz,說明辨識系統(tǒng)與原始系統(tǒng)達到了很好的逼近效果。為了說明補償后系統(tǒng)時間域改善效果,繪制圖6。

    圖6 補償前后傳感器階躍響應

    圖6從時間域表明,補償后傳感器的超調(diào)量和響應時間都得到很大程度上的改善。為了更加直觀的得到補償前后傳感器各動態(tài)參數(shù)的改善程度,繪制補償前后系統(tǒng)各項動態(tài)參數(shù)于表1。

    表1 傳感器補償前后各項動態(tài)參數(shù)

    表1中,ts為響應時間,表征的是傳感器響應速度性能;σp為超調(diào)量,表征的是傳感器穩(wěn)定性能;wq為工作頻帶,這里取幅值誤差在±5%范圍內(nèi)的頻帶。

    5.3 補償模型適應性驗證

    補償模型設計后,需要判斷該補償模型的適應性及抗干擾性。對于模擬的理想系統(tǒng),該環(huán)節(jié)是必不可少的,因為該環(huán)節(jié)可以驗證補償及建模方法。具體方法是:在理想系統(tǒng)的階躍響應輸出數(shù)據(jù)中添加不同類型,不同信噪比的噪聲信號構成混合信號,將這種混合信號通過補償系統(tǒng),以此檢驗補償模型的一般適應性,即補償模型對外界擾動保持自身穩(wěn)定的性能。

    為了驗證補償模型適應性,選取了信噪比分別為26 dB,30 dB和40 dB的均勻噪聲和隨機白噪聲構成混合信號,圖7是補償模型對不同混合信號的處理效果。

    圖7 補償模型對不同類型和信噪比噪聲響應曲線

    通過圖7可知,補償模型適應性能優(yōu)異,說明補償模型選取的適合,也進一步說明傳感器數(shù)據(jù)建模方法適用且建模精確。

    6 實際應用

    圖8是應用于爆炸沖擊波測試領域某型號壓阻傳感器階躍響應曲線(激波管標定實驗獲得)。運用矩陣QR分解方法判定該壓阻傳感器模型階次為2。

    圖9是辨識系統(tǒng)的階躍響應和激波標定響應,從圖上看出辨識系統(tǒng)模型階躍響應曲線與標定實驗數(shù)據(jù)曲線比較一致。

    圖8 某壓阻傳感器階躍響應曲線

    圖9 原始系統(tǒng)和辨識系統(tǒng)階躍響應曲線

    圖10為補償前后傳感器幅頻特性曲線,分析補償前系統(tǒng)幅頻特性曲線可知,該系統(tǒng)幅值誤差在±2%范圍內(nèi)的工作頻帶約為50 kHz,不能有效覆蓋爆炸沖擊波信號帶寬(約為70 kHz),因此需要進行補償,即拓寬工作頻帶。分析補償后系統(tǒng)幅頻特性曲線可知,系統(tǒng)幅值誤差在±2%范圍內(nèi)的工作頻帶約為180 kHz,完全覆蓋了沖擊波信號帶寬。

    圖11為補償前后傳感器對階躍信號響應曲線,從圖中可以看出補償后系統(tǒng)階躍響應的超調(diào)量降幅明顯,響應時間也明顯縮短。

    圖10 補償前后傳感器系統(tǒng)幅頻特性

    圖11 補償前后傳感器系統(tǒng)階躍響應

    7 結論

    矩陣QR分解方法用以判定傳感器模型階次有著簡單、準確和實用等優(yōu)點。傳感器模型結構確定之后,用改進的最小二乘算法優(yōu)化傳感器模型參數(shù),能夠尋求一組實用參數(shù),且該算法易收斂,建模精度高。基于零極點理論的傳感器動態(tài)補償方法對傳感器動態(tài)性能改善顯著,且方法簡便。

    [1]吳建.傳感器動態(tài)補償濾波器及其硬件實現(xiàn)方式研究[D].太原:中北大學,2012.

    [2]黃曉敏.傳感器動態(tài)特性補償方法研究及其FPGA實現(xiàn)[D].太原:中北大學,2015.

    [3]黃俊欽,季平.動態(tài)補償數(shù)字濾波器的設計方法[J].儀器儀表學報,1987,8(3):45-46.

    [4]張志杰.動態(tài)測試中單次性過程測試數(shù)據(jù)的置信度—環(huán)境因子和反濾波過程[D].北京:北京理工大學,1997.

    [5]劉清,曹國華.基于微粒子群算法優(yōu)化的微硅加速度傳感器動態(tài)補償研究[J].儀器儀表學報,2006,27(12):1707-1710.

    [6]戴先中,殷銘,王勤.傳感器動態(tài)補償?shù)纳窠?jīng)網(wǎng)絡逆系統(tǒng)方法[J].儀器儀表學報,2004,25(5):593-596.

    [7]俞阿龍.基于遺傳神經(jīng)網(wǎng)絡的加速度傳感器動態(tài)建模方法[J].儀器儀表學報,2006,27(3):315-318.

    [8]軒青春,軒志偉,陳保立.基于最小二乘與粒子群算法的壓力傳感器動態(tài)補償方法[J].傳感技術學報,2014,27(10):1363-1366.

    [9]軒青春,軒志偉,賴富文.壓力傳感器測試系統(tǒng)的動態(tài)校準及特性分析[J].傳感技術學報,2015,28(7):983-84.

    [10]黃俊欽.測試系統(tǒng)動力學[M].北京:國防工業(yè)出版社,2013.

    [11]賴富文,張志杰,張建宇,等.基于動態(tài)特性補償?shù)臎_擊波測試數(shù)據(jù)處理方法[J].爆炸與沖擊,2015,35(6):872-873.

    [12]楊文杰,張志杰,趙晨陽,等.基于零極點配置理論的壓力傳感器動態(tài)特性補償[J].科學技術與工程,2016,16(2),78-82.

    楊文杰(1990-),男,河南駐馬店人,碩士研究生,主要從事動態(tài)測試系統(tǒng)搭建,信號獲取,存儲與處理以及智能儀器開發(fā),1210921680@qq.com;

    張志杰(1965-),男,山西五臺人。現(xiàn)任中北大學教授、博導,儀器科學與技術學科帶頭人之一,主要從事動態(tài)測試理論、技術與應用,信號處理理論與技術的研究。教育部儀器科學教學指導委員會協(xié)作委員,zhangzhijie@nuc.edu.cn。

    The Application of QR Decomposition and GLS Improved Algorithm in Improving the Dynamic Performance of Pressure Sensor

    YANG Wenjie,ZHANG Zhijie*,LI Yanfeng,MU Xinrong
    (Key Laboratory for Instrumentation Science&Dynamic Measurement,Ministry of Education,North University of China,Taiyuan 030051,China)

    Combining rectangular matrix orthogonal(QR)decomposition with improved least squares(GLS)algo?rithm,the dynamic mathematical model of the sensor can be established and the dynamic performance of the sensor can be improved by utilizing zero-pole assignment theory.Firstly,using the dynamic calibration test data of the pressure sensor to construct a rectangular data matrix,the model order of sensor was determined by using the method of matrix orthogonal decomposition.Secondly,the mathematical model of the sensor was established by using improved least squares algorithm.Lastly,the performance of the sensor was improved by using the zero-pole assignment theory to reconfigure the zero-pole of the model.The experimental results show that the method is simple and the modeling precision is high.A mathematical model which reflecting the dynamic perfor?mance of the sensor can be obtained accurately and the dynamic performance of the sensor can be improved sig?nificantly.

    pressure sensor;dynamic modeling;orthogonal decomposition;least square;zero-pole

    TP274

    A

    1004-1699(2016)11-1698-07

    EEACC:7230 10.3969/j.issn.1004-1699.2016.11.012

    2016-05-11 修改日期:2016-07-02

    猜你喜歡
    階次階躍補償
    基于階躍雙包層光纖的螺旋型光纖傳感器
    階次分析在驅動橋異響中的應用
    無功補償電容器的應用
    山東冶金(2019年5期)2019-11-16 09:09:38
    基于Vold-Kalman濾波的階次分析系統(tǒng)設計與實現(xiàn)*
    探討單位階躍信號的教學
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    價值工程(2017年28期)2018-01-23 20:48:29
    解讀補償心理
    植物補償和超補償作用
    離婚時,能否要求家務補償
    海峽姐妹(2014年5期)2014-02-27 15:09:41
    一種階躍函數(shù)在矩形時間窗口頻域特性的分析方法
    曰老女人黄片| 亚洲欧美一区二区三区黑人| 老司机午夜十八禁免费视频| 国产精品电影一区二区三区| 亚洲精品在线美女| 国产91精品成人一区二区三区| av视频在线观看入口| 国产国语露脸激情在线看| 亚洲情色 制服丝袜| 国产精品98久久久久久宅男小说| 久久精品aⅴ一区二区三区四区| av在线天堂中文字幕| 香蕉国产在线看| 欧美国产日韩亚洲一区| 亚洲欧美日韩另类电影网站| 国产私拍福利视频在线观看| 午夜免费激情av| 午夜久久久在线观看| 国产亚洲精品av在线| 日韩免费av在线播放| 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 成人18禁在线播放| 亚洲欧美激情在线| 人人妻人人爽人人添夜夜欢视频| 久久久国产成人免费| 两个人看的免费小视频| 色哟哟哟哟哟哟| 国产91精品成人一区二区三区| 亚洲精华国产精华精| 亚洲免费av在线视频| 亚洲国产欧美网| 国产欧美日韩一区二区三区在线| 国产精品日韩av在线免费观看 | 在线观看66精品国产| 国产又爽黄色视频| 久久人妻熟女aⅴ| av片东京热男人的天堂| 欧美性长视频在线观看| 一区福利在线观看| 亚洲国产欧美一区二区综合| 国产精品一区二区在线不卡| 欧美日韩黄片免| 国产91精品成人一区二区三区| 又大又爽又粗| 久久伊人香网站| 又黄又爽又免费观看的视频| 国产高清激情床上av| 亚洲一码二码三码区别大吗| 成人18禁在线播放| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 午夜成年电影在线免费观看| 国产精品亚洲一级av第二区| 十八禁网站免费在线| 成人手机av| 夜夜躁狠狠躁天天躁| 黑人巨大精品欧美一区二区mp4| 国产高清激情床上av| 一级片免费观看大全| 国产激情久久老熟女| 757午夜福利合集在线观看| 欧美av亚洲av综合av国产av| 久久亚洲真实| 两个人看的免费小视频| 波多野结衣高清无吗| 一本久久中文字幕| 国产精品永久免费网站| 老司机在亚洲福利影院| 天堂动漫精品| 精品无人区乱码1区二区| 成人精品一区二区免费| 国产真人三级小视频在线观看| 黄频高清免费视频| 国内久久婷婷六月综合欲色啪| 国语自产精品视频在线第100页| av欧美777| 日韩av在线大香蕉| 9色porny在线观看| www.自偷自拍.com| 久久久精品欧美日韩精品| 国产区一区二久久| 国产精品亚洲美女久久久| 97人妻精品一区二区三区麻豆 | 天堂动漫精品| 99精品欧美一区二区三区四区| 91麻豆精品激情在线观看国产| 免费av毛片视频| 欧美黄色片欧美黄色片| 一本久久中文字幕| 亚洲国产欧美网| 国产高清视频在线播放一区| 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 亚洲成a人片在线一区二区| 国产亚洲精品一区二区www| 色在线成人网| 男女下面插进去视频免费观看| 少妇的丰满在线观看| 亚洲最大成人中文| 日韩精品免费视频一区二区三区| 亚洲人成电影观看| 一区二区三区国产精品乱码| 巨乳人妻的诱惑在线观看| 男女下面插进去视频免费观看| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸| 国产精品二区激情视频| 97超级碰碰碰精品色视频在线观看| 99国产精品99久久久久| 亚洲国产看品久久| 国产色视频综合| 日韩大尺度精品在线看网址 | 亚洲成国产人片在线观看| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 日韩大码丰满熟妇| 色综合婷婷激情| 亚洲美女黄片视频| 桃红色精品国产亚洲av| 91精品国产国语对白视频| 久久久久亚洲av毛片大全| 亚洲av第一区精品v没综合| svipshipincom国产片| 欧美日韩黄片免| 90打野战视频偷拍视频| 成人亚洲精品av一区二区| 黑人巨大精品欧美一区二区mp4| 日本精品一区二区三区蜜桃| 国产精品秋霞免费鲁丝片| 午夜福利影视在线免费观看| 美女高潮喷水抽搐中文字幕| 99久久国产精品久久久| 亚洲专区字幕在线| av在线播放免费不卡| 在线观看午夜福利视频| 韩国av一区二区三区四区| 涩涩av久久男人的天堂| 日韩大码丰满熟妇| 亚洲专区字幕在线| 免费高清视频大片| 中文字幕色久视频| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| а√天堂www在线а√下载| 婷婷丁香在线五月| 激情视频va一区二区三区| 国产一区二区在线av高清观看| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲人成电影观看| 老司机深夜福利视频在线观看| 免费av毛片视频| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 欧美激情 高清一区二区三区| 在线观看日韩欧美| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| 国产麻豆69| 亚洲va日本ⅴa欧美va伊人久久| 男女做爰动态图高潮gif福利片 | 亚洲中文字幕一区二区三区有码在线看 | 很黄的视频免费| 色婷婷久久久亚洲欧美| 久久久水蜜桃国产精品网| 性欧美人与动物交配| 侵犯人妻中文字幕一二三四区| 免费高清视频大片| 午夜成年电影在线免费观看| 久久久久久人人人人人| 日本免费a在线| 国产av一区在线观看免费| 黄色 视频免费看| 久久中文字幕人妻熟女| 在线播放国产精品三级| 大型黄色视频在线免费观看| 高清毛片免费观看视频网站| 久久久久久久久中文| 国产xxxxx性猛交| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| 亚洲欧美激情综合另类| 亚洲国产高清在线一区二区三 | 一个人观看的视频www高清免费观看 | 每晚都被弄得嗷嗷叫到高潮| a级毛片在线看网站| 一级毛片精品| 久久天堂一区二区三区四区| 精品欧美一区二区三区在线| 午夜视频精品福利| 国产高清有码在线观看视频 | 欧美人与性动交α欧美精品济南到| 中文字幕精品免费在线观看视频| 久久精品影院6| 欧美一级毛片孕妇| 日日爽夜夜爽网站| 少妇被粗大的猛进出69影院| e午夜精品久久久久久久| 黄色成人免费大全| 国产一区二区三区在线臀色熟女| 久久人妻av系列| 人人妻,人人澡人人爽秒播| 久久人妻熟女aⅴ| 我的亚洲天堂| 在线观看一区二区三区| 国产麻豆69| 侵犯人妻中文字幕一二三四区| 国产精品爽爽va在线观看网站 | 天堂影院成人在线观看| 久久久国产欧美日韩av| 操出白浆在线播放| 一边摸一边抽搐一进一出视频| 亚洲av成人av| 亚洲男人的天堂狠狠| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 亚洲国产欧美网| 久久天堂一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 亚洲伊人色综图| 免费看十八禁软件| 不卡av一区二区三区| 99热只有精品国产| 亚洲第一电影网av| 国产一区在线观看成人免费| 免费在线观看亚洲国产| 国产精品综合久久久久久久免费 | 亚洲国产欧美网| 国产91精品成人一区二区三区| 亚洲一区高清亚洲精品| 此物有八面人人有两片| 午夜精品久久久久久毛片777| 少妇熟女aⅴ在线视频| 亚洲av美国av| 欧美一级毛片孕妇| 亚洲 欧美 日韩 在线 免费| 99在线人妻在线中文字幕| 桃红色精品国产亚洲av| 校园春色视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 午夜影院日韩av| 搡老妇女老女人老熟妇| 精品午夜福利视频在线观看一区| 欧美 亚洲 国产 日韩一| 欧美午夜高清在线| 色综合欧美亚洲国产小说| 一区二区三区高清视频在线| 亚洲人成电影观看| 一边摸一边抽搐一进一出视频| 午夜福利免费观看在线| 日韩大尺度精品在线看网址 | 精品国内亚洲2022精品成人| 亚洲男人天堂网一区| 午夜福利成人在线免费观看| 国产精品久久久av美女十八| 老司机午夜十八禁免费视频| 久久人人97超碰香蕉20202| 欧美成人一区二区免费高清观看 | 午夜免费激情av| 天天一区二区日本电影三级 | 久久久久国产精品人妻aⅴ院| 国产午夜精品久久久久久| 级片在线观看| 国产免费av片在线观看野外av| 国产精品亚洲av一区麻豆| 亚洲精品中文字幕在线视频| 久久亚洲真实| 欧美激情久久久久久爽电影 | 久久久久国产精品人妻aⅴ院| 露出奶头的视频| 9色porny在线观看| 免费av毛片视频| 亚洲欧美精品综合久久99| 久久久久久人人人人人| 国产熟女午夜一区二区三区| 亚洲五月色婷婷综合| 中文字幕另类日韩欧美亚洲嫩草| 91精品国产国语对白视频| 黄片小视频在线播放| 91麻豆av在线| 欧美日本亚洲视频在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美精品综合一区二区三区| 久久亚洲真实| 国产精品乱码一区二三区的特点 | 曰老女人黄片| 国产成人精品无人区| 国产熟女xx| 亚洲精品美女久久av网站| 国产精品 欧美亚洲| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 最好的美女福利视频网| 国产成人av教育| 亚洲熟女毛片儿| 久久亚洲真实| 757午夜福利合集在线观看| 极品教师在线免费播放| 久久久久久久久中文| 99在线人妻在线中文字幕| 成人18禁高潮啪啪吃奶动态图| 日韩精品免费视频一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 午夜影院日韩av| 黄色片一级片一级黄色片| 亚洲电影在线观看av| 麻豆av在线久日| 老熟妇乱子伦视频在线观看| 欧美成人性av电影在线观看| 超碰成人久久| 国内毛片毛片毛片毛片毛片| 99久久99久久久精品蜜桃| 欧美激情极品国产一区二区三区| 日日爽夜夜爽网站| 琪琪午夜伦伦电影理论片6080| 欧美一级a爱片免费观看看 | 91麻豆精品激情在线观看国产| 熟女少妇亚洲综合色aaa.| 成人三级黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人色综图| 日本一区二区免费在线视频| 亚洲 国产 在线| av福利片在线| 国产精品综合久久久久久久免费 | 亚洲欧美日韩另类电影网站| 国产成人影院久久av| 亚洲专区国产一区二区| 又大又爽又粗| 精品国产亚洲在线| 男人操女人黄网站| 国产一区二区三区视频了| 夜夜爽天天搞| 9191精品国产免费久久| 精品一品国产午夜福利视频| 日本免费a在线| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 19禁男女啪啪无遮挡网站| 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 一二三四社区在线视频社区8| 可以免费在线观看a视频的电影网站| 妹子高潮喷水视频| 国产亚洲精品久久久久5区| xxx96com| 一级a爱片免费观看的视频| 91国产中文字幕| 9191精品国产免费久久| 久久天堂一区二区三区四区| 欧美丝袜亚洲另类 | 一边摸一边抽搐一进一出视频| 久久久水蜜桃国产精品网| cao死你这个sao货| 淫妇啪啪啪对白视频| 色在线成人网| 亚洲av熟女| 午夜久久久在线观看| 黑人巨大精品欧美一区二区蜜桃| 狠狠狠狠99中文字幕| 国产1区2区3区精品| 一级a爱片免费观看的视频| 精品熟女少妇八av免费久了| 久久精品亚洲熟妇少妇任你| 人人妻人人爽人人添夜夜欢视频| 国产成人一区二区三区免费视频网站| 又黄又爽又免费观看的视频| 国产欧美日韩精品亚洲av| e午夜精品久久久久久久| 亚洲国产高清在线一区二区三 | 欧美老熟妇乱子伦牲交| 日本vs欧美在线观看视频| 亚洲电影在线观看av| 久久精品国产亚洲av高清一级| 看免费av毛片| 亚洲国产看品久久| 国产成人一区二区三区免费视频网站| 侵犯人妻中文字幕一二三四区| 久久中文字幕一级| 亚洲 欧美一区二区三区| 国产av一区二区精品久久| 精品一区二区三区四区五区乱码| 久久香蕉精品热| 狠狠狠狠99中文字幕| 国产精品亚洲美女久久久| 黄频高清免费视频| 国产精品99久久99久久久不卡| 一级,二级,三级黄色视频| 黄色成人免费大全| 成人手机av| 亚洲成人久久性| 国产精品美女特级片免费视频播放器 | 母亲3免费完整高清在线观看| 曰老女人黄片| 一区二区三区精品91| 欧美亚洲日本最大视频资源| 久久久精品国产亚洲av高清涩受| 此物有八面人人有两片| 中文字幕色久视频| 久久精品国产亚洲av香蕉五月| svipshipincom国产片| 久久精品国产清高在天天线| 91老司机精品| 国产精品九九99| 91在线观看av| 日韩av在线大香蕉| 日本黄色视频三级网站网址| 欧美日韩一级在线毛片| 婷婷精品国产亚洲av在线| 久久精品国产99精品国产亚洲性色 | 精品免费久久久久久久清纯| 亚洲欧美精品综合一区二区三区| 午夜精品在线福利| 神马国产精品三级电影在线观看 | 国产精品美女特级片免费视频播放器 | tocl精华| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 精品久久久久久,| 欧美另类亚洲清纯唯美| 91麻豆精品激情在线观看国产| 黄色a级毛片大全视频| 亚洲中文av在线| 中文字幕最新亚洲高清| 精品乱码久久久久久99久播| 久久青草综合色| 久久精品国产清高在天天线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲第一青青草原| 国产精品野战在线观看| 亚洲精品美女久久av网站| 日韩欧美国产在线观看| 亚洲久久久国产精品| 精品久久久久久久人妻蜜臀av | 99国产精品99久久久久| 搡老妇女老女人老熟妇| 在线观看免费视频日本深夜| 一区在线观看完整版| 亚洲伊人色综图| 精品午夜福利视频在线观看一区| 国产精品免费一区二区三区在线| 91精品国产国语对白视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看66精品国产| 真人一进一出gif抽搐免费| 久久久久亚洲av毛片大全| 免费在线观看亚洲国产| avwww免费| 国产精品一区二区免费欧美| 国产亚洲欧美在线一区二区| 熟妇人妻久久中文字幕3abv| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 性少妇av在线| 久久婷婷人人爽人人干人人爱 | 搡老妇女老女人老熟妇| 亚洲欧美激情在线| 精品一区二区三区四区五区乱码| 男女下面进入的视频免费午夜 | 91在线观看av| 日本 欧美在线| 欧美亚洲日本最大视频资源| 成人国语在线视频| 18禁国产床啪视频网站| 亚洲国产欧美网| 欧美在线一区亚洲| 精品一区二区三区视频在线观看免费| 日本vs欧美在线观看视频| 日韩精品免费视频一区二区三区| 怎么达到女性高潮| 日本免费a在线| 精品不卡国产一区二区三区| 无人区码免费观看不卡| 精品一品国产午夜福利视频| 免费在线观看日本一区| 91av网站免费观看| 久久香蕉精品热| 久久久精品国产亚洲av高清涩受| av欧美777| 一级毛片女人18水好多| 香蕉丝袜av| av免费在线观看网站| 后天国语完整版免费观看| 亚洲伊人色综图| 午夜成年电影在线免费观看| 国产一区二区激情短视频| 久久精品亚洲熟妇少妇任你| 久久九九热精品免费| 此物有八面人人有两片| 久久伊人香网站| 国产精品九九99| 久久精品亚洲熟妇少妇任你| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看| 高清在线国产一区| 久久久久久国产a免费观看| 久久久久国产精品人妻aⅴ院| 十八禁人妻一区二区| 91字幕亚洲| 亚洲男人的天堂狠狠| 亚洲精品美女久久av网站| 不卡av一区二区三区| 久久九九热精品免费| 午夜影院日韩av| 日韩欧美国产在线观看| 亚洲性夜色夜夜综合| 久久国产精品影院| 国产极品粉嫩免费观看在线| 久久亚洲精品不卡| 麻豆一二三区av精品| 夜夜爽天天搞| 嫩草影视91久久| 日本在线视频免费播放| 久久久久久国产a免费观看| 激情在线观看视频在线高清| 成人手机av| 天天躁狠狠躁夜夜躁狠狠躁| 久久人妻av系列| 国内毛片毛片毛片毛片毛片| 久久精品91蜜桃| 性欧美人与动物交配| 88av欧美| 在线十欧美十亚洲十日本专区| 午夜精品在线福利| 免费少妇av软件| 亚洲国产精品久久男人天堂| 91麻豆精品激情在线观看国产| 中文字幕最新亚洲高清| 宅男免费午夜| 午夜成年电影在线免费观看| 99re在线观看精品视频| aaaaa片日本免费| 丝袜在线中文字幕| 欧美日韩一级在线毛片| 午夜福利18| 精品一区二区三区av网在线观看| 亚洲精品中文字幕一二三四区| 亚洲精品中文字幕在线视频| 成人三级做爰电影| 亚洲中文av在线| 国产精品免费视频内射| 日本vs欧美在线观看视频| 18禁国产床啪视频网站| 国产免费男女视频| 成年人黄色毛片网站| 啦啦啦免费观看视频1| 成年人黄色毛片网站| 亚洲精品一区av在线观看| 久久国产精品人妻蜜桃| 51午夜福利影视在线观看| 免费在线观看影片大全网站| 亚洲男人的天堂狠狠| 久久久国产精品麻豆| 麻豆成人av在线观看| 国产麻豆成人av免费视频| 色综合亚洲欧美另类图片| 91九色精品人成在线观看| 国产色视频综合| 中文字幕人成人乱码亚洲影| 精品久久久久久久人妻蜜臀av | 亚洲一区二区三区不卡视频| 亚洲色图 男人天堂 中文字幕| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久人妻蜜臀av | 成人国产综合亚洲| 美女 人体艺术 gogo| 久久中文看片网| 一本综合久久免费| 欧美人与性动交α欧美精品济南到| 男女下面进入的视频免费午夜 | 九色亚洲精品在线播放| 精品国产一区二区久久| 在线免费观看的www视频| 一本大道久久a久久精品| 欧美成狂野欧美在线观看| 国产精品香港三级国产av潘金莲| 色综合婷婷激情| 男人的好看免费观看在线视频 | 国产欧美日韩一区二区三区在线| www.熟女人妻精品国产| 欧美日韩乱码在线| 欧美一区二区精品小视频在线| 欧美另类亚洲清纯唯美| 精品久久久久久久毛片微露脸| 亚洲国产高清在线一区二区三 | 久久青草综合色| 中亚洲国语对白在线视频| 亚洲美女黄片视频| 国产激情欧美一区二区| 日本免费a在线| 欧美一区二区精品小视频在线| 国产片内射在线| 中文字幕另类日韩欧美亚洲嫩草| 一区二区三区高清视频在线| 亚洲精品国产色婷婷电影| 精品免费久久久久久久清纯| 伊人久久大香线蕉亚洲五| 久久热在线av| 免费人成视频x8x8入口观看| 男女午夜视频在线观看| 免费观看精品视频网站| 18美女黄网站色大片免费观看| 亚洲欧美日韩无卡精品| 久久国产亚洲av麻豆专区| 大香蕉久久成人网| 在线观看免费视频网站a站| 亚洲国产欧美网|