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

    含過程變量的K模型最優(yōu)正交區(qū)組設(shè)計(jì)

    2021-07-02 01:02:48朱小淵張崇岐
    關(guān)鍵詞:單純形區(qū)組混料

    龔 政,朱小淵,張崇岐

    (廣州大學(xué) 經(jīng)濟(jì)與統(tǒng)計(jì)學(xué)院,廣東廣州 510006)

    §1 引言

    隨著傳統(tǒng)醫(yī)學(xué)系統(tǒng)化與現(xiàn)代化進(jìn)程的推進(jìn),混料試驗(yàn)設(shè)計(jì)在臨床治療中的應(yīng)用越發(fā)普及.中醫(yī)理論來自于經(jīng)驗(yàn)醫(yī)學(xué),在進(jìn)行中藥試驗(yàn)研究時(shí),為減少受試人數(shù)是試驗(yàn)設(shè)計(jì)的主要研究目標(biāo),混料最優(yōu)設(shè)計(jì)能保證試驗(yàn)結(jié)果的同時(shí)顯著達(dá)到減小樣本量.近年來利用混料最優(yōu)設(shè)計(jì)的中醫(yī)研究越發(fā)普遍,部分研究者利用混料最優(yōu)設(shè)計(jì)優(yōu)化處方,例如多糖速溶片的處方[1]的優(yōu)化.傳統(tǒng)醫(yī)學(xué)需要進(jìn)步,應(yīng)當(dāng)接受循證醫(yī)學(xué)的檢驗(yàn),合理利用混料最優(yōu)設(shè)計(jì)在內(nèi)的現(xiàn)代統(tǒng)計(jì)手段來檢驗(yàn)其有效性.除此之外,混料試驗(yàn)設(shè)計(jì)在工業(yè)生產(chǎn)以及現(xiàn)代生物技術(shù)中,也有著大規(guī)模的應(yīng)用.

    混料試驗(yàn)設(shè)計(jì)主要研究各種混料成份與試驗(yàn)指標(biāo)之間的關(guān)系,試驗(yàn)指標(biāo)只與每種成份的比例有關(guān),而與總量無關(guān).在常用的最優(yōu)準(zhǔn)則中,D-最優(yōu)設(shè)計(jì)準(zhǔn)則以最大化信息矩陣行列式為目標(biāo),A-最優(yōu)準(zhǔn)則是需要滿足最小化試驗(yàn)域上設(shè)計(jì)的最大預(yù)測方差,而I-最優(yōu)準(zhǔn)則需要滿足試驗(yàn)域上設(shè)計(jì)的最小化平均預(yù)測方差,各類最優(yōu)準(zhǔn)則均有不同的優(yōu)點(diǎn),人們基于實(shí)際的需求,選擇最適合的模型與最優(yōu)準(zhǔn)則[2].本文主要研究滿足D-最優(yōu)準(zhǔn)則時(shí)的最優(yōu)設(shè)計(jì).

    在實(shí)際試驗(yàn)中,先有具體問題,再有模型,從規(guī)范多項(xiàng)式到中心多項(xiàng)式再到特殊模型,由于試驗(yàn)點(diǎn)隨著混料成分的增加而變多,因此需要適當(dāng)簡化模型,使得模型參數(shù)個(gè)數(shù)正好等于試驗(yàn)點(diǎn)數(shù),以達(dá)到飽和設(shè)計(jì).現(xiàn)今關(guān)于一般混料模型的最優(yōu)設(shè)計(jì)的研究已經(jīng)日趨完善,人們開始研究更多復(fù)雜模型的最優(yōu)設(shè)計(jì)[3].

    記混料模型中的q個(gè)分量為xi,i=1,2,...,q,其試驗(yàn)域?yàn)?/p>

    它定義了q分量混料試驗(yàn)所在的q-1維單純形試驗(yàn)區(qū)域,并對(duì)擬合的模型產(chǎn)生重大的影響.

    在混料模型中,變量是指各成份占混料總量的比例,一般混料模型可表示為

    其中y(x)是點(diǎn)x響應(yīng)的觀察值,f(x)=(f1(x),f2(x),...,fp(x))T是給定的關(guān)于x=(x1,x2,……,xq)的點(diǎn)函數(shù)向量,ε是隨機(jī)誤差項(xiàng),服從正態(tài)分布N(0,σ2),而β=(β1,β2,……,βp)T是待估參數(shù)向量,模型(2)常用η(x)=fT(x)β,其中η(x)表示響應(yīng)的均值.

    在混料試驗(yàn)中,影響試驗(yàn)指標(biāo)的變量除了混合物的比例外,還存在與成份的組成無關(guān)的變量,該變量統(tǒng)稱為過程變量,例如中藥藥材研磨的精細(xì)程度,藥品煎熬時(shí)的火候溫度,以及病人服用藥物的間隔時(shí)間等.一般的含過程變量混料模型可表示為η(x)=fT(x)β+αZ,這里的Z是過程變量,α為過程變量的待估參數(shù).

    關(guān)于混料模型的區(qū)組研究現(xiàn)如今的成果已有許多,為了保證各區(qū)組內(nèi)混合物不受區(qū)組效應(yīng)的影響,Nigam[4]首先提出需以混合物成份為自變量,從研究重心點(diǎn)測度轉(zhuǎn)向研究等測度的試驗(yàn)點(diǎn)的位置,Saxena[5]以此提出了對(duì)稱單純形區(qū)組設(shè)計(jì),John[6]對(duì)Nigam提出的混料正交條件進(jìn)行了簡化并給出了含過程變量的二階規(guī)范多項(xiàng)式模型的擬合,但是這樣設(shè)計(jì)雖然規(guī)避了區(qū)組效應(yīng)對(duì)試驗(yàn)的影響,卻使得試驗(yàn)點(diǎn)變多.

    因此Czitrom考慮以正交表來構(gòu)建混料正交設(shè)計(jì),對(duì)于含過程變量的3分量二階規(guī)范多項(xiàng)式模型,Czitrom[7]給出了正交區(qū)組設(shè)計(jì)并求出了滿足D-最優(yōu)準(zhǔn)則時(shí)的設(shè)計(jì)點(diǎn),隨后一些學(xué)者完善并拓展了正交區(qū)組設(shè)計(jì)的理論基礎(chǔ),成果突出的有Murthy[8]考慮了在混料限制下的區(qū)組設(shè)計(jì),Lewis[9]將規(guī)范多項(xiàng)式模型的正交區(qū)組設(shè)計(jì)推廣到了q分量的情形,Prescott[10]則是提出了各區(qū)組中設(shè)計(jì)點(diǎn)的利益區(qū)域應(yīng)是正規(guī)單純形在其內(nèi)部的投影,解決了方差函數(shù)驗(yàn)證的難題,現(xiàn)今已有研究者為提高模型參數(shù)估計(jì)精度,將研究擴(kuò)展到了缺失項(xiàng)模型的I-最優(yōu)設(shè)計(jì)[11].

    對(duì)于混料中單成份效應(yīng)強(qiáng)于交互效應(yīng)的情況,Draper等[12]基于向量與矩陣的Kronecker代數(shù),提出了K模型,它是一類齊次模型,它是用混料分量的基本約束條件,將一次項(xiàng)轉(zhuǎn)化為二次項(xiàng),統(tǒng)一去研究,無需分為線性和二次項(xiàng)分別考慮.

    關(guān)于K模型的研究十分廣泛,Manisha和Pal[13-14]將K模型轉(zhuǎn)化為矩陣形式研究,李俊鵬等[15]研究了混料K模型的R-最優(yōu)設(shè)計(jì).

    其中γi與γij是待估參數(shù),xi(i=1,2,……,q)為混料成份的比例.

    含過程變量Z的q分量K混料模型,表述如下

    §2 主要結(jié)果

    本文對(duì)含過程變量的K模型展開研究,并求出其q分量下滿足D-最優(yōu)準(zhǔn)則的正交區(qū)組設(shè)計(jì).

    混料模型的試驗(yàn)設(shè)計(jì)ξ是試驗(yàn)域(1)上的一個(gè)概率分布

    q分量混料模型的試驗(yàn)區(qū)域是一個(gè)q維空間上的q-1維正規(guī)單純形,xi(i=1,2,……,n) 為q-1維正規(guī)單純形上的任一試驗(yàn)點(diǎn),ωi表示各試驗(yàn)點(diǎn)的測度,且ωi滿足.

    得到設(shè)計(jì)后,可以得到設(shè)計(jì)矩陣

    由此可得到包含一個(gè)設(shè)計(jì)的所有信息的信息矩陣

    其中Λ為測度矩陣.

    同時(shí)可以得到設(shè)計(jì)矩陣對(duì)應(yīng)的方差函數(shù)

    在試驗(yàn)設(shè)計(jì)中,經(jīng)常會(huì)遇見一個(gè)可能對(duì)響應(yīng)產(chǎn)生影響的因子,對(duì)這個(gè)因子的效應(yīng)并不感興趣,將之稱為討厭因子,它的特性為未知的和不可控,具有隨機(jī)性.可以通過區(qū)組化設(shè)計(jì)來系統(tǒng)地消除它在處理統(tǒng)計(jì)比較中的效應(yīng),如果有a個(gè)待比較的處理和b個(gè)區(qū)組,它的區(qū)組設(shè)計(jì)如下.

    表1 a個(gè)處理和b 個(gè)區(qū)組的設(shè)計(jì)

    它的效應(yīng)模型為

    其中yij是在第j個(gè)區(qū)組中第i個(gè)處理的觀測數(shù)據(jù),μ是總均值,λi是第i個(gè)處理的效應(yīng),αj是第j個(gè)區(qū)組的效應(yīng),εij是通常的隨機(jī)誤差項(xiàng).

    為了系統(tǒng)性控制過程變量,可以使用正交拉丁方設(shè)計(jì)讓區(qū)組保持正交.由此,在處理含過程變量的混料模型時(shí),可將過程變量的效應(yīng)看作區(qū)組效應(yīng),將每個(gè)實(shí)驗(yàn)單元分為若干正交區(qū)組.

    當(dāng)q為素?cái)?shù)或者是素?cái)?shù)的方冪時(shí),必然存在q-1個(gè)正交的拉丁方[16].從常用的3分量混料模型開始討論,K模型的3分量形式為

    此時(shí)存在2個(gè)正交區(qū)組.若在混料試驗(yàn)過程中以zi(i=1,2,……,n)來表示n個(gè)過程變量,每個(gè)過程變量均為兩水平,其中的高水平取值為+1,低水平取值為-1,將之分為兩個(gè)區(qū)組,區(qū)組1滿足條件z1=z2=……·=zn=-1,區(qū)組2滿足z1=z2=……·=zn=+1[17].為使區(qū)組間具有正交性,將混合物的組成成分作為自變量進(jìn)行處理.

    由此可知K模型(4)在3分量情況下,作出如下二正交區(qū)組設(shè)計(jì).

    其設(shè)計(jì)矩陣為

    其信息矩陣表示為

    其中I3是3階單位陣,J3是各元素皆為1的3階方陣,并且

    當(dāng)滿足D-最優(yōu)準(zhǔn)則時(shí),在滿足正交性的各區(qū)組內(nèi),依舊保留了設(shè)計(jì)點(diǎn)在邊界處與中心處的特性,q分量K模型滿足D-最優(yōu)準(zhǔn)則的設(shè)計(jì)點(diǎn)必然在中心點(diǎn)與邊界處,此時(shí)信息矩陣行列式為

    其最大值取在邊界處,令c=0,b=1-a.經(jīng)計(jì)算,滿足D-最優(yōu)準(zhǔn)則時(shí)a=0.1685或0.8315,如圖1所示,可發(fā)現(xiàn)為雙峰圖形,每條邊界處各有兩個(gè)頂點(diǎn),共計(jì)6個(gè)邊界頂點(diǎn),再加上2個(gè)中心點(diǎn),共同構(gòu)成了試驗(yàn)的設(shè)計(jì)點(diǎn).

    圖1 信息矩陣行列式最大值

    含過程變量的3分量K模型的D-最優(yōu)正交區(qū)組設(shè)計(jì)如下所示

    為了驗(yàn)證其D-最優(yōu)設(shè)計(jì),需求解設(shè)計(jì)點(diǎn)以及各類重心點(diǎn)的方差函數(shù)值,三分量K模型方差函數(shù)公式為

    引理2.1[18]設(shè)f(x)為列向量,G為正定陣,則二次型fT(x)Gf(x)在正規(guī)單純形利益區(qū)域上的最大值必在其各類重心處到達(dá).

    各區(qū)組的利益區(qū)域?yàn)檎?guī)單純形在其內(nèi)部的投影,因此計(jì)算各重心點(diǎn),即

    各重心處方差函數(shù)為d(x1,ξ)=1,d(x2,ξ)=0.8864,d(x3,ξ)=0.5.

    由上所述,可以看到所求區(qū)域內(nèi)方差函數(shù)值小于等于1,滿足D-最優(yōu)準(zhǔn)則.同時(shí)可以觀察到方差最大值在頂點(diǎn)取得,得出定理2.1.

    定理2.1含過程變量的3分量K混料模型,其滿足D-最優(yōu)準(zhǔn)則的正交區(qū)組設(shè)計(jì)為

    其中P3表示單純形上滿足D-最優(yōu)準(zhǔn)則的單純形邊界上的6個(gè)點(diǎn).

    當(dāng)q=3時(shí),混料模型在2維單純形上,因此只會(huì)取到邊界線上的點(diǎn),當(dāng)q ≥4 時(shí),其所在的單純形區(qū)域是一個(gè)高維的凸多面體,而正交區(qū)組的利益區(qū)域?yàn)閱渭冃蝺?nèi)部的一個(gè)投影,其頂點(diǎn)依舊取在邊界線上.

    推廣到分量q為素?cái)?shù)或素?cái)?shù)冪方時(shí),共有q(q-1)個(gè)排列,可將試驗(yàn)分為q-1個(gè)正交區(qū)組,含過程變量的q分量K模型的信息矩陣為

    其中Iq是q階單位陣,Jq是各元素皆為1的q階方陣,并且

    根據(jù)D-最優(yōu)準(zhǔn)則,要求滿足max|XTX|,在此令x1=a,x2=1-a,xi≥3=0,此時(shí)(11)為

    其信息矩陣行列式經(jīng)計(jì)算為

    由于正交區(qū)組設(shè)計(jì)依賴于正交表的建立,需確定q分量,求(13)極大值來解a.對(duì)于q分量情況下的一個(gè)D-最優(yōu)設(shè)計(jì)點(diǎn)為(0,0,……,a,0,……,1-a,……,0,0).由此可以得出任意q為確定值時(shí)的設(shè)計(jì)點(diǎn)與正交區(qū)組.根據(jù)上述分析證明,可得出定理2.2.

    定理2.2含過程變量的q分量二階K模型,其滿足D-最優(yōu)正交區(qū)組設(shè)計(jì)共有q2-q個(gè)邊界頂點(diǎn),每個(gè)區(qū)組中的設(shè)計(jì)點(diǎn)由q個(gè)邊界點(diǎn)與一個(gè)中心點(diǎn)構(gòu)成,q分量時(shí)滿足D-最優(yōu)準(zhǔn)則的正交區(qū)組設(shè)計(jì)為

    其中Pq為滿足D-最優(yōu)準(zhǔn)則的單純形邊界上的點(diǎn).

    §3 結(jié)論

    本文對(duì)K模型構(gòu)建正交區(qū)組,并論證q種混合物時(shí)滿足D-最優(yōu)準(zhǔn)則的正交區(qū)組設(shè)計(jì).D-最優(yōu)準(zhǔn)則是應(yīng)用最廣泛的最優(yōu)準(zhǔn)則,且可以通過KWT等價(jià)定理來驗(yàn)證本設(shè)計(jì)同樣是G-最優(yōu)的,相較于規(guī)范多項(xiàng)式模型的擬合,精度有了進(jìn)一步提升.除此之外,還可以用此方法研究多項(xiàng)式模型,倒數(shù)模型以及齊次模型等滿足混料約束條件的模型.

    混料最優(yōu)設(shè)計(jì)在中醫(yī)領(lǐng)域的合理應(yīng)用,解決了以往試驗(yàn)設(shè)計(jì)中樣本量需求過大的不足,研究者可以根據(jù)需求來選擇各類最優(yōu)準(zhǔn)則,以提高模型參數(shù)估計(jì)的精度,但在中藥制造中,由于工藝精度,病人個(gè)體差異以及輔藥的存在,一般混料模型在擬合時(shí)存在一定的誤差,每一次建模時(shí)都需謹(jǐn)慎對(duì)待.

    隨著混料試驗(yàn)設(shè)計(jì)的推廣,更多的混料模型也會(huì)在中醫(yī)領(lǐng)域有著更多的應(yīng)用,并且針對(duì)更復(fù)雜的臨床試驗(yàn),研究領(lǐng)域?qū)?huì)逐步轉(zhuǎn)移到含特殊項(xiàng)的混料模型中,研究人員在進(jìn)行混料試驗(yàn)最優(yōu)設(shè)計(jì)時(shí)也會(huì)為提高參數(shù)估計(jì)效率,而從D-最優(yōu)準(zhǔn)則的理論推導(dǎo)發(fā)展推廣到E-最優(yōu)準(zhǔn)則與I-最優(yōu)準(zhǔn)則等滿足研究者需求的方向,可以預(yù)見混料模型的正交區(qū)組設(shè)計(jì)會(huì)在將來不斷推陳出新,有著長足的發(fā)展空間.

    猜你喜歡
    單純形區(qū)組混料
    變化區(qū)組隨機(jī)化及其SAS宏實(shí)現(xiàn)
    雙重稀疏約束優(yōu)化問題的一種貪婪單純形算法
    如何正確運(yùn)用方差分析
    ——平衡不完全區(qū)組設(shè)計(jì)定量資料一元方差分析
    一種具有混料功能的橡膠射出成型機(jī)
    中醫(yī)臨床研究中區(qū)組設(shè)計(jì)應(yīng)用現(xiàn)狀的計(jì)量學(xué)分析*
    基于改進(jìn)單純形算法的Topmodel參數(shù)優(yōu)化研究
    基于PLC的混料罐控制系統(tǒng)設(shè)計(jì)
    電子制作(2016年21期)2016-05-17 03:52:46
    聚丙烯/CaCO3共混料中CaCO3含量測定方法對(duì)比
    一類混料指數(shù)模型的D-最優(yōu)設(shè)計(jì)
    基于數(shù)據(jù)融合與單純形遺傳算法的管道損傷識(shí)別
    久久久国产精品麻豆| 亚洲五月婷婷丁香| 热99re8久久精品国产| 中国国产av一级| 在线 av 中文字幕| 国产一区有黄有色的免费视频| 亚洲久久久国产精品| 亚洲一区二区三区欧美精品| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁日日躁夜夜躁夜夜| 中文字幕制服av| 国产不卡av网站在线观看| 91av网站免费观看| 777米奇影视久久| 美国免费a级毛片| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 国产精品偷伦视频观看了| 男女高潮啪啪啪动态图| 欧美成人午夜精品| 伊人久久大香线蕉亚洲五| 国产日韩欧美在线精品| 午夜免费鲁丝| 欧美人与性动交α欧美软件| 每晚都被弄得嗷嗷叫到高潮| 日本欧美视频一区| 国产精品二区激情视频| 欧美精品一区二区大全| 精品国产一区二区久久| 欧美中文综合在线视频| 精品国内亚洲2022精品成人 | 久久精品aⅴ一区二区三区四区| 在线观看免费视频网站a站| 国产欧美亚洲国产| 老鸭窝网址在线观看| 18禁黄网站禁片午夜丰满| 青春草亚洲视频在线观看| 在线十欧美十亚洲十日本专区| 狠狠精品人妻久久久久久综合| 无遮挡黄片免费观看| 亚洲欧美色中文字幕在线| 丝袜喷水一区| 深夜精品福利| 51午夜福利影视在线观看| 免费看十八禁软件| 欧美精品一区二区免费开放| 午夜精品久久久久久毛片777| 国产男女内射视频| av在线播放精品| 美国免费a级毛片| 97人妻天天添夜夜摸| 动漫黄色视频在线观看| 亚洲全国av大片| 青春草视频在线免费观看| 另类亚洲欧美激情| 日韩欧美一区二区三区在线观看 | 国产黄色免费在线视频| 久久久久久久久久久久大奶| a级毛片在线看网站| 18在线观看网站| 国产成人影院久久av| 久久性视频一级片| 午夜精品久久久久久毛片777| 99精品久久久久人妻精品| 99精品久久久久人妻精品| 精品国产乱码久久久久久小说| 成人手机av| 一级片免费观看大全| 丝袜人妻中文字幕| 嫁个100分男人电影在线观看| 日韩中文字幕视频在线看片| 午夜日韩欧美国产| 在线观看www视频免费| 午夜影院在线不卡| 久久精品人人爽人人爽视色| 麻豆av在线久日| e午夜精品久久久久久久| 国产男女内射视频| 日日摸夜夜添夜夜添小说| 又黄又粗又硬又大视频| 国产精品免费视频内射| 每晚都被弄得嗷嗷叫到高潮| 欧美变态另类bdsm刘玥| 国产一区二区三区在线臀色熟女 | 日本av免费视频播放| 亚洲一卡2卡3卡4卡5卡精品中文| 免费不卡黄色视频| 日韩大片免费观看网站| 欧美大码av| 精品卡一卡二卡四卡免费| 久久久久国产精品人妻一区二区| 嫩草影视91久久| 午夜免费观看性视频| av有码第一页| 国产精品久久久久久人妻精品电影 | 老汉色∧v一级毛片| 精品亚洲成a人片在线观看| 热99久久久久精品小说推荐| 亚洲av日韩精品久久久久久密| 十八禁网站网址无遮挡| www.999成人在线观看| 午夜免费成人在线视频| 亚洲 国产 在线| 12—13女人毛片做爰片一| 亚洲av成人不卡在线观看播放网 | 黄片大片在线免费观看| 亚洲伊人久久精品综合| 真人做人爱边吃奶动态| 啦啦啦啦在线视频资源| 午夜福利一区二区在线看| 免费高清在线观看日韩| 高清黄色对白视频在线免费看| 国产成人免费观看mmmm| 啦啦啦视频在线资源免费观看| 啦啦啦在线免费观看视频4| 精品国产乱子伦一区二区三区 | 久久毛片免费看一区二区三区| 国产日韩欧美在线精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产淫语在线视频| 日韩电影二区| avwww免费| 日韩视频在线欧美| 欧美av亚洲av综合av国产av| 90打野战视频偷拍视频| 满18在线观看网站| 亚洲熟女精品中文字幕| 熟女少妇亚洲综合色aaa.| 搡老乐熟女国产| 一本综合久久免费| 在线观看www视频免费| 国产精品香港三级国产av潘金莲| 一本综合久久免费| 嫩草影视91久久| 久久av网站| 欧美xxⅹ黑人| 这个男人来自地球电影免费观看| 男女下面插进去视频免费观看| 最新的欧美精品一区二区| 国产欧美日韩一区二区精品| 一级片'在线观看视频| 中文字幕av电影在线播放| 国产1区2区3区精品| 亚洲国产欧美网| 少妇 在线观看| 国产深夜福利视频在线观看| 国产又爽黄色视频| 男女无遮挡免费网站观看| 啦啦啦免费观看视频1| 另类亚洲欧美激情| 深夜精品福利| 国产亚洲精品久久久久5区| 免费高清在线观看视频在线观看| 国产精品偷伦视频观看了| 99久久国产精品久久久| 成人国语在线视频| 午夜成年电影在线免费观看| svipshipincom国产片| 国产片内射在线| av国产精品久久久久影院| 国产免费一区二区三区四区乱码| 亚洲熟女毛片儿| 久久人人97超碰香蕉20202| 成年动漫av网址| 狠狠婷婷综合久久久久久88av| tocl精华| 精品一区在线观看国产| 久久久久网色| 久久久久久久国产电影| 18禁国产床啪视频网站| 人人妻人人澡人人爽人人夜夜| 亚洲人成77777在线视频| 少妇裸体淫交视频免费看高清 | 一级毛片女人18水好多| 女人爽到高潮嗷嗷叫在线视频| 免费日韩欧美在线观看| a级毛片黄视频| 狠狠精品人妻久久久久久综合| 久久久精品免费免费高清| 国产在线免费精品| 亚洲av日韩精品久久久久久密| 亚洲国产欧美一区二区综合| 午夜福利在线观看吧| 久久九九热精品免费| 国产精品亚洲av一区麻豆| 亚洲av电影在线观看一区二区三区| 国产免费一区二区三区四区乱码| 欧美成人午夜精品| 欧美av亚洲av综合av国产av| 在线天堂中文资源库| 亚洲av欧美aⅴ国产| 精品国产一区二区久久| 搡老乐熟女国产| 久久久久久久久免费视频了| av在线播放精品| av电影中文网址| 国产精品久久久久久精品电影小说| 水蜜桃什么品种好| 老熟妇仑乱视频hdxx| 久久久国产欧美日韩av| a级片在线免费高清观看视频| 五月开心婷婷网| 亚洲性夜色夜夜综合| 久久人妻熟女aⅴ| 久久久久国内视频| 操美女的视频在线观看| 精品国产乱子伦一区二区三区 | 国产精品av久久久久免费| av福利片在线| 亚洲国产精品一区三区| 欧美国产精品va在线观看不卡| 天堂俺去俺来也www色官网| 精品人妻1区二区| 国产又色又爽无遮挡免| 视频区欧美日本亚洲| 伊人亚洲综合成人网| tocl精华| 国产日韩一区二区三区精品不卡| 侵犯人妻中文字幕一二三四区| 午夜福利视频在线观看免费| 男女之事视频高清在线观看| 亚洲伊人久久精品综合| 日日夜夜操网爽| 97精品久久久久久久久久精品| 一个人免费看片子| 人人妻人人澡人人爽人人夜夜| 久久精品亚洲av国产电影网| 2018国产大陆天天弄谢| 热99国产精品久久久久久7| 悠悠久久av| 一本一本久久a久久精品综合妖精| 视频区欧美日本亚洲| 两人在一起打扑克的视频| 精品国产国语对白av| 老汉色∧v一级毛片| 久久国产精品男人的天堂亚洲| 十八禁网站免费在线| 免费观看av网站的网址| 黄片大片在线免费观看| 最新在线观看一区二区三区| 窝窝影院91人妻| 欧美人与性动交α欧美精品济南到| 日韩中文字幕视频在线看片| 女人爽到高潮嗷嗷叫在线视频| 亚洲av成人不卡在线观看播放网 | 一本色道久久久久久精品综合| 在线观看舔阴道视频| 国产精品影院久久| 久久人妻福利社区极品人妻图片| 精品国产乱码久久久久久小说| 脱女人内裤的视频| 丰满饥渴人妻一区二区三| kizo精华| 久久免费观看电影| 欧美另类亚洲清纯唯美| 久久精品亚洲熟妇少妇任你| videos熟女内射| 亚洲性夜色夜夜综合| 午夜老司机福利片| 久久99热这里只频精品6学生| 少妇 在线观看| 国产淫语在线视频| 国产有黄有色有爽视频| 亚洲av欧美aⅴ国产| 视频在线观看一区二区三区| 一区在线观看完整版| 国产高清国产精品国产三级| 亚洲成人免费电影在线观看| 欧美精品高潮呻吟av久久| a级毛片在线看网站| 少妇的丰满在线观看| 欧美性长视频在线观看| 国产野战对白在线观看| 黑丝袜美女国产一区| 日本撒尿小便嘘嘘汇集6| 色播在线永久视频| 欧美亚洲日本最大视频资源| 国产成人系列免费观看| 国产av精品麻豆| 制服诱惑二区| 久久午夜综合久久蜜桃| 欧美国产精品va在线观看不卡| 国产亚洲av片在线观看秒播厂| 国产一区二区三区在线臀色熟女 | 日韩欧美免费精品| 亚洲精品日韩在线中文字幕| 在线av久久热| 老司机亚洲免费影院| 97在线人人人人妻| 亚洲国产欧美一区二区综合| 丰满人妻熟妇乱又伦精品不卡| 久久精品亚洲熟妇少妇任你| 欧美亚洲日本最大视频资源| 深夜精品福利| 欧美成狂野欧美在线观看| 久久久久久久久免费视频了| 欧美日韩亚洲高清精品| 亚洲,欧美精品.| 国产黄频视频在线观看| 伦理电影免费视频| 亚洲五月色婷婷综合| 欧美变态另类bdsm刘玥| 狠狠婷婷综合久久久久久88av| 一区二区av电影网| 亚洲色图综合在线观看| 亚洲熟女毛片儿| 热99久久久久精品小说推荐| 伊人亚洲综合成人网| 欧美一级毛片孕妇| 一本大道久久a久久精品| 老熟女久久久| 欧美激情极品国产一区二区三区| a在线观看视频网站| 多毛熟女@视频| 欧美日韩精品网址| 巨乳人妻的诱惑在线观看| 久久久水蜜桃国产精品网| 黄片大片在线免费观看| 亚洲精品国产色婷婷电影| 又紧又爽又黄一区二区| 热99久久久久精品小说推荐| 日韩人妻精品一区2区三区| 夫妻午夜视频| 下体分泌物呈黄色| 亚洲av成人不卡在线观看播放网 | 精品一区在线观看国产| 又紧又爽又黄一区二区| 美女高潮喷水抽搐中文字幕| 国产欧美日韩综合在线一区二区| 少妇猛男粗大的猛烈进出视频| 久久久久久久大尺度免费视频| 国产欧美亚洲国产| 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三区在线| 色94色欧美一区二区| 极品人妻少妇av视频| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕| 成人手机av| 黄片播放在线免费| 不卡一级毛片| 亚洲,欧美精品.| 操出白浆在线播放| 欧美激情久久久久久爽电影 | 黑人猛操日本美女一级片| 久久久久久免费高清国产稀缺| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 国精品久久久久久国模美| 三上悠亚av全集在线观看| 日韩欧美免费精品| av又黄又爽大尺度在线免费看| 亚洲av成人一区二区三| 99久久人妻综合| 亚洲一码二码三码区别大吗| 国产亚洲欧美精品永久| 丁香六月欧美| 日韩 亚洲 欧美在线| 国产精品久久久人人做人人爽| 久久免费观看电影| av一本久久久久| 亚洲精品成人av观看孕妇| 蜜桃在线观看..| 久久久久国内视频| 久久热在线av| 欧美日韩黄片免| 日韩,欧美,国产一区二区三区| 中文字幕人妻熟女乱码| 久久国产亚洲av麻豆专区| 老熟妇仑乱视频hdxx| 国产成人精品在线电影| 精品久久久久久电影网| 精品久久久精品久久久| 欧美av亚洲av综合av国产av| 国产伦人伦偷精品视频| 成人三级做爰电影| 欧美精品亚洲一区二区| 日韩欧美一区二区三区在线观看 | 欧美国产精品一级二级三级| 一边摸一边做爽爽视频免费| 国产三级黄色录像| 久久久精品国产亚洲av高清涩受| 亚洲欧美一区二区三区久久| 国产精品久久久av美女十八| 最近最新中文字幕大全免费视频| 国产男人的电影天堂91| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精| 十分钟在线观看高清视频www| 欧美日韩精品网址| 精品少妇黑人巨大在线播放| 最新在线观看一区二区三区| 黄色视频在线播放观看不卡| 精品乱码久久久久久99久播| 日韩制服骚丝袜av| 不卡av一区二区三区| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 国产视频一区二区在线看| 女性生殖器流出的白浆| 欧美黄色淫秽网站| 老司机影院成人| 欧美黄色片欧美黄色片| 国产一区二区在线观看av| 欧美国产精品va在线观看不卡| 欧美黑人精品巨大| tocl精华| 亚洲欧美清纯卡通| 久久久欧美国产精品| 99热国产这里只有精品6| 丰满少妇做爰视频| 成年动漫av网址| 国产av又大| 中文字幕精品免费在线观看视频| 两个人免费观看高清视频| 香蕉国产在线看| 99九九在线精品视频| 国产主播在线观看一区二区| 国产精品偷伦视频观看了| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 男人操女人黄网站| 亚洲国产精品成人久久小说| 国产91精品成人一区二区三区 | 伊人亚洲综合成人网| 免费在线观看日本一区| 精品亚洲成a人片在线观看| 桃红色精品国产亚洲av| 国产成人a∨麻豆精品| 国产淫语在线视频| 男人操女人黄网站| 老司机深夜福利视频在线观看 | 欧美另类亚洲清纯唯美| 超碰97精品在线观看| 久久中文字幕一级| 欧美另类一区| 91字幕亚洲| 人成视频在线观看免费观看| 人人妻人人添人人爽欧美一区卜| videos熟女内射| 一区在线观看完整版| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 熟女少妇亚洲综合色aaa.| a级毛片在线看网站| 精品一区二区三区av网在线观看 | 亚洲欧美一区二区三区久久| 蜜桃在线观看..| 少妇的丰满在线观看| 少妇粗大呻吟视频| 一级片免费观看大全| 亚洲欧美一区二区三区久久| a级片在线免费高清观看视频| 日韩大片免费观看网站| 在线观看www视频免费| 国产成人精品无人区| 黑人操中国人逼视频| 中文字幕人妻丝袜一区二区| 多毛熟女@视频| 十八禁高潮呻吟视频| 欧美精品一区二区免费开放| 岛国在线观看网站| 五月开心婷婷网| 国产片内射在线| 一个人免费看片子| 三级毛片av免费| 999精品在线视频| av超薄肉色丝袜交足视频| 国产亚洲精品一区二区www | 啦啦啦中文免费视频观看日本| 国产又爽黄色视频| 国产国语露脸激情在线看| 国产日韩欧美视频二区| 亚洲精品中文字幕一二三四区 | 脱女人内裤的视频| 考比视频在线观看| 欧美一级毛片孕妇| www.自偷自拍.com| 久久性视频一级片| 一边摸一边抽搐一进一出视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产男女内射视频| 91字幕亚洲| 在线观看免费日韩欧美大片| 岛国在线观看网站| 亚洲九九香蕉| 在线观看免费高清a一片| 久久久久精品人妻al黑| 国产亚洲午夜精品一区二区久久| 在线观看人妻少妇| 别揉我奶头~嗯~啊~动态视频 | 国产免费福利视频在线观看| netflix在线观看网站| 捣出白浆h1v1| 青春草亚洲视频在线观看| 日韩视频在线欧美| 两人在一起打扑克的视频| 亚洲国产成人一精品久久久| 亚洲五月婷婷丁香| 久久久久国产精品人妻一区二区| 日本a在线网址| 欧美乱码精品一区二区三区| 国产在线视频一区二区| 亚洲情色 制服丝袜| 18禁国产床啪视频网站| 欧美少妇被猛烈插入视频| 久久国产精品男人的天堂亚洲| 日本wwww免费看| 成年av动漫网址| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品日韩在线中文字幕| 男人爽女人下面视频在线观看| 丝袜美腿诱惑在线| 午夜福利视频精品| 十八禁网站免费在线| 黄片播放在线免费| 亚洲精品一二三| 精品一区在线观看国产| av欧美777| 久久久久国产一级毛片高清牌| www.自偷自拍.com| 亚洲伊人久久精品综合| 乱人伦中国视频| 99久久国产精品久久久| 亚洲国产成人一精品久久久| 国产av又大| 婷婷成人精品国产| 桃红色精品国产亚洲av| 黄色视频在线播放观看不卡| 免费女性裸体啪啪无遮挡网站| 99久久99久久久精品蜜桃| 亚洲精品国产精品久久久不卡| 久久久久久久精品精品| 久久天堂一区二区三区四区| 91麻豆av在线| avwww免费| 国产男女内射视频| 精品一区二区三卡| 人妻久久中文字幕网| 精品一区二区三卡| 老司机午夜福利在线观看视频 | 女性被躁到高潮视频| 久久99一区二区三区| 无限看片的www在线观看| 成人亚洲精品一区在线观看| 精品少妇内射三级| 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 亚洲,欧美精品.| 伦理电影免费视频| 国产一卡二卡三卡精品| 久久热在线av| 精品乱码久久久久久99久播| 制服人妻中文乱码| 国产欧美日韩一区二区三 | 久久久久网色| 高清黄色对白视频在线免费看| 久久精品亚洲熟妇少妇任你| 肉色欧美久久久久久久蜜桃| av不卡在线播放| 51午夜福利影视在线观看| 国产老妇伦熟女老妇高清| 男女免费视频国产| 狂野欧美激情性bbbbbb| 1024视频免费在线观看| 又紧又爽又黄一区二区| 亚洲精品日韩在线中文字幕| 国产精品1区2区在线观看. | 国产高清videossex| 欧美激情 高清一区二区三区| 黄片播放在线免费| 国产欧美日韩精品亚洲av| 日日夜夜操网爽| 亚洲精品国产精品久久久不卡| 真人做人爱边吃奶动态| 91成年电影在线观看| 91九色精品人成在线观看| 国产av精品麻豆| 国产成+人综合+亚洲专区| 亚洲一码二码三码区别大吗| 国产老妇伦熟女老妇高清| 亚洲精品自拍成人| 久久人人爽人人片av| 亚洲美女黄色视频免费看| 少妇粗大呻吟视频| 亚洲免费av在线视频| 国产精品麻豆人妻色哟哟久久| 91麻豆av在线| 午夜福利乱码中文字幕| 久久午夜综合久久蜜桃| 在线天堂中文资源库| 中文字幕高清在线视频| 91av网站免费观看| 一二三四在线观看免费中文在| 午夜福利在线观看吧| 国产无遮挡羞羞视频在线观看| 国产黄频视频在线观看| 一二三四在线观看免费中文在| 男女之事视频高清在线观看| 老司机深夜福利视频在线观看 | 久久国产精品人妻蜜桃| 久久亚洲精品不卡| 黄片播放在线免费| 成人免费观看视频高清| 美国免费a级毛片| 天天操日日干夜夜撸| 亚洲熟女毛片儿| 日韩有码中文字幕| 天天躁夜夜躁狠狠躁躁| 可以免费在线观看a视频的电影网站| 精品一品国产午夜福利视频| 欧美久久黑人一区二区|