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

    一種改進(jìn)的LLL模糊度規(guī)約算法

    2017-12-02 03:02:02呂志平翟樹峰鄺英才王福林
    關(guān)鍵詞:規(guī)約尺度向量

    呂 浩,呂志平,翟樹峰,鄺英才,王福林

    (1.信息工程大學(xué) 地理空間信息學(xué)院,鄭州 450000;2.61287部隊(duì),昆明 650000)

    一種改進(jìn)的LLL模糊度規(guī)約算法

    呂 浩1,呂志平1,翟樹峰1,鄺英才1,王福林2

    (1.信息工程大學(xué) 地理空間信息學(xué)院,鄭州 450000;2.61287部隊(duì),昆明 650000)

    整周模糊度的高效解算是GNSS高精度數(shù)據(jù)處理中的關(guān)鍵,基于格論進(jìn)行GNSS模糊度估計(jì)時(shí)需要通過(guò)格基規(guī)約來(lái)實(shí)現(xiàn)最優(yōu)整周模糊度向量的快速搜索。針對(duì)高維情況下常規(guī)LLL規(guī)約算法輔助整周模糊度解算存在規(guī)約耗時(shí)較長(zhǎng)和規(guī)約性能有限的問題,引入最小列旋轉(zhuǎn)QR分解技術(shù)對(duì)基向量進(jìn)行預(yù)排序,采用延后尺度規(guī)約和部分尺度規(guī)約來(lái)減少規(guī)約過(guò)程中的冗余尺度規(guī)約,以改善LLL算法的執(zhí)行效果。分別通過(guò)模擬和實(shí)測(cè)數(shù)據(jù)進(jìn)行實(shí)驗(yàn),結(jié)果表明:改進(jìn)后的LLL算法可以明顯降低格基規(guī)約耗時(shí),實(shí)測(cè)環(huán)境下其規(guī)約效率相比于傳統(tǒng)方法提高了約 10倍,且能夠保證較好的規(guī)約性能,從而有效提升高維模糊度的解算效率。

    模糊度解算;格基規(guī)約;LLL規(guī)約;QR分解;尺度規(guī)約

    在GNSS精密定位數(shù)據(jù)處理中,相位整周模糊度的快速、精確解算是保證定位精度和可靠性的關(guān)鍵,也是GNSS研究領(lǐng)域中的熱點(diǎn)問題和難點(diǎn)問題。在眾多的模糊度解算方法中,基于模糊度域內(nèi)的LAMBDA算法應(yīng)用最為廣泛[1],該方法以整數(shù)最小二乘原理為基礎(chǔ),通過(guò)降低模糊度方差分量間的相關(guān)性來(lái)壓縮搜索空間,提高搜索效率。隨后,針對(duì)模糊度降相關(guān)問題的多種方法被相繼提出[2-5]。

    模糊度解算作為一個(gè)整數(shù)最小二乘問題,實(shí)際上等價(jià)于格理論中的最近向量問題(CVP)[6-7]。由于格的特殊性質(zhì),為保證高效、精確地求解CVP,首先要對(duì)格基進(jìn)行規(guī)約處理,其中以LLL規(guī)約算法最為流行[8-9]。文獻(xiàn)[10]最早將格基規(guī)約的思想引入到 GNSS整周模糊度解算當(dāng)中;文獻(xiàn)[11]從理論上證明了 LLL規(guī)約算法和LAMBDA降相關(guān)之間的等價(jià)性;文獻(xiàn)[12]指出單純的尺度規(guī)約過(guò)程對(duì)搜索空間和搜索效率沒有影響,并提出了部分尺度規(guī)約;文獻(xiàn)[13]通過(guò)改變尺度規(guī)約順序來(lái)降低LLL算法計(jì)算復(fù)雜度;文獻(xiàn)[14-15]分別采用最小列旋轉(zhuǎn)技術(shù)和貪心算法來(lái)減少規(guī)約過(guò)程中的基向量交換次數(shù)。近來(lái)的研究進(jìn)一步指出,為提高模糊度搜索效率,格基規(guī)約的本質(zhì)在于盡可能促使基向量按照一定方向排序[16-17]。已有的格基規(guī)約算法主要是圍繞基向量交換條件和尺度規(guī)約的有效性來(lái)尋求單一的規(guī)約計(jì)算復(fù)雜度的降低或規(guī)約性能的改善,而提高規(guī)約性能往往會(huì)導(dǎo)致計(jì)算復(fù)雜度的增大,因此需要對(duì)兩者進(jìn)行平衡來(lái)減少整體解算耗時(shí)。

    在解決高維模糊度參數(shù)時(shí),由于規(guī)約耗時(shí)通常遠(yuǎn)大于搜索耗時(shí),因而減少規(guī)約耗時(shí)對(duì)于提高模糊度整體解算效率影響顯著。鑒于此,本文在分析LLL規(guī)約算法的基礎(chǔ)上,對(duì)該算法進(jìn)行優(yōu)化改進(jìn),在降低規(guī)約復(fù)雜度的同時(shí)兼顧規(guī)約性能,并結(jié)合模擬和實(shí)測(cè)數(shù)據(jù)進(jìn)行了驗(yàn)證與分析。

    1 基本數(shù)學(xué)模型

    1.1 基于格的模糊度估計(jì)

    基于載波相位測(cè)量的GNSS定位模型是一種混合整數(shù)模型,通過(guò)經(jīng)典最小二乘可以得到模糊度的實(shí)數(shù)解及其方差協(xié)方差陣Q?a。考慮模糊度的整數(shù)約束,則對(duì)應(yīng)的整數(shù)最小二乘的估計(jì)準(zhǔn)則為:

    對(duì)Q?a進(jìn)行Cholesky分解:

    其中,B為上三角矩陣。則式(1)可以表示為:

    式中,y=B-T為常量。該式在格理論中被稱為最近向量問題,B為基于Q?a獲得的格基。

    于是,模糊度解算過(guò)程即可轉(zhuǎn)化為求解與目標(biāo)向量y最近的格上向量,常采用LLL規(guī)約算法來(lái)處理此類問題。

    1.2 LLL格基規(guī)約

    采用經(jīng)典 LLL規(guī)約算法對(duì)上述問題進(jìn)行處理時(shí)還需要對(duì)基矩陣B進(jìn)行分解:

    若矩陣B*和U中的元素滿足以下條件[8]:

    則稱B為L(zhǎng)LL規(guī)約基,其中第一項(xiàng)稱為尺度規(guī)約,第二項(xiàng)稱為基向量交換。

    在實(shí)際解算中,為了保證較好的數(shù)值穩(wěn)定性,常基于QR分解來(lái)實(shí)現(xiàn)LLL規(guī)約[18-19]。令基矩陣B=QR,其中,Q為正交矩陣,R為上三角矩陣,且對(duì)角線元素則有:

    于是,LLL算法的規(guī)約條件可以改寫為:

    為了滿足上述規(guī)約條件,通常需要構(gòu)造變換矩陣進(jìn)行規(guī)約操作。

    式中,[ ]int為取整操作,利用其右乘基矩陣B即可實(shí)現(xiàn)對(duì)應(yīng)元素的尺度規(guī)約,同時(shí)更新上三角矩陣R。

    2)基向量交換(Swap變換):若不滿足式(6)中條件2的要求,即時(shí),則構(gòu)造初等變換矩陣來(lái)調(diào)換的順序,并對(duì)上三角矩陣R進(jìn)行更新。

    基于QR分解的LLL規(guī)約同經(jīng)典LLL規(guī)約本質(zhì)上是相同的,但僅基于上三角陣進(jìn)行變換操作,較為簡(jiǎn)便,且擁有更高的浮點(diǎn)精度。

    2 LLL改進(jìn)算法

    根據(jù)LLL規(guī)約算法定義可知,其主要由尺度規(guī)約和基向量交換兩部分構(gòu)成,其中:尺度規(guī)約對(duì)模糊度搜索效率并無(wú)直接影響[14],而是為了促使基向量得以充分交換;LLL規(guī)約的核心正是通過(guò)基向量交換來(lái)實(shí)現(xiàn)模糊度快速搜索。因此,本文基于規(guī)約效率和規(guī)約性能的平衡這一目的,對(duì)LLL算法進(jìn)行改進(jìn)。

    2.1 延后尺度規(guī)約與部分尺度規(guī)約

    在引入新的尺度規(guī)約手段之前,首先對(duì)常規(guī)LLL算法中尺度規(guī)約的特點(diǎn)和有效性進(jìn)行分析。將LLL算法中的尺度規(guī)約分為兩類:第一類是次對(duì)角線元素的尺度規(guī)約;第二類是當(dāng)基向量交換條件(Swap條件)不滿足時(shí),對(duì)列向量中其余元素進(jìn)行尺度規(guī)約。這兩類尺度規(guī)約的計(jì)算過(guò)程如圖1所示。

    圖1 LLL算法尺度規(guī)約過(guò)程Fig.1 Process of LLL size reduction

    由圖1可知,LLL算法中尺度規(guī)約的執(zhí)行是根據(jù)是否滿足Swap條件,在第k層、第max(k-1, 2)層和第k+1層三個(gè)規(guī)約層之間循環(huán)往復(fù),因而相同的規(guī)約層由于基向量交換的影響,必然會(huì)導(dǎo)致部分第一類尺度規(guī)約和第二類尺度規(guī)約的重復(fù)執(zhí)行。

    1)延后尺度規(guī)約

    根據(jù)LLL算法尺度規(guī)約的特點(diǎn),由于同一規(guī)約層存在重復(fù)規(guī)約,所以只有當(dāng)規(guī)約后能夠進(jìn)行基向量交換時(shí),第一類規(guī)約才是必需的;除此以外的第一類和第二類規(guī)約均可在規(guī)約基B完成所有基向量交換后再執(zhí)行。

    將第一類尺度規(guī)約和基向量交換步驟進(jìn)行合并,此時(shí)的基向量交換條件如下:

    其中2維矩陣P表示如下:

    利用變換矩陣Zi對(duì)格基B和上三角矩陣R進(jìn)行更新,用Zi右乘R,同時(shí)消去R的對(duì)角線以下元素。

    該過(guò)程即為規(guī)約后基向量交換,在對(duì)B完成所有上述操作后可以將剩余的第一類和第二類尺度規(guī)約延遲到算法最后執(zhí)行。需要說(shuō)明的是,延后尺度規(guī)約僅消除LLL算法中基向量交換引起的冗余尺度規(guī)約,不改變基向量交換策略。

    2)部分尺度規(guī)約

    對(duì)于常規(guī)LLL算法,文獻(xiàn)[12]從幾何角度說(shuō)明了尺度規(guī)約不改變規(guī)約基的施密特正交化向量長(zhǎng)度,因而其對(duì)于搜索效率并不構(gòu)成影響??紤]到基向量交換對(duì)搜索效率改善的作用,一般僅對(duì)次對(duì)角線元素進(jìn)行尺度規(guī)約,亦即第一類尺度規(guī)約;但為了保證算法的有效性和穩(wěn)定性,還需要有選擇性地進(jìn)行第二類尺度規(guī)約。

    在延后尺度規(guī)約中,不影響基向量交換的第一類和第二類尺度規(guī)約被置于算法最后一步。因此,采用部分尺度規(guī)約,對(duì)于其中的第二類尺度規(guī)約只在滿足一定條件下進(jìn)行。

    部分尺度規(guī)約針對(duì)的是第二類尺度規(guī)約,為避免出現(xiàn)數(shù)值穩(wěn)定性問題而選擇部分執(zhí)行,在不影響基向量交換的基礎(chǔ)上進(jìn)一步減少了LLL算法的尺度規(guī)約數(shù)。

    2.2 最小列旋轉(zhuǎn)QR分解

    已有的研究表明,格基規(guī)約的目的是為了使基向量按照特定方向排序[17]。除了尺度規(guī)約的個(gè)數(shù),基向量交換的次數(shù)也會(huì)對(duì)格基規(guī)約的效率產(chǎn)生影響,且基向量的排序直接影響到規(guī)約基性能的好壞。因此,考慮在正式的格基規(guī)約前可以按照某種準(zhǔn)則對(duì)規(guī)約基向量進(jìn)行預(yù)排序,進(jìn)而減少規(guī)約過(guò)程中的基向量交換。

    通?;贖ouseholder變換實(shí)現(xiàn)QR分解,假設(shè)對(duì)格基B的前i-1個(gè)向量進(jìn)行 Householder變換后得到的矩陣為Ri-1。最小列旋轉(zhuǎn)QR分解是指在對(duì)格基B進(jìn)行第i步 QR分解之前,找到對(duì)應(yīng)當(dāng)前子矩陣中列向量長(zhǎng)度最小的那一列并構(gòu)造初等變換矩陣Ui交換Ri-1的第i列和第j列;然后再構(gòu)造Householder變換矩陣Hi進(jìn)行常規(guī)的第i步 QR分解,實(shí)現(xiàn)矩陣的上三角化最小列旋轉(zhuǎn)QR分解的整體變換過(guò)程可以表示為:

    顯然,在第i步最小列旋轉(zhuǎn)QR分解完成后,R中的元素滿足如下條件:

    根據(jù)Householder變換的性質(zhì),對(duì)矩陣Ri-1Ui進(jìn)行上三角化并不改變矩陣列向量長(zhǎng)度。式(12)表明,通過(guò)以上策略能夠保證每一步變換后得到的對(duì)角線元素ri,i小于子矩陣中其它列向量長(zhǎng)度。

    2.3 算法流程

    將上述基于LLL算法的改進(jìn)策略進(jìn)行整合,主要從尺度規(guī)約和基向量交換兩個(gè)方面對(duì)算法進(jìn)行優(yōu)化處理,具體流程見圖2。

    圖2 LLL改進(jìn)算法流程圖Fig.2 Flow chart of modified LLL algorithm

    如圖2所示,改進(jìn)后算法的處理步驟主要包括:基于模糊度方差協(xié)方差陣Q?a構(gòu)造格基B,對(duì)原始格基進(jìn)行最小列旋轉(zhuǎn)QR分解得到上三角陣R;設(shè)定初始的規(guī)約層k,根據(jù)式(8)的基向量交換條件進(jìn)行判斷,如果成立則進(jìn)行規(guī)約后基向量交換,即對(duì)應(yīng)式(9)的第一類尺度規(guī)約與基向量交換的合并操作,并更新規(guī)約層否則直接進(jìn)入規(guī)約層k=k+1;若規(guī)約層k≤n,則重復(fù)上步操作直到循環(huán)結(jié)束;最后在第一類尺度規(guī)約的基礎(chǔ)上,根據(jù)當(dāng)前次對(duì)角元素是否滿足進(jìn)行部分第二類尺度規(guī)約。

    3 算例分析

    文獻(xiàn)[7]提出的基于Householder變換的LLL規(guī)約(HLLL規(guī)約)其實(shí)就是基于QR分解的LLL規(guī)約;文獻(xiàn)[20]提出的 HELLL規(guī)約算法實(shí)際上是基于最小列旋轉(zhuǎn)和全部尺度規(guī)約的迭代操作。本文將基于延后尺度規(guī)約技術(shù)和QR分解的LLL規(guī)約稱為DLLL算法,將圖2中對(duì)應(yīng)的LLL改進(jìn)算法稱為DPLLL算法。

    為了評(píng)估 LLL改進(jìn)算法在模糊度解算中的應(yīng)用效果,分別基于模擬和實(shí)測(cè)數(shù)據(jù)結(jié)合不同指標(biāo)對(duì)LLL、HLLL、DLLL、HELLL和DPLLL算法進(jìn)行驗(yàn)證比較和統(tǒng)計(jì)分析。在模糊度解算過(guò)程中,模糊度的搜索部分統(tǒng)一采用目前廣泛應(yīng)用的SE-VB算法[5]。實(shí)驗(yàn)環(huán)境為:華碩PC機(jī)(Intel Core i7-6700 CPU,2.60GHz,內(nèi)存 8.0GB,64位 Windows10 操作系統(tǒng)),軟件MATLAB R2014a(8.3.0.532)。

    3.1 實(shí)驗(yàn)1

    利用文獻(xiàn)[2]的隨機(jī)模擬方法構(gòu)建模糊度協(xié)方差陣,具體構(gòu)造如下:

    式中:U為單位正交矩陣,通過(guò) Givens旋轉(zhuǎn)原理進(jìn)行構(gòu)造;Λ為對(duì)角陣,其元素為矩陣Q?a的特征值的升序排列。當(dāng)Q?a的維數(shù)大于20時(shí),條件數(shù)對(duì)數(shù)值取值范圍為[3, 4.5];低維情況下條件數(shù)對(duì)數(shù)值取值范圍為[0, 4.5]。條件數(shù)按照均勻分布在取值范圍內(nèi)隨機(jī)產(chǎn)生,特征值也按照均勻分布隨機(jī)產(chǎn)生,且最大值λmax與最小值λmin滿足與條件數(shù)c的關(guān)系:在實(shí)驗(yàn)1中,本文分兩組方案進(jìn)行模擬數(shù)據(jù)實(shí)驗(yàn)。

    1)第一組模擬實(shí)驗(yàn)

    由式(13)生成10、20、30、40維的仿真數(shù)據(jù),每一維構(gòu)造100組數(shù)據(jù),分別采用LLL、HLLL和DLLL算法對(duì)每一維數(shù)據(jù)進(jìn)行模糊度規(guī)約處理,并計(jì)算 100組數(shù)據(jù)的平均尺度規(guī)約個(gè)數(shù)和平均規(guī)約耗時(shí),結(jié)果如表1所示。由表1可以看出,隨著維數(shù)的增加,模糊度規(guī)約過(guò)程中的尺度規(guī)約個(gè)數(shù)和規(guī)約耗時(shí)均逐漸增加,HLLL與LLL相比僅基于上三角陣操作,其規(guī)約耗時(shí)更少但尺度規(guī)約數(shù)不變,而DLLL在HLLL基礎(chǔ)上引入延后尺度規(guī)約技術(shù),其規(guī)約耗時(shí)和尺度規(guī)約數(shù)均明顯小于 HLLL。此外,三種規(guī)約算法的基向量交換部分的步驟相同,因而其處理后的規(guī)約基是一致的。

    表1 不同維數(shù)下的尺度規(guī)約數(shù)和規(guī)約耗時(shí)Tab.1 Numbers of size reduction and the reduction time in different dimensions

    2)第二組模擬實(shí)驗(yàn)

    利用隨機(jī)模擬方法構(gòu)建 3~40維的模糊度方差協(xié)方差陣,同樣每一維構(gòu)造100組數(shù)據(jù),分別采用HLLL、HELLL和DPLLL算法對(duì)其進(jìn)行規(guī)約處理。圖3(a)和圖3(b)分別為三種算法對(duì)應(yīng)100組數(shù)據(jù)的平均尺度規(guī)約數(shù)和平均基向量交換數(shù)隨不同維數(shù)的變化趨勢(shì)。圖4為對(duì)應(yīng)不同維數(shù)下的平均規(guī)約耗時(shí)的情況。

    圖3 不同維數(shù)下的尺度規(guī)約數(shù)和基向量交換數(shù)Fig.3 Numbers of size reduction and base vectors swapping in different dimensions

    從圖3中可以看出,三種算法的尺度規(guī)約數(shù)和基向量交換數(shù)均與維數(shù)成正相關(guān)關(guān)系,即隨著維數(shù)的增大而呈整體上升趨勢(shì),其中改進(jìn)后的 DPLLL規(guī)約算法擁有最少的尺度規(guī)約數(shù)和基向量交換數(shù),且明顯小于HLLL和HELLL算法。

    從圖4中可以看出,三種算法的規(guī)約耗時(shí)均隨著維數(shù)的增加而變大,其中:HELLL在低維(10維以下)可能存在耗時(shí)大于HLLL的情況,隨著維數(shù)增加HELLL規(guī)約時(shí)間要小于HLLL算法;而DPLLL則由于大大減少了尺度規(guī)約和基向量交換的操作次數(shù),所需要的規(guī)約耗時(shí)最少。

    圖4 不同維數(shù)下的規(guī)約耗時(shí)Fig.4 Reduction time in different dimensions

    3.2 實(shí)驗(yàn)2

    為進(jìn)一步驗(yàn)證改進(jìn)后算法的有效性,采用HLLL、HELLL和DPLLL算法基于實(shí)測(cè)GNSS數(shù)據(jù)進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)選取基線長(zhǎng)約為10.5 km的GPS+BeiDou雙頻數(shù)據(jù)單歷元處理后的模糊度及其方差協(xié)方差陣,接收機(jī)型號(hào)為Trimble NETR9,采樣間隔為30 s,歷元數(shù)為200,模糊度維數(shù)均在30維以上。

    圖5為三種規(guī)約算法200個(gè)歷元下的規(guī)約耗時(shí)情況,從圖中可以看出,DPLLL算法的規(guī)約耗時(shí)最小,而HELLL與HLLL算法的規(guī)約時(shí)間在實(shí)測(cè)環(huán)境下無(wú)明顯大小差異規(guī)律,這可能與HELLL算法迭代操作過(guò)程有關(guān),其中DPLLL算法的平均規(guī)約耗時(shí)僅為傳統(tǒng)算法的1/10。

    圖5 不同歷元下的規(guī)約耗時(shí)Fig.5 Reduction time in different epochs

    進(jìn)一步地,鑒于格基規(guī)約的本質(zhì)在于促使基向量按特定方向排序,以規(guī)約得到的規(guī)約基的施密特正交化向量長(zhǎng)度排列趨勢(shì)來(lái)衡量不同規(guī)約算法的性能[7],如圖6所示。從圖6中可以看出:HLLL算法規(guī)約后的正交化向量長(zhǎng)度排列抖動(dòng)較大;DPLLL算法由于采用最小列旋轉(zhuǎn)QR分解進(jìn)行預(yù)排序,其排列趨勢(shì)明顯改善;而HELLL算法通過(guò)迭代處理可獲得滿足式(12)的更強(qiáng)的規(guī)約基,其正交化向量長(zhǎng)度排列最為平緩,從而說(shuō)明其規(guī)約后的性能要優(yōu)于HLLL和DPLLL算法。這里 HELLL算法可獲得性能最優(yōu)的規(guī)約基,但由圖4和圖5可知,相比于DPLLL算法,HELLL算法并不能有效提高規(guī)約效率。

    圖6 單歷元規(guī)約后正交化向量長(zhǎng)度變化趨勢(shì)Fig.6 Trend of the length of orthogonal vectors after reduction in single epoch

    圖7(a)和圖7(b)分別對(duì)應(yīng)三種算法不同歷元下的搜索耗時(shí)和模糊度整體解算耗時(shí)。由圖 7(a)可知,HELLL算法的搜索耗時(shí)小于其他兩種算法,這與圖6的分析是一致的,即HELLL算法獲得的規(guī)約基性能最優(yōu),DPLLL算法次之。但在高維模糊度解算中,搜索耗時(shí)占模糊度整體解算耗時(shí)的比例很小,減少模糊度整體解算耗時(shí)的關(guān)鍵在于提高規(guī)約效率,而不是單一地追求性能最優(yōu),因此考慮結(jié)合規(guī)約耗時(shí)來(lái)比較三種算法對(duì)模糊度解算效果的影響。根據(jù)圖7(b)可以發(fā)現(xiàn),同時(shí)顧及模糊度規(guī)約和搜索過(guò)程時(shí),本文提出的DPLLL算法的整體解算耗時(shí)最小,即模糊度處理整體效率最高。

    圖7 不同歷元下的搜索耗時(shí)和整體解算耗時(shí)Fig.7 Search time and entire computation time in different epochs

    4 結(jié) 論

    本文在分析 LLL規(guī)約算法的流程和特點(diǎn)的基礎(chǔ)上,通過(guò)引入最小列旋轉(zhuǎn)QR分解、延后尺度規(guī)約和部分尺度規(guī)約技術(shù)對(duì)常規(guī)LLL算法進(jìn)行改進(jìn),提出了DPLLL規(guī)約算法,并基于仿真和實(shí)測(cè)數(shù)據(jù)進(jìn)行了驗(yàn)證分析。從實(shí)驗(yàn)結(jié)果可以看出,改進(jìn)后的算法可以有效減少格基規(guī)約過(guò)程中的尺度規(guī)約數(shù)和基向量交換數(shù),并且能夠獲得較好的規(guī)約性能,顯著改善了LLL算法的規(guī)約效果。利用 DPLLL算法處理的規(guī)約效率和模糊度整體解算效率都得到了明顯提高,從而有利于模糊度解算的快速實(shí)現(xiàn)。同時(shí),隨著組合系統(tǒng)定位和觀測(cè)頻率的增加,如何更好地平衡高維模糊度解算過(guò)程中的規(guī)約效率和規(guī)約性能將是下一步研究的重點(diǎn)。

    (References):

    [1]Li B, Verhagen S, Teunissen P J G.GNSS integer ambiguity estimation and evaluation: LAMBDA and Ps-LAMBDA[C]//Proceedings of China Satellite Navigation Conference.Wuhan, 2013: 291-301.

    [2]Xu P L.Random simulation and GPS decorrelation[J].Journal of Geodesy, 2001, 75: 408-423.

    [3]Chang X, Yang X, Zhou T.MLAMBDA: a modified LAMBDA method for integer least-squares estimation[J].Journal of Geodesy, 2005, 79(9): 552-565.

    [4]Xu P L.Parallel cholesky-based reduction for the weighted integer least squares problem[J].Journal of Geodesy,2012, 86(1): 35-52.

    [5]李豹, 許江寧, 曹可勁, 等.改進(jìn) LAMBDA 算法實(shí)現(xiàn)單頻GPS 整周模糊度快速解算[J].中國(guó)慣性技術(shù)學(xué)報(bào),2013, 21(3): 365-368.Li B, Xu J N, Cao K J, et al.Fast resolution of single frequency GPS integer ambiguity realized by improved LAMBDA algorithm[J].Journal of Chinese Inertial Technology, 2013, 21(3): 365-368.

    [6]劉經(jīng)南, 于興旺, 張小紅.基于格論的GNSS模糊度解算[J].測(cè)繪學(xué)報(bào), 2012, 41(5): 636-645.Liu J N, Yu X W, Zhang X H.GNSS ambiguity resolution using lattice theory[J].Acta Geodaetica et Cartographica Sinica, 2012, 41(5): 636-645.

    [7]范龍.基于格理論的 GNSS模糊度估計(jì)方法研究[D].鄭州: 信息工程大學(xué), 2013.Fan L.Research on method of integer ambiguity estimation with lattice theory[D].Zhengzhou: Information Engineering University, 2013.

    [8]Aliev I, Henk M.LLL-reduction for integer knapsacks[J].Journal of Combinatorial Optimization, 2012, 24(4): 613-626.

    [9]Neumaier A, Stehle D.Faster LLL-type reduction of lattice bases[C]//Proceedings of the ACM on International Symposium on Symbolic and Algebraic Computation.Canada, 2016: 373-380.

    [10]Hassibi A, Boyed S.Integer parameter estimation in linear models with applications to GPS[J].IEEE Transactions on Signal Processing, 1998, 46(11): 2938-2952.

    [11]Lannes A.On the theoretical link between LLL-reduction and Lambda-decorrelation[J].Journal of Geodesy, 2013,87(4): 323-335.

    [12]Cong L, Howgrave-graham N.Effective LLL reduction for lattice decoding[C]//IEEE International Symposium on Information Theory.Nice, France, 2007: 196-200.

    [13]Luo Y X, Qiao S Z.A parallel LLL algorithm[C]//Proceedings of the Fourth International C* Conference on Computer Science and Software Engineering.Montreal.Quebec, Canada, 2011: 93-101.

    [14]Xie X, Chang X W.Borno M A.Partial LLL reduction[C]//Proceeding of IEEE GLOBECOM.Houston, USA, 2011.

    [15]盧立果, 劉萬(wàn)科, 李江衛(wèi).一種有效的LLL規(guī)約算法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2016, 41(8):1118-1124.Lu L, Liu W K, Li J W.An effective LLL reduction algorithm[J].Geomatics and Information Science of Wuhan University, 2016, 41(8): 1118-1124.

    [16]Borno M A, Chang X W, Xie X H.On decorrelation in solving integer least-squares for ambiguity determination[J].Survey Review, 2014, 46: 37-49.

    [17]盧立果,劉萬(wàn)科,李江衛(wèi).降相關(guān)對(duì)模糊度解算中搜索效率的影響分析[J].測(cè)繪學(xué)報(bào)2015, 44(5):481-487.Lu L G, Liu W K, Li J W.Impact of decorrelation on search efficiency of ambiguity resolution[J].Acta Geodaetica et Cartographica Sinica, 2015, 44(5): 481-487.

    [18]盧立果.基于球形搜索的模糊度格基規(guī)約方法研究與分析[D].武漢: 武漢大學(xué), 2013.Lu L G.The research and analysis based on sphere search for ambiguity on reduction[D].Wuhan: Wuhan University,2013.

    [19]Jazaerl S, Amiri-Simkooei A R, Sharifi M A.Fast integer least-squares estimation for GNSS high-dimensional ambiguity resolution using lattice theory[J].Journal of Geodesy,2012, 86: 123-136.

    [20]謝愷, 柴洪洲, 范龍, 等.一種改進(jìn)的LLL模糊度降相關(guān)算法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2014, 39(11):1363-1368.Xie K, Chai H Z, Fan L, et al.An improved LLL ambiguity decorrelation algorithm[J].Geomatics and Information Science of Wuhan University, 2014, 39(11): 1363-1368.

    [21]Wu M K, He J, Wang H C.Reliable single epoch ambiguity resolution for precise attitude determination using BeiDou triple-frequency observations[J].Journal of Chinese Inertial Technology, 2016, 24(5): 624-631.

    Improved LLL ambiguity reduction algorithm

    LYU Hao1, LYU Zhi-ping1, ZHAI Shu-feng1, KUANG Ying-cai1, WANG Fu-lin2
    (1.Institute of Geography and Spatial Information, Information Engineering University,Zhengzhou 450000, China; 2.Unit 61287 of PLA, Kunming 650000, China)

    Efficient ambiguity resolution is the key of GNSS data processing with high precision.Based on the lattice theory, the computational efficiency of searching the optimal integer ambiguity vector can be improved by using lattice reduction.According to the fact that the routine LLL reduction algorithm for ambiguity resolution under high-dimensional case has long reduction time and limited reduction performance, the QR factorization with minimum column pivoting was introduced to preorder the reduction base vectors, and the algorithms of delayed size reduction and partial size reduction were used to reduce the number of redundant size reduction.These methods were optimally combined to improve the implementation effects of LLL reduction algorithm.Experiments with simulated and real GNSS data were executed,respectively.The consequence shows that the modified LLL algorithm can effectively reduce the reduction time and provide better performance than those of current reduction methods.The computing efficiency of lattice reduction using the proposed method with real data is improved by nearly 10 times than those of traditional algorithms.Therefore the computational efficiency of ambiguity resolution under high-dimensional case can be improved significantly by using the proposed algorithm.

    ambiguity resolution; lattice reduction; LLL reduction; QR decomposition; size reduction

    P228.1

    A

    1005-6734(2017)05-0611-07

    10.13695/j.cnki.12-1222/o3.2017.05.010

    2017-07-25;

    2017-09-10

    國(guó)家自然科學(xué)基金(41674019);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFB0501701)

    呂浩(1989—),男,博士研究生,從事測(cè)量數(shù)據(jù)處理理論與方法研究。E-mail: lyangeo@126.com

    猜你喜歡
    規(guī)約尺度向量
    向量的分解
    聚焦“向量與三角”創(chuàng)新題
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    電力系統(tǒng)通信規(guī)約庫(kù)抽象設(shè)計(jì)與實(shí)現(xiàn)
    一種在復(fù)雜環(huán)境中支持容錯(cuò)的高性能規(guī)約框架
    基于XML的電力二次設(shè)備異構(gòu)規(guī)約建模與轉(zhuǎn)換*
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    向量垂直在解析幾何中的應(yīng)用
    向量五種“變身” 玩轉(zhuǎn)圓錐曲線
    9
    午夜免费成人在线视频| 中文字幕高清在线视频| 老司机影院毛片| 他把我摸到了高潮在线观看| 免费少妇av软件| 18禁观看日本| 午夜视频精品福利| 飞空精品影院首页| 午夜视频精品福利| 午夜精品国产一区二区电影| 麻豆乱淫一区二区| 黑人巨大精品欧美一区二区mp4| 国产成+人综合+亚洲专区| 久久久久久久久免费视频了| av天堂久久9| 免费av中文字幕在线| 中出人妻视频一区二区| 最新在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 久久精品成人免费网站| 中亚洲国语对白在线视频| 一夜夜www| 日韩视频一区二区在线观看| 大香蕉久久网| 亚洲精品国产精品久久久不卡| 亚洲av成人av| 人人妻人人爽人人添夜夜欢视频| 国产91精品成人一区二区三区| 宅男免费午夜| 一区福利在线观看| 久久亚洲精品不卡| 久久国产乱子伦精品免费另类| 69精品国产乱码久久久| 美女 人体艺术 gogo| 性少妇av在线| 亚洲人成伊人成综合网2020| 两人在一起打扑克的视频| 国产成人一区二区三区免费视频网站| 日本一区二区免费在线视频| 国产欧美日韩一区二区精品| 无人区码免费观看不卡| 成人黄色视频免费在线看| 亚洲熟女精品中文字幕| 久久久久国内视频| 免费女性裸体啪啪无遮挡网站| 中文字幕人妻丝袜制服| 天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 男女免费视频国产| 18禁美女被吸乳视频| av网站免费在线观看视频| 中国美女看黄片| 久久九九热精品免费| 在线观看免费日韩欧美大片| 后天国语完整版免费观看| 成年人黄色毛片网站| 少妇的丰满在线观看| 熟女少妇亚洲综合色aaa.| 成人影院久久| 亚洲精品久久午夜乱码| 国产成+人综合+亚洲专区| 老司机在亚洲福利影院| 一级黄色大片毛片| 婷婷成人精品国产| 欧美 日韩 精品 国产| 精品国产超薄肉色丝袜足j| 午夜两性在线视频| 又紧又爽又黄一区二区| 妹子高潮喷水视频| 久久人人97超碰香蕉20202| 大型黄色视频在线免费观看| 午夜影院日韩av| 在线观看免费视频网站a站| 亚洲成人国产一区在线观看| 啦啦啦 在线观看视频| 99国产综合亚洲精品| 国产精品自产拍在线观看55亚洲 | 日韩欧美免费精品| 18禁美女被吸乳视频| 男女免费视频国产| 亚洲精品自拍成人| 欧美日韩视频精品一区| 后天国语完整版免费观看| 成人黄色视频免费在线看| 很黄的视频免费| 十分钟在线观看高清视频www| 国产成人免费观看mmmm| 亚洲免费av在线视频| 一级片'在线观看视频| 777米奇影视久久| 欧美中文综合在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久电影中文字幕 | 超碰成人久久| 变态另类成人亚洲欧美熟女 | 老司机午夜十八禁免费视频| 18禁裸乳无遮挡动漫免费视频| 99久久综合精品五月天人人| 国产亚洲欧美精品永久| 制服人妻中文乱码| 热re99久久国产66热| av不卡在线播放| 一边摸一边做爽爽视频免费| 黄色 视频免费看| 亚洲国产精品sss在线观看 | 国产亚洲一区二区精品| 国产成人欧美| av天堂在线播放| 国产蜜桃级精品一区二区三区 | 亚洲成国产人片在线观看| 婷婷丁香在线五月| 久久国产精品人妻蜜桃| 精品国内亚洲2022精品成人 | 免费看a级黄色片| 一区福利在线观看| 老汉色av国产亚洲站长工具| 99国产综合亚洲精品| 亚洲精华国产精华精| 国产色视频综合| 久久久久久亚洲精品国产蜜桃av| 国产精品欧美亚洲77777| 窝窝影院91人妻| 国产在视频线精品| 老汉色∧v一级毛片| 国产真人三级小视频在线观看| 国产又色又爽无遮挡免费看| 久久精品国产清高在天天线| 天堂√8在线中文| 成人亚洲精品一区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产欧美日韩在线播放| 精品亚洲成a人片在线观看| 国产又色又爽无遮挡免费看| 91精品三级在线观看| 久久中文字幕一级| 交换朋友夫妻互换小说| 69精品国产乱码久久久| 国产亚洲欧美98| 我的亚洲天堂| 久久精品亚洲精品国产色婷小说| 亚洲午夜精品一区,二区,三区| 国产高清视频在线播放一区| 亚洲成人免费av在线播放| 国产高清激情床上av| 欧美日本中文国产一区发布| 久久99一区二区三区| 国产欧美日韩一区二区精品| 一进一出抽搐动态| 国产精品免费大片| 窝窝影院91人妻| 欧美不卡视频在线免费观看 | 午夜久久久在线观看| 高潮久久久久久久久久久不卡| 国产在视频线精品| 精品一区二区三区视频在线观看免费 | 亚洲欧美色中文字幕在线| 欧美黑人欧美精品刺激| 国内毛片毛片毛片毛片毛片| 久久99一区二区三区| 久久久国产成人免费| 99热网站在线观看| 高清在线国产一区| 天天躁日日躁夜夜躁夜夜| 黄色成人免费大全| 亚洲专区中文字幕在线| 国产精品香港三级国产av潘金莲| 国产精品九九99| 男女下面插进去视频免费观看| 欧美日韩福利视频一区二区| 亚洲第一欧美日韩一区二区三区| 91大片在线观看| 国产在视频线精品| 精品一区二区三区视频在线观看免费 | 一进一出抽搐动态| 免费av中文字幕在线| 好看av亚洲va欧美ⅴa在| 不卡一级毛片| 国产xxxxx性猛交| 久久久水蜜桃国产精品网| 中文字幕人妻丝袜制服| 99re6热这里在线精品视频| videosex国产| 亚洲熟妇熟女久久| 一夜夜www| 人人妻人人澡人人看| 制服人妻中文乱码| 免费日韩欧美在线观看| av天堂在线播放| 久久香蕉国产精品| 老鸭窝网址在线观看| 国产精品秋霞免费鲁丝片| 搡老乐熟女国产| 高清欧美精品videossex| 亚洲全国av大片| 一本一本久久a久久精品综合妖精| 亚洲av日韩在线播放| 国产欧美日韩综合在线一区二区| 国产成人啪精品午夜网站| 亚洲av片天天在线观看| 久久久久久免费高清国产稀缺| 一a级毛片在线观看| 啦啦啦 在线观看视频| 99riav亚洲国产免费| 国产精品欧美亚洲77777| 99国产精品一区二区三区| 欧美亚洲日本最大视频资源| 久久久精品免费免费高清| 午夜福利在线免费观看网站| 十八禁人妻一区二区| 欧美中文综合在线视频| 亚洲欧美一区二区三区久久| 欧美 日韩 精品 国产| 成年动漫av网址| 国产有黄有色有爽视频| 国产欧美日韩综合在线一区二区| 久久久久精品国产欧美久久久| 老司机深夜福利视频在线观看| 黄色 视频免费看| 亚洲片人在线观看| 99精国产麻豆久久婷婷| 久久精品亚洲av国产电影网| 欧美丝袜亚洲另类 | 操出白浆在线播放| 国产片内射在线| 国产精品亚洲一级av第二区| 日韩成人在线观看一区二区三区| 亚洲av成人av| 在线观看免费日韩欧美大片| 亚洲综合色网址| 亚洲午夜理论影院| 老熟妇仑乱视频hdxx| 欧美老熟妇乱子伦牲交| 国产精品久久视频播放| 午夜91福利影院| 曰老女人黄片| 亚洲美女黄片视频| 国产亚洲欧美在线一区二区| 男人操女人黄网站| 亚洲情色 制服丝袜| 国产精品一区二区精品视频观看| 久久精品熟女亚洲av麻豆精品| 亚洲一码二码三码区别大吗| 91麻豆av在线| 在线十欧美十亚洲十日本专区| 日韩免费av在线播放| 久久久久精品人妻al黑| 中文欧美无线码| 脱女人内裤的视频| 不卡一级毛片| 1024香蕉在线观看| 亚洲精品自拍成人| 91大片在线观看| 色老头精品视频在线观看| av欧美777| 交换朋友夫妻互换小说| 老司机午夜福利在线观看视频| 色婷婷久久久亚洲欧美| 搡老岳熟女国产| 性色av乱码一区二区三区2| 免费不卡黄色视频| 91精品国产国语对白视频| 每晚都被弄得嗷嗷叫到高潮| 老汉色av国产亚洲站长工具| 欧美精品一区二区免费开放| 久久热在线av| 国产亚洲av高清不卡| 免费在线观看完整版高清| 最新的欧美精品一区二区| 久久久国产一区二区| 无限看片的www在线观看| 新久久久久国产一级毛片| 老司机影院毛片| 成人国语在线视频| 色婷婷久久久亚洲欧美| 午夜福利乱码中文字幕| 91av网站免费观看| 欧美丝袜亚洲另类 | 精品一区二区三区四区五区乱码| 国产激情欧美一区二区| 黄色 视频免费看| 亚洲精品国产精品久久久不卡| 欧美日韩福利视频一区二区| 看片在线看免费视频| 国产精品影院久久| 天天添夜夜摸| 最新在线观看一区二区三区| 咕卡用的链子| 亚洲国产精品sss在线观看 | 久热这里只有精品99| 男女床上黄色一级片免费看| 欧美+亚洲+日韩+国产| 国产91精品成人一区二区三区| 18禁裸乳无遮挡免费网站照片 | avwww免费| 国产成人精品在线电影| 精品久久久久久电影网| 女人爽到高潮嗷嗷叫在线视频| 村上凉子中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 最新在线观看一区二区三区| 一本大道久久a久久精品| 男女午夜视频在线观看| 久久香蕉国产精品| 日本一区二区免费在线视频| 亚洲 欧美一区二区三区| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 首页视频小说图片口味搜索| 香蕉久久夜色| 搡老岳熟女国产| 免费人成视频x8x8入口观看| xxx96com| 亚洲全国av大片| 久久久国产精品麻豆| 精品无人区乱码1区二区| 欧美性长视频在线观看| 国产主播在线观看一区二区| 麻豆av在线久日| 国产在线精品亚洲第一网站| 国产一区二区三区视频了| 久久久久久久午夜电影 | 亚洲人成伊人成综合网2020| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美色中文字幕在线| 亚洲黑人精品在线| 国产精品 欧美亚洲| 欧美日韩黄片免| 亚洲精品粉嫩美女一区| 视频区图区小说| 欧美国产精品va在线观看不卡| 国产精品久久久av美女十八| 欧美在线一区亚洲| 最近最新中文字幕大全电影3 | 国产精品亚洲av一区麻豆| 久久精品成人免费网站| 啦啦啦 在线观看视频| 久久狼人影院| 亚洲精华国产精华精| 亚洲情色 制服丝袜| 日韩免费高清中文字幕av| 亚洲欧美一区二区三区久久| 免费久久久久久久精品成人欧美视频| 这个男人来自地球电影免费观看| 日韩视频一区二区在线观看| 欧美日韩国产mv在线观看视频| 丝瓜视频免费看黄片| 国产av精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 建设人人有责人人尽责人人享有的| 日本黄色日本黄色录像| 国产精品美女特级片免费视频播放器 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产精品合色在线| 午夜福利一区二区在线看| tocl精华| 亚洲欧美日韩高清在线视频| 色在线成人网| a级毛片在线看网站| 国产成人av激情在线播放| 大型黄色视频在线免费观看| 亚洲av日韩在线播放| 精品久久久久久久久久免费视频 | 丝袜美足系列| 一a级毛片在线观看| 大香蕉久久成人网| 午夜91福利影院| 午夜福利乱码中文字幕| 青草久久国产| 久久精品国产综合久久久| 久热这里只有精品99| 免费观看a级毛片全部| 国产成人精品久久二区二区免费| 亚洲欧美色中文字幕在线| 国产成人影院久久av| 中文字幕高清在线视频| 无遮挡黄片免费观看| 青草久久国产| 国产在视频线精品| 午夜福利一区二区在线看| 一区在线观看完整版| 黑人操中国人逼视频| 久久久精品免费免费高清| 男人操女人黄网站| 无限看片的www在线观看| 人人妻人人爽人人添夜夜欢视频| 成人国语在线视频| 国产亚洲一区二区精品| 成在线人永久免费视频| 伦理电影免费视频| 成人国产一区最新在线观看| 欧美在线一区亚洲| 丁香六月欧美| 在线观看一区二区三区激情| 国产黄色免费在线视频| 亚洲欧美日韩高清在线视频| 午夜福利视频在线观看免费| 久久天躁狠狠躁夜夜2o2o| 国产1区2区3区精品| av欧美777| 亚洲人成电影免费在线| 美女国产高潮福利片在线看| 亚洲男人天堂网一区| 18禁裸乳无遮挡动漫免费视频| 亚洲五月色婷婷综合| videosex国产| 欧美精品av麻豆av| av天堂在线播放| 视频在线观看一区二区三区| 一级毛片精品| 成年动漫av网址| 国产欧美日韩一区二区精品| 亚洲自偷自拍图片 自拍| 欧美黄色片欧美黄色片| 亚洲欧美激情综合另类| 岛国在线观看网站| 日韩三级视频一区二区三区| aaaaa片日本免费| 脱女人内裤的视频| videos熟女内射| 亚洲七黄色美女视频| 亚洲一区高清亚洲精品| 亚洲熟女毛片儿| 99精国产麻豆久久婷婷| 成人影院久久| 久久精品国产a三级三级三级| 在线观看午夜福利视频| 天天躁日日躁夜夜躁夜夜| 午夜福利乱码中文字幕| 在线观看免费午夜福利视频| 老熟妇乱子伦视频在线观看| 美女高潮到喷水免费观看| 男女高潮啪啪啪动态图| 每晚都被弄得嗷嗷叫到高潮| 黄频高清免费视频| 欧美久久黑人一区二区| 国产野战对白在线观看| 国产精品偷伦视频观看了| 精品卡一卡二卡四卡免费| 国产免费av片在线观看野外av| 操美女的视频在线观看| 亚洲五月色婷婷综合| 一夜夜www| 黄色毛片三级朝国网站| 国产蜜桃级精品一区二区三区 | 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 看免费av毛片| 午夜福利乱码中文字幕| 午夜老司机福利片| 在线观看www视频免费| a级毛片在线看网站| 午夜日韩欧美国产| 亚洲成av片中文字幕在线观看| 极品人妻少妇av视频| 午夜福利视频在线观看免费| 老汉色∧v一级毛片| 99国产精品一区二区蜜桃av | 久久中文看片网| 性少妇av在线| 飞空精品影院首页| 夜夜爽天天搞| 免费在线观看视频国产中文字幕亚洲| 99re6热这里在线精品视频| 成人影院久久| 自线自在国产av| 久久精品91无色码中文字幕| 成人永久免费在线观看视频| 中文字幕人妻丝袜一区二区| 婷婷丁香在线五月| www.自偷自拍.com| 成人手机av| 手机成人av网站| 黄色视频,在线免费观看| xxxhd国产人妻xxx| 久久亚洲精品不卡| av不卡在线播放| 十八禁网站免费在线| 丰满人妻熟妇乱又伦精品不卡| 精品电影一区二区在线| 成人手机av| 手机成人av网站| 丝袜美腿诱惑在线| 日日爽夜夜爽网站| 人妻一区二区av| √禁漫天堂资源中文www| 亚洲国产精品sss在线观看 | 欧美日韩乱码在线| 一区二区三区国产精品乱码| 国产精品二区激情视频| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 国产国语露脸激情在线看| 亚洲三区欧美一区| 男女下面插进去视频免费观看| 精品久久久精品久久久| 亚洲熟妇熟女久久| 久久久久精品国产欧美久久久| 久久国产乱子伦精品免费另类| 国产成人av教育| 国产99白浆流出| 欧美精品一区二区免费开放| 麻豆乱淫一区二区| 亚洲一码二码三码区别大吗| 精品国产国语对白av| 美女高潮到喷水免费观看| 人人妻人人添人人爽欧美一区卜| 91麻豆精品激情在线观看国产 | 国产精品欧美亚洲77777| 91成人精品电影| 操出白浆在线播放| 身体一侧抽搐| 欧美亚洲 丝袜 人妻 在线| 美女国产高潮福利片在线看| 中文字幕人妻丝袜制服| 午夜91福利影院| 日韩三级视频一区二区三区| 精品久久久久久久久久免费视频 | 动漫黄色视频在线观看| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 免费观看人在逋| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文av在线| 欧美最黄视频在线播放免费 | 久久久精品区二区三区| 这个男人来自地球电影免费观看| 岛国在线观看网站| 一夜夜www| 亚洲国产欧美网| 91成年电影在线观看| 色综合婷婷激情| av有码第一页| 精品国内亚洲2022精品成人 | 久久香蕉激情| 亚洲精品中文字幕在线视频| 男人操女人黄网站| 精品一区二区三区视频在线观看免费 | 欧美不卡视频在线免费观看 | 大型黄色视频在线免费观看| 后天国语完整版免费观看| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91| 两性午夜刺激爽爽歪歪视频在线观看 | 身体一侧抽搐| 久久久久久久精品吃奶| 电影成人av| 成人黄色视频免费在线看| 免费在线观看日本一区| 激情在线观看视频在线高清 | 欧美亚洲日本最大视频资源| 他把我摸到了高潮在线观看| 91国产中文字幕| 50天的宝宝边吃奶边哭怎么回事| 欧美成人午夜精品| videosex国产| 大片电影免费在线观看免费| 欧美日韩瑟瑟在线播放| 久久久精品免费免费高清| 精品欧美一区二区三区在线| 人人妻人人添人人爽欧美一区卜| 日韩制服丝袜自拍偷拍| 亚洲精品中文字幕一二三四区| 国产97色在线日韩免费| 色老头精品视频在线观看| 久久天堂一区二区三区四区| 亚洲人成电影观看| 亚洲欧美日韩高清在线视频| 日韩欧美国产一区二区入口| 国产欧美日韩精品亚洲av| 最新在线观看一区二区三区| av免费在线观看网站| 十八禁网站免费在线| 午夜老司机福利片| a级毛片黄视频| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 在线十欧美十亚洲十日本专区| 另类亚洲欧美激情| 国产欧美日韩综合在线一区二区| 免费在线观看视频国产中文字幕亚洲| 他把我摸到了高潮在线观看| 波多野结衣一区麻豆| 国产精品欧美亚洲77777| 亚洲精品在线美女| 亚洲一码二码三码区别大吗| 亚洲情色 制服丝袜| 亚洲国产精品sss在线观看 | 久久香蕉精品热| 欧美日韩中文字幕国产精品一区二区三区 | 一夜夜www| 精品国产国语对白av| 成年人免费黄色播放视频| 国产高清videossex| 国产真人三级小视频在线观看| 老司机午夜十八禁免费视频| 婷婷成人精品国产| 嫁个100分男人电影在线观看| 亚洲一区中文字幕在线| 一个人免费在线观看的高清视频| 国内久久婷婷六月综合欲色啪| 亚洲国产毛片av蜜桃av| bbb黄色大片| 一区二区三区激情视频| 色精品久久人妻99蜜桃| 欧美精品啪啪一区二区三区| 亚洲av电影在线进入| 亚洲色图综合在线观看| 久久中文看片网| 夫妻午夜视频| 国产一区在线观看成人免费| 美女 人体艺术 gogo| 欧美亚洲日本最大视频资源| 久久精品亚洲av国产电影网| 水蜜桃什么品种好| 交换朋友夫妻互换小说|