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

    多元總體最小二乘問題的牛頓解法

    2016-05-16 08:38:45王樂洋趙英文陳曉勇臧德彥
    測(cè)繪學(xué)報(bào) 2016年4期

    王樂洋,趙英文,陳曉勇,臧德彥

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013

    ?

    多元總體最小二乘問題的牛頓解法

    王樂洋1,2,3,趙英文1,陳曉勇1,2,3,臧德彥1,2,3

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013

    transformation

    Foundation support: The National Natural Science Foundation of China (Nos. 41204003; 41161069; 41304020; 41464001); The National Department Public Benefit Research Foundation (Surveying,Mapping and Geoinformation)(No. 201512026); The Natural Science Foundation of Jiangxi Province (No. 20151BAB203042); The Science and Technology Project of the Education Department of Jiangxi Province (Nos. GJJ150595; KJLD12077; KJLD14049); The Found of Key Laboratory of Watershed Ecology and Geographical Environment Monitoring (No. WE2015005); The Scientific Research Foundation of ECIT(No. DHBK201113); Innovation Fund Designated for Graduate Students of Jiangxi Province (Nos. YC2015-S266; YC2015-S267); Innovation Fund Designated for Graduate Students of ECIT(No. DHYC-2015005); The Found of the Key Laboratory of Mapping from Space, NASG (No.K201502)

    摘要:為提高多元總體最小二乘問題參數(shù)估值的解算效率,推導(dǎo)了基于牛頓法的多元加權(quán)總體最小二乘算法;分析比較了基于牛頓法的多元加權(quán)總體最小二乘解和基于拉格朗日乘數(shù)法多元加權(quán)總體最小二乘解之間的關(guān)系,根據(jù)協(xié)因數(shù)傳播律給出了多元總體最小二乘平差的16種協(xié)因數(shù)陣的近似計(jì)算公式。新算法能夠解決觀測(cè)矩陣和系數(shù)矩陣元素具有相關(guān)性的問題,并且可以把觀測(cè)矩陣和系數(shù)矩陣的隨機(jī)元素和常數(shù)元素納入到一個(gè)協(xié)因數(shù)陣中進(jìn)行處理。算例結(jié)果表明,本文提出的多元總體最小二乘問題的牛頓解法可行且收斂速度更快。

    關(guān)鍵詞:多元總體最小二乘;牛頓法;協(xié)因數(shù)傳播律;仿射變換

    總體最小二乘(total least squares,TLS)方法是近30多年來發(fā)展起來的一種能同時(shí)顧及觀測(cè)值誤差和模型系數(shù)矩陣誤差的數(shù)學(xué)方法,其理論及應(yīng)用研究是目前國(guó)內(nèi)外研究的熱點(diǎn)問題[1]。文獻(xiàn)[2]首次提出總體最小二乘概念,隨著變量誤差模型(errors-in-variables,EIV)由單變量模型擴(kuò)展到多元變量模型[3],總體最小二乘也相應(yīng)擴(kuò)展到多元總體最小二乘(multivariate TLS,MTLS)[4-5]。某些模型的系數(shù)矩陣會(huì)因未知參數(shù)寫成向量形式含有重復(fù)出現(xiàn)的隨機(jī)元素和大量的零元素;多元總體最小二乘將參數(shù)向量和觀測(cè)向量改寫成矩陣形式,避免了這種現(xiàn)象,進(jìn)而削弱了系數(shù)矩陣中行列之間的相關(guān)性,而且使得系數(shù)矩陣的構(gòu)造更加簡(jiǎn)單。

    在國(guó)外相關(guān)研究中,文獻(xiàn)[5—7]認(rèn)為觀測(cè)矩陣和系數(shù)矩陣中的隨機(jī)元素是等精度的,給出了基于奇異值分解法和拉格朗日乘數(shù)法的多元總體最小二乘算法;文獻(xiàn)[8]討論了基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘算法,將隨機(jī)元素的方差陣分解成兩個(gè)矩陣取直積的形式,但是當(dāng)觀測(cè)數(shù)據(jù)含有兩個(gè)以上的方差分量或是觀測(cè)數(shù)據(jù)含有多類觀測(cè)值時(shí),隨機(jī)變量的方差協(xié)方差陣就很難寫成或是不能寫成兩個(gè)矩陣取直積的形式;文獻(xiàn)[9]推導(dǎo)了基于奇異值分解法的多元加權(quán)總體最小二乘算法,但是這種算法只使用了行尺度權(quán)陣,其權(quán)陣不具有一般性。以上文獻(xiàn)的所有算法都不適用于系數(shù)矩陣與觀測(cè)矩陣中隨機(jī)元素具有相關(guān)性的情況。文獻(xiàn)[10]給出了基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘算法,其權(quán)陣具有一般性,但沒有具體給出系數(shù)矩陣含有常數(shù)元素的解法以及沒有解決算法的精度評(píng)定問題,這種算法的解算效率也還需要驗(yàn)證。在國(guó)內(nèi)相關(guān)研究中,文獻(xiàn)[11]將基于奇異值分解法的多元總體最小二乘算法應(yīng)用在三維坐標(biāo)轉(zhuǎn)換中;文獻(xiàn)[12]給出的多元總體最小二乘算法受到正交Procrustes算法制約僅適用于觀測(cè)值誤差對(duì)稱獨(dú)立的情況,隨機(jī)元素相關(guān)性問題也沒有解決;文獻(xiàn)[13]只給出了觀測(cè)值等精度同方差情況下的多元總體最小二乘算法。對(duì)于系數(shù)矩陣含有常量的情況,文獻(xiàn)[10,14—16]將系數(shù)矩陣拆分成分別包含隨機(jī)元素和常數(shù)元素的系數(shù)矩陣;文獻(xiàn)[11—13]采用消元法,分兩步計(jì)算隨機(jī)元素對(duì)應(yīng)的未知參數(shù)和常數(shù)元素對(duì)應(yīng)的未知參數(shù);文獻(xiàn)[8,17]使常數(shù)元素的方差和協(xié)方差為零,然后進(jìn)行加權(quán)解算;文獻(xiàn)[18]首次將EIV模型改寫成更具一般性的PEIV(partial EIV)模型,解決了系數(shù)矩陣任意位置含有常數(shù)元素的平差問題。但以上方法都沒有把系數(shù)矩陣和觀測(cè)向量(或矩陣)中的所有元素納入到一個(gè)協(xié)因數(shù)陣中進(jìn)行處理。

    基于以上問題,本文推導(dǎo)了顧及系數(shù)矩陣含有常數(shù)元素情況下的多元總體最小二乘的牛頓解法,并討論了多元總體最小二乘的牛頓法解和拉格朗日乘數(shù)法解之間的關(guān)系以及給出了多元總體最小二乘的精度評(píng)定公式,最后通過模擬數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)驗(yàn)證了本文算法的可行性和有效性。

    1多元總體最小二乘解法

    1.1多元變量EIV模型

    多元變量EIV模型[6-8]為

    Y=(A+EA)Ξ-EY

    (1)

    (2)

    顧及系數(shù)矩陣A中某列含有常數(shù)元素,在文獻(xiàn)[8,17—18]的基礎(chǔ)上,系數(shù)矩陣和相應(yīng)的協(xié)因數(shù)陣定義為

    (3)

    (4)

    式中,As是系數(shù)矩陣A的隨機(jī)元素部分;Ad是系數(shù)矩陣A的常數(shù)元素部分;Es為隨機(jī)元素的誤差矩陣;Qi(i=1,2,…,t)∈Rn×n(t+m)。

    此時(shí),協(xié)因數(shù)陣Q為奇異陣,但在計(jì)算時(shí)并沒有直接求取Q的逆,因此并不影響解算過程。

    1.2多元總體最小二乘問題的牛頓解法

    本文在文獻(xiàn)[10]推導(dǎo)牛頓法思路的基礎(chǔ)上,進(jìn)一步推導(dǎo)了一般情形下多元總體最小二乘問題的牛頓解法。根據(jù)總體最小二乘問題平差準(zhǔn)則VTPV=min,由求條件極值的拉格朗日乘數(shù)法構(gòu)造如下目標(biāo)函數(shù)

    Γ(V,Ξ,λ)=VTPV+2λT(vec(Y+EY)-

    (ΞT?In)vec(A+EA))

    (5)

    式中,λ∈Rnm×1為拉格朗日乘數(shù)向量;?為Kronecker-Zehfuss積。

    對(duì)式(5)中的元素分別求導(dǎo)得

    (6)

    (7)

    由式(6)和式(7),得到改正數(shù)最優(yōu)估值為[10]

    (8)

    進(jìn)而得到牛頓法求解無約束問題的目標(biāo)函數(shù)

    (9)

    式中

    式(9)對(duì)vec(Ξ)求一階導(dǎo)數(shù),得

    (10)

    經(jīng)推導(dǎo)可得式(10)第1部分

    (11)

    由逆矩陣求導(dǎo)原理,得式(10)第2部分

    (12)

    式中

    (13)

    (14)

    由式(11)和式(12)得

    (15)

    式(9)對(duì)vec(Ξ)求二階導(dǎo)數(shù),得

    (16)

    并令

    (17)

    (18)

    依據(jù)牛頓法原理得到未知參數(shù)的改正數(shù)

    (19)

    以及

    (20)

    多元總體最小二乘問題的牛頓解法的迭代步驟為(k為迭代次數(shù)):

    作為迭代初始值。

    第2步,由上一步獲得的數(shù)值,分別計(jì)算

    由式(15)和式(18),進(jìn)一步計(jì)算Gk+1、Hk+1和

    第3步,根據(jù)式(20)更新未知參數(shù)估值。

    作為未知參數(shù)的最佳估值。

    1.3牛頓法解和拉格朗日乘數(shù)法解之間的關(guān)系

    拉格朗日乘數(shù)法的多元總體最小二乘解為[10]

    (21)

    式(21)在近似值處展開,寫成式(20)的形式

    (22)

    式中,拉格朗日乘數(shù)法未知參數(shù)改正數(shù)為

    (23)

    (24)

    利用矩陣反演公式[10,19],得

    (25)

    又因?yàn)?/p>

    (26)

    聯(lián)合式(23)、式(25)和式(26),得

    (27)

    2多元總體最小二乘解法精度評(píng)定

    2.1單位權(quán)方差估值公式

    根據(jù)文獻(xiàn)[6—7],多元總體最小二乘解法的單位權(quán)方差估值公式為

    (28)

    式中,m(n-t)為自由度,mn為觀測(cè)矩陣Y元素的個(gè)數(shù),mt為參數(shù)矩陣Ξ元素的個(gè)數(shù)。

    2.2協(xié)因數(shù)陣計(jì)算公式

    多元總體最小二乘法中,數(shù)據(jù)多以矩陣形式存在。為計(jì)算協(xié)因數(shù)陣,把矩陣寫成向量形式,得到

    (29)

    表1多元總體最小二乘法協(xié)因數(shù)陣

    Tab.1Cofactor matrice of multivariate total least squares adjustment

    L^V^L^ξLQ-^X5Q-^X5-Q^ΣT^XT4^V-^X5^X50Q^ΣT^XT4^LQ-^X50Q-^X50^ξ-^X4^ΣQ^X4^ΣQ0^X3

    3算例分析

    3.1模擬算例

    根據(jù)文獻(xiàn)[17]中給出的源坐標(biāo)系的真值數(shù)據(jù),本文通過給定的二維仿射變換參數(shù)真值計(jì)算得到目標(biāo)坐標(biāo)系的真值數(shù)據(jù)。在得到兩組坐標(biāo)系的真值數(shù)據(jù)后,依據(jù)文獻(xiàn)[17],本文在模擬數(shù)據(jù)中源坐標(biāo)系所加誤差的協(xié)因數(shù)陣為

    Qxy=I2?(0.01diag([1,2,3,1,5,4,2,7,2,1,8,3,6]))

    模擬數(shù)據(jù)中目標(biāo)坐標(biāo)系所加誤差的協(xié)因數(shù)陣為

    QXY=I2?(0.01diag([1,3,6,1,1,8,4,3,6,5,4,5,2]))

    共模擬500次,分別使用本文提出的牛頓法(N法)和拉格朗日-牛頓法(L-N法),以及Fang(2011)法[10]、Schaffrin(2009)法[8]和加權(quán)最小二乘法(WLS)計(jì)算二維仿射變換參數(shù),迭代法的終止條件均為10-12,將獲得的二維仿射變換參數(shù)的平均值,以及二維仿射變換參數(shù)平均值與其真值的差值2范數(shù)和平均迭代次數(shù)列于表2中。

    參數(shù)為矩陣時(shí),文獻(xiàn)[10]并沒有給出在系數(shù)矩陣含有常數(shù)列時(shí)的具體公式算法。依據(jù)其提出的在參數(shù)為向量時(shí)系數(shù)矩陣含有常數(shù)列的加權(quán)總體最小二乘算法,本文編程實(shí)現(xiàn)了這種將系數(shù)矩陣分成兩部分的多元加權(quán)總體最小二乘算法,并簡(jiǎn)稱為Fang(2011)法。拉格朗日-牛頓法(L-N法)是指把Fang(2011)法(拉格朗日乘數(shù)法)獲得的第一次多元加權(quán)總體最小二乘參數(shù)解作為本文牛頓法迭代的初始值;在選取迭代初始值時(shí),如果使用多元最小二乘參數(shù)解作為初始值計(jì)算式(20),則為牛頓法(N法),如果使用多元加權(quán)總體最小二乘參數(shù)解作為初始值計(jì)算式(20),則為拉格朗日-牛頓法(L-N法);在能夠獲得多元加權(quán)總體最小二乘參數(shù)解的情況下,建議使用收斂速度更快的L-N法。

    以上算例沒有考慮系數(shù)矩陣元素與觀測(cè)矩陣元素的相關(guān)性。為了驗(yàn)證本文的牛頓法能夠解決更一般的多元總體最小二乘問題,在文獻(xiàn)[19]的基礎(chǔ)上,使用Matlab 7.11.0軟件,利用函數(shù)randn(100,13*4)T×randn(100,13*4)生成協(xié)因數(shù)陣Qr,利用函數(shù)mvnrnd生成均值均為0,方差-協(xié)方差陣為0.001 Qr的隨機(jī)誤差矩陣,并將其施加到真值數(shù)據(jù)中以作為計(jì)算的模擬數(shù)據(jù)。

    由于模擬數(shù)據(jù)使得兩套坐標(biāo)系具有了相關(guān)性,此時(shí)Schaffrin(2009)法不再適用,選用本文提出的牛頓法和拉格朗日-牛頓法,以及Fang(2011)法模擬計(jì)算10 000次,計(jì)算時(shí)取Q=Qr,迭代法的終止條件均為10-12,并由獲得的10 000次的參數(shù)估值繪制出的參數(shù)估值頻率直方圖及其相應(yīng)的擬合出的正態(tài)分布曲線圖如圖1所示。

    計(jì)算結(jié)果表明,牛頓法的平均迭代次數(shù)為6.972 6,拉格朗日-牛頓法的平均迭代次數(shù)為6.828 7,F(xiàn)ang(2011)法的平均迭代次數(shù)為12.188 2。牛頓法和拉格朗日-牛頓法使用的協(xié)因數(shù)陣Q為奇異矩陣,但這并不影響未知參數(shù)的解算。Fang(2011)法的迭代次數(shù)最多,拉格朗日-牛頓法迭代次數(shù)要少于牛頓法,這也與上文的結(jié)論相符合。參數(shù)估值的頻率直方圖近似于正態(tài)分布,相應(yīng)的正態(tài)分布擬合曲線大約也以參數(shù)的真值為對(duì)稱軸。以上的結(jié)論證明了多元總體最小二乘牛頓解法的可行性和有效性。

    表2 5種方法參數(shù)估計(jì)結(jié)果

    圖1 所估參數(shù)的直方圖及其正態(tài)分布擬合曲線Fig.1 Histogram and fitting curve of normal distribution of estimated parameters

    3.2實(shí)測(cè)算例

    選取文獻(xiàn)[19]中二維坐標(biāo)轉(zhuǎn)換數(shù)據(jù)作為計(jì)算數(shù)據(jù),共有5組在源坐標(biāo)系和目標(biāo)坐標(biāo)系下的公共點(diǎn)。

    源坐標(biāo)系和目標(biāo)坐標(biāo)系權(quán)陣為[19]

    Pxy=diag([18.781 7,6.377 4,12.648 9,

    17.476 9,22.272 6,23.982 3,13.680 4,

    3.465 6,3.732 4,6.437 7])

    PXY=diag([9.836 1,5.535 7,12.736 9,

    12.009 9,10.181,11.366 1,11.147,

    5.883 4,9.832 2,7.567 8])

    使用本文提出的牛頓法和拉格朗日-牛頓法,以及Fang(2011)法計(jì)算二維仿射變換參數(shù),迭代法的終止條件均為10-12。3種方法的未知參數(shù)估值列于表3中。利用表3中的未知參數(shù)估值結(jié)果和表1中的公式計(jì)算未知參數(shù)估值協(xié)因數(shù)陣,并將其列于表4中。

    表3 3種方法參數(shù)估計(jì)結(jié)果

    表4 參數(shù)估值協(xié)因數(shù)陣

    4結(jié)論

    本文推導(dǎo)了多元總體最小二乘問題的牛頓解法,并比較了基于牛頓法的多元加權(quán)總體最小二乘解和基于拉格朗日乘數(shù)法的多元加權(quán)總體最小二乘解之間的關(guān)系,給出了多元總體最小二乘算法精度評(píng)定的16種協(xié)因數(shù)陣的近似計(jì)算公式。通過算例驗(yàn)證,并與現(xiàn)有的解法作了比較,得到以下幾點(diǎn)結(jié)論:

    (1) 牛頓法有著更少的迭代次數(shù),更高的解算效率。

    (2) 在觀測(cè)矩陣和系數(shù)矩陣元素對(duì)應(yīng)的協(xié)因數(shù)陣中,本文取常數(shù)元素的協(xié)因數(shù)為0。這雖然導(dǎo)致協(xié)因數(shù)陣Q奇異,但并不影響參數(shù)的估計(jì)結(jié)果;本文的算法也能解決觀測(cè)矩陣與系數(shù)矩陣元素具有相關(guān)性的問題。

    (3) 牛頓法的解算效率會(huì)受到迭代初始值的影響,本文采用的拉格朗日-牛頓法因其迭代初始值更接近參數(shù)最優(yōu)估值,而使得解算效率比起牛頓法有所提高。

    本文給出的多元總體最小二乘解法精度評(píng)定的公式只是近似公式,如何推導(dǎo)更為嚴(yán)密的總體最小二乘精度評(píng)定公式還需進(jìn)一步研究。

    參考文獻(xiàn):

    [1]王樂洋.基于總體最小二乘的大地測(cè)量反演理論及應(yīng)用研究[J]. 測(cè)繪學(xué)報(bào), 2012, 41(4): 629.

    WANG Leyang. Research on Theory and Application of Total Least Squares in Geodetic Inversion[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 629.

    [2]GOLUB G H, VAN LOAN C F. An Analysis of the Total Least Squares Problem[J]. SIAM Journal on Numerical Analysis, 1980, 17(6): 883-893.

    [3]SPRENT P. Models in Regression and Related Topics[M]. London: Methuen & Co Ltd, 1969.

    [4]王樂洋, 許才軍. 總體最小二乘研究進(jìn)展[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2013, 38(7): 850-856, 878.

    WANG Leyang, XU Caijun. Progress in Total Least Squares[J]. Geomatics and Information Science of Wuhan University, 2013, 38(7): 850-856, 878.

    [5]VAN HUFFEL S, VANDEWALLE J. The Total Least Squares Problem: Computational Aspects and Analysis[D]. Philadelphia, Pennsylvania: SIAM, 1991.

    [6]SCHAFFRIN B, FELUS Y A. Multivariate Total Least-squares Adjustment for Empirical Affine Transformations[C]∥VI Hotine-Marussi Symposium on Theoretical and Computational Geodesy: Challenge and Role of Modern Geodesy. Berlin: Springer, 2008, 132: 238-242.

    [7]SCHAFFRIN B, FELUS Y A. On the Multivariate Total Least-squares Approach to Empirical Coordinate Transformations. Three Algorithms[J]. Journal of Geodesy, 2008, 82(6): 373-383.

    [8]SCHAFFRIN B, WIESER A. Empirical Affine Reference Frame Transformations by Weighted Multivariate TLS Adjustment[C]∥Geodetic Reference Frames. Berlin: Springer, 2009, 134: 213-218.

    [9]FELUS Y A, BURTCH R C. On Symmetrical Three-dimensional Datum Conversion[J]. GPS Solutions, 2009, 13(1): 65-74.

    [10]FANG X. Weighted Total Least Squares Solutions for Applications in Geodesy[D]. Hanover: Leibniz University of Hanover, 2011.

    [11]孫大雙, 張友陽, 黃令勇, 等. 多元總體最小二乘在大旋轉(zhuǎn)角三維坐標(biāo)轉(zhuǎn)換中的應(yīng)用[J]. 測(cè)繪科學(xué)技術(shù)學(xué)報(bào), 2014, 31(5): 481-485.

    SUN Dashuang, ZHANG Youyang, HUANG Lingyong, et al. Application of Multivariate Total Least Squares Method in Three-dimension Coordinate Conversion with Large Rotation Angle[J]. Journal of Geomatics Science and Technology, 2014, 31(5): 481-485.

    [12]黃令勇, 呂志平, 任雅奇, 等. 多元總體最小二乘在三維坐標(biāo)轉(zhuǎn)換中的應(yīng)用[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2014, 39(7): 793-798.

    HUANG Lingyong, Lü Zhiping, REN Yaqi, et al. Application of Multivariate Total Least Square in Three-dimensional Coordinate Transformation[J]. Geomatics and Information Science of Wuhan University, 2014, 39(7): 793-798.

    [13]錢承軍, 陳義. 多變量總體最小二乘在點(diǎn)云拼接中的應(yīng)用[J]. 測(cè)繪與空間地理信息, 2015, 38(1): 67-69, 76.

    QIAN Chengjun, CHEN Yi. Application of Multivariable Total Least Squares in the Registration of Point Clouds[J]. Geomatics & Spatial Information Technology, 2015, 38(1): 67-69, 76.

    [14]王樂洋, 許才軍, 魯鐵定. 邊長(zhǎng)變化反演應(yīng)變參數(shù)的總體最小二乘方法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(2): 181-184.

    WANG Leyang, XU Caijun, LU Tieding. Inversion of Strain Parameter Using Distance Changes Based on Total Least Squares[J]. Geomatics and Information Science of Wuhan University, 2010, 35(2): 181-184.

    [15]XU Caijun, WANG Leyang, WEN Yangmao, et al. Strain Rates in the Sichuan-Yunnan Region Based upon the Total Least Squares Heterogeneous Strain Model from GPS Data

    [J]. Terrestrial Atmospheric and Oceanic Sciences, 2011, 22(2): 133-147.

    [16]王樂洋, 于冬冬, 呂開云. 復(fù)數(shù)域總體最小二乘平差[J]. 測(cè)繪學(xué)報(bào), 2015, 44(8): 866-876. DOI: 10.11947/j.AGCS.2015.20130701.

    WANG Leyang, YU Dongdong, Lü Kaiyun. Complex Total Least Squares Adjustment[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(8): 866-876. DOI: 10.11947/j.AGCS.2015.20130701.

    [17]MAHBOUB V. On Weighted Total Least-squares for Geodetic Transformations[J]. Journal of Geodesy, 2012, 86(5): 359-367.

    [18]XU Peiliang,LIU Jingnan,SHI Chuang.Total Least Squares Adjustment in Partial Errors-in-variables Models: Algorithm and Statistical Analysis[J]. Journal of Geodesy, 2012, 86(8): 661-675.

    [19]JAZAERI S, AMIRI-SIMKOOEI A R, SHARIFI M A. Iterative Algorithm for Weighted Total Least Squares Adjustment[J]. Survey Review, 2014, 46(334): 19-27.

    [20]BOYD S, VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004.

    [21]王樂洋, 魯鐵定. 總體最小二乘平差法的誤差傳播定律[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2014, 34(2): 55-59.

    WANG Leyang, LU Tieding. Propagation Law of Errors in Total Least Squares Adjustment[J]. Journal of Geodesy and Geodynamics, 2014, 34(2): 55-59.

    [22]ZHOU Yongjun, KOU Xinjian, ZHU Jianjun, et al. A Newton Algorithm for Weighted Total Least-squares Solution to a Specific Errors-in-variables Model with Correlated Measurements[J]. Studia Geophysica et Geodaetica, 2014, 58(3): 349-375.

    [23]王樂洋, 許才軍, 魯鐵定. 病態(tài)加權(quán)總體最小二乘平差的嶺估計(jì)解法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(11): 1346-1350.

    WANG Leyang, XU Caijun, LU Tieding. Ridge Estimation Method in Ill-posed Weighted Total Least Squares Adjustment[J]. Geomatics and Information Science of Wuhan University, 2010, 35(11): 1346-1350.

    (責(zé)任編輯:叢樹平)

    修回日期: 2015-10-12

    First author: WANG Leyang (1983—), male, PhD, associate professor, majors in geodetic inversion and geodetic data processing.

    E-mail: wleyang@163.com

    A Newton Algorithm for Multivariate Total Least Squares Problems

    WANG Leyang1,2,3,ZHAO Yingwen1,CHEN Xiaoyong1,2,3,ZANG Deyan1,2,3

    1. Faculty of Geomatics, East China University of Technology, Nanchang 330013, China; 2. Key Laboratory of Watershed Ecology and Geographical Environment Monitoring, NASG,Nanchang 330013, China; 3. Jiangxi Province Key Lab for Digital Land, Nanchang 330013, China

    Abstract:In order to improve calculation efficiency of parameter estimation, an algorithm for multivariate weighted total least squares adjustment based on Newton method is derived. The relationship between the solution of this algorithm and that of multivariate weighted total least squares adjustment based on Lagrange multipliers method is analyzed. According to propagation of cofactor, 16 computational formulae of cofactor matrices of multivariate total least squares adjustment are also listed. The new algorithm could solve adjustment problems containing correlation between observation matrix and coefficient matrix. And it can also deal with their stochastic elements and deterministic elements with only one cofactor matrix. The results illustrate that the Newton algorithm for multivariate total least squares problems could be practiced and have higher convergence rate.

    Key words:multivariate total least squares; Newton method; propagation of cofactor; affine

    第一作者簡(jiǎn)介:王樂洋(1983—),男,博士,副教授,研究方向?yàn)榇蟮販y(cè)量反演及大地測(cè)量數(shù)據(jù)處理。

    收稿日期:2015-05-11

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(41204003;41161069;41304020;41464001);測(cè)繪地理信息公益性行業(yè)科研專項(xiàng)(201512026);江西省自然科學(xué)基金(20151BAB203042);江西省教育廳科技項(xiàng)目(GJJ150595;KJLD12077;KJLD14049);流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放基金(WE2015005);東華理工大學(xué)博士科研啟動(dòng)金(DHBK201113);江西省研究生創(chuàng)新專項(xiàng)資金(YC2015-S266;YC2015-S267);東華理工大學(xué)研究生創(chuàng)新專項(xiàng)資金(DHYC-2015005);對(duì)地觀測(cè)技術(shù)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放基金(K201502)

    中圖分類號(hào):P207

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1001-1595(2016)04-0411-07

    引文格式:王樂洋,趙英文,陳曉勇,等.多元總體最小二乘問題的牛頓解法[J].測(cè)繪學(xué)報(bào),2016,45(4):411-417,424. DOI:10.11947/j.AGCS.2016.20150246.

    WANG Leyang,ZHAO Yingwen,CHEN Xiaoyong,et al.A Newton Algorithm for Multivariate Total Least Squares Problems[J]. Acta Geodaetica et Cartographica Sinica,2016,45(4):411-417,424. DOI:10.11947/j.AGCS.2016.20150246.

    熟妇人妻久久中文字幕3abv| 亚洲自偷自拍三级| 国产精品.久久久| 国产精品伦人一区二区| 日日干狠狠操夜夜爽| 国产精品福利在线免费观看| 国产成人一区二区在线| 欧美高清性xxxxhd video| 插阴视频在线观看视频| 少妇被粗大猛烈的视频| ponron亚洲| 久久久欧美国产精品| 午夜久久久久精精品| 九九在线视频观看精品| 国产激情偷乱视频一区二区| 少妇人妻精品综合一区二区| 欧美xxxx性猛交bbbb| 水蜜桃什么品种好| 亚洲四区av| 观看美女的网站| 亚洲av福利一区| 国产精品一及| 亚洲美女视频黄频| 亚洲自拍偷在线| 美女大奶头视频| 亚洲人成网站在线播| 国产精品人妻久久久久久| 国产日韩欧美在线精品| 日本三级黄在线观看| 在线免费观看的www视频| 日韩欧美国产在线观看| 一级毛片我不卡| 中文精品一卡2卡3卡4更新| 97人妻精品一区二区三区麻豆| 精品酒店卫生间| 国产激情偷乱视频一区二区| 精品国产三级普通话版| 色综合亚洲欧美另类图片| 欧美日韩在线观看h| 大又大粗又爽又黄少妇毛片口| 国产精品一区www在线观看| 久久综合国产亚洲精品| 色综合亚洲欧美另类图片| 亚洲国产精品专区欧美| 嘟嘟电影网在线观看| 亚洲精品日韩av片在线观看| 美女大奶头视频| 日韩人妻高清精品专区| 蜜臀久久99精品久久宅男| 一区二区三区高清视频在线| 人妻夜夜爽99麻豆av| 亚洲成人中文字幕在线播放| 2018国产大陆天天弄谢| 国产黄色小视频在线观看| 欧美激情在线99| 国产精品久久久久久精品电影小说 | 久久6这里有精品| 亚洲欧美清纯卡通| 熟女人妻精品中文字幕| 亚洲精品成人av观看孕妇| 最后的刺客免费高清国语| 亚洲欧美一区二区三区国产| 成人亚洲精品av一区二区| 欧美极品一区二区三区四区| 18禁动态无遮挡网站| 干丝袜人妻中文字幕| 女的被弄到高潮叫床怎么办| 国产精品一二三区在线看| 麻豆av噜噜一区二区三区| 成年免费大片在线观看| 国产精品嫩草影院av在线观看| a级毛片免费高清观看在线播放| 国产精品一及| 精品久久久久久电影网| 草草在线视频免费看| 色综合亚洲欧美另类图片| 亚洲国产av新网站| 亚洲伊人久久精品综合| 亚洲成人精品中文字幕电影| 成人亚洲欧美一区二区av| 亚洲av中文字字幕乱码综合| 欧美日韩亚洲高清精品| 人妻系列 视频| 麻豆成人午夜福利视频| 在线观看美女被高潮喷水网站| 亚洲欧洲国产日韩| 久久久亚洲精品成人影院| 大陆偷拍与自拍| 搞女人的毛片| 国产日韩欧美在线精品| 九九爱精品视频在线观看| 少妇熟女欧美另类| 禁无遮挡网站| 久久97久久精品| 老司机影院成人| 日韩一区二区三区影片| 蜜桃亚洲精品一区二区三区| 午夜精品在线福利| 男人和女人高潮做爰伦理| 一本一本综合久久| 色综合站精品国产| 一级毛片久久久久久久久女| 国语对白做爰xxxⅹ性视频网站| 亚洲18禁久久av| 婷婷六月久久综合丁香| 欧美日韩精品成人综合77777| 婷婷色麻豆天堂久久| 国产午夜福利久久久久久| 尾随美女入室| 网址你懂的国产日韩在线| 人妻夜夜爽99麻豆av| 亚洲精品aⅴ在线观看| 国产免费一级a男人的天堂| 欧美日韩综合久久久久久| 欧美日韩一区二区视频在线观看视频在线 | 婷婷色麻豆天堂久久| 亚洲欧洲国产日韩| 国产三级在线视频| 久久久久国产网址| 日本免费在线观看一区| 国产精品国产三级专区第一集| 天天躁夜夜躁狠狠久久av| 亚洲精品视频女| 国产女主播在线喷水免费视频网站 | 免费观看在线日韩| 精品久久久久久久人妻蜜臀av| 久久久久性生活片| 毛片女人毛片| 亚州av有码| 国产黄色免费在线视频| 中国美白少妇内射xxxbb| 久久久久免费精品人妻一区二区| 中文资源天堂在线| 又黄又爽又刺激的免费视频.| av福利片在线观看| 午夜福利在线观看吧| av一本久久久久| 一级片'在线观看视频| 亚洲国产精品专区欧美| 亚洲欧美成人综合另类久久久| 两个人的视频大全免费| 国产精品.久久久| 黄色一级大片看看| 国产精品99久久久久久久久| 亚洲在久久综合| 一本一本综合久久| 天美传媒精品一区二区| 一区二区三区乱码不卡18| 三级毛片av免费| 日韩亚洲欧美综合| 观看美女的网站| 国产一区二区亚洲精品在线观看| 夫妻午夜视频| 丝袜美腿在线中文| 成人性生交大片免费视频hd| 久久草成人影院| 大又大粗又爽又黄少妇毛片口| 美女被艹到高潮喷水动态| 可以在线观看毛片的网站| 国产精品麻豆人妻色哟哟久久 | 久久热精品热| 久久99蜜桃精品久久| 亚洲av.av天堂| 欧美日韩视频高清一区二区三区二| 亚洲精品影视一区二区三区av| 日韩亚洲欧美综合| 一级黄片播放器| 久久精品久久久久久久性| 在线观看人妻少妇| 精品久久久久久久人妻蜜臀av| 天堂av国产一区二区熟女人妻| 国产成人午夜福利电影在线观看| 老司机影院毛片| 91aial.com中文字幕在线观看| 久久久久久久久大av| 搡老乐熟女国产| 青春草亚洲视频在线观看| 女人久久www免费人成看片| 伦精品一区二区三区| 少妇熟女aⅴ在线视频| 老司机影院毛片| 黄片无遮挡物在线观看| 国产av不卡久久| 日韩精品有码人妻一区| 日韩不卡一区二区三区视频在线| 国产精品无大码| 亚洲av成人av| 一本一本综合久久| 韩国高清视频一区二区三区| 欧美高清性xxxxhd video| 国产精品久久久久久久久免| 免费在线观看成人毛片| 亚洲人与动物交配视频| 性色avwww在线观看| 国内精品宾馆在线| 成人av在线播放网站| 亚洲国产色片| 非洲黑人性xxxx精品又粗又长| 大香蕉97超碰在线| 人妻少妇偷人精品九色| 麻豆乱淫一区二区| 女人被狂操c到高潮| 亚洲av电影不卡..在线观看| 国产乱来视频区| 国产亚洲一区二区精品| 如何舔出高潮| 高清在线视频一区二区三区| 男女那种视频在线观看| 亚洲av在线观看美女高潮| 免费黄频网站在线观看国产| 乱码一卡2卡4卡精品| 亚洲美女搞黄在线观看| 久久精品熟女亚洲av麻豆精品 | 男女啪啪激烈高潮av片| 国产成年人精品一区二区| 美女国产视频在线观看| 午夜福利高清视频| 久久久久久久久中文| 久久久久免费精品人妻一区二区| 国产精品一区二区在线观看99 | 免费av毛片视频| 亚洲熟妇中文字幕五十中出| 热99在线观看视频| 久久人人爽人人爽人人片va| 国产高清国产精品国产三级 | 日日啪夜夜爽| 国产精品99久久久久久久久| 欧美精品一区二区大全| 国内揄拍国产精品人妻在线| 免费av毛片视频| 深爱激情五月婷婷| 久99久视频精品免费| 在线 av 中文字幕| 亚洲婷婷狠狠爱综合网| 亚洲精品成人av观看孕妇| 亚洲人成网站高清观看| 国产伦理片在线播放av一区| 22中文网久久字幕| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 一本久久精品| 最近2019中文字幕mv第一页| 99久国产av精品| 男女那种视频在线观看| 国产淫语在线视频| 国产午夜精品论理片| 亚洲最大成人中文| 三级男女做爰猛烈吃奶摸视频| 精品国产一区二区三区久久久樱花 | 老师上课跳d突然被开到最大视频| 秋霞伦理黄片| 国产精品一区二区三区四区免费观看| 亚洲av二区三区四区| 午夜日本视频在线| 中文精品一卡2卡3卡4更新| 日本一二三区视频观看| 亚洲国产成人一精品久久久| 联通29元200g的流量卡| 十八禁网站网址无遮挡 | 男女那种视频在线观看| 看十八女毛片水多多多| 少妇被粗大猛烈的视频| 麻豆国产97在线/欧美| 高清av免费在线| 好男人在线观看高清免费视频| 国产精品久久久久久久久免| 日韩,欧美,国产一区二区三区| 日本免费在线观看一区| 久久精品国产自在天天线| 色吧在线观看| 国产黄a三级三级三级人| 国产亚洲午夜精品一区二区久久 | 国产欧美另类精品又又久久亚洲欧美| 亚洲av免费在线观看| 国产精品av视频在线免费观看| 亚洲av成人av| 嫩草影院入口| 啦啦啦啦在线视频资源| 天堂网av新在线| 欧美不卡视频在线免费观看| 亚洲精品日韩在线中文字幕| 少妇被粗大猛烈的视频| 最近中文字幕高清免费大全6| 亚洲欧美中文字幕日韩二区| 国产91av在线免费观看| 成年女人在线观看亚洲视频 | 国产视频内射| 亚洲美女搞黄在线观看| 久久久久久久久中文| 岛国毛片在线播放| 久久久亚洲精品成人影院| 亚洲伊人久久精品综合| 青春草亚洲视频在线观看| 毛片女人毛片| 成人二区视频| 国产伦精品一区二区三区四那| 国产成人午夜福利电影在线观看| 少妇高潮的动态图| 亚洲精品日韩在线中文字幕| av免费在线看不卡| 综合色丁香网| av女优亚洲男人天堂| 直男gayav资源| 国产精品99久久久久久久久| 两个人的视频大全免费| 精华霜和精华液先用哪个| 国产精品人妻久久久影院| 国产黄片视频在线免费观看| 国产成人aa在线观看| 日本三级黄在线观看| 亚洲欧美成人精品一区二区| 人体艺术视频欧美日本| 亚洲真实伦在线观看| 午夜激情久久久久久久| 精品欧美国产一区二区三| 亚洲天堂国产精品一区在线| 少妇熟女aⅴ在线视频| 国内精品美女久久久久久| 精品久久久久久久久亚洲| 青青草视频在线视频观看| 免费看日本二区| 国产伦在线观看视频一区| 汤姆久久久久久久影院中文字幕 | 亚洲欧美清纯卡通| 国产v大片淫在线免费观看| 亚洲国产成人一精品久久久| 欧美+日韩+精品| 亚洲真实伦在线观看| 国产成人免费观看mmmm| 成人无遮挡网站| 国产乱人偷精品视频| 美女脱内裤让男人舔精品视频| 久久久久精品久久久久真实原创| 在线观看美女被高潮喷水网站| 国产亚洲一区二区精品| 看免费成人av毛片| 2021天堂中文幕一二区在线观| 色哟哟·www| a级毛色黄片| 大香蕉久久网| 久久人人爽人人片av| 赤兔流量卡办理| 神马国产精品三级电影在线观看| 亚洲av.av天堂| 插阴视频在线观看视频| 精品久久久噜噜| 亚洲精品视频女| 免费大片18禁| 亚洲av不卡在线观看| 看十八女毛片水多多多| 国产伦精品一区二区三区视频9| 亚洲av电影在线观看一区二区三区 | 大又大粗又爽又黄少妇毛片口| 极品少妇高潮喷水抽搐| 国产免费一级a男人的天堂| 91久久精品国产一区二区三区| 99热全是精品| 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 亚洲精品国产成人久久av| 免费黄色在线免费观看| 麻豆成人av视频| 国产精品av视频在线免费观看| 精品久久久久久久末码| 99久国产av精品| 波多野结衣巨乳人妻| 插逼视频在线观看| 成人国产麻豆网| 一二三四中文在线观看免费高清| 免费黄色在线免费观看| 久久久精品免费免费高清| 国产高潮美女av| 九色成人免费人妻av| 插阴视频在线观看视频| 久久久久久久久久久丰满| 白带黄色成豆腐渣| 国产精品福利在线免费观看| 国产高清国产精品国产三级 | 特级一级黄色大片| 日韩av在线免费看完整版不卡| 精品久久久久久电影网| 久久午夜福利片| 青春草视频在线免费观看| 国产极品天堂在线| 一个人免费在线观看电影| 91精品国产九色| 中国美白少妇内射xxxbb| 亚洲精品,欧美精品| 国产黄色小视频在线观看| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 日韩av在线免费看完整版不卡| 综合色av麻豆| 一级av片app| 精品国产露脸久久av麻豆 | 2021天堂中文幕一二区在线观| 一级毛片电影观看| 一级av片app| 岛国毛片在线播放| 爱豆传媒免费全集在线观看| 五月天丁香电影| 高清视频免费观看一区二区 | 啦啦啦啦在线视频资源| 老司机影院毛片| 搡老乐熟女国产| 日韩欧美国产在线观看| 成年女人看的毛片在线观看| 亚洲av成人av| 成年免费大片在线观看| 欧美不卡视频在线免费观看| 一边亲一边摸免费视频| 久久久久性生活片| 黄色欧美视频在线观看| 黄片无遮挡物在线观看| 日本一二三区视频观看| 成人午夜精彩视频在线观看| 亚洲美女视频黄频| 五月天丁香电影| 免费av不卡在线播放| 床上黄色一级片| av女优亚洲男人天堂| 日韩成人伦理影院| 黄片无遮挡物在线观看| 在线观看免费高清a一片| 午夜久久久久精精品| 人体艺术视频欧美日本| 亚洲18禁久久av| 简卡轻食公司| 中文资源天堂在线| av在线观看视频网站免费| 嫩草影院新地址| 亚洲图色成人| 人体艺术视频欧美日本| 五月玫瑰六月丁香| 亚洲va在线va天堂va国产| 精品久久久噜噜| 亚洲精品视频女| 69人妻影院| 亚洲天堂国产精品一区在线| 亚洲电影在线观看av| 中文精品一卡2卡3卡4更新| 午夜激情欧美在线| 亚洲最大成人中文| 国内揄拍国产精品人妻在线| 午夜福利在线观看吧| 80岁老熟妇乱子伦牲交| 国产精品福利在线免费观看| 免费观看无遮挡的男女| 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| 久久韩国三级中文字幕| 国产伦精品一区二区三区四那| 婷婷色综合大香蕉| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 免费高清在线观看视频在线观看| av在线播放精品| 国产激情偷乱视频一区二区| 国产一区二区在线观看日韩| 国产精品国产三级国产专区5o| 日韩 亚洲 欧美在线| 国产探花在线观看一区二区| 精品熟女少妇av免费看| 国产精品av视频在线免费观看| 亚洲av免费高清在线观看| 伦精品一区二区三区| 99久久中文字幕三级久久日本| 一级毛片我不卡| 18+在线观看网站| 国产精品一及| 免费黄色在线免费观看| 网址你懂的国产日韩在线| 激情五月婷婷亚洲| av在线老鸭窝| 国产黄频视频在线观看| 好男人视频免费观看在线| 99久久精品一区二区三区| 人人妻人人澡欧美一区二区| 综合色丁香网| 一级爰片在线观看| 亚洲乱码一区二区免费版| 丝袜美腿在线中文| 大又大粗又爽又黄少妇毛片口| 日本与韩国留学比较| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 久久精品国产自在天天线| 人妻系列 视频| 三级国产精品欧美在线观看| 色播亚洲综合网| 久久亚洲国产成人精品v| 天堂√8在线中文| 欧美 日韩 精品 国产| 精品一区二区三区人妻视频| 亚洲av一区综合| 国产成人精品一,二区| 精华霜和精华液先用哪个| 久久国内精品自在自线图片| 免费看光身美女| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 神马国产精品三级电影在线观看| 亚洲精品影视一区二区三区av| 黄色欧美视频在线观看| 深夜a级毛片| 亚洲成人久久爱视频| 午夜视频国产福利| 国产精品久久久久久精品电影| 国产老妇伦熟女老妇高清| 精品国产露脸久久av麻豆 | 激情五月婷婷亚洲| 国产 一区 欧美 日韩| 男人舔女人下体高潮全视频| 亚洲成人久久爱视频| 最近中文字幕2019免费版| 亚洲精品国产成人久久av| 又爽又黄a免费视频| 国产av在哪里看| 美女cb高潮喷水在线观看| 国产精品1区2区在线观看.| 男女国产视频网站| av黄色大香蕉| 亚洲熟女精品中文字幕| 人妻少妇偷人精品九色| 高清午夜精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 日韩欧美一区视频在线观看 | 亚洲成人精品中文字幕电影| 毛片一级片免费看久久久久| 国产精品不卡视频一区二区| 国产一区亚洲一区在线观看| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 国产精品久久久久久久久免| 91av网一区二区| 99热6这里只有精品| 麻豆久久精品国产亚洲av| videossex国产| 网址你懂的国产日韩在线| 丝袜美腿在线中文| 黄色一级大片看看| 日日干狠狠操夜夜爽| 国产精品不卡视频一区二区| 欧美另类一区| 精品一区二区三卡| 国产精品一及| 国产精品蜜桃在线观看| 久久99蜜桃精品久久| 久久99热6这里只有精品| 午夜福利在线在线| 国产成人aa在线观看| 寂寞人妻少妇视频99o| 99久国产av精品国产电影| 九九在线视频观看精品| 日韩av在线大香蕉| 日韩在线高清观看一区二区三区| 久久久久久久久大av| 国产91av在线免费观看| 国产精品无大码| 久久久久性生活片| 寂寞人妻少妇视频99o| 午夜老司机福利剧场| 免费大片黄手机在线观看| 深爱激情五月婷婷| 欧美变态另类bdsm刘玥| 色吧在线观看| 日日啪夜夜爽| 久久精品国产自在天天线| 女人被狂操c到高潮| 蜜臀久久99精品久久宅男| 国产成人aa在线观看| 国产黄片美女视频| 美女xxoo啪啪120秒动态图| 国产午夜精品一二区理论片| 高清av免费在线| 黄色欧美视频在线观看| 日韩欧美三级三区| 国产亚洲av片在线观看秒播厂 | 97超视频在线观看视频| 日本三级黄在线观看| 国产三级在线视频| 国产黄色视频一区二区在线观看| 全区人妻精品视频| 91精品伊人久久大香线蕉| 大香蕉久久网| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 精品一区在线观看国产| 国产又色又爽无遮挡免| 久久久精品欧美日韩精品| 午夜视频国产福利| 日韩强制内射视频| 91狼人影院| 亚洲aⅴ乱码一区二区在线播放| 国产伦一二天堂av在线观看| 国国产精品蜜臀av免费| 免费少妇av软件| 亚洲18禁久久av| 免费黄频网站在线观看国产| 99久久中文字幕三级久久日本| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| 在线观看人妻少妇| 中文精品一卡2卡3卡4更新| 秋霞伦理黄片| 成人漫画全彩无遮挡| 高清视频免费观看一区二区 | 免费黄频网站在线观看国产| 亚洲内射少妇av| 亚洲成人久久爱视频| 两个人的视频大全免费| 成人无遮挡网站| 国产成人精品久久久久久| 简卡轻食公司| 蜜臀久久99精品久久宅男|