• <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.

    免费少妇av软件| 一个人免费看片子| 久久人人爽人人爽人人片va| tube8黄色片| 欧美一级a爱片免费观看看| 国产色爽女视频免费观看| 尾随美女入室| 精品一品国产午夜福利视频| 中国美白少妇内射xxxbb| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久久久免| 亚洲精品国产av蜜桃| 免费在线观看成人毛片| 国产 精品1| 国产成人aa在线观看| 2021少妇久久久久久久久久久| 能在线免费看毛片的网站| 午夜91福利影院| 嫩草影院新地址| 最近中文字幕高清免费大全6| 久久综合国产亚洲精品| 中国国产av一级| a级毛片免费高清观看在线播放| 午夜久久久在线观看| 国产亚洲欧美精品永久| 亚洲精品色激情综合| 在线观看三级黄色| 国产男人的电影天堂91| 又粗又硬又长又爽又黄的视频| 丰满人妻一区二区三区视频av| 99热网站在线观看| 亚洲国产av新网站| 亚洲伊人久久精品综合| 一区在线观看完整版| 少妇丰满av| 久久精品久久久久久噜噜老黄| 欧美日韩视频精品一区| 国产一区二区三区av在线| 日韩av不卡免费在线播放| 国产欧美亚洲国产| 亚洲经典国产精华液单| 九草在线视频观看| 亚洲中文av在线| 日韩,欧美,国产一区二区三区| 国产精品一区二区三区四区免费观看| 十八禁网站网址无遮挡 | 国国产精品蜜臀av免费| 九色成人免费人妻av| 日韩三级伦理在线观看| 最新中文字幕久久久久| 纵有疾风起免费观看全集完整版| 寂寞人妻少妇视频99o| 亚洲av欧美aⅴ国产| 91午夜精品亚洲一区二区三区| av播播在线观看一区| 3wmmmm亚洲av在线观看| 女人精品久久久久毛片| 国产精品免费大片| 我要看日韩黄色一级片| 天堂8中文在线网| 久久精品久久精品一区二区三区| 亚洲婷婷狠狠爱综合网| 国产免费一区二区三区四区乱码| 久久国产亚洲av麻豆专区| 中文字幕av电影在线播放| 日本爱情动作片www.在线观看| 国产精品女同一区二区软件| 美女xxoo啪啪120秒动态图| av国产精品久久久久影院| 日日摸夜夜添夜夜爱| 国产高清有码在线观看视频| 日韩强制内射视频| 人人澡人人妻人| 久久影院123| 成人18禁高潮啪啪吃奶动态图 | 这个男人来自地球电影免费观看 | 久久精品久久久久久久性| 曰老女人黄片| 免费黄色在线免费观看| 99久国产av精品国产电影| 色婷婷久久久亚洲欧美| 老司机亚洲免费影院| 欧美日韩av久久| 人妻系列 视频| 91久久精品国产一区二区成人| 精品少妇久久久久久888优播| 一区二区三区乱码不卡18| 91精品国产国语对白视频| 91精品国产国语对白视频| 亚洲综合精品二区| 在线 av 中文字幕| 国产淫片久久久久久久久| 国产亚洲一区二区精品| 老熟女久久久| 91成人精品电影| 两个人免费观看高清视频 | 交换朋友夫妻互换小说| 女人久久www免费人成看片| 80岁老熟妇乱子伦牲交| 欧美高清成人免费视频www| 在线播放无遮挡| 欧美变态另类bdsm刘玥| 亚洲在久久综合| 精品国产一区二区三区久久久樱花| 精品熟女少妇av免费看| 天天躁夜夜躁狠狠久久av| 日韩一本色道免费dvd| 观看av在线不卡| 哪个播放器可以免费观看大片| 亚洲伊人久久精品综合| 色哟哟·www| 狂野欧美激情性xxxx在线观看| 欧美日韩视频高清一区二区三区二| 好男人视频免费观看在线| 日本欧美国产在线视频| 天堂俺去俺来也www色官网| 亚洲成人一二三区av| 亚洲国产精品国产精品| 青春草视频在线免费观看| 99久久精品一区二区三区| 欧美激情国产日韩精品一区| 国产精品.久久久| 亚洲电影在线观看av| 亚洲激情五月婷婷啪啪| 成人黄色视频免费在线看| 亚洲av成人精品一二三区| 国产亚洲欧美精品永久| 综合色丁香网| 精品久久久精品久久久| 亚洲精品第二区| 久久精品国产鲁丝片午夜精品| 80岁老熟妇乱子伦牲交| 黄色一级大片看看| 精品熟女少妇av免费看| 欧美日韩视频精品一区| 中国三级夫妇交换| 国产男女内射视频| 久久青草综合色| 免费观看的影片在线观看| 人妻系列 视频| 十八禁高潮呻吟视频 | 国产成人91sexporn| 九色成人免费人妻av| 免费观看无遮挡的男女| 亚洲成人手机| 国产淫片久久久久久久久| 国产精品蜜桃在线观看| 国产成人精品一,二区| 欧美激情国产日韩精品一区| 永久网站在线| 国产真实伦视频高清在线观看| 亚洲综合色惰| 久久久精品94久久精品| 卡戴珊不雅视频在线播放| 老司机亚洲免费影院| 久久99热6这里只有精品| 国产亚洲av片在线观看秒播厂| 亚洲久久久国产精品| 欧美精品人与动牲交sv欧美| 日韩精品有码人妻一区| 国产男女内射视频| 夜夜骑夜夜射夜夜干| 水蜜桃什么品种好| 免费看av在线观看网站| 亚洲三级黄色毛片| 亚洲欧美日韩东京热| 欧美精品一区二区大全| 国产极品粉嫩免费观看在线 | 国产成人精品久久久久久| 日本与韩国留学比较| 777米奇影视久久| 久久久久久久精品精品| 日韩亚洲欧美综合| 国产精品久久久久久精品古装| 久久国产精品男人的天堂亚洲 | 久久青草综合色| av女优亚洲男人天堂| 精品99又大又爽又粗少妇毛片| 成人午夜精彩视频在线观看| 少妇高潮的动态图| 久久热精品热| 99热网站在线观看| 国产精品一区二区在线不卡| 最近手机中文字幕大全| 欧美日韩亚洲高清精品| 日韩电影二区| 少妇高潮的动态图| 国产极品天堂在线| 777米奇影视久久| 国产高清有码在线观看视频| 高清欧美精品videossex| 亚洲欧洲精品一区二区精品久久久 | 三级经典国产精品| 久久久久久久久久久久大奶| 观看av在线不卡| √禁漫天堂资源中文www| 亚洲,一卡二卡三卡| 日韩 亚洲 欧美在线| 久久精品国产鲁丝片午夜精品| 久久精品国产亚洲av涩爱| 日韩av不卡免费在线播放| 久久精品久久久久久噜噜老黄| 午夜激情久久久久久久| 日日爽夜夜爽网站| 国产精品一二三区在线看| 国产高清三级在线| 五月开心婷婷网| 日韩av在线免费看完整版不卡| 日韩精品有码人妻一区| 91久久精品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 人妻系列 视频| 婷婷色av中文字幕| 欧美精品亚洲一区二区| 一级a做视频免费观看| 精品一品国产午夜福利视频| 极品人妻少妇av视频| 51国产日韩欧美| 久久精品久久精品一区二区三区| 精品国产国语对白av| 亚洲精品456在线播放app| 国产在线男女| 大话2 男鬼变身卡| 我的女老师完整版在线观看| 天堂8中文在线网| 久久久久网色| 涩涩av久久男人的天堂| 久久精品熟女亚洲av麻豆精品| 欧美区成人在线视频| 亚洲国产av新网站| 中国国产av一级| 亚洲精品日本国产第一区| 自线自在国产av| 麻豆成人av视频| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| 亚洲,一卡二卡三卡| 另类精品久久| 尾随美女入室| 亚洲av中文av极速乱| 女人精品久久久久毛片| 国产成人免费无遮挡视频| 亚洲av综合色区一区| 夜夜爽夜夜爽视频| av福利片在线观看| 亚洲天堂av无毛| 一级毛片久久久久久久久女| 少妇 在线观看| 少妇被粗大的猛进出69影院 | 精品久久久久久久久av| 国产成人a∨麻豆精品| 香蕉精品网在线| 午夜视频国产福利| 国产成人精品无人区| 欧美精品国产亚洲| 卡戴珊不雅视频在线播放| 国产精品蜜桃在线观看| 伊人久久精品亚洲午夜| 免费高清在线观看视频在线观看| 午夜久久久在线观看| 亚洲欧美清纯卡通| 黑人巨大精品欧美一区二区蜜桃 | 欧美另类一区| 插阴视频在线观看视频| 新久久久久国产一级毛片| 99re6热这里在线精品视频| 婷婷色麻豆天堂久久| 九色成人免费人妻av| 最近手机中文字幕大全| 免费黄网站久久成人精品| 伦精品一区二区三区| 国产黄色免费在线视频| 精品久久国产蜜桃| 下体分泌物呈黄色| tube8黄色片| 精品国产露脸久久av麻豆| 99热网站在线观看| 精品少妇内射三级| 国产爽快片一区二区三区| 亚洲国产日韩一区二区| 国产精品成人在线| 色视频在线一区二区三区| 国产无遮挡羞羞视频在线观看| 在线播放无遮挡| 欧美精品一区二区免费开放| 在线观看国产h片| 国产视频内射| 一级毛片aaaaaa免费看小| 国产成人精品久久久久久| 日韩av免费高清视频| 国产淫片久久久久久久久| 日韩一区二区视频免费看| 日韩一本色道免费dvd| 一级二级三级毛片免费看| 91精品国产九色| 日本黄色日本黄色录像| 看非洲黑人一级黄片| 久久99热这里只频精品6学生| 国产精品蜜桃在线观看| www.av在线官网国产| 国内揄拍国产精品人妻在线| 卡戴珊不雅视频在线播放| 国产深夜福利视频在线观看| 人体艺术视频欧美日本| 亚洲精品乱码久久久久久按摩| 大码成人一级视频| 曰老女人黄片| h视频一区二区三区| 亚洲一区二区三区欧美精品| 在线观看美女被高潮喷水网站| 成人二区视频| 久久av网站| 免费人妻精品一区二区三区视频| 99九九线精品视频在线观看视频| 另类精品久久| 伦理电影大哥的女人| 99久久精品热视频| 插阴视频在线观看视频| 91成人精品电影| 91午夜精品亚洲一区二区三区| 97超视频在线观看视频| 2021少妇久久久久久久久久久| 亚洲丝袜综合中文字幕| 免费高清在线观看视频在线观看| 18禁在线播放成人免费| 免费看日本二区| 精品酒店卫生间| 一本—道久久a久久精品蜜桃钙片| 高清欧美精品videossex| 日韩电影二区| 免费人妻精品一区二区三区视频| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 一级毛片黄色毛片免费观看视频| 国产亚洲91精品色在线| 99热6这里只有精品| 日日摸夜夜添夜夜添av毛片| av福利片在线| 成人无遮挡网站| 国产一区有黄有色的免费视频| 久久久亚洲精品成人影院| 嫩草影院新地址| 亚洲精品国产av蜜桃| 不卡视频在线观看欧美| 一区二区av电影网| 亚洲国产精品成人久久小说| a级一级毛片免费在线观看| 久久99蜜桃精品久久| 一二三四中文在线观看免费高清| 亚洲不卡免费看| 我的老师免费观看完整版| 亚洲av电影在线观看一区二区三区| 日本黄色片子视频| 国产精品.久久久| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 国产白丝娇喘喷水9色精品| 日韩中文字幕视频在线看片| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 国产成人免费观看mmmm| 狂野欧美激情性bbbbbb| 麻豆成人av视频| 日韩大片免费观看网站| 国产中年淑女户外野战色| 性高湖久久久久久久久免费观看| 国产亚洲欧美精品永久| 如何舔出高潮| 亚洲国产精品一区二区三区在线| 亚洲av电影在线观看一区二区三区| 久久精品国产a三级三级三级| 在线观看av片永久免费下载| 亚洲高清免费不卡视频| 久久这里有精品视频免费| 免费观看a级毛片全部| 免费观看的影片在线观看| 另类亚洲欧美激情| 人人妻人人澡人人看| 免费在线观看成人毛片| 人体艺术视频欧美日本| 男女国产视频网站| 亚洲精品一区蜜桃| 免费人妻精品一区二区三区视频| 欧美日韩国产mv在线观看视频| 少妇的逼水好多| 亚洲国产av新网站| 麻豆乱淫一区二区| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久 | 麻豆乱淫一区二区| 色94色欧美一区二区| 伦理电影免费视频| 涩涩av久久男人的天堂| 自拍偷自拍亚洲精品老妇| 少妇裸体淫交视频免费看高清| 欧美97在线视频| 视频区图区小说| 天堂俺去俺来也www色官网| 国产熟女欧美一区二区| 黑人猛操日本美女一级片| 精品亚洲乱码少妇综合久久| 一区二区三区精品91| av不卡在线播放| 欧美日本中文国产一区发布| 国产男人的电影天堂91| 晚上一个人看的免费电影| 成人毛片a级毛片在线播放| 免费观看无遮挡的男女| 午夜免费观看性视频| 欧美bdsm另类| 国产淫片久久久久久久久| 亚洲久久久国产精品| 丰满人妻一区二区三区视频av| 狂野欧美激情性bbbbbb| 国产高清不卡午夜福利| 久久久久人妻精品一区果冻| 亚洲av在线观看美女高潮| 日本黄大片高清| 校园人妻丝袜中文字幕| 高清在线视频一区二区三区| 男人舔奶头视频| 午夜福利视频精品| 午夜免费观看性视频| 欧美bdsm另类| 51国产日韩欧美| 亚洲欧美日韩东京热| 久久国产乱子免费精品| 一区二区三区免费毛片| 国产色爽女视频免费观看| 欧美三级亚洲精品| 嫩草影院新地址| 午夜视频国产福利| 一本色道久久久久久精品综合| 妹子高潮喷水视频| 久久99一区二区三区| 久久国产亚洲av麻豆专区| 少妇被粗大的猛进出69影院 | 99九九在线精品视频 | 亚洲精品国产成人久久av| 久久精品国产亚洲av天美| 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 久久av网站| 高清av免费在线| 久久精品国产a三级三级三级| 最新的欧美精品一区二区| 亚洲四区av| 久久鲁丝午夜福利片| 久久ye,这里只有精品| 国产日韩欧美在线精品| 午夜久久久在线观看| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 美女福利国产在线| 国产伦在线观看视频一区| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 大香蕉久久网| 亚洲自偷自拍三级| 欧美日韩亚洲高清精品| 国产淫片久久久久久久久| 日韩欧美精品免费久久| 国产探花极品一区二区| 国产精品三级大全| 亚洲av不卡在线观看| 免费人妻精品一区二区三区视频| 亚洲国产最新在线播放| 国产成人freesex在线| 成人18禁高潮啪啪吃奶动态图 | 日本欧美国产在线视频| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看| 午夜视频国产福利| 精品午夜福利在线看| 自线自在国产av| 一区在线观看完整版| 建设人人有责人人尽责人人享有的| 日韩av不卡免费在线播放| 亚洲四区av| 国产精品一区二区在线观看99| 色视频在线一区二区三区| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 特大巨黑吊av在线直播| 日本与韩国留学比较| 51国产日韩欧美| 国产精品一区二区性色av| 国模一区二区三区四区视频| 国产成人精品福利久久| 视频中文字幕在线观看| 18+在线观看网站| √禁漫天堂资源中文www| 久久国产精品男人的天堂亚洲 | 在线观看一区二区三区激情| 少妇人妻久久综合中文| 中文字幕亚洲精品专区| 黄片无遮挡物在线观看| 青春草国产在线视频| 男女边摸边吃奶| av女优亚洲男人天堂| 久久99一区二区三区| 春色校园在线视频观看| 久久久国产精品麻豆| 嫩草影院入口| 十八禁网站网址无遮挡 | 两个人免费观看高清视频 | 中文字幕久久专区| 97超视频在线观看视频| 新久久久久国产一级毛片| 亚洲欧洲国产日韩| 在线亚洲精品国产二区图片欧美 | av福利片在线| 免费看av在线观看网站| 街头女战士在线观看网站| 国产精品国产三级国产专区5o| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 永久免费av网站大全| 欧美精品亚洲一区二区| 2022亚洲国产成人精品| 国产69精品久久久久777片| 国产极品天堂在线| 嘟嘟电影网在线观看| 又大又黄又爽视频免费| 久久久久久久国产电影| 国产欧美日韩综合在线一区二区 | www.av在线官网国产| 国产成人a∨麻豆精品| 性高湖久久久久久久久免费观看| 日本av手机在线免费观看| 久久av网站| 色视频在线一区二区三区| 十八禁网站网址无遮挡 | 亚洲,欧美,日韩| 一级毛片久久久久久久久女| 色5月婷婷丁香| 午夜免费鲁丝| 国产在线视频一区二区| 天堂8中文在线网| 精品人妻熟女毛片av久久网站| 校园人妻丝袜中文字幕| 一个人看视频在线观看www免费| 天堂8中文在线网| 欧美日本中文国产一区发布| 色视频在线一区二区三区| 亚洲久久久国产精品| 偷拍熟女少妇极品色| 久久毛片免费看一区二区三区| 噜噜噜噜噜久久久久久91| 亚洲国产精品999| 午夜免费鲁丝| 大香蕉久久网| av卡一久久| 免费av不卡在线播放| 性色avwww在线观看| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 免费观看av网站的网址| 成年人午夜在线观看视频| 国产免费一级a男人的天堂| a 毛片基地| 2022亚洲国产成人精品| 国产69精品久久久久777片| 如日韩欧美国产精品一区二区三区 | 你懂的网址亚洲精品在线观看| 一本久久精品| 夜夜看夜夜爽夜夜摸| 最后的刺客免费高清国语| 夫妻性生交免费视频一级片| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 国产精品久久久久久精品电影小说| 性高湖久久久久久久久免费观看| 成人特级av手机在线观看| 国产又色又爽无遮挡免| 久久久国产一区二区| 亚洲精品乱码久久久v下载方式| 自拍偷自拍亚洲精品老妇| 色视频www国产| 免费观看a级毛片全部| 免费黄网站久久成人精品| 七月丁香在线播放| 成人综合一区亚洲| 黑人猛操日本美女一级片| 中文资源天堂在线| 国产高清有码在线观看视频| 国产黄片美女视频| 国产精品久久久久成人av| 22中文网久久字幕| 少妇 在线观看| 日产精品乱码卡一卡2卡三| 国产精品蜜桃在线观看| 汤姆久久久久久久影院中文字幕| av免费在线看不卡| 日韩精品免费视频一区二区三区 | 人妻人人澡人人爽人人| 久久久亚洲精品成人影院| 中文字幕免费在线视频6| 一级片'在线观看视频| 日韩一区二区视频免费看| 亚洲精品久久午夜乱码| 精品国产露脸久久av麻豆| 99久国产av精品国产电影| 男女边摸边吃奶| 在线免费观看不下载黄p国产| 日韩中字成人| 亚洲av成人精品一区久久| 成人18禁高潮啪啪吃奶动态图 | 久久亚洲国产成人精品v| 久久韩国三级中文字幕| 午夜福利网站1000一区二区三区| 欧美丝袜亚洲另类| 国产成人免费观看mmmm|