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

    基于BLUP的單性狀動物育種模型的研究及實現(xiàn)

    2016-12-02 01:45:43胡世洋楊慧贊張永德林勇孫全亮
    賀州學(xué)院學(xué)報 2016年1期
    關(guān)鍵詞:關(guān)聯(lián)矩陣遺傳力賀州

    胡世洋,楊慧贊,張永德,林勇,孫全亮

    (1.賀州學(xué)院,計算機(jī)科學(xué)與信息工程學(xué)院,廣西賀州542899;2.廣西水產(chǎn)科學(xué)研究院,廣西南寧530021)

    基于BLUP的單性狀動物育種模型的研究及實現(xiàn)

    胡世洋1,楊慧贊2,張永德2,林勇2,孫全亮1

    (1.賀州學(xué)院,計算機(jī)科學(xué)與信息工程學(xué)院,廣西賀州542899;2.廣西水產(chǎn)科學(xué)研究院,廣西南寧530021)

    BLUP算法可以顯著提高動物遺傳進(jìn)展速度,但是它的計算過程比較復(fù)雜,數(shù)據(jù)量也相當(dāng)浩繁,若將這種復(fù)雜、抽象的數(shù)學(xué)算法轉(zhuǎn)換成可視化的計算機(jī)系統(tǒng)將大大有利于生產(chǎn)力的提高,針對這一問題,以羅非魚的遺傳育種為研究切入點,將BLUP算法應(yīng)用于羅非魚的育種值估計當(dāng)中,構(gòu)建出系統(tǒng)的模型,并采用面向?qū)ο蟮乃季S方法將BLUP的算法轉(zhuǎn)化為可執(zhí)行程序。實驗結(jié)果表明,開發(fā)出的軟件可以方便、高效地進(jìn)行羅非魚育種值的估計,當(dāng)影響羅非魚的外部因素發(fā)生改變時,最終的計算結(jié)果都能較準(zhǔn)確地反應(yīng)出相關(guān)因素對羅非魚的影響效應(yīng)。

    BLUP;算法;羅非魚;軟件

    1 引言

    動物的選擇育種問題一直是養(yǎng)殖業(yè)的重要科題之一,隨著數(shù)理統(tǒng)計學(xué)與線性模型理論、計算機(jī)科學(xué)與互聯(lián)網(wǎng)絡(luò)技術(shù)的迅速發(fā)展,動物的育種值估計的方法發(fā)生了根本的變化[1]。最佳線性無偏預(yù)測法(Best Linear Unbiased Prediction,BLUP)由于能顯著提高遺傳進(jìn)展速度,特別是對于低遺傳力性狀和顯性性狀,因而被廣泛運用到動物的遺傳改良當(dāng)中。BLUP法可以通過對稱矩陣求解的方式把動物的種類效應(yīng)、動物自身效應(yīng)和機(jī)誤誤差區(qū)分開來,可以有效提高動物遺傳值的準(zhǔn)確性[2]。由于BLUP法的過程比較復(fù)雜,計算量也相當(dāng)浩繁,如何將BLUP算法的數(shù)學(xué)思想轉(zhuǎn)化成計算機(jī)思維,用程序設(shè)計語言將其運行并且實現(xiàn)代碼的可重用性是一個值得深入研究的問題[3]。文章以羅非魚的養(yǎng)殖育種為例研究BLUP算法的應(yīng)用及實現(xiàn)方法。

    2 BLUP算法的工作原理

    一般混合模型可表示為:

    y是所有觀察值構(gòu)成的向量;b是所有固定效應(yīng)(包括)構(gòu)成的向量;X是固定效應(yīng)的關(guān)聯(lián)矩陣;u是所有隨機(jī)效應(yīng)構(gòu)成的向量;Z是隨機(jī)效應(yīng)的關(guān)聯(lián)矩陣;e是隨機(jī)殘差向量[4]。

    隨機(jī)變量的數(shù)學(xué)期望:

    方差-協(xié)方差矩陣結(jié)構(gòu):

    BLUP估計一般方程:

    混合模型方程組的簡化形式:

    3 BLUP算法在羅非魚研究過程中的具體研究及分析

    某種魚場有如下性能測定資料,測定性狀為達(dá)500g日齡,已知該性狀的遺傳力為0.33,已知某群達(dá)500g的羅非魚相關(guān)記錄數(shù)據(jù)如表1所示,下面將對該群羅非魚的性狀資料進(jìn)行個體育種值估計。

    表1 羅非魚達(dá)500g日齡記錄

    根據(jù)資料性質(zhì),可對羅非魚達(dá)500g日齡寫出如下動物模型:

    用矩陣形式表示,則對于該資料有:

    構(gòu)建混合模型方程組:

    求解混合模型方程組有:

    4 Blup算法在計算機(jī)中的算法實現(xiàn)方案

    4.1系統(tǒng)框架的構(gòu)建

    系統(tǒng)采用MVC的設(shè)計框架,由ODBC接口、相關(guān)數(shù)據(jù)庫構(gòu)件以及文件流操作模塊組成Model層[5],由矩陣生成模塊、矩陣運算模塊、方程求解模塊以及相關(guān)賦值運算模塊組成Controller層,由于相關(guān)顯示窗體模塊組成View層[6]。系統(tǒng)可獲取的外部數(shù)據(jù)類型包括數(shù)據(jù)庫文件,XML文件、txt文件或xls文件,系統(tǒng)框架如果圖1所示。

    圖1 系統(tǒng)框架圖

    4.2系統(tǒng)相關(guān)操作類的實現(xiàn)

    本系統(tǒng)是采用C#程序語言進(jìn)行開發(fā),在系統(tǒng)實現(xiàn)過程中需要大量的矩陣運算,因此在系統(tǒng)中必須建立Matrix類來管理矩陣數(shù)據(jù),同時需要有一個MatrixOperator類來操作Matrix類,從而實現(xiàn)矩陣的混合運算及相關(guān)信息的輸出[7]。系統(tǒng)要從數(shù)據(jù)庫文件中獲取外部數(shù)據(jù)文件并存取在內(nèi)存當(dāng)中,因此需要有DataModel類來負(fù)責(zé)外部數(shù)據(jù)的獲取來組織[8]。最后內(nèi)存中的數(shù)據(jù)需要在AnalysisOperator類中完成矩陣的構(gòu)建和運算,并將最終結(jié)果輸出到顯示窗體當(dāng)中。系統(tǒng)的相關(guān)類型如圖2所示。

    圖2 系統(tǒng)類圖

    根據(jù)上一節(jié)基于BLUP算法羅非魚育種值計算規(guī)則,可以構(gòu)建出其在計算機(jī)中的實現(xiàn)算法為:

    步驟1:讀取數(shù)據(jù),基礎(chǔ)數(shù)據(jù)可以是來自數(shù)據(jù),或者是其它外部文件。

    步驟2:構(gòu)建系數(shù)矩陣。針對羅非魚的數(shù)據(jù)特征,可以對體重達(dá)500g的日齡構(gòu)建動物模型,并由此得出固定因素和非固定因素的系數(shù)據(jù)矩陣。

    步驟3:構(gòu)建親緣矩陣。根據(jù)羅非魚的個體系譜列表來構(gòu)建分子親緣矩陣。

    步驟4:構(gòu)建混合模型矩陣。根據(jù)BLUP估計一般方程,混合模型方程組的簡化形式。

    步驟5:混合矩陣方程的求解。得到了矩陣向量元素將是各個影響因素的育種估計值。

    4.3算法實現(xiàn)的幾個具體問題

    4.3.1矩陣類的定義

    矩陣的存儲方式是采用二維數(shù)組存儲。數(shù)組類型為double型,Matrix類中不僅有行和列的相關(guān)屬性,而且有其對應(yīng)的構(gòu)造方法以及實例化矩陣的相關(guān)方法。

    4.3.2矩陣操作類的定義算法

    BLUP算法在運算過程中涉及矩陣相乘、矩陣數(shù)乘、矩陣轉(zhuǎn)置和矩陣求逆等相關(guān)操作。因此我們在MatrixOperator類必須要聲明這些方法,其中MatrixPrint()可以實現(xiàn)矩陣內(nèi)部元素的輸出。

    4.3.3系數(shù)矩陣的構(gòu)建算法

    由上一節(jié)的動物模型可以得知,在單性狀遺傳值估計中,固定效應(yīng)的關(guān)聯(lián)矩陣X的維數(shù)是由測量記錄數(shù)和魚池這一固定效應(yīng)的數(shù)量決定的,而矩陣X的元素值是因固定效應(yīng)的量值決定。關(guān)聯(lián)矩陣X的構(gòu)建算法則是先確定魚場的個數(shù),再遍歷測量記錄條目,若當(dāng)前記錄的魚場數(shù)與先前設(shè)置的魚池編號相同則將其在X矩陣中所對應(yīng)的無素值置為1,其它元素置為0。而隨機(jī)效應(yīng)的關(guān)聯(lián)矩陣Z則是與測量記錄數(shù)相關(guān)的單位矩陣。顯性向量矩陣y則是由各條記錄數(shù)中的日齡字段決定。這兩個矩陣的構(gòu)建較為簡單,在此不再贅述。

    4.3.4親緣矩陣的構(gòu)建算法

    構(gòu)建親緣矩陣的基本規(guī)則:(1)當(dāng)編號為t的個體的父母均未知時,att=1,ati=ait=0,i=1、2…t-1。(2)當(dāng)編號為t的個體的父母中的一個已知為p,而另一個未知時,att=1,ati=ait=0.5aip,i=1、2…t-1。(3)當(dāng)編號為t的個體的父母均已知(假定為p和q)時,att=1+0.5apq,ati=ait=0.5(aip+aiq),i=1、2…t-1。因此確定親緣矩陣的構(gòu)建算法為:(1)將羅非魚的觀測記錄讀入到數(shù)組arr[][]當(dāng)中,創(chuàng)建空矩陣A,并使用數(shù)組a[][]來存儲矩陣的元素。(2)循環(huán)數(shù)組arr[][]中行記錄即逐個確定觀測記錄中的每一條記錄。在第t(1<=t<=arr.length)條記錄中,假定這條記錄的第a列是父本編號,第b列是其母本編號,若arr[t,a]和arr[t,b]的值均為空,那么判定t是不是第一條記錄,如果是那么a[t,t]=1。否則建立t的循環(huán)變量i(1<=i<=t),a[t,t]=1;a[t,i]=a[i,t]=0。(3)若arr[t,a]和arr[t,b]的值只有一個不為空記為p,那么判定t是不是第一條記錄,如果是那么a[t,t]= 1。否則建立t的循環(huán)變量i(1<=i<=t),a[t,t]=1,a[t,i]=a[i,t]=0.5*a[i,p]。(4)若arr[t,a]和arr[t,b]的值分別為p和q時,那么判定t是不是第一條記錄,如果是那么a[t,t]=1。否則建立t的循環(huán)變量i(1<=i<=t),a[t,t]=1+0.5*a[p,q],a[t,i]=a[i,t]= 0.5*(a[i,p]+a[i,q])。

    4.3.5混合矩陣的構(gòu)建算法

    5 系統(tǒng)的運行結(jié)果及應(yīng)用分析

    系統(tǒng)測試環(huán)境:硬件環(huán)境:(1)硬件環(huán)境:PC機(jī)一臺。(2)軟件環(huán)境:Window7操作系統(tǒng),Microsoft. NET Framework 4.0組件,SQL server2005數(shù)據(jù)庫。

    實例測試1:提供2個魚場的5條記錄,并輸入性狀遺傳力值0.33,最終計算的育種值向量值為:(h1 h2 a1 a2 a3 a4 a5)=(114.44 112.72-0.46 2.57 -3.404.70-3.47)。測量記錄數(shù)據(jù)如圖3所示,系統(tǒng)輸出結(jié)果如圖4所示。

    圖3 第一次測量記錄數(shù)據(jù)

    圖4 第一次實驗輸出結(jié)果

    實例測試2:提供3個魚場的6條記錄,并輸入性狀遺傳力值0.33,最終計算的育種值向量值為:(h1h2h3a1a2a3a4a5)=(142.84152.12169.62-2. 463.05-2.12-1.392.150.38)。測量記錄數(shù)據(jù)如圖5所示,實現(xiàn)輸出結(jié)果如圖6所示。

    圖5 第二次測量記錄數(shù)據(jù)

    圖6 第二次實驗輸出結(jié)果

    6 結(jié)語

    文章分析了BLUP算法的基本原理,并研究了其在羅非魚單性狀遺傳育種研究中的實現(xiàn)方法。在系統(tǒng)建模型的基礎(chǔ)上,開發(fā)了一套與之相配套的軟件系統(tǒng),并對軟件實現(xiàn)過程的具體算法進(jìn)行了研究和描述。經(jīng)過軟件的測試和運行表明,對于單性態(tài)的育種案例,本軟件可以順利計算出羅非魚的育種值,因此本軟件有著一定的現(xiàn)實意義。

    BLUP算法不僅可以用于單性狀育種值估計,而且同樣可以應(yīng)用于多性狀的育種值估計?,F(xiàn)實生活中影響魚類遺傳育種的性狀常常不僅僅一個。因此多性狀的育種值估計的研究和實現(xiàn)將是下一步工作的主要內(nèi)容。

    [1]曾聰,曹小娟,高澤霞,等.團(tuán)頭魴生長性狀的遺傳力和育種值估計[J].華中農(nóng)業(yè)大學(xué)學(xué)報,2014,33(2):89-95.

    [2]楊瑞飛,張凱,汪平,等.獺兔育種與生產(chǎn)信息管理系統(tǒng)的研發(fā)[J].四川農(nóng)業(yè)大學(xué)學(xué)報,2014,32(1):82-86.

    [3]劉峰,陳松林,王磊,等.不同牙鲆群體遺傳力和育種值分析[J].中國水產(chǎn)科學(xué),2013,20(4):691-697.

    [4]劉瑩瑩,肖兵南,燕海峰,等.利用BLUP法估計荷斯坦種公牛的抗熱應(yīng)激育種值[J].家畜生態(tài)學(xué)報,2012,33(1):29-32.

    [5]王光瓊,杜天行,未永慶,等.基于SOA構(gòu)架的分布式租車公司管理系統(tǒng)設(shè)計與實現(xiàn)[J].軟件導(dǎo)刊,2014,13(6):52-54.

    [6]岳紹敏,王國春,董亞則,等.基于Spring MVC框架的電商平臺的研究與設(shè)計[J].軟件工程師,2014,17(7):3-7.

    [7]鄒燕燕,安虹,崔慧敏,等.矩陣形狀無關(guān)的并行編譯優(yōu)化[J].小型微型計算機(jī)系統(tǒng),2014,35(7):1518-1522.

    [8]Seo-Young NOH,Shashi K.GADIA,Haengjin JANG.Comparisons of three data storage models in parametric temporal databases[J].Journal of Central South University,2013,20(7):1919-1927.

    [責(zé)任編輯]張琴芳

    Research and Implementation of Single Trait Animal Breeding Models Based BLUP

    HuShi-yang1,YANG Hui-zan2,ZHANG Yong-de2,LINYong2,SUNQuan-liang1
    (1.School ofComputerScienceandInformationEngineering,Hezhouuniversity,HezhouGuangxi 542899;2.Guangxi AcademyofFisherySciences,NanningGuangxi 530021)

    BLUP algorithm can significantly improve the speed of animal genetic progress,but the calculations are more complex and the data size is quite voluminous.If this complex abstract mathematical algorithm changes into a visualization computer system,it will greatly help to improve productivity for this problem.This paper applies the BLUP algorithm in tilapia breeding value estimation to build the model of the system and use an object-oriented way of thinking to turn the BLUP into executable program.Experimental results show that the development of the software can easily and efficiently estimate breeding values of tilapia.When external factors change,the final result can more accurately reflect the effect of related factors affecting tilapia.

    BLUP;algorithm;tilapia;software

    TP311.11

    A

    1673—8861(2016)01—0144—05

    2015-02-23

    胡世洋(1987-),男,河南南陽人,賀州學(xué)院講師,碩士。主要研究方向:計算機(jī)網(wǎng)絡(luò)與信息管理系統(tǒng)。

    廣西水產(chǎn)遺傳育種與健康養(yǎng)殖重點實驗室資助項目(GXKL-AQUA-2013-C-05)、賀州學(xué)院2014年度校級科研項目(2014YBZK09)、賀州學(xué)院2013年度校級科研項目(2013ZXZK14)。

    猜你喜歡
    關(guān)聯(lián)矩陣遺傳力賀州
    某國家生豬核心育種場加系大白母豬繁殖性狀遺傳參數(shù)估計
    n階圈圖關(guān)聯(lián)矩陣的特征值
    美系大白種豬生長性狀和繁殖性狀遺傳力估計
    養(yǎng)豬(2022年4期)2022-08-17 07:07:02
    游賀州紫云洞
    單圈圖關(guān)聯(lián)矩陣的特征值
    印象賀州
    黃河之聲(2021年20期)2021-02-24 02:48:52
    基于關(guān)聯(lián)矩陣主對角線譜理論的歐拉圖研究
    n階圈圖的一些代數(shù)性質(zhì)
    我的家鄉(xiāng)最美之賀州
    吾城·吾鄉(xiāng)
    成在线人永久免费视频| 欧美绝顶高潮抽搐喷水| 亚洲中文字幕日韩| 色综合亚洲欧美另类图片| videosex国产| 久久午夜亚洲精品久久| 色尼玛亚洲综合影院| 久久精品国产综合久久久| 亚洲五月色婷婷综合| 亚洲美女黄片视频| or卡值多少钱| 99在线视频只有这里精品首页| 国产精品永久免费网站| 午夜激情av网站| 亚洲美女黄片视频| 男女之事视频高清在线观看| 精品一区二区三区四区五区乱码| 久久国产精品人妻蜜桃| 国产乱人伦免费视频| 亚洲欧美日韩无卡精品| 国产亚洲精品久久久久5区| 波多野结衣一区麻豆| 伦理电影免费视频| 亚洲五月色婷婷综合| 99re在线观看精品视频| 国产成人精品久久二区二区91| 午夜老司机福利片| 在线国产一区二区在线| 人人妻,人人澡人人爽秒播| 亚洲精品粉嫩美女一区| 国产欧美日韩综合在线一区二区| 91老司机精品| 色哟哟哟哟哟哟| 国产精品免费一区二区三区在线| 亚洲人成电影观看| 国产在线精品亚洲第一网站| 亚洲 欧美 日韩 在线 免费| 国产精品一区二区三区四区久久 | 纯流量卡能插随身wifi吗| 美女午夜性视频免费| 男人舔女人下体高潮全视频| 满18在线观看网站| 久久伊人香网站| 亚洲av五月六月丁香网| 欧美午夜高清在线| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 老熟妇仑乱视频hdxx| 国产精品一区二区三区四区久久 | 日本一区二区免费在线视频| 色综合婷婷激情| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 亚洲人成电影观看| 免费观看精品视频网站| 国产欧美日韩一区二区三区在线| 国语自产精品视频在线第100页| 午夜免费成人在线视频| 久久人妻熟女aⅴ| 在线av久久热| 久久久久国产精品人妻aⅴ院| 无遮挡黄片免费观看| 黄色视频,在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产激情欧美一区二区| 国产男靠女视频免费网站| 淫秽高清视频在线观看| 少妇的丰满在线观看| 一本大道久久a久久精品| 欧美日本亚洲视频在线播放| 一夜夜www| 女人精品久久久久毛片| 欧美激情久久久久久爽电影 | 亚洲第一电影网av| 国产精品1区2区在线观看.| 啦啦啦免费观看视频1| 午夜亚洲福利在线播放| 国产激情久久老熟女| 岛国在线观看网站| 一区二区日韩欧美中文字幕| 国产亚洲av嫩草精品影院| 久久久久久久久免费视频了| 日韩中文字幕欧美一区二区| 啦啦啦 在线观看视频| 69精品国产乱码久久久| 黄色毛片三级朝国网站| 亚洲美女黄片视频| 久久久久久久久久久久大奶| 法律面前人人平等表现在哪些方面| 日韩精品中文字幕看吧| 国产精品久久久人人做人人爽| 久久亚洲真实| 中文字幕最新亚洲高清| 动漫黄色视频在线观看| 欧美黄色片欧美黄色片| 欧美日韩福利视频一区二区| 久久精品亚洲精品国产色婷小说| 久久精品aⅴ一区二区三区四区| 麻豆av在线久日| 中文字幕最新亚洲高清| 日本免费一区二区三区高清不卡 | 又黄又粗又硬又大视频| 99re在线观看精品视频| 丝袜人妻中文字幕| 黄色毛片三级朝国网站| 免费少妇av软件| 亚洲av熟女| 狂野欧美激情性xxxx| www日本在线高清视频| 国产又色又爽无遮挡免费看| 又大又爽又粗| 91老司机精品| 婷婷精品国产亚洲av在线| 露出奶头的视频| 亚洲人成伊人成综合网2020| 不卡av一区二区三区| 欧美在线黄色| 1024香蕉在线观看| 久久久久久大精品| 少妇裸体淫交视频免费看高清 | 婷婷丁香在线五月| 亚洲在线自拍视频| av电影中文网址| 一级黄色大片毛片| 亚洲成人久久性| 精品国产一区二区三区四区第35| 亚洲精品粉嫩美女一区| 亚洲国产欧美一区二区综合| 国产成人一区二区三区免费视频网站| 精品久久蜜臀av无| 淫妇啪啪啪对白视频| 国产精品国产高清国产av| 精品一区二区三区视频在线观看免费| 美女大奶头视频| 国产三级黄色录像| 久久久久久久久久久久大奶| aaaaa片日本免费| 91大片在线观看| 午夜福利影视在线免费观看| 欧美日韩黄片免| 久久亚洲精品不卡| 亚洲九九香蕉| 99香蕉大伊视频| 国产亚洲欧美精品永久| 热re99久久国产66热| 国产午夜福利久久久久久| 日韩精品免费视频一区二区三区| 国产精品爽爽va在线观看网站 | 亚洲九九香蕉| 99久久国产精品久久久| 精品久久久久久久久久免费视频| 老熟妇乱子伦视频在线观看| 亚洲精品一区av在线观看| 亚洲熟妇熟女久久| 电影成人av| 中文字幕色久视频| 非洲黑人性xxxx精品又粗又长| 国产av一区在线观看免费| 亚洲国产日韩欧美精品在线观看 | www日本在线高清视频| 97碰自拍视频| 久久精品国产清高在天天线| 国产成人精品久久二区二区91| 国产成人一区二区三区免费视频网站| 天堂影院成人在线观看| 国产伦一二天堂av在线观看| 午夜福利高清视频| 88av欧美| 最近最新免费中文字幕在线| 亚洲人成电影免费在线| 日韩大尺度精品在线看网址 | 亚洲熟妇熟女久久| 日韩欧美免费精品| 亚洲 欧美 日韩 在线 免费| 91字幕亚洲| 色综合亚洲欧美另类图片| 日日干狠狠操夜夜爽| 亚洲人成电影免费在线| 色尼玛亚洲综合影院| 亚洲国产日韩欧美精品在线观看 | 国产aⅴ精品一区二区三区波| 日韩国内少妇激情av| 精品电影一区二区在线| 欧美激情 高清一区二区三区| 大型av网站在线播放| 三级毛片av免费| 国产精品一区二区免费欧美| 岛国在线观看网站| 日韩成人在线观看一区二区三区| 国产三级黄色录像| 午夜精品在线福利| 国产aⅴ精品一区二区三区波| 97人妻天天添夜夜摸| 久久国产乱子伦精品免费另类| 视频区欧美日本亚洲| 麻豆国产av国片精品| 久久人妻av系列| 国语自产精品视频在线第100页| 久久人人精品亚洲av| 亚洲男人天堂网一区| 69av精品久久久久久| 久久久精品国产亚洲av高清涩受| 人人澡人人妻人| 在线观看www视频免费| 精品国产美女av久久久久小说| 亚洲人成77777在线视频| 在线观看免费日韩欧美大片| 精品午夜福利视频在线观看一区| 韩国av一区二区三区四区| 999精品在线视频| 精品人妻在线不人妻| 制服丝袜大香蕉在线| 午夜老司机福利片| 久久国产精品影院| 亚洲av日韩精品久久久久久密| 大码成人一级视频| 亚洲精品国产区一区二| 高清黄色对白视频在线免费看| 精品第一国产精品| 搡老岳熟女国产| 青草久久国产| 少妇熟女aⅴ在线视频| www.999成人在线观看| 真人一进一出gif抽搐免费| 国内精品久久久久精免费| 国产亚洲av高清不卡| 欧美午夜高清在线| 性色av乱码一区二区三区2| 中文字幕人妻熟女乱码| 久久久久久久久免费视频了| 精品不卡国产一区二区三区| 黄色女人牲交| 视频在线观看一区二区三区| 午夜成年电影在线免费观看| 9热在线视频观看99| 国产xxxxx性猛交| 18美女黄网站色大片免费观看| 亚洲欧美精品综合久久99| 欧美乱色亚洲激情| 人人妻,人人澡人人爽秒播| 久久久精品欧美日韩精品| 欧美日韩福利视频一区二区| 中文字幕av电影在线播放| 国产亚洲精品久久久久久毛片| 国产男靠女视频免费网站| 国产精品九九99| 久9热在线精品视频| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 波多野结衣av一区二区av| 久久中文字幕一级| 欧美丝袜亚洲另类 | 免费搜索国产男女视频| 中文字幕高清在线视频| 国产成年人精品一区二区| 老司机福利观看| 少妇粗大呻吟视频| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| 欧美激情极品国产一区二区三区| 国产在线精品亚洲第一网站| 男女下面进入的视频免费午夜 | 欧美一级毛片孕妇| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 亚洲欧美日韩另类电影网站| 午夜亚洲福利在线播放| 中文字幕av电影在线播放| 色综合亚洲欧美另类图片| 成年女人毛片免费观看观看9| 女人被躁到高潮嗷嗷叫费观| 午夜免费激情av| 免费在线观看视频国产中文字幕亚洲| 大陆偷拍与自拍| 午夜久久久久精精品| 91成人精品电影| 久久香蕉激情| 欧美成人午夜精品| 免费人成视频x8x8入口观看| 一区二区三区国产精品乱码| 88av欧美| 最近最新中文字幕大全电影3 | 美女高潮喷水抽搐中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 村上凉子中文字幕在线| 国产成人一区二区三区免费视频网站| 一进一出抽搐gif免费好疼| 电影成人av| 一本大道久久a久久精品| 国产精品九九99| 国产日韩一区二区三区精品不卡| 久久久国产欧美日韩av| 99久久精品国产亚洲精品| 999精品在线视频| 欧美成人一区二区免费高清观看 | 亚洲五月婷婷丁香| 成人国产综合亚洲| 日日干狠狠操夜夜爽| 男男h啪啪无遮挡| 亚洲,欧美精品.| 欧美日韩亚洲国产一区二区在线观看| 色在线成人网| 十分钟在线观看高清视频www| 国产成人欧美| 成年人黄色毛片网站| 一二三四在线观看免费中文在| 国产麻豆69| 久久久久国产精品人妻aⅴ院| 国产精品精品国产色婷婷| 成人欧美大片| 99在线人妻在线中文字幕| 日韩欧美在线二视频| 亚洲精品在线观看二区| 免费在线观看亚洲国产| 日韩免费av在线播放| 国产色视频综合| 黄色成人免费大全| 大型av网站在线播放| 午夜福利影视在线免费观看| 久久精品人人爽人人爽视色| 亚洲第一电影网av| 国产蜜桃级精品一区二区三区| 精品久久蜜臀av无| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码| 亚洲国产欧美日韩在线播放| 久久香蕉激情| 欧美一区二区精品小视频在线| 91国产中文字幕| 少妇 在线观看| 亚洲国产精品成人综合色| 一区二区三区国产精品乱码| 久久久国产精品麻豆| 免费在线观看日本一区| 99在线人妻在线中文字幕| 久久人人爽av亚洲精品天堂| 岛国视频午夜一区免费看| 在线播放国产精品三级| 久久婷婷人人爽人人干人人爱 | 国产成年人精品一区二区| 级片在线观看| 精品人妻在线不人妻| 人妻久久中文字幕网| 在线视频色国产色| 久久草成人影院| 午夜影院日韩av| 日韩一卡2卡3卡4卡2021年| 亚洲成人久久性| 久久九九热精品免费| 亚洲 欧美 日韩 在线 免费| 男人操女人黄网站| 日本免费一区二区三区高清不卡 | 美女午夜性视频免费| 中文字幕色久视频| 在线观看一区二区三区| 亚洲第一电影网av| 宅男免费午夜| 亚洲欧美日韩高清在线视频| 欧美激情极品国产一区二区三区| 久久 成人 亚洲| 国产黄a三级三级三级人| 怎么达到女性高潮| 成人免费观看视频高清| 波多野结衣巨乳人妻| 久久香蕉国产精品| 国产成人精品久久二区二区免费| 国产精品精品国产色婷婷| 国产亚洲精品av在线| 淫秽高清视频在线观看| 亚洲专区字幕在线| 51午夜福利影视在线观看| 啦啦啦韩国在线观看视频| av欧美777| 欧美日韩福利视频一区二区| 色婷婷久久久亚洲欧美| 男女做爰动态图高潮gif福利片 | www.精华液| 自线自在国产av| 这个男人来自地球电影免费观看| 欧美日本视频| 免费在线观看完整版高清| 精品久久蜜臀av无| 久久人人97超碰香蕉20202| av网站免费在线观看视频| 精品国产国语对白av| 宅男免费午夜| 欧美黄色淫秽网站| 午夜激情av网站| 国内毛片毛片毛片毛片毛片| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 国产精品久久久久久人妻精品电影| 真人做人爱边吃奶动态| 男女做爰动态图高潮gif福利片 | 人人澡人人妻人| 国产不卡一卡二| 亚洲欧美激情综合另类| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 久久人妻av系列| 国产男靠女视频免费网站| 精品欧美国产一区二区三| 欧美绝顶高潮抽搐喷水| 高潮久久久久久久久久久不卡| 色播在线永久视频| 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 国产亚洲欧美精品永久| 午夜视频精品福利| 视频区欧美日本亚洲| 成人国语在线视频| 999久久久精品免费观看国产| 啦啦啦 在线观看视频| 国产97色在线日韩免费| 国产精品免费视频内射| 国产精品av久久久久免费| 老司机福利观看| 激情视频va一区二区三区| 中文字幕人成人乱码亚洲影| 色精品久久人妻99蜜桃| 我的亚洲天堂| 日韩成人在线观看一区二区三区| 久久精品成人免费网站| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成av片中文字幕在线观看| 国产av一区二区精品久久| 午夜福利视频1000在线观看 | 波多野结衣巨乳人妻| 国产成人欧美| 一级,二级,三级黄色视频| 成人精品一区二区免费| 91麻豆精品激情在线观看国产| 亚洲人成电影免费在线| 国产视频一区二区在线看| 亚洲第一av免费看| 美女 人体艺术 gogo| 变态另类成人亚洲欧美熟女 | 国产成人av激情在线播放| 村上凉子中文字幕在线| 热99re8久久精品国产| 又黄又爽又免费观看的视频| 免费av毛片视频| 91麻豆av在线| 免费在线观看视频国产中文字幕亚洲| 琪琪午夜伦伦电影理论片6080| 一本综合久久免费| 男女之事视频高清在线观看| 九色亚洲精品在线播放| 午夜福利在线观看吧| 一二三四社区在线视频社区8| 丁香六月欧美| 黄频高清免费视频| 国产免费男女视频| 亚洲人成电影观看| ponron亚洲| 久久香蕉激情| 无遮挡黄片免费观看| 欧美丝袜亚洲另类 | 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品久久久久5区| 在线国产一区二区在线| 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 18禁黄网站禁片午夜丰满| 亚洲国产毛片av蜜桃av| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 国产麻豆成人av免费视频| 18禁观看日本| 国产成人av教育| 免费在线观看影片大全网站| 一进一出抽搐动态| 亚洲 国产 在线| 这个男人来自地球电影免费观看| 欧美日韩精品网址| 亚洲中文av在线| 黄色丝袜av网址大全| 国产又色又爽无遮挡免费看| 亚洲人成伊人成综合网2020| 999久久久国产精品视频| 色综合站精品国产| 制服人妻中文乱码| 精品久久久久久,| 久久九九热精品免费| 757午夜福利合集在线观看| 欧美日韩亚洲综合一区二区三区_| 性欧美人与动物交配| 99久久久亚洲精品蜜臀av| 午夜免费激情av| 咕卡用的链子| 亚洲视频免费观看视频| 亚洲国产精品久久男人天堂| 嫩草影视91久久| 国产成人一区二区三区免费视频网站| 侵犯人妻中文字幕一二三四区| 啪啪无遮挡十八禁网站| 国产成人精品久久二区二区免费| 黄色视频不卡| 亚洲第一电影网av| 精品人妻在线不人妻| 精品一区二区三区四区五区乱码| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| 久久久久久亚洲精品国产蜜桃av| 国产男靠女视频免费网站| 成年人黄色毛片网站| 一区二区三区激情视频| 夜夜爽天天搞| 老司机在亚洲福利影院| 欧美日韩亚洲综合一区二区三区_| 91在线观看av| 老汉色av国产亚洲站长工具| 可以在线观看的亚洲视频| 日韩免费av在线播放| 欧美在线一区亚洲| 不卡av一区二区三区| 黄色视频不卡| av天堂久久9| 亚洲国产毛片av蜜桃av| 99久久综合精品五月天人人| 日韩免费av在线播放| 操出白浆在线播放| www日本在线高清视频| 国产亚洲精品久久久久5区| 亚洲色图综合在线观看| 日韩精品免费视频一区二区三区| 禁无遮挡网站| 亚洲av片天天在线观看| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 好男人电影高清在线观看| 久久久精品欧美日韩精品| 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| 大型av网站在线播放| 美女国产高潮福利片在线看| 俄罗斯特黄特色一大片| 又黄又粗又硬又大视频| 国产精品日韩av在线免费观看 | 99re在线观看精品视频| 极品人妻少妇av视频| av中文乱码字幕在线| 久99久视频精品免费| 婷婷精品国产亚洲av在线| 99精品久久久久人妻精品| 天天躁夜夜躁狠狠躁躁| av福利片在线| ponron亚洲| 高清毛片免费观看视频网站| 亚洲精品在线美女| 日韩大码丰满熟妇| 身体一侧抽搐| 日韩大码丰满熟妇| 国产成+人综合+亚洲专区| 国产91精品成人一区二区三区| 99久久综合精品五月天人人| 亚洲中文字幕日韩| 亚洲精华国产精华精| 欧美乱色亚洲激情| 91成年电影在线观看| 日日夜夜操网爽| 999精品在线视频| 久久久久久人人人人人| 国产av又大| 久久久久久人人人人人| 给我免费播放毛片高清在线观看| 精品少妇一区二区三区视频日本电影| 久热爱精品视频在线9| 国产高清有码在线观看视频 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品一区二区三区四区久久 | 久久香蕉激情| 亚洲伊人色综图| 久99久视频精品免费| 黑人巨大精品欧美一区二区蜜桃| 99精品久久久久人妻精品| 乱人伦中国视频| 欧美一级毛片孕妇| 人人妻,人人澡人人爽秒播| 久久久国产成人精品二区| 国产三级在线视频| 亚洲av电影不卡..在线观看| 亚洲国产精品sss在线观看| www.精华液| 精品久久久久久,| 亚洲av电影不卡..在线观看| 美女大奶头视频| 久久伊人香网站| 巨乳人妻的诱惑在线观看| 两个人免费观看高清视频| 午夜福利成人在线免费观看| 亚洲全国av大片| 三级毛片av免费| 国产在线观看jvid| 成年女人毛片免费观看观看9| 午夜a级毛片| 精品欧美一区二区三区在线| 久久精品国产亚洲av香蕉五月| 日本在线视频免费播放| 中文字幕人妻丝袜一区二区| 夜夜躁狠狠躁天天躁| videosex国产| 亚洲国产欧美网| 亚洲一码二码三码区别大吗| 欧美日韩瑟瑟在线播放| 嫁个100分男人电影在线观看| 精品一区二区三区视频在线观看免费| 18禁黄网站禁片午夜丰满| 久久人人97超碰香蕉20202| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三区在线| 成人国语在线视频| 亚洲伊人色综图|