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

    先導(dǎo)式定壓減壓閥CAT系統(tǒng)的數(shù)據(jù)處理與誤差分析

    2014-11-10 14:41:29彭明哲李群松
    科技創(chuàng)新導(dǎo)報(bào) 2014年12期
    關(guān)鍵詞:誤差分析減壓閥數(shù)據(jù)處理

    彭明哲 李群松

    摘 要:以先導(dǎo)式減壓閥壓力傳感器的標(biāo)定、減壓閥的卸壓-建壓特性試驗(yàn)為例,對(duì)先導(dǎo)式定壓減壓閥CAT系統(tǒng)的數(shù)據(jù)處理與誤差分析做了詳細(xì)的闡述。給出了試驗(yàn)中傳感器標(biāo)定方法及數(shù)據(jù)預(yù)處理的一種具體算法。結(jié)果表明:將正交多項(xiàng)式回歸應(yīng)用于動(dòng)、靜態(tài)特性試驗(yàn)的數(shù)據(jù)處理中,回歸次數(shù)越高,回歸曲線越逼近真實(shí)曲線,正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差就越小,數(shù)據(jù)處理的精度就越高。該法可以應(yīng)用于相關(guān)的數(shù)據(jù)處理中。

    關(guān)鍵詞:減壓閥 傳感器 數(shù)據(jù)處理 誤差分析

    中圖分類(lèi)號(hào):TH137 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1674-098X(2014)04(c)-0049-04

    壓力控制閥幾乎在每個(gè)液壓系統(tǒng)中都要用到,其性能的好壞對(duì)整個(gè)系統(tǒng)的正常工作有很大影響。定壓減壓閥應(yīng)用最廣,它是用來(lái)控制出口壓力為定值,使液壓系統(tǒng)中某一部分得到較供油壓力低的穩(wěn)定壓力。而在壓力測(cè)定的過(guò)程中,不可能獲得絕對(duì)準(zhǔn)確的測(cè)量結(jié)果。因此我們通過(guò)測(cè)量所獲得的原始數(shù)據(jù),都含有誤差。只有將這些含有誤差的數(shù)據(jù)進(jìn)行數(shù)學(xué)加工處理,盡量消除誤差,使測(cè)量結(jié)果更加可靠,并估計(jì)其精確程度[1]。文章就先導(dǎo)式壓力減壓閥的數(shù)據(jù)處理提出了回歸等數(shù)據(jù)預(yù)處理和處理等辦法,以盡量消除測(cè)試誤差,使采集的數(shù)據(jù)更真實(shí)有效,以提高系統(tǒng)的測(cè)試精度。

    1 先導(dǎo)式定壓減壓閥工作原理

    圖1是先導(dǎo)式定壓減壓閥的工作原理圖。進(jìn)口壓力p經(jīng)減壓縫隙減壓后,壓力變?yōu)閜1(即出口壓力),出口壓力油經(jīng)主閥芯滑閥之軸向槽a進(jìn)入滑閥底部,然后經(jīng)滑閥中部軸向孔進(jìn)入滑閥上端油腔,進(jìn)入上腔的壓力油再經(jīng)閥蓋上通孔b、先導(dǎo)閥閥座孔c,作用在先導(dǎo)閥錐閥上。當(dāng)出口壓力低于調(diào)定壓力時(shí),先導(dǎo)閥在調(diào)壓彈簧的作用下關(guān)閉閥口,滑閥上下腔的油壓均等于出口壓力,滑閥在彈簧力的作用下處于最下端位置,滑閥中間凸間與閥體之間構(gòu)成的減壓閥閥口全開(kāi),不起減壓作用。當(dāng)出口壓力上升至調(diào)定壓力時(shí),先導(dǎo)閥錐閥在油壓作用下被打開(kāi)溢流,流入先導(dǎo)閥彈簧腔的液體經(jīng)上蓋上的泄油孔(圖中未畫(huà)出)流回油箱,出口壓力油經(jīng)上述通道不斷補(bǔ)充、最后經(jīng)泄油孔流回油箱。由于滑閥中部軸向孔為阻尼孔,因此油液經(jīng)該孔時(shí)有壓力損失存在,滑閥上腔壓力p2低于滑閥下腔壓力p1。于是滑閥在上下兩端壓力差的作用下克服彈簧力向上運(yùn)動(dòng),使得減壓閥閥口減小,運(yùn)動(dòng)到某一位置,作用在滑閥上的所有力處于平衡,記此平衡位置的減壓閥閥口減壓縫隙長(zhǎng)度為y。因?yàn)殚y口的減壓作用,使減壓閥的出口壓力低于進(jìn)口壓力。此時(shí),先導(dǎo)閥錐閥在油壓p2和先導(dǎo)閥調(diào)壓彈簧的彈簧力以及閥口液動(dòng)力的共同作用下處于某一平衡位置,閥口開(kāi)度為x。如果先導(dǎo)閥調(diào)壓彈簧的彈簧力調(diào)定,則先導(dǎo)閥前的壓力p2為定值,相應(yīng)閥的出口壓力亦為定值,此值即為減壓閥的調(diào)定壓力。各參數(shù)的測(cè)定均通過(guò)計(jì)算機(jī)輔助測(cè)試系統(tǒng)來(lái)完成。

    2 CAT系統(tǒng)中傳感器的數(shù)字標(biāo)定

    為便于計(jì)算機(jī)數(shù)據(jù)處理,可通過(guò)標(biāo)定幾組有效的值來(lái)對(duì)傳感器進(jìn)行標(biāo)定,并把標(biāo)定值保存到文件,以便后繼試驗(yàn)時(shí)調(diào)用。傳感器的數(shù)字標(biāo)定可以避免傳感器零點(diǎn)的校正、便于上、下位機(jī)間的數(shù)據(jù)傳輸。在本CAT系統(tǒng)中的實(shí)際值與傳感器所測(cè)量的上、下位機(jī)中的標(biāo)示值(也稱為測(cè)量值)符合線性關(guān)系,可以用回歸分

    析法來(lái)進(jìn)行傳感器的標(biāo)定。假設(shè)實(shí)際值為Y,標(biāo)示值為X,這樣X(jué)和Y就可用式1來(lái)表達(dá):

    (1)

    為使式2的E值達(dá)到最小,根據(jù)最小二乘法,可以確定式1的系數(shù)a,b的值

    (2)

    其中(xi,yi),i=0,1,…n-1,為(X,Y)的n組有效標(biāo)定值。

    根據(jù)極值原理,a 與b應(yīng)滿足下列等式(3):

    (3)

    從式(3)解得:

    ,

    (4)

    式4中:,

    根據(jù)以上的數(shù)據(jù)處理方法,下面給出了本CAT系統(tǒng)的測(cè)試回路中A口的壓力傳感器標(biāo)定的演示,如圖2所示。圖3給出液壓元件CAT系統(tǒng)各壓力傳感器的標(biāo)定的曲線,并計(jì)算得出相應(yīng)的標(biāo)定誤差的標(biāo)準(zhǔn)差。

    ,,, 。

    3 數(shù)據(jù)處理

    在壓力控制閥的靜特性試驗(yàn)中,其測(cè)試的數(shù)據(jù)是序列觀測(cè)值,要從實(shí)際的序列觀測(cè)數(shù)據(jù)中盡可能排除噪音等成分而讓真實(shí)信號(hào)保留下來(lái),在此采用數(shù)字濾波方法。當(dāng)測(cè)試數(shù)據(jù)足夠致密,且每一測(cè)試點(diǎn)附近二階導(dǎo)數(shù)為常數(shù)時(shí),即幾何上反映為曲線的凹凸方向一致時(shí),數(shù)字濾波可取得較好的平滑效果[2-4]。本系統(tǒng)采用非線性滑動(dòng)平滑方法,具體實(shí)現(xiàn)如下。

    設(shè)在某一點(diǎn)的左右兩側(cè)各取兩個(gè)相鄰數(shù)據(jù)點(diǎn),橫坐標(biāo)值分別是x=-2,-1,0,1,2,設(shè)定函數(shù)對(duì)自變量

    x的二階導(dǎo)數(shù)為常數(shù),可取擬合函數(shù):作為被測(cè)信號(hào)在〔 -2,2〕上的近似。

    據(jù)最小二乘原理 ,誤差的平方和為:

    由,得:

    ,

    我們的目的是校正x=0處的y0的實(shí)際值,得y0的估計(jì)值為:

    (5)

    式5即為五點(diǎn)加權(quán)平滑濾波公式,它的實(shí)質(zhì)是在小區(qū)間[-2,2]上的曲線擬合問(wèn)題,顯然具有局部性的特點(diǎn),這一特點(diǎn)決定了它不能適用于大幅值干擾的場(chǎng)合;另外,由于推導(dǎo)的前提條件是二階導(dǎo)數(shù)為常數(shù),對(duì)有周期性脈動(dòng)干擾的處理,效果必然不佳。它適合于噪聲較微弱,并且不存在高頻周期性干擾的緩變信號(hào),故我們用于靜態(tài)試驗(yàn)數(shù)據(jù)預(yù)處理的濾波中[5-6]。

    運(yùn)用此平滑濾波公式對(duì)先導(dǎo)式定壓減壓閥的穩(wěn)態(tài)壓力-流量特性試驗(yàn)數(shù)據(jù)進(jìn)行處理,其處理前后試驗(yàn)曲線如圖4所示。圖4(a)為濾波前的原始數(shù)據(jù)曲線,圖4(b)為濾波后的數(shù)據(jù)曲線。處理后的曲線明顯光滑,去除了原始曲線中的毛刺。

    在測(cè)試數(shù)據(jù)的回歸分析中,為了求得回歸參數(shù)的大小,最常用的是采用最小二乘法。這是因?yàn)榫€性回歸模型在符合Gauss-Markov假定的條件下,采用最小二乘法估計(jì)其回歸參數(shù)具有良好的統(tǒng)計(jì)性質(zhì),如無(wú)偏性、一致性等[7]。endprint

    為不失一般性,假設(shè)要擬合的動(dòng)態(tài)測(cè)試數(shù)據(jù)對(duì)為(i=0,1,…,N-1),擬合函數(shù)為f(x,),它是的非線性函數(shù);,為待求的參數(shù)向量。定義殘差向量e = [e,e,…,e]。其中y,i = 0,1,…,N-1,通常N>m。要求選擇使得擬合函數(shù)f(x,)在某種準(zhǔn)則函數(shù)意義上盡可能好地?cái)M合數(shù)據(jù)。定義準(zhǔn)則函數(shù)為殘差向量e的l范數(shù),即:

    (6)

    極小化式6求得參數(shù),就得到所謂的數(shù)據(jù)擬合問(wèn)題。為了求解的方便,極小化式6等價(jià)于極小化如下準(zhǔn)則函數(shù):

    ,p=1,2, (7)

    當(dāng)p=2 時(shí),就是我們要求的二次準(zhǔn)則函數(shù):

    對(duì)于測(cè)試定壓減壓閥的靜特性和動(dòng)特試驗(yàn)曲線,很難找到線性化的函數(shù)形式來(lái)表達(dá),而用多項(xiàng)式逼近是一種有力的工具。然而用一般多項(xiàng)式回歸,存在著諸如計(jì)算趨于復(fù)雜等一些缺點(diǎn)[8]。一種有效的克服辦法是采用正交多項(xiàng)式的回歸分析。它的法方程組系數(shù)矩陣為對(duì)角陣,所以不會(huì)出現(xiàn)病態(tài)問(wèn)題,其逆陣也是對(duì)角陣,且求逆非常方便。由正交多項(xiàng)式組成的回歸方程,除了計(jì)算回歸系數(shù)方便外,還有一個(gè)重要優(yōu)點(diǎn),從已有的回歸方程中刪去或增加一個(gè)自變量的多項(xiàng)式函數(shù)后,其余的自變量多項(xiàng)式函數(shù)的回歸系數(shù)都沒(méi)有變化。下面是基于二次準(zhǔn)則下的正交多項(xiàng)式的回歸,其基本實(shí)現(xiàn)方法如下:

    設(shè)已知N數(shù)據(jù)點(diǎn)(i=0,1,…,N-1),求(m-1)次最小二乘擬合多項(xiàng)式

    (8)

    這一多項(xiàng)式函數(shù)組的生成步驟為:

    可以證明,由上述遞推構(gòu)造的多項(xiàng)式函數(shù)組(j=0,1,…,m-1)是相互正交的[9]。根據(jù)最小二乘原理,可得的表達(dá)式為:,j=0,1,…,m-1

    最后可以化成一般的m-1次多項(xiàng)式為:

    以上算法便于計(jì)算機(jī)實(shí)現(xiàn)。在實(shí)際計(jì)算過(guò)程中,為了防止運(yùn)算溢出,用代替,其中, i=0,1,..,N-1。

    擬合多項(xiàng)式如下:

    在數(shù)據(jù)處理的軟件模塊中,設(shè)置了交互式的數(shù)據(jù)處理模式,即對(duì)采集的試驗(yàn)數(shù)據(jù)進(jìn)行多項(xiàng)式回歸時(shí),根據(jù)用戶輸入要求進(jìn)行m次的回歸,軟件并將所得的各多項(xiàng)式要進(jìn)行顯著水平的檢驗(yàn)和正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差計(jì)算,直至達(dá)到試驗(yàn)要求為止。

    運(yùn)用此正交多項(xiàng)式回歸分析,對(duì)減壓閥出口流量階躍壓力響應(yīng)特性試驗(yàn)數(shù)據(jù)進(jìn)行處理,經(jīng)過(guò)6次回歸,得到其處理前后的曲線如圖5所示。

    4 CAT系統(tǒng)的誤差分析

    一個(gè)CAT系統(tǒng)的測(cè)試誤差來(lái)源是多方面的,實(shí)際情況很復(fù)雜。主要包括硬件、軟件、及數(shù)據(jù)處理幾個(gè)方面。對(duì)于傳感器、A/D轉(zhuǎn)換引起的誤差、上、下位機(jī)間數(shù)碼轉(zhuǎn)換產(chǎn)生的誤差,其相對(duì)誤差很小,在此不再討論。主要對(duì)應(yīng)用正交多項(xiàng)式回歸分析處理數(shù)據(jù)時(shí)產(chǎn)生的誤差進(jìn)行分析。

    以減壓閥的卸壓-建壓特性試驗(yàn)為例,具體分析正交多項(xiàng)式回歸的誤差分析。在卸壓-建壓特性曲線中,右邊的建壓特性曲線是用正交多項(xiàng)式最小二乘擬合的,曲線如圖6所示。圖6(a)為擬合前的曲線,圖6(b)為擬合后的曲線。

    在此試驗(yàn)中,右邊建壓曲線所采集的點(diǎn)數(shù)N=625,回歸次數(shù)為m-1=5次。則其[10]:

    總離差平方和為:,回歸差平方和為:jy

    第j個(gè)多項(xiàng)式的偏回歸平方和:,剩余平方和:

    其中為正交多項(xiàng)式回歸的系數(shù)(前已推導(dǎo));為:

    方差比:

    其中:的相應(yīng)自由度為,所以偏回歸方差估計(jì)值為:

    剩余平方和的相應(yīng)自由度為,所以剩余的方差估計(jì)值為:

    根據(jù)以上計(jì)算,可以得到以下正交多項(xiàng)式回歸的方差分析,如表1所示:

    經(jīng)查F檢驗(yàn)的臨界值(Fa)表,取a=0.01,f1=1,f2=(因?yàn)閒2=618>500)。

    在方差分析表1中,F(xiàn)計(jì)算值全都大于,故認(rèn)為都是高度顯著的,所以表中注以“**”,即表示1-5次項(xiàng)都是高度顯著,則所配制的多項(xiàng)式能夠正確地代表試驗(yàn)結(jié)果。其正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差為

    這樣的精度是比較滿意的。如果對(duì)數(shù)據(jù)處理的結(jié)果不滿意,還可以提高正交多項(xiàng)式的回歸次數(shù)來(lái)更加逼近測(cè)試數(shù)據(jù),獲得更高的精度。

    5 結(jié)語(yǔ)

    (1)在CAT系統(tǒng)中,數(shù)據(jù)測(cè)量的誤差來(lái)自多個(gè)方面。其中,由A/D轉(zhuǎn)換引起的誤差以及上、下位機(jī)間數(shù)據(jù)的數(shù)碼轉(zhuǎn)換引起的誤差極其微小,可以忽略。

    (2)在壓力控制閥的靜特性試驗(yàn)中,當(dāng)測(cè)試數(shù)據(jù)足夠致密,且每一測(cè)試點(diǎn)附近二階導(dǎo)數(shù)為常數(shù)時(shí),即幾何上反映為曲線的凹凸方向一致時(shí),對(duì)采集的數(shù)據(jù)采用非線性滑動(dòng)平滑方法,可以取得比較好的平滑效果。

    (3)在用正交多項(xiàng)式回歸試驗(yàn)曲線時(shí),回歸次數(shù)越高,回歸曲線越逼近真實(shí)曲線,正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差就越小,數(shù)據(jù)處理的精度就越高,可以根據(jù)需要選擇合適的回歸次數(shù),以獲得所需的測(cè)量精度。

    參考文獻(xiàn)

    [1] Malrino A. Digital principles and applications[M]. New York: McGraw-Hill, 1981.

    [2] 段長(zhǎng)寶.液壓測(cè)試[M].國(guó)防工業(yè)出版社,1984.

    [3] 朱堅(jiān)民,周福章.溢流閥計(jì)算機(jī)輔助測(cè)試及數(shù)據(jù)處理[J].機(jī)床與液壓,1998(3):56-58.

    [4] Khalatova L.V., Accuracy Estimation for Automated Computer- measurement Process[J].Meas Tech,1985,1985,123(3).

    [5] L.R. Rabiner. Application of digital signal processing[M]. New Jersey: Prentice-Hall Inc.,1975.

    [6] R.F.Coughlin and F.F.Driscoll: Operational Amplifiers and Linear Integrated Circuits[M].Prentice Hall Inc.,New Jersey,1977.

    [7] D.E.Johnson and J.L.Hilburn: Rapid Practical Designs of Active Filters[M]. John Wiley & Sons, Inc., New York,1975.

    [8] 周江文,黃幼才.抗差最小二乘法[M].華中理工大學(xué)出版社,1997.

    [9] P.J. Torpey. Minicomputericing analog data collection[M].USA: McGraw-Hill, 1980.

    [10] 胡上序,陳德釗.觀測(cè)數(shù)據(jù)的分析與處理[M].浙江:浙江大學(xué)出版社,1996.endprint

    為不失一般性,假設(shè)要擬合的動(dòng)態(tài)測(cè)試數(shù)據(jù)對(duì)為(i=0,1,…,N-1),擬合函數(shù)為f(x,),它是的非線性函數(shù);,為待求的參數(shù)向量。定義殘差向量e = [e,e,…,e]。其中y,i = 0,1,…,N-1,通常N>m。要求選擇使得擬合函數(shù)f(x,)在某種準(zhǔn)則函數(shù)意義上盡可能好地?cái)M合數(shù)據(jù)。定義準(zhǔn)則函數(shù)為殘差向量e的l范數(shù),即:

    (6)

    極小化式6求得參數(shù),就得到所謂的數(shù)據(jù)擬合問(wèn)題。為了求解的方便,極小化式6等價(jià)于極小化如下準(zhǔn)則函數(shù):

    ,p=1,2, (7)

    當(dāng)p=2 時(shí),就是我們要求的二次準(zhǔn)則函數(shù):

    對(duì)于測(cè)試定壓減壓閥的靜特性和動(dòng)特試驗(yàn)曲線,很難找到線性化的函數(shù)形式來(lái)表達(dá),而用多項(xiàng)式逼近是一種有力的工具。然而用一般多項(xiàng)式回歸,存在著諸如計(jì)算趨于復(fù)雜等一些缺點(diǎn)[8]。一種有效的克服辦法是采用正交多項(xiàng)式的回歸分析。它的法方程組系數(shù)矩陣為對(duì)角陣,所以不會(huì)出現(xiàn)病態(tài)問(wèn)題,其逆陣也是對(duì)角陣,且求逆非常方便。由正交多項(xiàng)式組成的回歸方程,除了計(jì)算回歸系數(shù)方便外,還有一個(gè)重要優(yōu)點(diǎn),從已有的回歸方程中刪去或增加一個(gè)自變量的多項(xiàng)式函數(shù)后,其余的自變量多項(xiàng)式函數(shù)的回歸系數(shù)都沒(méi)有變化。下面是基于二次準(zhǔn)則下的正交多項(xiàng)式的回歸,其基本實(shí)現(xiàn)方法如下:

    設(shè)已知N數(shù)據(jù)點(diǎn)(i=0,1,…,N-1),求(m-1)次最小二乘擬合多項(xiàng)式

    (8)

    這一多項(xiàng)式函數(shù)組的生成步驟為:

    可以證明,由上述遞推構(gòu)造的多項(xiàng)式函數(shù)組(j=0,1,…,m-1)是相互正交的[9]。根據(jù)最小二乘原理,可得的表達(dá)式為:,j=0,1,…,m-1

    最后可以化成一般的m-1次多項(xiàng)式為:

    以上算法便于計(jì)算機(jī)實(shí)現(xiàn)。在實(shí)際計(jì)算過(guò)程中,為了防止運(yùn)算溢出,用代替,其中, i=0,1,..,N-1。

    擬合多項(xiàng)式如下:

    在數(shù)據(jù)處理的軟件模塊中,設(shè)置了交互式的數(shù)據(jù)處理模式,即對(duì)采集的試驗(yàn)數(shù)據(jù)進(jìn)行多項(xiàng)式回歸時(shí),根據(jù)用戶輸入要求進(jìn)行m次的回歸,軟件并將所得的各多項(xiàng)式要進(jìn)行顯著水平的檢驗(yàn)和正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差計(jì)算,直至達(dá)到試驗(yàn)要求為止。

    運(yùn)用此正交多項(xiàng)式回歸分析,對(duì)減壓閥出口流量階躍壓力響應(yīng)特性試驗(yàn)數(shù)據(jù)進(jìn)行處理,經(jīng)過(guò)6次回歸,得到其處理前后的曲線如圖5所示。

    4 CAT系統(tǒng)的誤差分析

    一個(gè)CAT系統(tǒng)的測(cè)試誤差來(lái)源是多方面的,實(shí)際情況很復(fù)雜。主要包括硬件、軟件、及數(shù)據(jù)處理幾個(gè)方面。對(duì)于傳感器、A/D轉(zhuǎn)換引起的誤差、上、下位機(jī)間數(shù)碼轉(zhuǎn)換產(chǎn)生的誤差,其相對(duì)誤差很小,在此不再討論。主要對(duì)應(yīng)用正交多項(xiàng)式回歸分析處理數(shù)據(jù)時(shí)產(chǎn)生的誤差進(jìn)行分析。

    以減壓閥的卸壓-建壓特性試驗(yàn)為例,具體分析正交多項(xiàng)式回歸的誤差分析。在卸壓-建壓特性曲線中,右邊的建壓特性曲線是用正交多項(xiàng)式最小二乘擬合的,曲線如圖6所示。圖6(a)為擬合前的曲線,圖6(b)為擬合后的曲線。

    在此試驗(yàn)中,右邊建壓曲線所采集的點(diǎn)數(shù)N=625,回歸次數(shù)為m-1=5次。則其[10]:

    總離差平方和為:,回歸差平方和為:jy

    第j個(gè)多項(xiàng)式的偏回歸平方和:,剩余平方和:

    其中為正交多項(xiàng)式回歸的系數(shù)(前已推導(dǎo));為:

    方差比:

    其中:的相應(yīng)自由度為,所以偏回歸方差估計(jì)值為:

    剩余平方和的相應(yīng)自由度為,所以剩余的方差估計(jì)值為:

    根據(jù)以上計(jì)算,可以得到以下正交多項(xiàng)式回歸的方差分析,如表1所示:

    經(jīng)查F檢驗(yàn)的臨界值(Fa)表,取a=0.01,f1=1,f2=(因?yàn)閒2=618>500)。

    在方差分析表1中,F(xiàn)計(jì)算值全都大于,故認(rèn)為都是高度顯著的,所以表中注以“**”,即表示1-5次項(xiàng)都是高度顯著,則所配制的多項(xiàng)式能夠正確地代表試驗(yàn)結(jié)果。其正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差為

    這樣的精度是比較滿意的。如果對(duì)數(shù)據(jù)處理的結(jié)果不滿意,還可以提高正交多項(xiàng)式的回歸次數(shù)來(lái)更加逼近測(cè)試數(shù)據(jù),獲得更高的精度。

    5 結(jié)語(yǔ)

    (1)在CAT系統(tǒng)中,數(shù)據(jù)測(cè)量的誤差來(lái)自多個(gè)方面。其中,由A/D轉(zhuǎn)換引起的誤差以及上、下位機(jī)間數(shù)據(jù)的數(shù)碼轉(zhuǎn)換引起的誤差極其微小,可以忽略。

    (2)在壓力控制閥的靜特性試驗(yàn)中,當(dāng)測(cè)試數(shù)據(jù)足夠致密,且每一測(cè)試點(diǎn)附近二階導(dǎo)數(shù)為常數(shù)時(shí),即幾何上反映為曲線的凹凸方向一致時(shí),對(duì)采集的數(shù)據(jù)采用非線性滑動(dòng)平滑方法,可以取得比較好的平滑效果。

    (3)在用正交多項(xiàng)式回歸試驗(yàn)曲線時(shí),回歸次數(shù)越高,回歸曲線越逼近真實(shí)曲線,正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差就越小,數(shù)據(jù)處理的精度就越高,可以根據(jù)需要選擇合適的回歸次數(shù),以獲得所需的測(cè)量精度。

    參考文獻(xiàn)

    [1] Malrino A. Digital principles and applications[M]. New York: McGraw-Hill, 1981.

    [2] 段長(zhǎng)寶.液壓測(cè)試[M].國(guó)防工業(yè)出版社,1984.

    [3] 朱堅(jiān)民,周福章.溢流閥計(jì)算機(jī)輔助測(cè)試及數(shù)據(jù)處理[J].機(jī)床與液壓,1998(3):56-58.

    [4] Khalatova L.V., Accuracy Estimation for Automated Computer- measurement Process[J].Meas Tech,1985,1985,123(3).

    [5] L.R. Rabiner. Application of digital signal processing[M]. New Jersey: Prentice-Hall Inc.,1975.

    [6] R.F.Coughlin and F.F.Driscoll: Operational Amplifiers and Linear Integrated Circuits[M].Prentice Hall Inc.,New Jersey,1977.

    [7] D.E.Johnson and J.L.Hilburn: Rapid Practical Designs of Active Filters[M]. John Wiley & Sons, Inc., New York,1975.

    [8] 周江文,黃幼才.抗差最小二乘法[M].華中理工大學(xué)出版社,1997.

    [9] P.J. Torpey. Minicomputericing analog data collection[M].USA: McGraw-Hill, 1980.

    [10] 胡上序,陳德釗.觀測(cè)數(shù)據(jù)的分析與處理[M].浙江:浙江大學(xué)出版社,1996.endprint

    為不失一般性,假設(shè)要擬合的動(dòng)態(tài)測(cè)試數(shù)據(jù)對(duì)為(i=0,1,…,N-1),擬合函數(shù)為f(x,),它是的非線性函數(shù);,為待求的參數(shù)向量。定義殘差向量e = [e,e,…,e]。其中y,i = 0,1,…,N-1,通常N>m。要求選擇使得擬合函數(shù)f(x,)在某種準(zhǔn)則函數(shù)意義上盡可能好地?cái)M合數(shù)據(jù)。定義準(zhǔn)則函數(shù)為殘差向量e的l范數(shù),即:

    (6)

    極小化式6求得參數(shù),就得到所謂的數(shù)據(jù)擬合問(wèn)題。為了求解的方便,極小化式6等價(jià)于極小化如下準(zhǔn)則函數(shù):

    ,p=1,2, (7)

    當(dāng)p=2 時(shí),就是我們要求的二次準(zhǔn)則函數(shù):

    對(duì)于測(cè)試定壓減壓閥的靜特性和動(dòng)特試驗(yàn)曲線,很難找到線性化的函數(shù)形式來(lái)表達(dá),而用多項(xiàng)式逼近是一種有力的工具。然而用一般多項(xiàng)式回歸,存在著諸如計(jì)算趨于復(fù)雜等一些缺點(diǎn)[8]。一種有效的克服辦法是采用正交多項(xiàng)式的回歸分析。它的法方程組系數(shù)矩陣為對(duì)角陣,所以不會(huì)出現(xiàn)病態(tài)問(wèn)題,其逆陣也是對(duì)角陣,且求逆非常方便。由正交多項(xiàng)式組成的回歸方程,除了計(jì)算回歸系數(shù)方便外,還有一個(gè)重要優(yōu)點(diǎn),從已有的回歸方程中刪去或增加一個(gè)自變量的多項(xiàng)式函數(shù)后,其余的自變量多項(xiàng)式函數(shù)的回歸系數(shù)都沒(méi)有變化。下面是基于二次準(zhǔn)則下的正交多項(xiàng)式的回歸,其基本實(shí)現(xiàn)方法如下:

    設(shè)已知N數(shù)據(jù)點(diǎn)(i=0,1,…,N-1),求(m-1)次最小二乘擬合多項(xiàng)式

    (8)

    這一多項(xiàng)式函數(shù)組的生成步驟為:

    可以證明,由上述遞推構(gòu)造的多項(xiàng)式函數(shù)組(j=0,1,…,m-1)是相互正交的[9]。根據(jù)最小二乘原理,可得的表達(dá)式為:,j=0,1,…,m-1

    最后可以化成一般的m-1次多項(xiàng)式為:

    以上算法便于計(jì)算機(jī)實(shí)現(xiàn)。在實(shí)際計(jì)算過(guò)程中,為了防止運(yùn)算溢出,用代替,其中, i=0,1,..,N-1。

    擬合多項(xiàng)式如下:

    在數(shù)據(jù)處理的軟件模塊中,設(shè)置了交互式的數(shù)據(jù)處理模式,即對(duì)采集的試驗(yàn)數(shù)據(jù)進(jìn)行多項(xiàng)式回歸時(shí),根據(jù)用戶輸入要求進(jìn)行m次的回歸,軟件并將所得的各多項(xiàng)式要進(jìn)行顯著水平的檢驗(yàn)和正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差計(jì)算,直至達(dá)到試驗(yàn)要求為止。

    運(yùn)用此正交多項(xiàng)式回歸分析,對(duì)減壓閥出口流量階躍壓力響應(yīng)特性試驗(yàn)數(shù)據(jù)進(jìn)行處理,經(jīng)過(guò)6次回歸,得到其處理前后的曲線如圖5所示。

    4 CAT系統(tǒng)的誤差分析

    一個(gè)CAT系統(tǒng)的測(cè)試誤差來(lái)源是多方面的,實(shí)際情況很復(fù)雜。主要包括硬件、軟件、及數(shù)據(jù)處理幾個(gè)方面。對(duì)于傳感器、A/D轉(zhuǎn)換引起的誤差、上、下位機(jī)間數(shù)碼轉(zhuǎn)換產(chǎn)生的誤差,其相對(duì)誤差很小,在此不再討論。主要對(duì)應(yīng)用正交多項(xiàng)式回歸分析處理數(shù)據(jù)時(shí)產(chǎn)生的誤差進(jìn)行分析。

    以減壓閥的卸壓-建壓特性試驗(yàn)為例,具體分析正交多項(xiàng)式回歸的誤差分析。在卸壓-建壓特性曲線中,右邊的建壓特性曲線是用正交多項(xiàng)式最小二乘擬合的,曲線如圖6所示。圖6(a)為擬合前的曲線,圖6(b)為擬合后的曲線。

    在此試驗(yàn)中,右邊建壓曲線所采集的點(diǎn)數(shù)N=625,回歸次數(shù)為m-1=5次。則其[10]:

    總離差平方和為:,回歸差平方和為:jy

    第j個(gè)多項(xiàng)式的偏回歸平方和:,剩余平方和:

    其中為正交多項(xiàng)式回歸的系數(shù)(前已推導(dǎo));為:

    方差比:

    其中:的相應(yīng)自由度為,所以偏回歸方差估計(jì)值為:

    剩余平方和的相應(yīng)自由度為,所以剩余的方差估計(jì)值為:

    根據(jù)以上計(jì)算,可以得到以下正交多項(xiàng)式回歸的方差分析,如表1所示:

    經(jīng)查F檢驗(yàn)的臨界值(Fa)表,取a=0.01,f1=1,f2=(因?yàn)閒2=618>500)。

    在方差分析表1中,F(xiàn)計(jì)算值全都大于,故認(rèn)為都是高度顯著的,所以表中注以“**”,即表示1-5次項(xiàng)都是高度顯著,則所配制的多項(xiàng)式能夠正確地代表試驗(yàn)結(jié)果。其正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差為

    這樣的精度是比較滿意的。如果對(duì)數(shù)據(jù)處理的結(jié)果不滿意,還可以提高正交多項(xiàng)式的回歸次數(shù)來(lái)更加逼近測(cè)試數(shù)據(jù),獲得更高的精度。

    5 結(jié)語(yǔ)

    (1)在CAT系統(tǒng)中,數(shù)據(jù)測(cè)量的誤差來(lái)自多個(gè)方面。其中,由A/D轉(zhuǎn)換引起的誤差以及上、下位機(jī)間數(shù)據(jù)的數(shù)碼轉(zhuǎn)換引起的誤差極其微小,可以忽略。

    (2)在壓力控制閥的靜特性試驗(yàn)中,當(dāng)測(cè)試數(shù)據(jù)足夠致密,且每一測(cè)試點(diǎn)附近二階導(dǎo)數(shù)為常數(shù)時(shí),即幾何上反映為曲線的凹凸方向一致時(shí),對(duì)采集的數(shù)據(jù)采用非線性滑動(dòng)平滑方法,可以取得比較好的平滑效果。

    (3)在用正交多項(xiàng)式回歸試驗(yàn)曲線時(shí),回歸次數(shù)越高,回歸曲線越逼近真實(shí)曲線,正交多項(xiàng)式回歸的標(biāo)準(zhǔn)差就越小,數(shù)據(jù)處理的精度就越高,可以根據(jù)需要選擇合適的回歸次數(shù),以獲得所需的測(cè)量精度。

    參考文獻(xiàn)

    [1] Malrino A. Digital principles and applications[M]. New York: McGraw-Hill, 1981.

    [2] 段長(zhǎng)寶.液壓測(cè)試[M].國(guó)防工業(yè)出版社,1984.

    [3] 朱堅(jiān)民,周福章.溢流閥計(jì)算機(jī)輔助測(cè)試及數(shù)據(jù)處理[J].機(jī)床與液壓,1998(3):56-58.

    [4] Khalatova L.V., Accuracy Estimation for Automated Computer- measurement Process[J].Meas Tech,1985,1985,123(3).

    [5] L.R. Rabiner. Application of digital signal processing[M]. New Jersey: Prentice-Hall Inc.,1975.

    [6] R.F.Coughlin and F.F.Driscoll: Operational Amplifiers and Linear Integrated Circuits[M].Prentice Hall Inc.,New Jersey,1977.

    [7] D.E.Johnson and J.L.Hilburn: Rapid Practical Designs of Active Filters[M]. John Wiley & Sons, Inc., New York,1975.

    [8] 周江文,黃幼才.抗差最小二乘法[M].華中理工大學(xué)出版社,1997.

    [9] P.J. Torpey. Minicomputericing analog data collection[M].USA: McGraw-Hill, 1980.

    [10] 胡上序,陳德釗.觀測(cè)數(shù)據(jù)的分析與處理[M].浙江:浙江大學(xué)出版社,1996.endprint

    猜你喜歡
    誤差分析減壓閥數(shù)據(jù)處理
    氣動(dòng)舵機(jī)高壓反向直動(dòng)式減壓閥的設(shè)計(jì)及特性
    認(rèn)知診斷缺失數(shù)據(jù)處理方法的比較:零替換、多重插補(bǔ)與極大似然估計(jì)法*
    ILWT-EEMD數(shù)據(jù)處理的ELM滾動(dòng)軸承故障診斷
    氣壓高度計(jì)的測(cè)量誤差分析及修正方法
    基于“比值比較法”確定最佳實(shí)驗(yàn)方案的研究
    電學(xué)計(jì)量的誤差分析及不確定度理論探究分析
    車(chē)載雷達(dá)陣面舉升裝置的位置誤差
    基于希爾伯特- 黃變換的去噪法在外測(cè)數(shù)據(jù)處理中的應(yīng)用
    氣動(dòng)用電氣比例/伺服減壓閥的合理調(diào)壓方式探討
    氧氣呼吸器減壓閥的設(shè)計(jì)及試驗(yàn)研究
    河南科技(2015年7期)2015-03-11 16:23:18
    久久国产精品影院| 女人十人毛片免费观看3o分钟| 国产午夜福利久久久久久| 人妻久久中文字幕网| 香蕉av资源在线| 国内精品一区二区在线观看| 无人区码免费观看不卡| 国产成人福利小说| av中文乱码字幕在线| 精品久久久久久久人妻蜜臀av| aaaaa片日本免费| 国产精品三级大全| 欧美日韩黄片免| 国产精品不卡视频一区二区 | 校园春色视频在线观看| АⅤ资源中文在线天堂| 国产极品精品免费视频能看的| 内地一区二区视频在线| 夜夜爽天天搞| 亚洲欧美日韩无卡精品| 美女被艹到高潮喷水动态| 性色av乱码一区二区三区2| 婷婷色综合大香蕉| 国产免费男女视频| 欧美另类亚洲清纯唯美| 日本免费a在线| 日韩免费av在线播放| 国产高清三级在线| 高清毛片免费观看视频网站| 3wmmmm亚洲av在线观看| 国语自产精品视频在线第100页| 校园春色视频在线观看| 十八禁人妻一区二区| 亚洲国产欧洲综合997久久,| 精品国产三级普通话版| 国产免费一级a男人的天堂| 欧美日韩亚洲国产一区二区在线观看| 有码 亚洲区| 动漫黄色视频在线观看| 欧美区成人在线视频| 国产伦一二天堂av在线观看| 国产亚洲精品av在线| 亚洲成人中文字幕在线播放| 亚洲av成人精品一区久久| 日韩欧美精品免费久久 | 夜夜爽天天搞| 国产伦精品一区二区三区视频9| 伊人久久精品亚洲午夜| 精品日产1卡2卡| 毛片一级片免费看久久久久 | 国产综合懂色| 毛片女人毛片| 美女高潮喷水抽搐中文字幕| 久久久久久九九精品二区国产| 精品久久久久久久人妻蜜臀av| 日韩有码中文字幕| 亚洲国产欧洲综合997久久,| 日本五十路高清| 亚洲va日本ⅴa欧美va伊人久久| 亚洲久久久久久中文字幕| 国产精品久久视频播放| 少妇被粗大猛烈的视频| 亚洲专区国产一区二区| 激情在线观看视频在线高清| 欧美潮喷喷水| 很黄的视频免费| 亚洲精品影视一区二区三区av| 88av欧美| 亚洲成人精品中文字幕电影| av天堂中文字幕网| 免费在线观看日本一区| 美女 人体艺术 gogo| 欧美一区二区国产精品久久精品| 99国产精品一区二区三区| 如何舔出高潮| 欧美日韩福利视频一区二区| 国产高清有码在线观看视频| 精品一区二区免费观看| 久久精品91蜜桃| 精品一区二区三区av网在线观看| 五月玫瑰六月丁香| 老司机午夜福利在线观看视频| 嫩草影视91久久| 熟女人妻精品中文字幕| 国产三级在线视频| 99热这里只有是精品50| 欧美午夜高清在线| 可以在线观看毛片的网站| 一二三四社区在线视频社区8| 免费在线观看日本一区| 精品久久国产蜜桃| 免费人成在线观看视频色| 日韩欧美免费精品| 亚洲最大成人中文| 免费搜索国产男女视频| 午夜两性在线视频| 别揉我奶头~嗯~啊~动态视频| 少妇的逼好多水| 国产三级中文精品| 国产成+人综合+亚洲专区| 亚洲avbb在线观看| www.www免费av| 久久久久久久久久黄片| 2021天堂中文幕一二区在线观| 成年版毛片免费区| 午夜免费男女啪啪视频观看 | 亚洲中文字幕日韩| av福利片在线观看| 成人国产综合亚洲| 窝窝影院91人妻| 国产成人a区在线观看| 亚洲人与动物交配视频| 狠狠狠狠99中文字幕| 亚洲国产欧洲综合997久久,| 精华霜和精华液先用哪个| 在线a可以看的网站| 天堂影院成人在线观看| 两人在一起打扑克的视频| 观看美女的网站| 少妇被粗大猛烈的视频| 久久香蕉精品热| 国产色爽女视频免费观看| 国产淫片久久久久久久久 | 男女视频在线观看网站免费| 哪里可以看免费的av片| 美女大奶头视频| 欧美激情国产日韩精品一区| 一本精品99久久精品77| 欧美色视频一区免费| 十八禁网站免费在线| 亚洲精品日韩av片在线观看| 91字幕亚洲| h日本视频在线播放| 国产精品野战在线观看| 两人在一起打扑克的视频| 毛片一级片免费看久久久久 | 99久久99久久久精品蜜桃| 蜜桃久久精品国产亚洲av| 久99久视频精品免费| 乱人视频在线观看| 久久婷婷人人爽人人干人人爱| 我的老师免费观看完整版| netflix在线观看网站| 国产在视频线在精品| 成年免费大片在线观看| 一本一本综合久久| 欧美一区二区亚洲| 亚洲熟妇中文字幕五十中出| 丁香欧美五月| 精品一区二区三区人妻视频| 18禁裸乳无遮挡免费网站照片| 男人舔女人下体高潮全视频| 波多野结衣高清无吗| 网址你懂的国产日韩在线| 成人一区二区视频在线观看| 久久久久精品国产欧美久久久| 在线播放国产精品三级| 欧美中文日本在线观看视频| 给我免费播放毛片高清在线观看| 日日干狠狠操夜夜爽| 色哟哟·www| 国产精品,欧美在线| 精品国产三级普通话版| 国产乱人视频| 夜夜看夜夜爽夜夜摸| 又黄又爽又免费观看的视频| 欧美日本视频| 18+在线观看网站| 男女那种视频在线观看| 校园春色视频在线观看| 久久精品国产99精品国产亚洲性色| 无人区码免费观看不卡| 别揉我奶头~嗯~啊~动态视频| 毛片女人毛片| 男人的好看免费观看在线视频| 色播亚洲综合网| 免费观看人在逋| 欧美一区二区国产精品久久精品| 国产高清视频在线观看网站| 国产一区二区在线观看日韩| 亚洲精品一区av在线观看| 国产三级中文精品| 国内少妇人妻偷人精品xxx网站| 91午夜精品亚洲一区二区三区 | 又爽又黄无遮挡网站| 国产麻豆成人av免费视频| 亚洲国产色片| 午夜老司机福利剧场| 日韩av在线大香蕉| 免费大片18禁| 国产一区二区三区在线臀色熟女| 天堂动漫精品| 一二三四社区在线视频社区8| 久久久久久久久久黄片| 黄色女人牲交| 国产一区二区三区在线臀色熟女| 免费看美女性在线毛片视频| 亚洲国产欧美人成| 色综合亚洲欧美另类图片| 一区福利在线观看| 人人妻人人澡欧美一区二区| ponron亚洲| 91久久精品电影网| 国内精品美女久久久久久| 日本黄色视频三级网站网址| 悠悠久久av| 国产三级中文精品| 在线观看av片永久免费下载| 三级国产精品欧美在线观看| 两个人的视频大全免费| 女同久久另类99精品国产91| 看免费av毛片| 一级作爱视频免费观看| 黄片小视频在线播放| 伦理电影大哥的女人| 精品无人区乱码1区二区| 波多野结衣高清作品| 国产视频一区二区在线看| 久久国产乱子伦精品免费另类| 成熟少妇高潮喷水视频| 免费观看人在逋| 国产一区二区三区视频了| 黄色丝袜av网址大全| 欧美成人免费av一区二区三区| 午夜福利成人在线免费观看| 精品日产1卡2卡| 亚洲人成网站在线播| 久久伊人香网站| 国产又黄又爽又无遮挡在线| 精品99又大又爽又粗少妇毛片 | 69av精品久久久久久| 男女视频在线观看网站免费| 两个人视频免费观看高清| 淫秽高清视频在线观看| 又紧又爽又黄一区二区| 好男人在线观看高清免费视频| 午夜两性在线视频| www.www免费av| 长腿黑丝高跟| 国产伦在线观看视频一区| 免费人成在线观看视频色| 国产在线男女| 看黄色毛片网站| 天天一区二区日本电影三级| 夜夜看夜夜爽夜夜摸| 永久网站在线| ponron亚洲| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av香蕉五月| 丰满的人妻完整版| 国产欧美日韩精品亚洲av| 九色成人免费人妻av| 女生性感内裤真人,穿戴方法视频| 国产一区二区亚洲精品在线观看| 国产91精品成人一区二区三区| www.色视频.com| av欧美777| 国内精品美女久久久久久| 中文字幕av在线有码专区| 女人被狂操c到高潮| 男女视频在线观看网站免费| 床上黄色一级片| 日韩中字成人| 精品一区二区三区视频在线观看免费| 人妻丰满熟妇av一区二区三区| 少妇人妻精品综合一区二区 | 午夜福利免费观看在线| 精品久久久久久,| 三级男女做爰猛烈吃奶摸视频| 成人欧美大片| 亚洲avbb在线观看| av女优亚洲男人天堂| 18禁在线播放成人免费| 一进一出好大好爽视频| 成人国产一区最新在线观看| 精品国产亚洲在线| 久久欧美精品欧美久久欧美| 欧美bdsm另类| 国产av不卡久久| 免费观看人在逋| 午夜福利视频1000在线观看| 99热这里只有是精品在线观看 | 国产私拍福利视频在线观看| 亚洲国产欧美人成| 天美传媒精品一区二区| 久久天躁狠狠躁夜夜2o2o| 成人av在线播放网站| 欧美黄色淫秽网站| 中文字幕人妻熟人妻熟丝袜美| 国产欧美日韩精品一区二区| 欧美激情久久久久久爽电影| 欧美日韩黄片免| 无遮挡黄片免费观看| 久久精品国产自在天天线| 变态另类丝袜制服| 亚洲最大成人av| 成年人黄色毛片网站| 国产精品久久电影中文字幕| 怎么达到女性高潮| 在线观看舔阴道视频| 精品福利观看| 日韩免费av在线播放| 五月伊人婷婷丁香| 波多野结衣高清作品| 亚洲av美国av| 国产精品精品国产色婷婷| 午夜福利在线观看吧| 一进一出抽搐动态| 国产v大片淫在线免费观看| 午夜精品一区二区三区免费看| 小蜜桃在线观看免费完整版高清| 欧美最新免费一区二区三区 | 国产三级在线视频| 亚洲自偷自拍三级| 亚洲国产精品成人综合色| 久久性视频一级片| 又爽又黄无遮挡网站| 脱女人内裤的视频| 好看av亚洲va欧美ⅴa在| 在线天堂最新版资源| 免费看a级黄色片| 欧美黄色片欧美黄色片| 一夜夜www| 成年人黄色毛片网站| 校园春色视频在线观看| 搡老熟女国产l中国老女人| 不卡一级毛片| 久久久成人免费电影| 狠狠狠狠99中文字幕| 狂野欧美白嫩少妇大欣赏| 日本一本二区三区精品| 国产精品国产高清国产av| 国产在线精品亚洲第一网站| 久久久精品大字幕| 免费无遮挡裸体视频| 无遮挡黄片免费观看| 成人特级av手机在线观看| 久久精品综合一区二区三区| 日日干狠狠操夜夜爽| 老司机午夜十八禁免费视频| 级片在线观看| 激情在线观看视频在线高清| 亚洲乱码一区二区免费版| 免费高清视频大片| 啪啪无遮挡十八禁网站| 日本黄色片子视频| 波多野结衣巨乳人妻| 国内精品美女久久久久久| 精品国产三级普通话版| 五月玫瑰六月丁香| 欧美性猛交黑人性爽| 高清毛片免费观看视频网站| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点| 日韩精品青青久久久久久| 免费高清视频大片| 一本久久中文字幕| 国产久久久一区二区三区| 熟女电影av网| 色噜噜av男人的天堂激情| 国产淫片久久久久久久久 | 真人做人爱边吃奶动态| 麻豆成人av在线观看| 免费av毛片视频| 欧美成人a在线观看| 色尼玛亚洲综合影院| 婷婷亚洲欧美| 少妇丰满av| 久久亚洲真实| 国产精品永久免费网站| 欧美最新免费一区二区三区 | 国产一区二区在线av高清观看| 99在线视频只有这里精品首页| 国产色婷婷99| 12—13女人毛片做爰片一| 亚洲自偷自拍三级| 亚洲欧美日韩高清专用| 欧美最黄视频在线播放免费| 91久久精品电影网| 99在线人妻在线中文字幕| 国产精品嫩草影院av在线观看 | 欧美黑人巨大hd| 亚洲av不卡在线观看| 人妻制服诱惑在线中文字幕| 日韩欧美精品v在线| 欧美成人性av电影在线观看| 久久精品国产99精品国产亚洲性色| 久久久久久久久久成人| www日本黄色视频网| or卡值多少钱| 成人亚洲精品av一区二区| 色5月婷婷丁香| 国产一区二区在线观看日韩| 精品久久久久久久末码| 级片在线观看| 国产乱人伦免费视频| 国产一区二区在线观看日韩| 少妇裸体淫交视频免费看高清| 国内精品久久久久久久电影| 国产亚洲精品av在线| 国产男靠女视频免费网站| 免费观看人在逋| 中出人妻视频一区二区| 国产探花极品一区二区| eeuss影院久久| 欧美黑人巨大hd| 精品久久久久久久久久久久久| 亚洲成人免费电影在线观看| 男女床上黄色一级片免费看| 能在线免费观看的黄片| 中文字幕av在线有码专区| 成人av一区二区三区在线看| 99在线人妻在线中文字幕| 能在线免费观看的黄片| 少妇的逼水好多| 午夜影院日韩av| 亚洲aⅴ乱码一区二区在线播放| 国产精品伦人一区二区| 欧美+日韩+精品| 少妇被粗大猛烈的视频| 日本一二三区视频观看| 99热这里只有是精品50| 亚洲成人久久性| 综合色av麻豆| 亚洲精品456在线播放app | 久久久精品大字幕| 全区人妻精品视频| 国产av不卡久久| 老司机福利观看| 熟女电影av网| 一二三四社区在线视频社区8| 精品人妻熟女av久视频| 国产精品不卡视频一区二区 | av女优亚洲男人天堂| 欧美潮喷喷水| 精品人妻1区二区| av视频在线观看入口| 又爽又黄a免费视频| 免费av不卡在线播放| 亚洲av成人av| 精品无人区乱码1区二区| 国产免费男女视频| 久久99热这里只有精品18| 午夜免费激情av| 人妻丰满熟妇av一区二区三区| 动漫黄色视频在线观看| a级毛片免费高清观看在线播放| 国产精品亚洲av一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 国产蜜桃级精品一区二区三区| 俺也久久电影网| 麻豆一二三区av精品| 久久久久亚洲av毛片大全| 69人妻影院| 午夜久久久久精精品| 午夜福利高清视频| 18禁在线播放成人免费| 国产精品久久电影中文字幕| 国产极品精品免费视频能看的| 九九热线精品视视频播放| 99久国产av精品| 久久精品91蜜桃| 偷拍熟女少妇极品色| 久久伊人香网站| 黄色女人牲交| 毛片女人毛片| 看黄色毛片网站| 男女下面进入的视频免费午夜| 亚洲五月婷婷丁香| 在线播放无遮挡| 国产欧美日韩一区二区三| 亚洲国产色片| 99国产精品一区二区三区| 亚洲自拍偷在线| 在线观看美女被高潮喷水网站 | 国产伦在线观看视频一区| 香蕉av资源在线| 在线免费观看不下载黄p国产 | 国产一区二区三区在线臀色熟女| 我要看日韩黄色一级片| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美 国产精品| 亚洲国产精品成人综合色| 欧美最黄视频在线播放免费| 精品久久久久久久久亚洲 | 狠狠狠狠99中文字幕| 国产高清三级在线| av在线观看视频网站免费| 免费看美女性在线毛片视频| x7x7x7水蜜桃| 国产大屁股一区二区在线视频| 欧美日韩国产亚洲二区| 99久久久亚洲精品蜜臀av| 久久国产乱子伦精品免费另类| 美女高潮喷水抽搐中文字幕| 欧美成人性av电影在线观看| 国产综合懂色| 久久欧美精品欧美久久欧美| 99riav亚洲国产免费| 午夜免费男女啪啪视频观看 | 成人鲁丝片一二三区免费| 国产精品一区二区性色av| а√天堂www在线а√下载| 国产探花在线观看一区二区| 一区福利在线观看| 久久久久性生活片| 国语自产精品视频在线第100页| 午夜福利免费观看在线| 啦啦啦观看免费观看视频高清| 久久精品夜夜夜夜夜久久蜜豆| 国产成人av教育| 一级黄色大片毛片| 18禁在线播放成人免费| 精品久久久久久久久久免费视频| 村上凉子中文字幕在线| 国产精品嫩草影院av在线观看 | 免费在线观看亚洲国产| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区四区激情视频 | 亚洲在线观看片| 亚洲自拍偷在线| 亚洲成人久久爱视频| 又紧又爽又黄一区二区| 最近中文字幕高清免费大全6 | 毛片女人毛片| 亚洲国产精品合色在线| 国产一区二区三区在线臀色熟女| 最近最新中文字幕大全电影3| 色噜噜av男人的天堂激情| 很黄的视频免费| 亚洲成人久久性| 俺也久久电影网| 亚洲av二区三区四区| 嫩草影院新地址| 国产欧美日韩精品一区二区| 欧美黄色片欧美黄色片| 免费黄网站久久成人精品 | 欧美乱妇无乱码| .国产精品久久| 日韩av在线大香蕉| 中文资源天堂在线| 色哟哟·www| 夜夜躁狠狠躁天天躁| 男女做爰动态图高潮gif福利片| 久久久久性生活片| 久久久久久九九精品二区国产| a在线观看视频网站| 国产精品国产高清国产av| 免费看美女性在线毛片视频| 国产伦精品一区二区三区四那| 国产高清视频在线观看网站| 亚洲不卡免费看| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 97碰自拍视频| 午夜精品久久久久久毛片777| 亚洲成人免费电影在线观看| 一个人观看的视频www高清免费观看| 特级一级黄色大片| 99国产综合亚洲精品| 99久久99久久久精品蜜桃| 精品久久国产蜜桃| 国产探花在线观看一区二区| 深夜精品福利| 亚洲第一电影网av| 亚洲片人在线观看| 色av中文字幕| 少妇被粗大猛烈的视频| 老司机午夜福利在线观看视频| 亚洲精品一区av在线观看| 亚洲人与动物交配视频| 91久久精品电影网| 精品一区二区三区视频在线| 精品一区二区三区视频在线观看免费| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| а√天堂www在线а√下载| 亚洲va日本ⅴa欧美va伊人久久| 观看美女的网站| 亚洲一区二区三区不卡视频| 欧美性猛交╳xxx乱大交人| 国产精品人妻久久久久久| 久久精品综合一区二区三区| 99精品在免费线老司机午夜| 少妇的逼好多水| 免费在线观看日本一区| 免费高清视频大片| 怎么达到女性高潮| 88av欧美| 成人美女网站在线观看视频| 天堂动漫精品| 午夜精品一区二区三区免费看| www.999成人在线观看| 久久人妻av系列| 国产成年人精品一区二区| 看黄色毛片网站| 亚洲国产欧美人成| 99久久精品一区二区三区| 国产精品三级大全| 波多野结衣巨乳人妻| 露出奶头的视频| 国产精品电影一区二区三区| 亚洲片人在线观看| 色视频www国产| 少妇人妻一区二区三区视频| 色视频www国产| 亚洲精品成人久久久久久| 婷婷丁香在线五月| ponron亚洲| 欧美xxxx黑人xx丫x性爽| 亚洲专区国产一区二区| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看日本一区| 久久久久久久久大av| 亚洲人成网站高清观看|