漢語連續語音識別系統的制作方法

            文檔序號:2836098閱讀:312來源:國知局
            專利名稱:漢語連續語音識別系統的制作方法
            技術領域
            本發明涉及一種漢語連續語音識別系統,尤其是一種采用基于段長分布的隱含馬爾可夫模型的漢語連續語音識別系統。
            背景技術
            語音識別就是從人的語音波形中提取信息并確定其語言含義的過程。語音識別作為一門交叉學科,涉及到了信號處理、統計模式識別、人工智能、計算機科學、語言學和認知科學等眾多學科。語音識別技術大致可以分為以下幾類語音信號檢測和特征提取技術,語音信號發生機理和數學模型建立,語音識別單元選擇,訓練和識別搜索算法,語言模型處理技術,語法和語義分析技術等等。
            隱含馬爾可夫模型(Hidden Markov Model-HMM)是目前非特定人語音識別的最重要的方法,HMM模型是語音產生機理的一種符合自然的描述,早在1913年A.A.Markov就提出了將Markov模型用于字符串的匹配,HMM數學理論方面的工作主要由Baum等人所完善,70年代中期Baker和Jelinik等人首次把HMM用于語音識別,當時他們研究的是連續語音的識別問題。隨著語音識別的研究由連續語音轉向孤立詞和70年代末80年代初矢量量化(Vector quantification-VQ)技術的出現,Rabiner和Levinson等人把HMM用于與人無關的孤立詞識別獲得了成功,隨著語音識別研究工作的深入開展,HMM語音識別方法愈來愈受到人們的重視,但與此同時,人們也愈來愈認識到經典HMM語音識別模型在一些重要方面存在嚴重缺陷,這就是第一、經典HMM是一個齊次的Markov模型,狀態轉移概率與狀態駐留長度無關,與語音的實際過程不符;第二、經典HMM模型用于大詞匯表的識別系統時,其模型的訓練量是災難性的;第三、無法自適應,從而對于進一步提高系統的識別率將是困難的;第四、模型的存儲量太大。
            由于上述方面的缺陷,導致采用經典HMM模型的語音識別系統的識別率低下,同時因為系統運算量和存儲量巨大使得系統工作效率不高,自我學習能力降低,系統開發成本偏高,難以普及到廣泛的應用當中去。

            發明內容
            本發明的目的在于提出了一種采用基于段長分布的隱含馬爾可夫模型的漢語連續語音識別系統,該語音識別系統比現有語音識別系統的識別率有顯著提高,更重要的是該語音識別系統的訓練算法比經典HMM的Baum-Welch訓練算法要快近二個數量級,存儲量也小得多,解決了所謂的訓練災難,同時提高了系統學習能力,從而可以在使用過程當中不斷地提高語音識別率。
            為達到上述目的,本發明是這樣實現的本發明包括語音采集裝置、前端處理模塊、特征提取模塊、基音特征提取模塊、聲調識別模塊、訓練模塊、聲學層識別模塊、拼音文法理解模塊、語言理解模塊九個部分,該聲學層識別模塊采用以狀態的段長概率來表示的隱含馬爾可夫模型,即基于段長分布的隱含馬爾可夫模型,簡稱DDBHMM模型,該模型的段長服從具有凸性的分布,其中,該模型中的狀態與語音中的音素或音節等語音單元相對應,而這些語音單元讀音的語音信號特征作為對應語音單元的觀測量;該訓練模塊的訓練方法步驟如下對訓練模塊輸入一個或多個包含有若干特定詞的讀音樣本文件,該文件中的每句話的一系列幀的特征向量以幀為單位,通過訓練模塊中的訓練搜索算法,對每幀語音信號中的每個詞的每個狀態進行搜索比較,得到語音信號特征矢量的矢量量化碼本(VQ碼本)和該特定詞的DDBHMM模型參數,輸入給聲學層識別模塊;所述的漢語連續語音識別系統的語音識別方法步驟如下通過語音采集裝置接收人們的語音信號,對于輸入的語音信號進行前端處理,并進行MFCC語音特征(基于Mel倒譜系數的語音特征,Mel-Frequency CepstralCoefficients)序列的提取,得到的該MFCC語音特征序列被送入聲學層識別模塊,通過聲學層識別模塊的搜索算法,產生拼音格形式的識別結果,同時,語音信號的基音特征矢量也被提取出來,送入聲調識別模塊,聲調識別模塊利用基音特征信息和拼音的分割點信息,得到拼音的聲調信息并加入到拼音格中,接著,通過拼音文法理解模塊對拼音格進行修剪,精簡后的音節格被送入語言理解模塊,轉化為拼音圖和詞圖,并在詞圖中進行搜索,得到最后的理解結果。
            對語音信號的前端處理包括先進行語音信號檢測,濾去確認為靜音或雜聲的信號,然后進行包括高頻預加重、加窗、分幀的預處理操作。
            該基于段長分布的隱含馬爾可夫模型是一個六狀態的模型。
            該基于段長分布的隱含馬爾可夫模型是一個段間不相關的模型。
            所述的具有凸性的分布是均勻分布或者高斯分布。
            所說的VQ碼本數量為857,音節總數為1254。
            訓練模塊中的訓練搜索算法步驟如下首先獲得該狀態對應的初始碼本,計算出該特征向量對狀態的概率,接著計算上一時刻該狀態的路徑在原地駐留的似然距離,然后對于新加入的此幀信號的特征向量進行最佳路徑的匹配,即找到它到原駐地的最佳狀態路徑,此后進入下一幀的搜索,直到一句話的訓練結束,然后統計和修改初始碼本的DDBHMM模型參數,完成本句話的訓練和碼本參數統計后,接著再進行輸入文件的下一句話的訓練和碼本參數統計,累加各句話訓練后得到的概率,以得到當前文件經訓練得到的概率,并將本次訓練過的每個文件得到的概率累加到本次訓練得到的總概率中,判斷本次訓練與上一次訓練后得到的總概率是否相等,不相等則進入下一次的循環訓練,相等則訓練結束。
            計算上一時刻該狀態的路徑在原地駐留的似然距離的方法是,對該狀態上一時刻到原駐地的最佳路徑隊列中的每個單元,計算相對于本碼本的狀態駐留概率和狀態轉移概率,然后累加即可。
            訓練時對于新加入的此幀信號的特征向量進行最佳路徑的匹配的方法分為兩種情況一是該特征向量所對應的幀所處的狀態與上一個特征向量所對應的幀所處的狀態屬于同一個字,處理的方法為復制前一時刻本狀態到原駐地的最佳路徑,并計算這條路徑到原駐留地的總的概率,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中,具體加入的操作過程為首先看此幸存狀態路徑隊列是否為空,若是,則令幸存路徑雙向列表的指針都指向此路徑;否則,則將此新加的路徑加到幸存路徑隊列的隊尾,接著對似然距離進行排序,比較各條在不同的時刻進入此狀態的路徑,若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該先進入本狀態的路徑被刪除,按此方法沿著幸存路徑隊列向前推進,直到隊頭;二是該特征向量所對應的幀處在某一個字的第一個狀態,這又分為兩種情況此狀態對應的字不是第一個字,則也應該復制以前的最佳路徑,分配新的路徑,鏈接最佳路徑并計算出該狀態累加似然距離,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;此狀態對應的是第一個字,則直接分配路徑,記錄進入當前路徑的時刻,計算該路徑在原駐地的似然距離,然后將此路徑加入到當前狀態的幸存路徑的雙向列表中。
            所說的初始碼本的參數的統計,其步驟為在統計的時候是每個狀態對應一個碼本,對于給定的本句話的一系列特征向量,從這句話對應的幸存路徑的最后一個字開始,對每個字的每個狀態分別進行統計,包括獲取每個狀態的段長、計算每個狀態的持續時間,然后統計出該字對應碼本的段長的均值和方差,用與當前的狀態相對應的特征向量的值統計出該狀態的碼本各參數的均值,對于方差的統計,根據事先定義的標志來確定,這個標志表明系統對速度和精度的偏好,當它的值為1時,表明偏好速度,只統計碼本的方差;否則,表明偏好精度,則要計算各個特征之間的協方差矩陣,做完上述統計以后,則進入下一個字的各狀態的碼本的特征的統計,直到本句話統計完成。
            聲學層識別模塊工作方式如下首先調入經過訓練得到的碼本及其DDBHMM模型參數,對于輸入的MFCC語音特征向量,獲得其幀數,計算出該特征向量對所有碼本的概率,對每一幀輸入連續語音識別器進行幀同步識別,在使用幀同步識別方法得到一系列的路徑以后,獲取當前幀的最佳字串路徑;累加每句話的路徑概率,得到整句話的累加概率;直到整句話結束,通過回溯得到最優路徑以及路徑上各個音節的候選,并保存為識別結果。
            所述的幀同步識別方法為,在此語音識別器中,對每一幀到來的特征向量采用搜索算法進行識別,對于到來的這幀信號,根據位置的不同分為兩種情況一、此特征向量是某一句話中的第一幀信號,則它對應的狀態對任何字而言都是第一個狀態,因此計算出該特征向量對所有的候選字的第一個狀態對應的碼本的概率,為此特征向量分配一條新的狀態路徑,然后設置進入當前路徑的時刻為當前幀的標號、設置當前路徑的概率為該特征向量對所有的候選字的第一個狀態對應的碼本的概率、記錄當前路徑的字的標號,將段長值設為1,計算出當前路徑的狀態駐留概率,將兩個概率,即該特征向量對所有的候選字的第一個狀態對應的碼本的概率和當前路徑的狀態駐留概率求和,即得到該路徑的完全概率,最后將此路徑加到當前狀態的幸存路徑的雙向列表中,返回;二、此特征向量對應的不是一句話的第一幀信號,處理的方法為將最優狀態路徑的指針設為空,若此次識別是連續語音的識別,則需要考慮當前幀產生的最優字輸出結果,搜索出最佳的狀態路徑,搜索方法為從每個字對應的最后一個狀態開始搜索,找到該字的最佳幸存狀態路徑隊列的頭指針,若頭指針為空,則表示沒有與該字對應的最佳狀態路徑序列,直接返回;若不為空,則比較各個字的識別器的輸出概率,得到最大的那一個,然后記錄字間跳轉的位置,即當前幀的標號,返回最佳字的狀態路徑隊列的頭指針,這就是搜索到的最佳路徑,在得到最佳的字的狀態路徑以后,若不為空,則復制該最優路徑,若該路徑對應的字不是第一個字,則要增加對前導字的引用,并使最優路徑指針指向該新分配的狀態路徑;處理完每一幀的特征向量,則對當前所有可能的路徑進行剪枝。
            當對情況二,即此特征向量對應的不是一句話的第一幀信號,對于狀態間的跳轉,分為兩種情況,即字內狀態跳轉和字間狀態跳轉,處理的方法分別為對字內狀態跳轉的情況,則為當前狀態分配一個新的路經,具體方法是取得前一狀態路徑,為本幀特征向量創建新的路徑并復制前一狀態路徑、設置進入當前路徑的時刻,計算新的狀態路徑的完全概率,最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;字間跳轉的情況字間的跳轉只有在連續語音識別時才存在,對孤立字的識別時,此時不用執行任何操作;對連續語音識別的情況,不用獲取前一狀態的最佳路徑而直接創建新的路徑并鏈接該路徑、增加前一字的最佳路徑的引用,計算出新的狀態路徑到原駐留地的完全概率,即狀態駐留概率和路徑概率之和,最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中。
            獲取當前幀的最佳字串路徑的方法為對每一句話產生的各條幸存的路徑,對所有的字逐一比較找出其中似然距離最大的那一個,就得到最佳字串狀態路徑。
            對當前所有可能的路徑進行剪枝的步驟為從當前狀態的狀態幸存路徑隊列的隊尾開始,對似然距離進行排序,若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該先進入本狀態的路徑被刪除,采用此種方法向前推進,刪除其中所有不可能成為最佳的路徑,直到幸存狀態路徑隊列的隊頭。
            在本發明中,搜索最優路徑的過程是一個幀同步算法,在對每一幀進行搜索時,進行剪枝和排序,最后獲得最優路徑。在連續語音識別中采用此方法進行最優路徑搜索的計算量可以比全搜索的方法下降2個數量級以上,存儲量也大大減少。對本發明的非特定人連續語音實驗結果表明,本發明的DDBHMM語音識別模型比經典的HMM模型的誤識率降低了12.9%,展示了DDBHMM的良好的性能;另外,本發明提出的訓練模塊很容易自適應,方法是由系統使用人讀預先設計好的文本,以其特征修改HMM模型參數后再進行識別。這種自適應能力還體現在其實時性,在使用時并不一定需要專門念自適應文本,而是在使用過程中自動進行的,這就使得系統在使用中能自我學習、自我完善。因此,本發明與現有漢語連續語音識別技術相比,速度快,存儲量小是本發明的極大優點,因此不但識別率有顯著提高,而且還降低了本發明的開發成本,完全可以在一般的微機(如PC機)上應用。


            圖1為基于段長分布的隱含馬爾可夫模型示意圖;圖2為基于段長分布的隱含馬爾可夫模型狀態序列示意圖;圖3為本發明實施例的總體流程框圖;圖4為圖1所示實施例的訓練模塊流程框圖;圖5為圖4所示“訓練模塊中對每個文件進行訓練”模塊的進一步描述流程框圖;圖6為圖5所示“依次對文件中的每句話進行訓練”模塊的進一步描述流程框圖;圖7為圖6所示“統計碼本的參數”模塊的進一步描述流程框圖;圖8為圖6所示“對相應的語音段進行訓練”模塊的進一步描述流程框圖;圖9為圖1所示實施例的聲學層模塊流程框圖;圖10為圖9所示的聲學層模塊中“幀同步語音識別”模塊的進一步描述流程框圖;圖11為圖10所示的流程框圖中“對每個字進行幀同步識別”模塊的進一步描述流程框圖。
            具體實施例方式
            以下結合附圖和具體實施例對本發明做進一步的闡述如圖1、2所示,本發明提出了一種用段長概率來表示的HMM模型(隱含馬爾可夫模型),顯然這是一個非齊次馬爾可夫模型。將這個非齊次馬爾可夫模型用于語音識別,并且把模型中的狀態與語音中的音素或音節等語音單元對應起來,而把這些語音單元讀音的語音信號特征作為對應語音單元的觀測量,我們就得到了一個基于段長分布的HMM模型(Duration DistributionBased Hidden Markov Models),簡稱為DDBHMM。在DDBHMM中只需要兩個參數λ=(B,D)。這里矩陣B是特征觀測概率矩陣B=[b1(o),b2(o),…bN(o)],矩陣D是狀態的段長概率矩陣D=[d1(τ),d2(τ),…dN(τ)],這里τ表示段長,服從概率分布di(τ)(0≤τ≤Dmax,i=1,2…N),N是模型狀態的個數。
            為了計算DDBHMM模型產生觀測矢量O=[o1,o2,…,oT]的概率,我們設S=[s1,s2,…,sT]為觀測矢量O所對應的系統的狀態序列,Ti為系統在狀態i的段長,用ti=Σk=liTk(1≤i≤N)]]>表示狀態序列s1…sT的分割點(t0=0,tN=T)。于是對于給定的DDBHMM模型λ有p(O,S/λ)=p(s1···sT)·b(o1···oT/s1···sT)]]>=a1,1(1)Πi=lNΠk=2Tiai,i(k)·ai,i+1(Ti+1)·Πi=lNb(oTi-1+1···oTi/o1···oTi-1)]]>=Πk=lNPk(Tk/T1···Tk-1)b(oTk-1+1···oTk/o1···oTk-l)---(1)]]>這是一個相當一般化的描述它既包容了狀態段長之間的相關,也包含了幀間特征之間的相關,而且段長分布函數Pk(τ)(k=1,2,…,N)可以取任何形式。因此DDBHMM模型具備了很強的對語音的描述能力。
            假設觀測矢量是段(狀態)間不相關的,則觀測矢量概率就可以寫成P(O,S/λ)=Πk=lNPk(Tk)bk(otk-1+1···otk)---(2)]]>DDBHMM模型是非齊次的HMM模型,由于非齊次馬爾可夫鏈是有后效過程,不能使用諸如Viterbi、Baum-Welch等基于Bellman動態規劃理論的解碼算法和參數重估算法,需要有完全新的、高效的模型訓練和解碼算法。DDBHMM模型應用的關鍵是構造高效的最優分割算法,即在某一模型λ下快速地找到一種分割(t1,t2,…,tN)使函數式(2)取值最大。函數式(3)所示的算法是一種基于最優狀態路徑的識別算法。在模型訓練時,則是對指定的訓練模型λ求出最佳的狀態路徑,然后再對模型參數進行分段訓練。
            λ^=argmax(λ;S)P(O,S/λ)]]>=atgmax(λ;t1,t2,···,tN)Πk=lNPk(Tk)bk(otk-1+1···otk)]]>=argmaxλmax(t1,t2,···,tN)Πk=lNPk(Tk)bk(otk-1+1···otk)---(3)]]>為了獲得函數式(3)中的最佳分割點t1,t2,…,tN,可以構造如下的模型訓練和識別算法。
            定義Lk(t,n)為k時刻由狀態(n-1)進入狀態n的路徑在t時刻(t≥k)的累積似然比值。
            Lk(t,n)=Σi=0n-1{lnPi(Ti)+Σm=ti-l+1tilnbi(Oti)}+Σm=ktlnbn(Om),]]>I(t,n)=Lk(t,n)+lnPn(t-tn-1)=Lk(t,n)+lnPn(t-k+1)]]>n=1,···N;t=1,···T---(4)]]>則λ^=arg{maxλmax(t1,···,tN)I(T,N)}---(5)]]>為了求解(3)式的最優路徑問題,需求出每一時刻進入每個狀態的最佳路徑。可以證明引理1設lgPn(τ)是定義域D(即τ∈D)中的上凸函數,若j≤k≤t,Lj(t,n)+lnPn(t+1-j)≤Lk(t,n)+lnPn(t+1-k)(6)則t’>t有Lj(t’,n)+lnPn(t’+1-j)≤Lk(t’,n)+lnPn(t’+1-k)(7)引理1成立的唯一條件是關于段長概率分布的凸性假設,它包含了一大類有用的分布(例如均勻分布、三角分布、指數分布、正態分布、Reilygh分布、余弦分布等均滿足引理假設),從而這一條件在實際應用中總是成立的,所以引理1的假設不會成為算法應用的限制。
            引理1表明,如果進入狀態n的某條路徑在(t+1)時刻向(n+1)狀態轉移時不比其它更晚進入的路徑占優,則在以后任何時刻(t’+1)(t’>t)從n狀態向(n+1)狀態轉移時都不可能成為占優者,從而在t時刻就可以把它刪除而不會影響算法的最優性。在搜索算法中,我們通過只保留幸存路徑,來保證最優的路徑不被錯誤地剪枝掉。幸存路徑的定義如下對于任意的k,當j<k≤t,在j時刻進入狀態n的路徑在t(t>j)時刻均有Lj(t,n)+lgPn(t+1-j)>Lk(t,n)+lgPn(t+1-k) (8)則稱該路徑是幸存路徑。
            根據引理1,如果我們在每一幀時刻,對狀態路徑都進行上式的比較并且只保留每個狀態下的幸存路徑,則可以大大減少最優搜索的路徑數目。據此容易推出引理2引理2設t時刻狀態n有m條幸存路徑,對應的進入n狀態的時間為tn(0)<tn(1)<…<tn(m-1),累積似然比為Ltn(0)(t,n),Ltn(1)(t,n),···,Ltn(m-1)(t,n),]]>則有Ltn(j)(t,n)+lgPn(t+1-tn(j))>Ltn(j+1)(t,n)++lgPn(t+1-tn(j+1)),]]>j=0,1,…,m-2 (9)這就是說,越早進入的幸存路徑在下一個時刻的轉移中具有越大的優先權,從而在(t+1)時刻由狀態n向狀態(n+1)轉移的最佳路徑總是第0條,對應的累加似然比值為Ltn(0)(t,n)。引理2保證了最佳路徑總是最早進入當前狀態的幸存路徑。因此可以推出如下定理定理1在T時刻滿足式(3)的最佳路徑是T時刻狀態N所對應的所有幸存路徑中最早進入狀態N的路徑。
            在DDBHMM模型中,搜索最優路徑的過程是一個幀同步算法。在對每一幀進行搜索時,用引理1和引理2進行剪枝和排序,最后由定理1保證最優路徑的獲得。在連續語音識別中采用此方法進行最優路徑搜索的計算量可以比全搜索的方法下降2個數量級以上。
            圖3是本發明的一個實施例,它包括有語音采集裝置、前端處理模塊、特征提取模塊、基音特征提取模塊、聲調識別模塊、訓練模塊、聲學層識別模塊、拼音文法理解模塊、語言理解模塊九個部分,該訓練模塊和聲學層識別模塊采用基于段長分布的隱含馬爾可夫模型(DDBHMM模型),該段長服從高斯分布,并且段間不相關。其中,該模型中的狀態與語音中的音素或音節等語音單元相對應,并且每個字的狀態為六個,而這些語音單元讀音的語音信號特征作為對應語音單元的觀測量。對訓練模塊輸入一個或多個包含有若干特定詞的讀音樣本文件,該文件中的每句話的一系列幀的特征向量以幀為單位,通過訓練模塊中的訓練搜索算法,對每幀語音信號中的每個詞的每個狀態進行搜索比較,得到一個語音信號特征矢量的VQ碼本和該特定詞的DDBHMM模型參數,輸入給聲學層識別模塊;本發明的語音識別方法步驟如下通過語音采集裝置接收人們的語音信號,對于輸入的語音信號進行前端處理,并進行MFCC語音特征(基于Mel倒譜系數的語音特征,Mel-Frequency Cepstral Coefficients)序列的提取,得到的該MFCC語音特征序列被送入聲學層識別模塊,通過聲學層識別模塊的搜索算法,產生拼音格形式的識別結果,同時,語音信號的基音特征矢量也被提取出來,送入聲調識別模塊,聲調識別模塊利用基音特征信息和拼音的分割點信息,得到拼音的聲調信息并加入到拼音格中,接著,通過拼音文法理解模塊對拼音格進行修剪,精簡后的音節格被送入語言理解模塊,轉化為拼音圖和詞圖,并在詞圖中進行搜索,得到最后的理解結果。
            上述訓練模塊中的訓練搜索算法是這樣實現的,如圖4至圖8所示首先,初始化訓練模塊,包括打開系統要使用的批處理命令文件;讀取初始語音訓練碼本文件名;讀取輸出語音訓練碼本的文件名;讀取段長加權系數,并設置段長加權系數;從批處理命令文件中讀取本程序使用的訓練的文件數量,此值應該在0到200之間,并檢驗訓練文件的數量是否符合要求;從上述批處理命令文件中讀取用來訓練的語音文件、語音描述文件和語音端點文件的文件名,到此,批處理命令文件的命令參數讀取完畢,關閉該文件;打開語音描述文件,并為每個語音描述文件在內存中分配相應的內存鏡像,具體的步驟為打開文件后,獲取文件的長度,然后在內存中分配相應大小的緩沖區,最后移動文件指針到數據的起始地點,讀入文件中的數據到內存緩沖區中,關閉該文件。進入下一個文件數據到內存緩沖區鏡像的建立,直到所有文件執行完畢;顯示訓練文件的數量和文件名;顯示初始語音訓練碼本文件名和輸出語音碼本文件名;設置保存上一輪訓練結束時的概率的變量為0;加載原始訓練碼本,即按照碼本數量等于857和單詞總數等于1254生成碼本對象,創建對象的過程中,完成的操作包括設置碼本的數量、設置詞的總數、在內存中分配碼本的緩沖區和碼本的臨時緩沖區、分配保存碼本的段長值的數組、分配存儲每個特征對碼本概率的數組、分配記錄每個狀態所對應的碼本號的二維數組和碼本所對應的特征向量的數量的數組(各碼本對應的特征向量對此數組中相應于該碼本的值的平均即為碼本的各參數);從記錄狀態到碼本的映射的文件中讀入每個狀態對應的在該文件中的碼本號,存儲到二維數組中;從初始化語音訓練文件中讀入各碼本參數的初始值;設置段長加權系數;其次,創建DDBHMM語音訓練器,在創建過程中,完成的操作包括設置使用的碼本的指針,設置用來訓練的詞的數量和詞表中每部分的詞的數量的指針。同時初始化幸存路徑的雙向隊列的指針為空和本狀態中幸存路徑的數量(為0值);然后,開始循環訓練,次數最多為200次,訓練時對每個文件一次進行,具體的步驟如下獲取每一輪訓練開始的時刻,將其顯示出來,并保存在文件中;設置每一輪訓練的概率Probability的初值,為0;清除碼本緩沖區,開始對每個語音文件進行分別的訓練,步驟為獲取當前數據文件在內存鏡像中的地址;計算與該文件對應的語音描述文件中拼音的數量,讀取的方法是用文件的有用的數據區的大小除以每個結構體的大小;創建特征管理器,用于特征讀取管理,創建的過程中,完成的操作為在內存中完成語音特征文件的內存鏡像,同時計算語音特征文件中存儲的語音的個數,并創建語音特征緩沖區。創建語音端點文件的內存鏡像,同時獲取語音端點文件的內容;獲取當前文件中語音的數量,并檢驗拼音的數量和語音的數量是否相等,不相等則返回,相等則繼續執行;設置每個文件訓練的初始概率為0;依次對文件中的每一個語音(每句話)進行訓練,訓練的方法為獲取與當前這句話所對應的特征向量的幀數;設置訓練器中訓練的碼本號參數,其中有詞表索引中詞的數量,內存鏡像中當前語音存儲區的地址;初始化各個字的狀態,其設置的值與創建語音存儲訓練器的類的構造函數相同;對相應的語音段進行訓練,并獲得該句話的經訓練后得到的概率。訓練中,對傳入的與某句話相應的一系列幀的特征向量進行訓練。在此過程中,以幀為單位進行搜索比較。具體的方法為
            對此段語音中的每個詞的每個狀態進行搜索比較。在獲取當前狀態對應的碼本號后,計算出該特征向量對狀態的概率StateProb。
            計算上一時刻本狀態的路徑在原地駐留的似然距離,計算方法為找到本狀態上一時刻到原駐地的最佳路徑的隊列的尾指針,對隊列中的每個單元計算相對于本碼本的狀態駐留概率CurTraceDurationProb和狀態轉移概率PathProb然后累加即可。
            訓練時對于新加入的此幀信號的特征向量進行最佳路徑的匹配,即找到它到原駐地的最佳狀態路徑。分為兩種情況此特征向量所對應的幀所處的狀態與上一個特征向量所對應的幀所處的狀態屬于同一個字,即此時狀態號HmmStateNo>0。處理的方法為復制前一時刻本狀態到原駐地的最佳路徑,并計算這條路徑到原駐留地的總的概率,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中。具體加入的操作過程為首先看此幸存狀態路徑隊列是否為空,若是,則令幸存路徑雙向列表的指針都指向此路徑;否則,則將此新加的路徑加到幸存路徑隊列的隊尾。接著對似然距離進行排序,考慮到引理1的內容,比較各條在不同的時刻進入此狀態的路徑,若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該先進入本狀態的路徑被刪除,按此方法沿著幸存路徑隊列向前推進,直到到隊頭;此特征向量所對應的幀處在某一個字的第一個狀態即此時狀態號HmmStateNo=0。這又分為兩種情況第一、此狀態對應的字不是第一個字,則也應該復制以前的最佳路徑,分配新的路徑,鏈接最佳路徑并計算出該狀態累加似然距離,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中。加入的方法同于上述加入方法。
            第二、此狀態對應的是第一個字,則直接分配路徑,記錄進入當前路徑的時刻,計算該路徑在原駐地的似然距離。然后將此路徑加入到當前狀態的幸存路徑的雙向列表中,加入的方法同上。至此,對當前幀的搜索結束,進入下一幀的搜索,直到一句話的訓練結束;累加各句話訓練后得到的概率WordProb,以得到整個文件經訓練得到的概率FileProbability;統計碼本的參數,統計的過程中,傳遞的參數為一整句話的所有的特征向量。在統計的時候則是每個狀態對應一個碼本。對于給定的一系列特征向量,首先找到這句話對應的幸存路徑的頭指針,然后從幸存路徑的最后一個字開始逐個的進行統計。統計的步驟為獲取該字的標號及進入某個狀態的時刻,對每個字的6個狀態分別進行統計。包括獲取碼本號和每個狀態的段長、計算每個狀態的持續時間。
            計算出對應碼本的臨時緩沖區地址,然后統計其段長的均值和方差。
            用與當前的狀態相對應的特征向量的值統計出該狀態的碼本各參數的均值。對于方差的統計,根據事先定義的標志DIAG_CODE_BOOK_FLAG來確定,這個標志表明系統對速度和精度的偏好,當它的值為1時,表明偏好系統計算速度,只統計碼本的方差;否則,表明偏好系統計算精度,則要計算各個特征之間的協方差矩陣。做完上述工作以后,則進入下一個字的各狀態的碼本的特征的統計,直到本句話統計完成。最后,返回本條路徑的到原駐地的似然距離,進入下一句話的訓練搜索和個狀態碼本參數的統計,到本文件中所有語音進行完畢;累加當前文件經訓練后得到的概率FileProbability到Probability中;釋放特征管理器占用的資源;顯示和保存該文件的循環選優的次數,文件編號,訓練后得到該文件的概率FileProbability和訓練過的文件得到的概率的總和Probability;修改碼本的參數和保存碼本到碼本輸出參數文件中;獲得此輪訓練結束的時刻,并顯示該輪訓練結束的時刻和保存該輪訓練結束的時刻到文件中;判斷兩次訓練后得到的總的概率是否相等,不相等則進入下一次的循環訓練,相等則訓練結束;最后,釋放碼本對象,訓練算法結束。
            本實施例中聲學層識別模塊是這樣實現的,如圖9至圖11所示首先是初始化聲學層識別模塊打開語音識別批處理命令文件,從該批處理命令文件中讀取詞的狀態的文件名;讀取碼本文件名;讀取待識別的文件的數量;檢驗待識別的語音文件的數量是否在0和本系統最大識別文件數之間,是則繼續;否則退出系統;獲取聲學層識別模塊中使用的識別文件名、識別拼音文件名、端點文件名和識別結果輸出文件名,到此,批處理命令文件讀取完畢,關閉批處理命令文件;顯示待識別的文件的數量,顯示識別文件的文件名、識別拼音文件的文件名、端點文件的文件名和識別結果輸出文件的文件名;打開用于將1254個有調音節轉化成408個無調音節的文件,然后讀取該文件的內容,確定有調音節到無調音節的轉化關系;讀取拼音文件的內容,并將其賦給字符串數組;初始化每個文件的字的識別率、音節的識別率、總的字的識別率和總的音節的識別率;設置識別字的數量的初值為0;加載語音識別碼本,初始化碼本,其碼本的數量和字的數量分別為857和1254;從字的狀態集文件中讀入每個字的每個狀態對應的碼本編號;從碼本文件中加載語音識別碼本,即從碼本參數文件中讀入每個碼本的參數,包括特征向量的均值、方差或協方差矩陣、段長的均值和方差;獲取段加權系數,該系數來自于碼本參數文件;顯示碼本文件名和段長加權系數的值;創建連續語音識別器,碼本為上面創建的碼本的指針指向的碼本;接著,在得到需要的相應的參數和建立好上述語音識別器后,對每個待識別文件逐一進行識別,對每個文件識別的步驟為創建與待識別文件組相對應的語音識別結果輸出文件;創建特征管理器,同時從待識別文件和端點文件中讀取相應的語音特征;獲取當前語音文件中的句子的數量,并將此句子的數量寫入保存結果的文件中;獲取該文件訓練開始的時刻,設置文件的初始概率為0,以一句話為單位開始進行識別,具體過程為獲取該句話中特征向量的幀數,初始化語音識別器,計算出該特征向量對所有碼本的概率,對每一幀進行幀同步識別,識別方法為,在此幀同步語音識別器中,對每一幀到來的特征向量采用搜索方法進行識別,對于到來的這幀信號,根據位置的不同分為兩種情況一、此特征向量是某一句話中的第一幀信號,則它對應的狀態對任何字而言都是第一個狀態,因此計算出此特征向量對所有的候選字的第一個狀態對應的碼本的概率StateProb。為此特征向量分配一條新的狀態路徑,然后設置進入當前路徑的時刻為當前幀的標號、設置路徑的概率PathProb為StateProb、記錄當前路徑的字的標號。將段長值設為1,計算出路徑的狀態駐留概率CurTraceDurationProb,將兩個概率StateProb和CurTraceDurationProb求和,即得到該路徑的完全概率。最后將此路徑加到當前狀態的幸存路徑的雙向列表中,返回;二、特征向量對應的不是一句話的第一幀信號,處理的方法為將最優狀態路徑的指針設為空,若此次識別是連續語音的識別,則需要考慮當前幀產生的最優字輸出結果。搜索出最佳的狀態路徑,搜索方法為從每個字對應的最后一個狀態開始搜索,找到該字的最佳幸存狀態路徑隊列的頭指針,若頭指針為空,則表示沒有與該字對應的最佳狀態路徑序列,直接返回;若不為空,則比較各個字的識別器的輸出概率,得到最大的那一個。然后記錄字間跳轉的位置,即當前幀的標號,返回最佳字的狀態路徑隊列的頭指針,這就是我們搜索到的最佳路徑。在得到最佳的字的狀態路徑以后,若不為空,則復制該最優路徑。若該路徑對應的字不是第一個字,則要增加對前導字的引用,并使最優路徑指針指向該新分配的狀態路徑。對于狀態間的跳轉,分為兩種情況,即字內狀態跳轉和字間狀態跳轉,處理的方法分別為對字內狀態跳轉即狀態號HmmStateNo>0的情況,則為當前狀態分配一個新的路經,具體方法是獲得前一狀態路徑的隊列的頭指針,為本幀特征向量創建新的路徑并復制前一狀態路徑、設置進入當前路徑的時刻,計算新的狀態路徑的完全概率,最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;字間跳轉即狀態號HmmStateNo=0的情況。因為字間的跳轉只有在連續語音識別時才存在,對孤立字的識別時,不用執行任何操作。對連續語音識別的情況,不用獲取前一狀態的最佳路徑而直接創建新的路徑并鏈接該路徑、增加前一字的最佳路徑的引用,計算出新的狀態路徑到原駐留地的完全概率(為狀態駐留概率和路徑概率之和),最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;處理完每一幀的特征向量,則對當前所有可能的路徑進行剪枝。剪枝函數執行的操作為在獲得當前狀態的狀態幸存路徑隊列的指針后,從此隊列的尾部開始,對似然距離進行排序,由前述的引理1可知若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該(先進入本狀態的)路徑被刪除,采用此種方法向前推進,刪除其中所有不可能成為最佳的路徑,直到幸存狀態路徑隊列的隊頭,達到剪枝的目的;在使用幀同步識別器得到一系列的路徑以后,獲取當前幀的最佳字串路徑,獲取的方法為對每一句話產生的各條幸存的路徑,由于它的最后一個字肯定是所有可能的候選字中的某一個,故對所有的字逐一比較找出其中似然距離最大的那一個,就得到我們所說的最佳狀態路徑;累加每句話的路徑概率PathProb,得到整句話的累加概率,直到整句話結束,通過回溯得到最優路徑以及路徑上各個音節的候選,并保存為識別結果。
            顯示這句話識別出的結果,并將此結果保存到存儲結果的文件中;最后,釋放碼本和語音識別訓練器,語音識別過程結束。
            以上所述的是本發明的一個精度優選實施例,如果系統對速度要求較高,可以將上述實施例中的段長分布由高斯分布改為均勻分布,這樣系統的計算量大大降低,而識別精度雖然有所降低,但還是比現有技術要優越。
            從上面的說明可以看出,本發明不但提出了具有國際領先水平的基于段長分布的HMM(DDBHMM)模型,而且給出了比經典HMM模型傳統算法更優越、而且具有實用性的訓練算法和識別算法。可以廣泛應用到計算機語音識別輸入、電話語音合成、連續語音聲控電話交換機、機器翻譯等技術項目和產品中去。
            權利要求
            1.一種漢語連續語音識別系統,其特征在于所述的漢語連續語音識別系統包括語音采集裝置、前端處理模塊、特征提取模塊、基音特征提取模塊、聲調識別模塊、訓練模塊、聲學層識別模塊、拼音文法理解模塊、語言理解模塊九個部分,該聲學層識別模塊采用以狀態的段長概率來表示的隱含馬爾可夫模型,即基于段長分布的隱含馬爾可夫模型,簡稱DDBHMM模型,該模型的段長服從具有凸性的分布,其中,該模型中的狀態與語音中的音素或音節等語音單元相對應,而這些語音單元讀音的語音信號特征作為對應語音單元的觀測量;該訓練模塊的訓練方法步驟如下對訓練模塊輸入一個或多個包含有若干特定詞的讀音樣本文件,該文件中的每句話的一系列幀的特征向量以幀為單位,通過訓練模塊中的訓練搜索算法,對每幀語音信號中的每個詞的每個狀態進行搜索比較,得到語音信號特征矢量的矢量量化碼本(VQ碼本)和該特定詞的DDBHMM模型參數,輸入給聲學層識別模塊;所述的漢語連續語音識別系統的語音識別方法步驟如下通過語音采集裝置接收人們的語音信號,對于輸入的語音信號進行前端處理,并進行MFCC語音特征(基于Mel倒譜系數的語音特征,Mel-Frequency CepstralCoefficients)序列的提取,得到的該MFCC語音特征序列被送入聲學層識別模塊,通過聲學層識別模塊的搜索算法,產生拼音格形式的識別結果,同時,語音信號的基音特征矢量也被提取出來,送入聲調識別模塊,聲調識別模塊利用基音特征信息和拼音的分割點信息,得到拼音的聲調信息并加入到拼音格中,接著,通過拼音文法理解模塊對拼音格進行修剪,精簡后的音節格被送入語言理解模塊,轉化為拼音圖和詞圖,并在詞圖中進行搜索,得到最后的理解結果。
            2.如權利要求1所述的漢語連續語音識別系統,其特征在于對語音信號的前端處理包括先進行語音信號檢測,濾去確認為靜音或雜聲的信號,然后進行包括高頻預加重、加窗、分幀的預處理操作。
            3.如權利要求1所述的漢語連續語音識別系統,其特征在于該基于段長分布的隱含馬爾可夫模型是一個六狀態的模型。
            4.如權利要求1所述的漢語連續語音識別系統,其特征在于該基于段長分布的隱含馬爾可夫模型是一個段間不相關的模型。
            5.如權利要求1所述的漢語連續語音識別系統,其特征在于所述的具有凸性的分布是均勻分布或者高斯分布。
            6.如權利要求1所述的漢語連續語音識別系統,其特征在于所說的VQ碼本數量為857,音節總數為1254。
            7.如權利要求1所述的漢語連續語音識別系統,其特征在于訓練模塊中的訓練搜索算法步驟如下首先獲得該狀態對應的初始碼本,計算出該特征向量對狀態的概率,接著計算上一時刻該狀態的路徑在原地駐留的似然距離,然后對于新加入的此幀信號的特征向量進行最佳路徑的匹配,即找到它到原駐地的最佳狀態路徑,此后進入下一幀的搜索,直到一句話的訓練結束,然后統計和修改初始碼本的DDBHMM模型參數,完成本句話的訓練和碼本參數統計后,接著再進行輸入文件的下一句話的訓練和碼本參數統計,累加各句話訓練后得到的概率,以得到當前文件經訓練得到的概率,并將本次訓練過的每個文件得到的概率累加到本次訓練得到的總概率中,判斷本次訓練與上一次訓練后得到的總概率是否相等,不相等則進入下一次的循環訓練,相等則訓練結束。
            8.如權利要求7所述的漢語連續語音識別系統,其特征在于計算上一時刻該狀態的路徑在原地駐留的似然距離的方法是,對該狀態上一時刻到原駐地的最佳路徑隊列中的每個單元,計算相對于本碼本的狀態駐留概率和狀態轉移概率,然后累加即可。
            9.如權利要求7所述的漢語連續語音識別系統,其特征在于訓練時對于新加入的此幀信號的特征向量進行最佳路徑的匹配的方法分為兩種情況一是該特征向量所對應的幀所處的狀態與上一個特征向量所對應的幀所處的狀態屬于同一個字,處理的方法為復制前一時刻本狀態到原駐地的最佳路徑,并計算這條路徑到原駐留地的總的概率,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中,具體加入的操作過程為首先看此幸存狀態路徑隊列是否為空,若是,則令幸存路徑雙向列表的指針都指向此路徑;否則,則將此新加的路徑加到幸存路徑隊列的隊尾,接著對似然距離進行排序,比較各條在不同的時刻進入此狀態的路徑,若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該先進入本狀態的路徑被刪除,按此方法沿著幸存路徑隊列向前推進,直到隊頭;二是該特征向量所對應的幀處在某一個字的第一個狀態,這又分為兩種情況此狀態對應的字不是第一個字,則也應該復制以前的最佳路徑,分配新的路徑,鏈接最佳路徑并計算出該狀態累加似然距離,然后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;此狀態對應的是第一個字,則直接分配路徑,記錄進入當前路徑的時刻,計算該路徑在原駐地的似然距離,然后將此路徑加入到當前狀態的幸存路徑的雙向列表中。
            10.如權利要求7所述的漢語連續語音識別系統,其特征在于所說的初始碼本的參數的統計,其步驟為在統計的時候是每個狀態對應一個碼本,對于給定的本句話的一系列特征向量,從這句話對應的幸存路徑的最后一個字開始,對每個字的每個狀態分別進行統計,包括獲取每個狀態的段長、計算每個狀態的持續時間,然后統計出該字對應碼本的段長的均值和方差,用與當前的狀態相對應的特征向量的值統計出該狀態的碼本各參數的均值,對于方差的統計,根據事先定義的標志來確定,這個標志表明系統對速度和精度的偏好,當它的值為1時,表明偏好速度,只統計碼本的方差;否則,表明偏好精度,則要計算各個特征之間的協方差矩陣,做完上述統計以后,則進入下一個字的各狀態的碼本的特征的統計,直到本句話統計完成。
            11.如權利要求1所述的漢語連續語音識別系統,其特征在于聲學層識別模塊工作方式如下首先調入經過訓練得到的碼本及其DDBHMM模型參數,對于輸入的MFCC語音特征向量,獲得其幀數,計算出該特征向量對所有碼本的概率,對每一幀輸入連續語音識別器進行幀同步識別,在使用幀同步識別方法得到一系列的路徑以后,獲取當前幀的最佳字串路徑;累加每句話的路徑概率,得到整句話的累加概率;直到整句話結束,通過回溯得到最優路徑以及路徑上各個音節的候選,并保存為識別結果。
            12.如權利要求11所述的漢語連續語音識別系統,其特征在于所述的幀同步識別方法為,在此語音識別器中,對每一幀到來的特征向量采用搜索算法進行識別,對于到來的這幀信號,根據位置的不同分為兩種情況一、此特征向量是某一句話中的第一幀信號,則它對應的狀態對任何字而言都是第一個狀態,因此計算出該特征向量對所有的候選字的第一個狀態對應的碼本的概率,為此特征向量分配一條新的狀態路徑,然后設置進入當前路徑的時刻為當前幀的標號、設置當前路徑的概率為該特征向量對所有的候選字的第一個狀態對應的碼本的概率、記錄當前路徑的字的標號,將段長值設為1,計算出當前路徑的狀態駐留概率,將兩個概率,即該特征向量對所有的候選字的第一個狀態對應的碼本的概率和當前路徑的狀態駐留概率求和,即得到該路徑的完全概率,最后將此路徑加到當前狀態的幸存路徑的雙向列表中,返回;二、此特征向量對應的不是一句話的第一幀信號,處理的方法為將最優狀態路徑的指針設為空,若此次識別是連續語音的識別,則需要考慮當前幀產生的最優字輸出結果,搜索出最佳的狀態路徑,搜索方法為從每個字對應的最后一個狀態開始搜索,找到該字的最佳幸存狀態路徑隊列的頭指針,若頭指針為空,則表示沒有與該字對應的最佳狀態路徑序列,直接返回;若不為空,則比較各個字的識別器的輸出概率,得到最大的那一個,然后記錄字間跳轉的位置,即當前幀的標號,返回最佳字的狀態路徑隊列的頭指針,這就是搜索到的最佳路徑,在得到最佳的字的狀態路徑以后,若不為空,則復制該最優路徑,若該路徑對應的字不是第一個字,則要增加對前導字的引用,并使最優路徑指針指向該新分配的狀態路徑;處理完每一幀的特征向量,則對當前所有可能的路徑進行剪枝。
            13.如權利要求12所述的漢語連續語音識別系統,其特征在于當對情況二,即此特征向量對應的不是一句話的第一幀信號,對于狀態間的跳轉,分為兩種情況,即字內狀態跳轉和字間狀態跳轉,處理的方法分別為對字內狀態跳轉的情況,則為當前狀態分配一個新的路經,具體方法是取得前一狀態路徑,為本幀特征向量創建新的路徑并復制前一狀態路徑、設置進入當前路徑的時刻,計算新的狀態路徑的完全概率,最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中;字間跳轉的情況字間的跳轉只有在連續語音識別時才存在,對孤立字的識別時,此時不用執行任何操作;對連續語音識別的情況,不用獲取前一狀態的最佳路徑而直接創建新的路徑并鏈接該路徑、增加前一字的最佳路徑的引用,計算出新的狀態路徑到原駐留地的完全概率,即狀態駐留概率和路徑概率之和,最后將新增加的路徑加入到當前狀態的幸存路徑的雙向列表中。
            14.如權利要求11所述的漢語連續語音識別系統,其特征在于獲取當前幀的最佳字串路徑的方法為對每一句話產生的各條幸存的路徑,對所有的字逐一比較找出其中似然距離最大的那一個,就得到最佳字串狀態路徑。
            15.如權利要求12所述的漢語連續語音識別系統,其特征在于對當前所有可能的路徑進行剪枝的步驟為從當前狀態的狀態幸存路徑隊列的隊尾開始,對似然距離進行排序,若先進入本狀態的路徑在概率上若不比后進入本狀態的路徑占優,則該先進入本狀態的路徑被刪除,采用此種方法向前推進,刪除其中所有不可能成為最佳的路徑,直到幸存狀態路徑隊列的隊頭。
            全文摘要
            本發明提供了一種采用基于段長分布的隱含馬爾可夫模型的漢語連續語音識別系統,包括語音采集裝置、前端處理模塊、特征提取模塊、基音特征提取模塊、聲調識別模塊、訓練模塊、聲學層識別模塊、拼音文法理解模塊、語言理解模塊;通過訓練模塊得到碼本和模型參數;通過聲學層識別模塊對MFCC語音特征序列的搜索算法,產生拼音格形式的識別結果,同時,語音信號的基音特征矢量也被提取出來,送入聲調識別模塊,聲調識別模塊利用基音特征信息和拼音的分割點信息,得到拼音的聲調信息并加入到拼音格中,接著,通過拼音文法理解模塊對拼音格進行修剪,送入語言理解模塊,轉化為拼音圖和詞圖,并在詞圖中進行搜索,得到最后的理解結果。
            文檔編號G10L15/00GK1499484SQ02146100
            公開日2004年5月26日 申請日期2002年11月6日 優先權日2002年11月6日
            發明者王作英, 吳及, 肖熙, 孫甲松, 王俠 申請人:北京天朗語音科技有限公司
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品