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

    亚洲在久久综合| 国产精品久久久久久精品古装| 久久久久久久大尺度免费视频| 国产午夜精品一二区理论片| 国产成人精品久久久久久| 日韩伦理黄色片| 久久久精品欧美日韩精品| 一区二区三区乱码不卡18| 特大巨黑吊av在线直播| 亚洲国产日韩一区二区| 五月伊人婷婷丁香| 69人妻影院| 国产视频内射| 免费av观看视频| 国产午夜福利久久久久久| 五月伊人婷婷丁香| 国产成人a∨麻豆精品| 99热全是精品| 人人妻人人爽人人添夜夜欢视频 | 一级黄片播放器| 成人一区二区视频在线观看| 制服丝袜香蕉在线| 一边亲一边摸免费视频| 久久久久网色| 成人一区二区视频在线观看| 国产黄色免费在线视频| 久久ye,这里只有精品| 色婷婷久久久亚洲欧美| 美女被艹到高潮喷水动态| 大话2 男鬼变身卡| 丰满乱子伦码专区| 日韩 亚洲 欧美在线| 久久精品熟女亚洲av麻豆精品| 亚洲第一区二区三区不卡| 日日撸夜夜添| 久久久色成人| 欧美高清成人免费视频www| 久久久亚洲精品成人影院| 日韩av免费高清视频| 国内精品宾馆在线| av专区在线播放| 纵有疾风起免费观看全集完整版| 97在线视频观看| 精品一区二区三区视频在线| 久热这里只有精品99| 午夜福利网站1000一区二区三区| 日本猛色少妇xxxxx猛交久久| 日韩亚洲欧美综合| 成人午夜精彩视频在线观看| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 国产 精品1| a级一级毛片免费在线观看| 亚洲精品aⅴ在线观看| 亚洲av免费在线观看| 国产亚洲精品久久久com| 老女人水多毛片| 精品国产一区二区三区久久久樱花 | 在线精品无人区一区二区三 | 综合色丁香网| 成年av动漫网址| 亚洲经典国产精华液单| 亚洲自偷自拍三级| 精品一区二区三区视频在线| 精品久久久久久久人妻蜜臀av| 亚洲精品亚洲一区二区| 校园人妻丝袜中文字幕| 极品少妇高潮喷水抽搐| 欧美日韩一区二区视频在线观看视频在线 | 91精品伊人久久大香线蕉| 精品99又大又爽又粗少妇毛片| 黄色一级大片看看| 亚洲国产av新网站| 久久久久网色| 国产伦精品一区二区三区四那| 国产永久视频网站| 少妇裸体淫交视频免费看高清| 一级av片app| 九草在线视频观看| 一本色道久久久久久精品综合| 美女cb高潮喷水在线观看| 国产日韩欧美在线精品| 亚洲av中文av极速乱| 狂野欧美激情性xxxx在线观看| 麻豆精品久久久久久蜜桃| 国产成人一区二区在线| 婷婷色综合大香蕉| 婷婷色综合www| 成人美女网站在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 少妇猛男粗大的猛烈进出视频 | 久久久久久国产a免费观看| 最近的中文字幕免费完整| 国产精品国产av在线观看| 婷婷色麻豆天堂久久| 听说在线观看完整版免费高清| 三级国产精品片| 精品少妇黑人巨大在线播放| 国产成人91sexporn| 黄色怎么调成土黄色| 热99国产精品久久久久久7| 国产精品99久久99久久久不卡 | 97精品久久久久久久久久精品| 香蕉精品网在线| 免费人成在线观看视频色| 日韩国内少妇激情av| 最近中文字幕2019免费版| 婷婷色麻豆天堂久久| 美女国产视频在线观看| 亚洲精品国产色婷婷电影| 18禁在线无遮挡免费观看视频| 少妇 在线观看| 欧美老熟妇乱子伦牲交| 精品久久久久久久久亚洲| 午夜福利视频精品| 欧美潮喷喷水| 3wmmmm亚洲av在线观看| 在现免费观看毛片| 人妻 亚洲 视频| 亚洲精品久久久久久婷婷小说| 少妇裸体淫交视频免费看高清| 熟妇人妻不卡中文字幕| 精品人妻偷拍中文字幕| 99九九线精品视频在线观看视频| 亚洲av在线观看美女高潮| 精品熟女少妇av免费看| 国产精品女同一区二区软件| av黄色大香蕉| 欧美亚洲 丝袜 人妻 在线| 69人妻影院| 极品少妇高潮喷水抽搐| 国内少妇人妻偷人精品xxx网站| 综合色av麻豆| 97在线人人人人妻| 国产亚洲5aaaaa淫片| 欧美激情在线99| 免费黄频网站在线观看国产| 自拍偷自拍亚洲精品老妇| 啦啦啦中文免费视频观看日本| 别揉我奶头 嗯啊视频| 亚洲欧美精品自产自拍| 黄色怎么调成土黄色| 高清在线视频一区二区三区| 国产高潮美女av| 在线观看美女被高潮喷水网站| 亚洲av不卡在线观看| 99精国产麻豆久久婷婷| 听说在线观看完整版免费高清| 国产91av在线免费观看| 边亲边吃奶的免费视频| 女人十人毛片免费观看3o分钟| videos熟女内射| 国产高潮美女av| 在线观看美女被高潮喷水网站| 久久精品国产a三级三级三级| 亚洲电影在线观看av| 欧美国产精品一级二级三级 | 国产黄色视频一区二区在线观看| 白带黄色成豆腐渣| 欧美国产精品一级二级三级 | 天堂俺去俺来也www色官网| 欧美激情久久久久久爽电影| 国产男人的电影天堂91| 国产爽快片一区二区三区| 欧美三级亚洲精品| 91久久精品国产一区二区三区| 18禁动态无遮挡网站| 国产有黄有色有爽视频| 国产黄片视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 男的添女的下面高潮视频| 综合色丁香网| 秋霞伦理黄片| 99热网站在线观看| 精品人妻熟女av久视频| 大陆偷拍与自拍| 激情五月婷婷亚洲| 在线观看人妻少妇| 精品人妻视频免费看| 久久久久九九精品影院| 久久鲁丝午夜福利片| 最近最新中文字幕免费大全7| 熟女人妻精品中文字幕| 黄色配什么色好看| 亚洲自偷自拍三级| 久久精品国产亚洲av天美| 欧美精品人与动牲交sv欧美| 午夜福利视频精品| 色5月婷婷丁香| 日日撸夜夜添| 晚上一个人看的免费电影| 国产午夜福利久久久久久| 听说在线观看完整版免费高清| av线在线观看网站| 国产一区二区三区av在线| 国产在视频线精品| 精品少妇黑人巨大在线播放| 97在线视频观看| 99热这里只有精品一区| 汤姆久久久久久久影院中文字幕| 亚洲最大成人手机在线| 国产毛片a区久久久久| 99热网站在线观看| 久久久久久久大尺度免费视频| 国内精品美女久久久久久| 综合色丁香网| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 亚洲成人一二三区av| 国产午夜福利久久久久久| 在线观看美女被高潮喷水网站| 免费黄频网站在线观看国产| 99热国产这里只有精品6| 国产淫片久久久久久久久| 内射极品少妇av片p| 各种免费的搞黄视频| 干丝袜人妻中文字幕| 熟女电影av网| 69av精品久久久久久| 联通29元200g的流量卡| 美女国产视频在线观看| 久久国内精品自在自线图片| 老司机影院成人| 高清日韩中文字幕在线| 69av精品久久久久久| 亚洲国产成人一精品久久久| 伦精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 最近2019中文字幕mv第一页| 午夜视频国产福利| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美在线一区| videossex国产| 中文在线观看免费www的网站| 少妇高潮的动态图| 汤姆久久久久久久影院中文字幕| 久久久久久久午夜电影| 在线 av 中文字幕| 免费人成在线观看视频色| 国内精品美女久久久久久| 国产av不卡久久| 午夜爱爱视频在线播放| 天堂网av新在线| 精品国产一区二区三区久久久樱花 | 国产爽快片一区二区三区| 亚洲精品一二三| 国产精品一二三区在线看| 日日啪夜夜爽| 国内精品美女久久久久久| 亚洲欧洲日产国产| 国产精品99久久99久久久不卡 | 久久精品夜色国产| 国产伦精品一区二区三区四那| 中文字幕av成人在线电影| av一本久久久久| 国产乱来视频区| 中国美白少妇内射xxxbb| 亚洲真实伦在线观看| 深爱激情五月婷婷| 亚洲国产成人一精品久久久| 美女脱内裤让男人舔精品视频| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻少妇偷人精品九色| 天堂中文最新版在线下载 | 亚洲精华国产精华液的使用体验| 色吧在线观看| 精品人妻一区二区三区麻豆| 精品国产三级普通话版| 欧美日韩亚洲高清精品| 亚洲精品日韩av片在线观看| 国产69精品久久久久777片| 日本-黄色视频高清免费观看| 日韩成人av中文字幕在线观看| 国产欧美日韩精品一区二区| 又粗又硬又长又爽又黄的视频| 日本三级黄在线观看| 91aial.com中文字幕在线观看| 在线免费十八禁| 精品国产露脸久久av麻豆| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 国产黄a三级三级三级人| videos熟女内射| 男女下面进入的视频免费午夜| 三级国产精品片| 国产高潮美女av| 国产在视频线精品| 国产av码专区亚洲av| 久久精品久久精品一区二区三区| 在线精品无人区一区二区三 | 三级国产精品片| 交换朋友夫妻互换小说| 亚洲精品日本国产第一区| 三级男女做爰猛烈吃奶摸视频| 青青草视频在线视频观看| 久热这里只有精品99| av又黄又爽大尺度在线免费看| 国产久久久一区二区三区| 欧美日韩视频精品一区| 又爽又黄无遮挡网站| 亚洲一区二区三区欧美精品 | 亚洲欧美精品专区久久| 国产精品一二三区在线看| 国产永久视频网站| 黄色视频在线播放观看不卡| 亚洲精品aⅴ在线观看| 春色校园在线视频观看| 国产有黄有色有爽视频| 日韩免费高清中文字幕av| 白带黄色成豆腐渣| 一个人观看的视频www高清免费观看| 国产成人一区二区在线| 久久人人爽av亚洲精品天堂 | 亚洲真实伦在线观看| 国产精品久久久久久精品电影小说 | av在线蜜桃| 性色avwww在线观看| 国产黄片美女视频| 精品久久久久久久末码| 色网站视频免费| 一区二区三区四区激情视频| 久久鲁丝午夜福利片| 在线观看美女被高潮喷水网站| 国产精品人妻久久久久久| 麻豆国产97在线/欧美| 午夜精品一区二区三区免费看| 色吧在线观看| 国产日韩欧美在线精品| 国产成人精品婷婷| av专区在线播放| 大片电影免费在线观看免费| 男女无遮挡免费网站观看| 日韩制服骚丝袜av| 天堂俺去俺来也www色官网| 天堂中文最新版在线下载 | 成年av动漫网址| 街头女战士在线观看网站| 秋霞在线观看毛片| 2021天堂中文幕一二区在线观| 亚洲精品成人av观看孕妇| 日韩一本色道免费dvd| av女优亚洲男人天堂| 大话2 男鬼变身卡| 国产人妻一区二区三区在| 国产成人精品久久久久久| 日本一二三区视频观看| 欧美日韩一区二区视频在线观看视频在线 | 日本三级黄在线观看| 欧美bdsm另类| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区视频9| 国产精品蜜桃在线观看| 成人国产麻豆网| 成年av动漫网址| 又大又黄又爽视频免费| 日日啪夜夜爽| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 一本色道久久久久久精品综合| 老女人水多毛片| 中国三级夫妇交换| 真实男女啪啪啪动态图| 亚洲自偷自拍三级| 日韩电影二区| 亚洲综合色惰| 蜜桃亚洲精品一区二区三区| 边亲边吃奶的免费视频| 18禁在线无遮挡免费观看视频| 麻豆成人av视频| 国产国拍精品亚洲av在线观看| 女的被弄到高潮叫床怎么办| 你懂的网址亚洲精品在线观看| 色视频在线一区二区三区| 亚洲av中文字字幕乱码综合| 在线观看国产h片| 欧美高清成人免费视频www| 激情 狠狠 欧美| 日日摸夜夜添夜夜添av毛片| 综合色丁香网| 少妇高潮的动态图| 日韩欧美精品免费久久| 身体一侧抽搐| 99热这里只有是精品在线观看| 亚洲一区二区三区欧美精品 | 草草在线视频免费看| 成年av动漫网址| 亚洲天堂av无毛| 晚上一个人看的免费电影| 男女下面进入的视频免费午夜| 亚洲美女搞黄在线观看| 99久久中文字幕三级久久日本| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲丝袜综合中文字幕| 综合色av麻豆| 免费看光身美女| 一二三四中文在线观看免费高清| 国产精品一区二区在线观看99| 日韩欧美精品免费久久| 两个人的视频大全免费| 97精品久久久久久久久久精品| 免费高清在线观看视频在线观看| 欧美成人精品欧美一级黄| 亚洲精品成人久久久久久| 狠狠精品人妻久久久久久综合| 视频区图区小说| 伊人久久国产一区二区| 国产亚洲91精品色在线| 97超碰精品成人国产| 欧美性猛交╳xxx乱大交人| 99re6热这里在线精品视频| 久久99热这里只频精品6学生| 欧美激情久久久久久爽电影| 一区二区三区乱码不卡18| 永久网站在线| 在线免费观看不下载黄p国产| freevideosex欧美| 成人亚洲精品av一区二区| 亚洲美女视频黄频| 国产高清三级在线| 色哟哟·www| 国产成人午夜福利电影在线观看| 免费观看性生交大片5| h日本视频在线播放| 一级毛片我不卡| 午夜福利在线观看免费完整高清在| 色婷婷久久久亚洲欧美| 蜜臀久久99精品久久宅男| 热re99久久精品国产66热6| 亚洲熟女精品中文字幕| 97在线人人人人妻| 有码 亚洲区| 日日摸夜夜添夜夜添av毛片| 一个人观看的视频www高清免费观看| 黄色日韩在线| 精品一区二区三区视频在线| 国产精品偷伦视频观看了| 极品教师在线视频| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 晚上一个人看的免费电影| www.色视频.com| 国产一区有黄有色的免费视频| 自拍欧美九色日韩亚洲蝌蚪91 | 中文在线观看免费www的网站| 又粗又硬又长又爽又黄的视频| 人人妻人人看人人澡| 嫩草影院新地址| 亚洲经典国产精华液单| 国产精品嫩草影院av在线观看| 亚洲一区二区三区欧美精品 | 久久久久久伊人网av| 男人添女人高潮全过程视频| 狂野欧美激情性xxxx在线观看| 国产探花极品一区二区| 在线观看免费高清a一片| 国产精品久久久久久久电影| 水蜜桃什么品种好| 国产 一区精品| 别揉我奶头 嗯啊视频| 五月玫瑰六月丁香| 成年女人在线观看亚洲视频 | 亚洲经典国产精华液单| 亚洲最大成人手机在线| 汤姆久久久久久久影院中文字幕| 18禁动态无遮挡网站| av免费在线看不卡| 可以在线观看毛片的网站| av在线天堂中文字幕| www.av在线官网国产| 日本一二三区视频观看| 在线观看av片永久免费下载| 久热久热在线精品观看| 日本wwww免费看| 日本一本二区三区精品| av一本久久久久| 岛国毛片在线播放| 日韩免费高清中文字幕av| 男女边摸边吃奶| 婷婷色av中文字幕| 免费电影在线观看免费观看| 少妇的逼水好多| 最近中文字幕高清免费大全6| 99精国产麻豆久久婷婷| 国产黄色视频一区二区在线观看| 婷婷色综合大香蕉| 国产色爽女视频免费观看| 亚洲最大成人av| 白带黄色成豆腐渣| 亚洲国产精品成人久久小说| 国产真实伦视频高清在线观看| 99久久人妻综合| 亚洲欧洲日产国产| 黄色一级大片看看| 亚洲精品日本国产第一区| 街头女战士在线观看网站| 97在线人人人人妻| 91久久精品国产一区二区三区| 99久久九九国产精品国产免费| 久久久色成人| 国产黄频视频在线观看| av黄色大香蕉| 噜噜噜噜噜久久久久久91| 亚洲国产日韩一区二区| 男人和女人高潮做爰伦理| 亚洲一区二区三区欧美精品 | 爱豆传媒免费全集在线观看| 别揉我奶头 嗯啊视频| 亚洲国产精品专区欧美| 精品亚洲乱码少妇综合久久| 久久久成人免费电影| 极品少妇高潮喷水抽搐| 亚洲国产欧美在线一区| 特级一级黄色大片| 亚洲va在线va天堂va国产| 日韩av不卡免费在线播放| 欧美 日韩 精品 国产| 亚洲国产欧美在线一区| 国产探花极品一区二区| 国产精品熟女久久久久浪| 国产真实伦视频高清在线观看| 免费黄频网站在线观看国产| 美女xxoo啪啪120秒动态图| 三级男女做爰猛烈吃奶摸视频| 97精品久久久久久久久久精品| 国产精品av视频在线免费观看| 国产伦精品一区二区三区四那| 亚洲国产精品成人综合色| 亚洲av欧美aⅴ国产| 中文字幕免费在线视频6| 日韩av免费高清视频| 国产毛片a区久久久久| 午夜爱爱视频在线播放| 国内揄拍国产精品人妻在线| 成人欧美大片| 欧美成人a在线观看| 22中文网久久字幕| 精品99又大又爽又粗少妇毛片| 国产男女内射视频| 国产白丝娇喘喷水9色精品| 水蜜桃什么品种好| 99久久精品热视频| 91精品一卡2卡3卡4卡| 亚洲综合精品二区| 深夜a级毛片| 欧美97在线视频| 日日摸夜夜添夜夜爱| 麻豆久久精品国产亚洲av| 国产v大片淫在线免费观看| 大片电影免费在线观看免费| 99久久九九国产精品国产免费| 日韩 亚洲 欧美在线| 久久久成人免费电影| 欧美日本视频| 高清欧美精品videossex| 久久国产乱子免费精品| 国产探花在线观看一区二区| 久久韩国三级中文字幕| 2022亚洲国产成人精品| 又黄又爽又刺激的免费视频.| 国产精品爽爽va在线观看网站| 99视频精品全部免费 在线| 青春草国产在线视频| 黑人高潮一二区| 伊人久久精品亚洲午夜| 亚洲欧美中文字幕日韩二区| 少妇丰满av| 人人妻人人看人人澡| 最近手机中文字幕大全| 亚洲精品日韩av片在线观看| 一级毛片aaaaaa免费看小| 亚洲精品影视一区二区三区av| 午夜福利在线在线| tube8黄色片| 天美传媒精品一区二区| 亚洲精品aⅴ在线观看| 别揉我奶头 嗯啊视频| 深爱激情五月婷婷| 成人毛片60女人毛片免费| 中文字幕免费在线视频6| 国产亚洲最大av| 免费观看av网站的网址| 欧美 日韩 精品 国产| 久久精品综合一区二区三区| 九九在线视频观看精品| 香蕉精品网在线| 一级毛片电影观看| 夜夜看夜夜爽夜夜摸| 18禁动态无遮挡网站| 日韩强制内射视频| 免费观看在线日韩| 狠狠精品人妻久久久久久综合| 激情 狠狠 欧美| 国产免费福利视频在线观看| 欧美高清性xxxxhd video| 欧美性感艳星| 亚洲精品aⅴ在线观看| 伊人久久国产一区二区| 校园人妻丝袜中文字幕| 国内少妇人妻偷人精品xxx网站| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频| 永久网站在线| 国产久久久一区二区三区| 超碰av人人做人人爽久久| 男的添女的下面高潮视频| videos熟女内射| 欧美激情在线99| 国产亚洲5aaaaa淫片| 联通29元200g的流量卡| 午夜老司机福利剧场| eeuss影院久久| av专区在线播放|