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

    依不同緯度變量的子午線弧長正反解公式的級數(shù)展開

    2016-09-02 06:05:45過家春李厚樸莊云玲李大軍吳艷蘭
    測繪學報 2016年5期

    過家春,李厚樸,莊云玲,李大軍,吳艷蘭

    1. 安徽農業(yè)大學理學院,安徽 合肥 230036; 2. 東華理工大學江西省數(shù)字國土重點實驗室,江西 南昌 330013; 3. 安徽大學資源與環(huán)境工程學院,安徽 合肥 230601; 4. 海軍工程大學導航工程系,湖北 武漢 430033

    ?

    依不同緯度變量的子午線弧長正反解公式的級數(shù)展開

    過家春1, 2, 3,李厚樸4,莊云玲1,李大軍2,吳艷蘭3

    1. 安徽農業(yè)大學理學院,安徽 合肥 230036; 2. 東華理工大學江西省數(shù)字國土重點實驗室,江西 南昌 330013; 3. 安徽大學資源與環(huán)境工程學院,安徽 合肥 230601; 4. 海軍工程大學導航工程系,湖北 武漢 430033

    Foundationsupport:TheNationalNaturalScienceFoundationofChina(Nos. 41504031; 41571441);FoundationofJiangxiProvinceKeyLabforDigitalLand(No.DLLJ201507)

    推導了以歸化緯度、地心緯度解算子午線弧長的展開公式,同時又根據(jù)拉格朗日反演定理,得到了由子午線弧長反解歸化緯度、地心緯度的直接公式。該組公式與子午線弧長正反解公式的大地緯度表達在結構形式上保持一致,進一步揭示了子午線弧長同3種緯度變量之間的內在聯(lián)系。分析表明,基于歸化緯度的子午線弧長解算與大地主題解算方法具有理論上的統(tǒng)一性,正反解精度均高于傳統(tǒng)基于大地緯度的展開。

    子午線弧長;級數(shù)展開;大地緯度;歸化緯度;地心緯度

    子午線弧長計算是經(jīng)典大地測量問題之一[1-12],圍繞這一問題的計算和應用,近年來各國學者提出了許多新的方法和見解[13-16]。過去,人們主要關注子午線弧長與大地緯度的關系,但在測量和地圖投影的相關理論中,也常用到歸化緯度、地心緯度等輔助緯度,以使問題得到簡化,尤其是應用歸化緯度。例如,文獻[13]曾給出分別基于大地緯度和歸化緯度的子午線弧長公式正反解的復變函數(shù)表達,并分析指出以歸化緯度為自變量的子午線弧長公式更為簡潔,且精度更高;文獻[14—16]中給出的子午線弧長公式與第二類橢圓積分的關系表明以歸化緯度表示的子午線弧長公式表達更為簡潔;文獻[17—18]分析了多種輔助緯度之間的關系,并指出其實用意義和價值。在高斯投影、UTM投影等正形投影中,應用復變函數(shù)和橢圓積分的表達與解算方法,也常引入歸化緯度[19-22],在空間直角坐標與大地坐標的轉換中的一些方法中也常使用到歸化緯度[19-22]。同時,在貝塞爾大地主題解算方法及近年來提出的許多新的方法中,均以歸化緯度作為過渡或直接解算結果[23-24]。事實上,子午線是特殊的大地線,分析子午線弧長與歸化緯度、地心緯度的正反解關系,有助于建立子午線弧長與大地主題解算理論上的聯(lián)系。

    另外,子午線弧長問題的本質是橢圓弧長問題。而在其他學科領域,一般更多關注的是橢圓弧長與其離心角、極角的關系[25-26],即對應于大地測量中的歸化緯度和球心緯度。

    鑒于此,本文將子午線弧長正、反解公式分別以大地緯度、歸化緯度和地心緯度表示,實現(xiàn)3種緯度變量與子午線弧長的直接轉換,并給出理論與精度分析,以資參考。

    1 子午線弧長正解公式的不同緯度變量表達

    1.1以大地緯度為自變量的子午線弧長公式

    以大地緯度為自變量的子午線弧長正解公式為

    B′sin 2B+C′sin 4B-D′sin 6B+E′sin 8B)

    (1)

    式中,各系數(shù)參考文獻[5]。文獻[16]通過引入第三扁率n和高斯超幾何函數(shù)F將其化為

    (2)

    該公式使子午線弧長公式得到簡化,且精度提高顯著。式中

    (3)

    1.2以歸化緯度為自變量的子午線弧長公式

    子午線弧長公式以歸化緯度u為自變量的公式為

    (4)

    按文獻[16]的思路,先將被積函數(shù)按二項式定理展開,再積分,然后引入第三扁率n和高斯超幾何函數(shù)F,可得

    (5)

    式(5)的推導過程本文從略。為避免復雜的推導過程,可在數(shù)學軟件Mathematica中輸入如下命令一步實現(xiàn)

    (6)

    式(6)中各命令功能參見Mathematica手冊[27]。

    1.3以地心緯度為自變量的子午線弧長公式

    以地心緯度Φ為自變量的子午線弧長公式為

    (7)

    按同樣思路展開可得

    (8)

    2 由子午線弧長反解不同的緯度變量

    2.1子午線弧長反解大地緯度的Helmert公式

    Helmert于1880年應用拉格朗日反演定理[28-29]給出子午弧長反解的直接公式,也常稱之為三角級數(shù)回代方法,本文直接引用如下

    (9)

    式中,β=S/aF。

    2.2子午線弧長反解歸化緯度的直接公式

    (10)

    根據(jù)拉格朗日反演定理,有

    (11)

    用β替換式(10)中的u,根據(jù)式(11)可得

    (12)

    (13)

    (14)

    依次類推,求至5階導,并代入式(11),可得

    (15)

    此即為由子午線弧長反解歸化緯度的直接公式。

    2.3子午線弧長反解地心緯度的直接公式

    類似的,略去推導過程,可以得到由子午線弧長反解地心緯度的直接公式

    (16)

    至此,得到了分別基于大地緯度、歸化緯度和地心緯度的子午線弧長正解公式,式(2)、式(5)和式(8),以及由子午線弧長反解各緯度變量的直接公式,式(9)、式(15)和式(16)。3組公式在結構形式上保持較高的一致性,可實現(xiàn)子午線弧長與3種不同緯度相互間直接的轉換計算。

    3 基于歸化緯度的子午線弧長解算與大地主題解算方法的統(tǒng)一性分析

    上述3種緯度變量的子午線弧長公式的表達中,以大地緯度最為常用,其他兩種表達應用較少。但基于歸化緯度的子午線弧長理論應用上有許多優(yōu)點,主要表現(xiàn)在:①形式上更為簡潔,本文式(1)、式(4)、式(7)表明以歸化緯度的積分表達更為簡潔,與第二類橢圓積分的關系也更為明顯;②精度更高,文獻[13]以復變函數(shù)理論解決子午線弧長問題,分析表明,基于大地緯度的子午線弧長計算最弱精度在7.0×10-6m,而以歸化緯度的最弱精度約為5.0×10-7m,且公式更為簡潔,本文在下一節(jié)內容中分析本文方法的精度比較。

    除以上兩方面優(yōu)勢以外,基于歸化緯度的子午線弧長與大地主題解算相關方法具有理論上的統(tǒng)一性。

    在文獻[22—24]中,學者給出了新的大地主題解算方法,用歸化緯度表達大地線的微分方程作為解算的基礎方程

    (17)

    式中,C=cosumax=cosuisinAi,稱為大地線常數(shù)或克萊羅常數(shù)(Clairault’s constant)。

    由式(17)可得

    (18)

    特別的,當大地方位角A=0°時,則C=0,該大地線為子午線,式(18)即化為式(4)。在此種情況下,文獻[24]給出按下式以迭代算法計算歸化緯度

    (19)

    式中

    (20)

    文獻[22—24]給出的以上過程表明了歸化緯度在大地主題解算中的應用價值,而子午線是大地線的一種特殊形式。本文分析的基于歸化緯度的子午線弧長正反解建立了子午線弧長理論與大地主題解算相關方法的聯(lián)系,表明二者具有理論上的統(tǒng)一性。

    4 精度驗算分析

    為比較3組公式的精度,筆者以WGS-84橢球參數(shù)為例,在Mathematica 8.0中調用第2類橢圓積分函數(shù)得到各公式分別展開至n2、n3、n4項的最大誤差,計算結果如表1所示,正反解誤差曲線如圖1所示。

    圖1 基于不同緯度變量的子午線弧長正反解誤差曲線Fig.1 Error curves of direct and inverse solutions for meridian by using different latitude variables

    驗算表明:

    (1) 正解方面,展開至n2項以大地緯度為自變量的解算精度略低,地心緯度次之;展開至n3、n4項,基于地心緯度的解算精度略低,大地緯度次之;展開至n2~n4各階次項均以歸化緯度的解算精度最好。

    (2) 反解方面,展開至n2項反解歸化緯度的精度略高于地心緯度;展開至n3、n4項反解地心緯度的精度略高于歸化緯度;展開至n2~n4各階次項反解大地緯度的精度略低;

    (3) 實用上來講,各正反解公式展開至n4項能夠滿足大地測量領域的應用精度要求,并可根據(jù)需要展開至更高階以獲得更高精度。

    (4) 對比3組公式可發(fā)現(xiàn),以歸化緯度為變量的正解公式最為簡潔。事實上,如前所述,在大地測量領域里的許多方面,歸化緯度都有類似的性質,常使復雜的問題得到簡化。

    5 結 論

    本文按以大地緯度為自變量的子午線弧長公式思路,推導了以歸化緯度、地心緯度解算子午線弧長的展開公式,同時又按拉格朗日反演定理,得到了由子午線弧長反解歸化緯度、地心緯度的直接公式,并分析了子午線弧長理論與相關大地主題解算方法的統(tǒng)一性?;?種不同緯度變量的子午線弧長正反解公式結構、形式上一致,收斂速度也基本相當。該組公式實現(xiàn)了子午線弧長同3種緯度變量的直接正反解,揭示了子午線弧長同3種緯度變量之間的關系。子午線弧長理論與大地主題解算理論的統(tǒng)一性分析表明,歸化緯度在大地主題解算方面有其應用價值,而子午線是大地線的一種特殊形式,二者表現(xiàn)為一般與特殊的關系,具有理論上的統(tǒng)一性。

    [1]TORGE W. Geodesy[M]. 3rd ed. Berlin: Walter de Gruyter, 2001: 91-98.

    [2]孔祥元, 郭際明, 劉宗泉. 大地測量學基礎[M]. 武漢: 武漢大學出版社, 2001: 64-73.

    KONG Xiangyuan, GUO Jiming, LIU Zongquan. Foundation of Geodesy[M]. Wuhan: Wuhan University Press, 2001: 64-73.

    [3]HELMERT F R. Die Mathematischen und Physikalischen Theorien der H?heren Geod?sie [M]. Leipzig: Druck und Verlag von B.G. Teubner, 1880: 46-48.

    [4]DEAKIN R E, HUNTER M N. Geometric Geodesy: Part A[R]. Melbourne: School of Mathematical & Geospatial Science, RMIT University, 2010: 60-77.

    [5]程鵬飛, 成英燕, 文漢江, 等. 2000國家大地坐標系實用寶典[M]. 北京: 測繪出版社, 2008: 147-148. CHENG Pengfei, CHENG Yingyan, WEN Hanjiang, et al. Practical Manual on CGCS 2000[M]. Beijing: Surveying and Mapping Press, 2008: 147-148.

    [6]劉正才. 子午線弧長公式的簡化及通用高斯投影計算程序介紹[J]. 測繪工程, 2001, 10(1): 55-56, 62.

    LIU Zhengcai. Simplification of Formula of Meridian Arc Length & Program of General Gauss Projection[J]. Engineering of Surveying and Mapping, 2001, 10(1): 55-56, 62.

    [7]易維勇, 邊少鋒, 朱漢泉. 子午線弧長的解析型冪級數(shù)確定[J]. 測繪學院學報, 2000, 17(3): 167-171. YI Weiyong, BIAN Shaofeng, ZHU Hanquan. Determination of Foot Point Latitude by Analytic Positive Serires[J]. Journal of Institute of Surveying and Mapping, 2000, 17(3): 167-171.

    [8]劉仁釗, 伍吉倉. 任意精度的子午線弧長遞歸計算[J]. 大地測量與地球動力學, 2007, 27(5): 59-62. LIU Renzhao, WU Jicang. Recursive Computation of Meridian Arc Length with Discretionary Precision[J]. Journal of Geodesy and Geodynamics, 2007, 27(5): 59-62.

    [9]BIAN S F, CHEN Y B. Solving an Inverse Problem of a Meridian Arc in Terms of Computer Algebra System[J]. Journal of Surveying Engineering, 2006, 132(1): 7-10.

    [10]牛卓立. 以空間直角坐標為參數(shù)的子午線弧長計算公式[J]. 測繪通報, 2001(11): 14-15. DOI:10.3969/j.issn.0494-0911.2001.11.006.

    NIU Zhuoli. Formulae for Calculation of Meridian Arc Length by the Parameters of Space Rectangular Coordinates[J]. Bulletin of Surveying and Mapping, 2001(11): 14-15. DOI: 10.3969/j.issn.0494-0911.2001.11.006.

    [11]劉修善. 計算子午線弧長的數(shù)值積分法[J]. 測繪通報, 2006(5): 4-6. DOI:10.3969/j.issn.0494-0911.2006.05.002.

    LIU Xiushan. Numerical Integral Method of Calculating Meridian Arc Length[J]. Bulletin of Surveying and Mapping, 2006(5):4-6. DOI:10.3969/j.issn.0494-0911.2006.05.002.

    [12]BOWRING B R, New Equations for Meridional Distance[J]. Bulletin Géodésique, 1983, 57(1-4): 374-381.

    [13]KAWASE K. A General Formula for Calculating Meridian Arc Length and Its Application to Coordinate Conversion in the Gauss-Krüger Projection[J]. Bulletin of the Geospatial Information Authority of Japan, 2011, 59: 1-13.

    [14]過家春, 趙秀俠, 徐麗, 等. 基于第二類橢圓積分的子午線弧長公式變換及解算[J]. 大地測量與地球動力學, 2011, 31(4): 94-98.

    GUO Jiachun, ZHAO Xiuxia, XU Li, et al. Calculating Meridian Arc Length by Transforming Its Formulae into Elliptic Integral of Second Kind[J]. Journal of Geodesy and Geodynamics, 2011, 31(4): 94-98.

    [15]過家春. 基于第二類橢圓積分的子午線弧長反解新方法[J]. 大地測量與地球動力學, 2012, 32(3): 116-120. GUO Jiachun. New Method for Inverse Solution of Meridian Based on Elliptic Integral of Second Kind[J]. Journal of Geodesy and Geodynamics, 2012, 32(3): 116-120.

    [16]過家春. 子午線弧長公式的簡化及其泰勒級數(shù)解釋[J]. 測繪學報, 2014, 43(2): 125-130. DOI: 10.13485/j.cnki.11-2089.2014.0017.

    GUO Jiachun. A Simplification of the Meridian Formula and Its Taylor-series Interpretation[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(2): 125-130. DOI: 10.13485/j.cnki.11-2089.2014.0017.

    [17]李忠美, 李厚樸, 邊少鋒. 常用緯度差異極值符號表達式[J]. 測繪學報, 2014, 43(2): 214-220. DOI: 10.13485/j.cnki.11-2089.2014.0031.LI Zhongmei, LI Houpu, BIAN Shaofeng. Symbolic Expressions of Difference Extrema between Regular Latitudes[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(2): 214-220. DOI: 10.13485/j.cnki.11-2089.2014.0031.

    [18]孫群, 楊啟和. 底點緯度解算以及等量緯度和面積函數(shù)反解問題的探討[J]. 解放軍測繪學院學報, 1985(2): 64-75. SUN Qun, YANG Qihe. The Research on the Computation of the Foot-point Latitude and the Inverse Solution of Isometric Latitude and Area Function[J]. Journal of PLA of Surveying and Mapping, 1985(2): 64-75.

    [19]李厚樸, 邊少鋒. 高斯投影的復變函數(shù)表示[J]. 測繪學報, 2008, 37(1): 5-9. DOI: 10.3321/j.issn:1001-1595.2008.01.002.

    LI Houpu, BIAN Shaofeng. The Expressions of Gauss Projection by Complex Numbers[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(1): 5-9. DOI: 10.3321/j.issn:1001-1595.2008.01.002.

    [20]DORRER E. From Elliptic Arc Length to Gauss-Krüger Coordinates by Analytical Continuation[M]∥GRAFAREND E W, KRUMM F W, SCHWARZE V S. Geodesy-The Challenge of the 3rd Millennium. Berlin: Springer, 2003: 293-298.

    [21]LEE L P. Some Conformal Projections Based on Elliptic Functions[J]. Geographical Review, 1965, 55(4): 563-580.

    [22]BERMEJO-SOLERA M, OTERO J. Simple and Highly Accurate Formulas for the Computation of Transverse Mercator Coordinates from Longitude and Isometric Latitude[J]. Journal of Geodesy, 2009, 83(1): 1-12.

    [23]SJ?BERG L E. New Solutions to the Direct and Indirect Geodetic Problems on the Ellipsoid[J]. Zeitschriftfuer Vermessungswesen, 2006, 131: 35-39.

    [24]SJ?BERG L E, SHIRAZIAN M. Solving the Direct and Inverse Geodetic Problems on the Ellipsoid by Numerical Integration[J]. Journal of Surveying Engineering, 2012, 138(1): 9-16.

    [25]張彥博. 等弧長橢圓時間分割插補算法[J]. 機床與液壓, 2005(7): 41, 110. ZHANG Yanbo. Equal Arc Length Time Dividing Interpolation of Ellipse[J]. Machine Tool & Hydraulics, 2005(7): 41, 110.

    [26]劉有余, 韓江, 夏鏈 等. 外嚙合橢圓齒輪多方案插齒三維仿真與分析[J]. 機械科學與技術, 2014, 33(7): 1031-1035.

    LIU Youyu, HAN Jiang, XIA Lian, et al. 3D-simulation and Analysis of Multiple Slotting External Elliptic Gears[J]. Mechanical Science and Technology for Aerospace Engineering, 2014, 33(7): 1031-1035.

    [27]WOLFRAM S. The Mathematica Book[M]. 5th ed. Champaign: Wolfram Media Inc., 2003.

    [28]ABRAMOWITZ M, STEGUN I A. Handbook of Mathematical Functions, with Formulas, Graphs, and Mathematical Tables[M]. New York: Dover Publications Inc., 1972: 14-16.

    [29]WEISSTEIN E W, “Series Reversion” from MathWorld: A Wolfram Web Resource[DB/OL]. [2015-05-14]. http:∥mathworld.wolfram.com/SeriesReversion.html.

    (責任編輯:叢樹平)

    GUO Jiachun(1981—), male, PhD candidate, associate professor, majors in geodesy.

    SeriesExpansionforDirectandInverseSolutionsofMeridianinTermsofDifferentLatitudeVariables

    GUOJiachun1, 2, 3,LIHoupu4,ZHUANGYunling1,LIDajun2,WUYanlan3

    1.SchoolofScience,AnhuiAgriculturalUniversity,Hefei230036,China; 2.JiangxiProvinceKeyLabforDigitalLand,EastChinaInstituteofTechnology,Nanchang330013,China; 3.SchoolofResourcesandEnvironmentalEngineering,AnhuiUniversity,Hefei230601,China; 4.DepartmentofNavigation,NavalUniversityofEngineering,Wuhan430033,China

    Formulasfordirectsolutionsofmeridianwrittenbythereducedandgeocentriclatitudesrespectivelywerederivedbyseriesexpansion.Meanwhile,accordingtoLagrangeinversiontheorem,formulasforinversesolutionsoftheissuewerealsoexpressedintermsofthesamelatitudes.Thesetwoformulaswerestructurallyconsistentwiththatexpressedbygeodeticlatitudeones.Inthesesetsofformulas,internalconnectionbetweenmeridianandthreedifferenttypesoflatitudewererealized.Analysisandnumericalcalculationshowedthatthedirectandinversemeridionalsolutionwithreducedlatitudewasofhigherprecisionthanthatwithgeodeticlatitude,andfurthermore,therehadaunifiedtheorybetweenmeridiantheoryandclassicalgeodeticproblemsexpressedbyreducedlatitude.

    meridianarclength;seriesexpansion;geodeticlatitude;reducedlatitude;geocentriclatitude

    2014-11-05

    2015-05-15

    過家春(1981—),男,博士生,副教授,研究方向為大地測量學。

    E-mail: guojiachun@ahau.edu.cn

    吳艷蘭

    WU Yanlan

    E-mail: wylmq@sina.com

    GUOJiachun,LIHoupu,ZHUANGYunling,etal.SeriesExpansionforDirectandInverseSolutionsofMeridianinTermsofDifferentLatitudeVariables[J].ActaGeodaeticaetCartographicaSinica,2016,45(5):560-565.DOI:10.11947/j.AGCS.2016.20140575.

    P226

    A

    1001-1595(2016)05-0560-06

    國家自然科學基金(41504031; 41571441);東華理工大學江西省數(shù)字國土重點實驗室開放研究基金(DLLJ201507)

    引文格式:過家春,李厚樸,莊云玲,等.依不同緯度變量的子午線弧長正反解公式的級數(shù)展開[J].測繪學報,2016,45(5):560-565.DOI:10.11947/j.AGCS.2016.20140575.

    日韩中文字幕欧美一区二区 | 久久精品aⅴ一区二区三区四区| 十八禁高潮呻吟视频| 老熟女久久久| 精品亚洲乱码少妇综合久久| 精品亚洲乱码少妇综合久久| 国产精品香港三级国产av潘金莲 | 国产成人一区二区在线| 在线天堂中文资源库| 肉色欧美久久久久久久蜜桃| 午夜福利免费观看在线| 国产一区二区三区综合在线观看| 欧美精品亚洲一区二区| 欧美国产精品一级二级三级| 午夜激情久久久久久久| 99精品久久久久人妻精品| 五月开心婷婷网| 免费在线观看完整版高清| 啦啦啦 在线观看视频| 老司机影院毛片| 欧美激情高清一区二区三区 | xxxhd国产人妻xxx| 久久精品国产综合久久久| 老司机影院成人| 一区二区三区激情视频| 国产深夜福利视频在线观看| 成年动漫av网址| 国产精品av久久久久免费| 观看av在线不卡| 黄片播放在线免费| 国产精品久久久久久久久免| 国产不卡av网站在线观看| 中文天堂在线官网| 最近最新中文字幕免费大全7| 日本av手机在线免费观看| 一级黄片播放器| 国产淫语在线视频| 国产精品 国内视频| 国产精品人妻久久久影院| netflix在线观看网站| 中文字幕亚洲精品专区| 看免费成人av毛片| 精品一区二区三区四区五区乱码 | 午夜91福利影院| 嫩草影视91久久| av有码第一页| 操出白浆在线播放| 亚洲一级一片aⅴ在线观看| 亚洲精品一区蜜桃| 国产片特级美女逼逼视频| 免费在线观看视频国产中文字幕亚洲 | 久久久欧美国产精品| 99国产精品免费福利视频| 老司机影院毛片| 又大又爽又粗| 极品人妻少妇av视频| 免费av中文字幕在线| 美女福利国产在线| 夜夜骑夜夜射夜夜干| 嫩草影视91久久| 日韩熟女老妇一区二区性免费视频| 哪个播放器可以免费观看大片| 人人妻人人添人人爽欧美一区卜| 老汉色∧v一级毛片| 另类亚洲欧美激情| 人体艺术视频欧美日本| 99久久精品国产亚洲精品| 亚洲精品av麻豆狂野| 如何舔出高潮| 亚洲精品久久成人aⅴ小说| 久久久久久久久久久免费av| 视频在线观看一区二区三区| 操美女的视频在线观看| 丝袜脚勾引网站| 一区福利在线观看| 亚洲精品久久久久久婷婷小说| www.熟女人妻精品国产| 老司机靠b影院| 捣出白浆h1v1| 久久久精品国产亚洲av高清涩受| 亚洲欧美清纯卡通| 欧美黑人精品巨大| 日韩av不卡免费在线播放| 国产亚洲午夜精品一区二区久久| 国产精品一区二区精品视频观看| 国产精品久久久人人做人人爽| 亚洲激情五月婷婷啪啪| 久久精品久久精品一区二区三区| 久久久久久人妻| videosex国产| 免费不卡黄色视频| 美女国产高潮福利片在线看| 国产午夜精品一二区理论片| 制服诱惑二区| 美女福利国产在线| 在线观看免费视频网站a站| av网站免费在线观看视频| 欧美中文综合在线视频| 亚洲精品中文字幕在线视频| av有码第一页| 久久久国产一区二区| 9热在线视频观看99| 久久精品aⅴ一区二区三区四区| 视频在线观看一区二区三区| 日本av手机在线免费观看| 三上悠亚av全集在线观看| 久久久久久久久久久免费av| 狂野欧美激情性xxxx| 久久久久久久久免费视频了| 免费久久久久久久精品成人欧美视频| 精品国产一区二区三区四区第35| 女人高潮潮喷娇喘18禁视频| 欧美激情 高清一区二区三区| 亚洲欧美一区二区三区黑人| 最黄视频免费看| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜爱| 国产精品99久久99久久久不卡 | 亚洲精品视频女| 一级毛片 在线播放| 国产日韩一区二区三区精品不卡| 亚洲欧美清纯卡通| 亚洲av成人不卡在线观看播放网 | 欧美日韩国产mv在线观看视频| 在线观看人妻少妇| 色精品久久人妻99蜜桃| 最新在线观看一区二区三区 | 飞空精品影院首页| 精品人妻一区二区三区麻豆| 久久天堂一区二区三区四区| 中文字幕人妻丝袜制服| 日韩av在线免费看完整版不卡| 一区二区av电影网| 亚洲情色 制服丝袜| 老司机深夜福利视频在线观看 | av网站在线播放免费| 亚洲色图综合在线观看| 少妇 在线观看| 美女午夜性视频免费| 波野结衣二区三区在线| 在线观看免费日韩欧美大片| 午夜激情久久久久久久| 亚洲精品国产av成人精品| 国产爽快片一区二区三区| 亚洲,欧美,日韩| 又大又爽又粗| 午夜福利影视在线免费观看| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 18在线观看网站| 日韩制服骚丝袜av| 久久久久国产一级毛片高清牌| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品电影小说| 精品亚洲乱码少妇综合久久| 岛国毛片在线播放| 中文欧美无线码| www.精华液| 免费日韩欧美在线观看| 999久久久国产精品视频| 欧美黑人欧美精品刺激| 亚洲一区二区三区欧美精品| 色综合欧美亚洲国产小说| 精品久久久精品久久久| 国产欧美日韩一区二区三区在线| 亚洲精华国产精华液的使用体验| 欧美变态另类bdsm刘玥| 丰满饥渴人妻一区二区三| 激情视频va一区二区三区| 新久久久久国产一级毛片| 18禁国产床啪视频网站| 精品国产一区二区三区久久久樱花| av免费观看日本| 91老司机精品| 最近的中文字幕免费完整| 丝袜喷水一区| 国产有黄有色有爽视频| 极品人妻少妇av视频| 在线观看免费高清a一片| 啦啦啦啦在线视频资源| 国产精品国产av在线观看| 欧美精品一区二区大全| 日韩大片免费观看网站| 大片免费播放器 马上看| 亚洲av在线观看美女高潮| 国产一区二区激情短视频 | 欧美精品人与动牲交sv欧美| 精品一区二区三区av网在线观看 | 欧美 日韩 精品 国产| 91精品伊人久久大香线蕉| 老汉色∧v一级毛片| 五月开心婷婷网| 亚洲美女黄色视频免费看| 天堂8中文在线网| 久久久精品区二区三区| 日韩欧美一区视频在线观看| 国产成人精品无人区| 伦理电影免费视频| 女的被弄到高潮叫床怎么办| 午夜福利免费观看在线| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 美女视频免费永久观看网站| 国产在线免费精品| 性少妇av在线| a级毛片黄视频| 免费看不卡的av| 亚洲精品,欧美精品| 成人国语在线视频| 亚洲一区中文字幕在线| kizo精华| 美女主播在线视频| 精品少妇内射三级| 国产色婷婷99| 日韩大码丰满熟妇| 国产亚洲午夜精品一区二区久久| 在线观看www视频免费| 伦理电影大哥的女人| 最新在线观看一区二区三区 | 国产精品成人在线| 国产精品欧美亚洲77777| 国产日韩欧美视频二区| 亚洲精品美女久久久久99蜜臀 | √禁漫天堂资源中文www| 国产精品一二三区在线看| 国产一区二区激情短视频 | 美女午夜性视频免费| 亚洲欧美精品自产自拍| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美国产在线视频| 在线观看一区二区三区激情| 欧美日韩亚洲综合一区二区三区_| 免费观看a级毛片全部| 久久青草综合色| 久久久久久人妻| 国产高清不卡午夜福利| 99久久精品国产亚洲精品| 高清欧美精品videossex| 亚洲国产毛片av蜜桃av| 男女免费视频国产| 2018国产大陆天天弄谢| 在线免费观看不下载黄p国产| 黄色一级大片看看| 免费不卡黄色视频| 亚洲自偷自拍图片 自拍| 99re6热这里在线精品视频| 大码成人一级视频| 18禁观看日本| 亚洲熟女精品中文字幕| 亚洲欧美精品综合一区二区三区| av.在线天堂| 午夜免费鲁丝| 宅男免费午夜| 国产精品 欧美亚洲| 高清在线视频一区二区三区| 国产 一区精品| 亚洲国产精品成人久久小说| 一级毛片黄色毛片免费观看视频| 国产成人av激情在线播放| 亚洲国产看品久久| 国产成人免费观看mmmm| 国产深夜福利视频在线观看| 一区在线观看完整版| 777久久人妻少妇嫩草av网站| 黄色视频不卡| 国产精品女同一区二区软件| 赤兔流量卡办理| 少妇的丰满在线观看| 狠狠精品人妻久久久久久综合| 日韩av在线免费看完整版不卡| 日韩免费高清中文字幕av| 午夜福利,免费看| 亚洲av福利一区| 国产精品蜜桃在线观看| 99精国产麻豆久久婷婷| 视频区图区小说| 日日摸夜夜添夜夜爱| 激情视频va一区二区三区| 看免费成人av毛片| 国产极品粉嫩免费观看在线| 丝袜美足系列| 久久国产亚洲av麻豆专区| 久久 成人 亚洲| 十八禁高潮呻吟视频| 成人亚洲欧美一区二区av| 狂野欧美激情性xxxx| 我要看黄色一级片免费的| 亚洲精品第二区| 免费在线观看完整版高清| 母亲3免费完整高清在线观看| 久久久久精品国产欧美久久久 | 啦啦啦 在线观看视频| 亚洲免费av在线视频| 女性被躁到高潮视频| 少妇被粗大的猛进出69影院| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 国产麻豆69| 9热在线视频观看99| 另类亚洲欧美激情| 两性夫妻黄色片| 国产99久久九九免费精品| 精品少妇一区二区三区视频日本电影 | 91aial.com中文字幕在线观看| 一级片'在线观看视频| 日韩一区二区三区影片| 19禁男女啪啪无遮挡网站| 日本vs欧美在线观看视频| 老司机在亚洲福利影院| 一本一本久久a久久精品综合妖精| 欧美日韩福利视频一区二区| 日韩欧美精品免费久久| 一区福利在线观看| 日本av免费视频播放| 午夜日韩欧美国产| 丝袜喷水一区| 只有这里有精品99| av电影中文网址| 热99久久久久精品小说推荐| 久久这里只有精品19| 久久久久精品性色| 免费女性裸体啪啪无遮挡网站| 美女福利国产在线| 国产一区二区激情短视频 | 在线观看三级黄色| av在线app专区| 日本av免费视频播放| 侵犯人妻中文字幕一二三四区| 亚洲欧美一区二区三区黑人| 高清不卡的av网站| 国产成人午夜福利电影在线观看| 中文字幕色久视频| 伊人久久国产一区二区| 国产一区二区三区av在线| 在线免费观看不下载黄p国产| 最近最新中文字幕大全免费视频 | 人人妻,人人澡人人爽秒播 | 女人精品久久久久毛片| 亚洲少妇的诱惑av| 日本一区二区免费在线视频| 亚洲国产欧美日韩在线播放| 纵有疾风起免费观看全集完整版| 2021少妇久久久久久久久久久| 亚洲中文av在线| 中国国产av一级| 欧美变态另类bdsm刘玥| www.精华液| 美女午夜性视频免费| 亚洲av男天堂| 国产乱人偷精品视频| 亚洲色图 男人天堂 中文字幕| 九色亚洲精品在线播放| 1024视频免费在线观看| 一级a爱视频在线免费观看| 叶爱在线成人免费视频播放| 久久 成人 亚洲| 高清视频免费观看一区二区| 成人国语在线视频| 亚洲精品中文字幕在线视频| 国产精品无大码| 一级黄片播放器| 欧美精品一区二区大全| 黄色视频在线播放观看不卡| 久久天躁狠狠躁夜夜2o2o | 久久精品aⅴ一区二区三区四区| 日韩av免费高清视频| 精品亚洲成a人片在线观看| 免费在线观看黄色视频的| 久久久久久免费高清国产稀缺| 伊人亚洲综合成人网| 国产精品熟女久久久久浪| 免费观看人在逋| 国产在视频线精品| 免费观看人在逋| 亚洲综合精品二区| 亚洲成人av在线免费| 午夜日韩欧美国产| 精品国产露脸久久av麻豆| 亚洲国产欧美网| 人人妻人人澡人人爽人人夜夜| 啦啦啦中文免费视频观看日本| 捣出白浆h1v1| 久久韩国三级中文字幕| 色94色欧美一区二区| 男男h啪啪无遮挡| 性高湖久久久久久久久免费观看| 久久青草综合色| 老汉色av国产亚洲站长工具| av国产精品久久久久影院| 欧美黑人欧美精品刺激| 欧美精品av麻豆av| 熟女av电影| 成人影院久久| 免费在线观看完整版高清| 亚洲av电影在线观看一区二区三区| 久久久欧美国产精品| 国产视频首页在线观看| 久久久久精品人妻al黑| 99久久99久久久精品蜜桃| 亚洲av电影在线进入| 王馨瑶露胸无遮挡在线观看| 制服丝袜香蕉在线| 97精品久久久久久久久久精品| 日韩av在线免费看完整版不卡| 一边摸一边做爽爽视频免费| 日本午夜av视频| 天天躁夜夜躁狠狠躁躁| e午夜精品久久久久久久| 欧美精品高潮呻吟av久久| 久久久国产一区二区| av视频免费观看在线观看| 亚洲av电影在线观看一区二区三区| 亚洲成av片中文字幕在线观看| 一区二区三区四区激情视频| 国产av精品麻豆| 又黄又粗又硬又大视频| 精品福利永久在线观看| 亚洲熟女毛片儿| 91aial.com中文字幕在线观看| 亚洲国产成人一精品久久久| 最新在线观看一区二区三区 | 国产成人一区二区在线| 卡戴珊不雅视频在线播放| 人体艺术视频欧美日本| 丝袜脚勾引网站| 欧美亚洲 丝袜 人妻 在线| 老汉色∧v一级毛片| 少妇人妻 视频| 午夜福利视频在线观看免费| 一本—道久久a久久精品蜜桃钙片| 国产精品麻豆人妻色哟哟久久| 美女扒开内裤让男人捅视频| 欧美黄色片欧美黄色片| 国产1区2区3区精品| 色网站视频免费| 日本欧美国产在线视频| 亚洲国产精品一区二区三区在线| 国产乱来视频区| 69精品国产乱码久久久| 亚洲少妇的诱惑av| 一级毛片 在线播放| 乱人伦中国视频| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲av高清不卡| 久久99精品国语久久久| 91老司机精品| av视频免费观看在线观看| 汤姆久久久久久久影院中文字幕| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产欧美一区二区综合| 免费在线观看视频国产中文字幕亚洲 | 国产人伦9x9x在线观看| 一级爰片在线观看| 久久久精品国产亚洲av高清涩受| 国产精品一二三区在线看| 成人毛片60女人毛片免费| 99九九在线精品视频| 在线天堂中文资源库| 日韩欧美精品免费久久| 日韩av不卡免费在线播放| 欧美亚洲日本最大视频资源| 七月丁香在线播放| 国产午夜精品一二区理论片| 一二三四中文在线观看免费高清| 制服丝袜香蕉在线| 亚洲精品成人av观看孕妇| 中文字幕人妻熟女乱码| 在线观看人妻少妇| 欧美精品av麻豆av| 女性生殖器流出的白浆| 久久久久网色| 亚洲美女黄色视频免费看| 亚洲av在线观看美女高潮| 欧美日韩亚洲国产一区二区在线观看 | 午夜影院在线不卡| 国产日韩一区二区三区精品不卡| 亚洲成色77777| 亚洲精品日本国产第一区| 国产精品偷伦视频观看了| 九九爱精品视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久亚洲国产成人精品v| 91精品伊人久久大香线蕉| 最近最新中文字幕大全免费视频 | 久久精品国产亚洲av涩爱| 午夜91福利影院| 中文字幕另类日韩欧美亚洲嫩草| 最近中文字幕高清免费大全6| 欧美黄色片欧美黄色片| 久久久久人妻精品一区果冻| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| 国产无遮挡羞羞视频在线观看| 超色免费av| 少妇 在线观看| 美女大奶头黄色视频| 国产 精品1| 午夜福利网站1000一区二区三区| 王馨瑶露胸无遮挡在线观看| 波野结衣二区三区在线| 午夜日韩欧美国产| 亚洲成色77777| 午夜日韩欧美国产| 亚洲av电影在线进入| 精品国产国语对白av| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜爱| 性高湖久久久久久久久免费观看| 久久精品久久精品一区二区三区| 中文字幕人妻丝袜一区二区 | 91精品三级在线观看| 日本av手机在线免费观看| 国产精品二区激情视频| 9色porny在线观看| 国产熟女欧美一区二区| 一边摸一边做爽爽视频免费| 色婷婷av一区二区三区视频| 亚洲视频免费观看视频| 精品国产国语对白av| 电影成人av| 亚洲av日韩精品久久久久久密 | 国产精品欧美亚洲77777| 亚洲一区中文字幕在线| 在线精品无人区一区二区三| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品成人久久小说| 亚洲男人天堂网一区| 观看av在线不卡| 2021少妇久久久久久久久久久| 午夜福利一区二区在线看| 欧美亚洲 丝袜 人妻 在线| 女性被躁到高潮视频| 色94色欧美一区二区| 两个人看的免费小视频| 日韩一区二区视频免费看| 爱豆传媒免费全集在线观看| 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 老汉色∧v一级毛片| 国产精品无大码| 热re99久久精品国产66热6| 亚洲精品,欧美精品| 亚洲熟女毛片儿| 欧美日韩成人在线一区二区| 亚洲五月色婷婷综合| 一级毛片黄色毛片免费观看视频| 咕卡用的链子| 亚洲av电影在线进入| 男男h啪啪无遮挡| 女性生殖器流出的白浆| av在线播放精品| 亚洲成av片中文字幕在线观看| 成人漫画全彩无遮挡| 国产午夜精品一二区理论片| 午夜久久久在线观看| 成年动漫av网址| 亚洲欧洲日产国产| 精品视频人人做人人爽| 各种免费的搞黄视频| 老司机靠b影院| 伦理电影大哥的女人| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频| 欧美人与善性xxx| 国产精品亚洲av一区麻豆 | 久久久久久久久久久久大奶| 亚洲欧美精品综合一区二区三区| 亚洲精品久久久久久婷婷小说| 国产无遮挡羞羞视频在线观看| av福利片在线| 乱人伦中国视频| 桃花免费在线播放| 97在线人人人人妻| 成年女人毛片免费观看观看9 | 看免费av毛片| 美女高潮到喷水免费观看| av国产久精品久网站免费入址| 肉色欧美久久久久久久蜜桃| 老司机在亚洲福利影院| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 香蕉国产在线看| 国产精品久久久久久久久免| 少妇被粗大猛烈的视频| 欧美成人精品欧美一级黄| 最黄视频免费看| 新久久久久国产一级毛片| 可以免费在线观看a视频的电影网站 | 国产成人午夜福利电影在线观看| √禁漫天堂资源中文www| 亚洲第一青青草原| xxxhd国产人妻xxx| 亚洲av日韩精品久久久久久密 | 国产精品成人在线| 天天躁日日躁夜夜躁夜夜| 国产在线视频一区二区| 青青草视频在线视频观看| 久久毛片免费看一区二区三区| 亚洲国产中文字幕在线视频| 成人午夜精彩视频在线观看| 欧美另类一区| 久久毛片免费看一区二区三区| 亚洲国产精品一区三区| 精品人妻熟女毛片av久久网站| 两个人看的免费小视频| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 老汉色av国产亚洲站长工具| 免费观看性生交大片5| 久久鲁丝午夜福利片| 精品卡一卡二卡四卡免费| 你懂的网址亚洲精品在线观看| 老司机影院成人| 亚洲一区二区三区欧美精品| 一区二区三区激情视频|