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

    小變量情況下第一類整數(shù)階Bessel函數(shù)的計(jì)算

    2016-12-09 06:34:14馮佳計(jì)賈曉偉沈建琪
    電子學(xué)報(bào) 2016年11期
    關(guān)鍵詞:展開式相位角階數(shù)

    馮佳計(jì),賈曉偉,沈建琪

    (上海理工大學(xué)理學(xué)院,上海200093)

    ?

    小變量情況下第一類整數(shù)階Bessel函數(shù)的計(jì)算

    馮佳計(jì),賈曉偉,沈建琪

    (上海理工大學(xué)理學(xué)院,上海200093)

    在計(jì)算第一類整數(shù)階Bessel函數(shù)時(shí),后向遞推算法穩(wěn)定高效.然而,起始點(diǎn)的選取必須有足夠高的階數(shù),并且需要進(jìn)行歸一化處理.本文對(duì)Taylor級(jí)數(shù)展開算法進(jìn)行研究,并對(duì)其級(jí)數(shù)展開規(guī)律、計(jì)算精度,以及求和項(xiàng)與參數(shù)間的關(guān)系進(jìn)行了討論.此外,本文利用指數(shù)形式,極大擴(kuò)展了該算法的可計(jì)算范圍.與du Toit算法、MATLAB和Mathematica應(yīng)用軟件的計(jì)算結(jié)果比較顯示,本文的算法具有較高的準(zhǔn)確性和穩(wěn)定性.

    Bessel函數(shù);Taylor級(jí)數(shù)展開;指數(shù)擴(kuò)展

    1 引言

    第一類整數(shù)階Bessel函數(shù)Jn(z)在物理學(xué)、力學(xué)、電磁學(xué)、環(huán)境科學(xué)等領(lǐng)域都有著重要的用途[1,2].由于應(yīng)用廣泛,大量文獻(xiàn)致力于討論它的性質(zhì)并給出了多種不同的算法[3~16].其中后向遞推算法選取兩個(gè)連續(xù)的高階函數(shù)作為起始點(diǎn)賦任意值并后向遞推計(jì)算,最后做歸一化處理[9~11].這種算法穩(wěn)定且高效.然而,為了確保計(jì)算的順利進(jìn)行,后向遞推算法起始點(diǎn)的階數(shù)必須足夠高,需要更多的計(jì)算機(jī)資源和更大的計(jì)算量.在采用后向遞推算法計(jì)算Bessel函數(shù)時(shí),可考慮采用Taylor級(jí)數(shù)展開法精確計(jì)算出起始階數(shù)的函數(shù)值以提高計(jì)算效率.原則上,Bessel函數(shù)的Taylor級(jí)數(shù)展開對(duì)所有變量均收斂.然而,其收斂速度卻依賴于階數(shù)n和變量z的大小.在階數(shù)n很小、參數(shù)z很大的情況下,需要計(jì)算的級(jí)數(shù)展開項(xiàng)增多,且在計(jì)算過程中有可能會(huì)丟失有效數(shù)字,導(dǎo)致計(jì)算效率降低、計(jì)算結(jié)果出現(xiàn)偏差.

    2 遞推算法

    整數(shù)階Bessel函數(shù)滿足前向遞推關(guān)系

    (1)

    或者后向遞推關(guān)系

    (2)

    對(duì)于給定參數(shù)z,前向遞推從兩個(gè)連續(xù)低階級(jí)數(shù)開始計(jì)算,而后向遞推從兩個(gè)連續(xù)高階級(jí)數(shù)開始計(jì)算.像大多數(shù)遞推過程一樣,上述公式易受傳播誤差的影響[7].

    當(dāng)z為實(shí)數(shù)或者虛部很小(|lmz|≤1.0)的復(fù)數(shù)時(shí),Bessel函數(shù)的絕對(duì)值|Jn(z)|在n<0.9|z|范圍變化不大,而在n≥0.9|z|范圍隨階數(shù)的增大快速衰減.對(duì)于n<0.9|z|的范圍,前向遞推或者后向遞推都能用來計(jì)算Jn(z);但當(dāng)n>0.9|z|時(shí),則只能采用后向遞推法[5,9].如果z的虛部較大,Bessel函數(shù)的絕對(duì)值|Jn(z)|在n≥0范圍內(nèi)隨階數(shù)的增大迅速衰減,只能采用后向遞推法計(jì)算.

    用后向遞推方法計(jì)算Bessel函數(shù)Jn(z)需從兩個(gè)連續(xù)的高階級(jí)數(shù)Bq(z)和Bq+1(z)開始,Bq(z)和Bq+1(z)的值可以任意選取,du.Toit建議Bq(z)=1和Bq+1(z)=0[9].然而,階數(shù)q必須同時(shí)滿足env|Jq+1(z)|≤10-p/2env|JN(z)|和env|2Jq(z)|≤10-penv|J0(z)|以確保收斂,其中p為所需的有效位數(shù)、N是所需計(jì)算的最高階數(shù),env|x|表示函數(shù)的包絡(luò)線.如果N≤|z|,q的最小值需滿足

    (3)

    若N>|z|,q值的選取則相對(duì)繁瑣[5,11].完成Bn(z)的計(jì)算后,需進(jìn)行歸一化處理得到相應(yīng)的函數(shù)值Jn(z).針對(duì)不同的變量值,采用不同的歸一化方法.當(dāng)|z|≤25時(shí),Bn(z)歸一化處理為

    (4)

    當(dāng)|z|>25時(shí),使用S=B0(z)/J0(z)或S=B1(z)/J1(z)歸一化效率更高.使用時(shí)應(yīng)避開J0(z)和J1(z)的零點(diǎn).

    從以上討論可知,后向遞推法從兩個(gè)任意取值Bq(z)和Bq+1(z)開始計(jì)算,不僅需要計(jì)算起始階數(shù)q以及從階數(shù)q到階數(shù)N的函數(shù)Bn(z),還需對(duì)Bn(z)做歸一化處理.因此,為了使計(jì)算更加高效,在使用后向遞推方法時(shí),建議采用Taylor級(jí)數(shù)展開法預(yù)先計(jì)算Bq(z)和Bq+1(z)的精確值.

    3 Bessel函數(shù)的Taylor級(jí)數(shù)展開

    整數(shù)階Bessel函數(shù)Jn(z)的Taylor級(jí)數(shù)展開式為:

    (5)

    (6)

    (7)

    則有遞推關(guān)系:

    (8)

    (9)

    (10)

    對(duì)于復(fù)變量的一般情況z=|z|eiδz,Taylor級(jí)數(shù)展開形式寫成如下形式:

    (11)

    (12)

    (13)

    (14)

    公式(13)中的求和項(xiàng)在計(jì)算過程中出現(xiàn)正負(fù)交替的情況,如果kpeak比較大,必然導(dǎo)致數(shù)值比較大的數(shù)字相減并導(dǎo)致有效數(shù)字丟失從而計(jì)算結(jié)果變差.反之,公式(14)中的求和項(xiàng)一直為正,不會(huì)出現(xiàn)公式(13)的情況.因此,在這種情況下,即使kpeak比較大,也能得到合理的計(jì)算結(jié)果,但求和項(xiàng)比較多、收斂較慢.

    (15)

    在采用Taylor級(jí)數(shù)展開式時(shí),并不一定要求嚴(yán)格滿足(|z|/2)2

    相位角δz=0時(shí)有效數(shù)字丟失最嚴(yán)重,因此在該相位角情況下滿足S-T<18的所有階數(shù)ncal對(duì)所有的相位角計(jì)算均是可信的(見圖5).可以看出,ncal≈|z|2/40-10表示滿足n>ncal時(shí),對(duì)于任意相位角δz的計(jì)算結(jié)果保證至少7位有效數(shù),這與n>|z|2/4-1的要求降低了大約10倍.

    圖6給出了|z|=600時(shí)階數(shù)n和相位角δz的關(guān)系,以曲線為界,上部區(qū)域表示Taylor級(jí)數(shù)展開式適用的范圍,下部區(qū)域中計(jì)算得到的結(jié)果其有效位數(shù)低于7位.顯然,如果提高對(duì)有效位數(shù)的要求,則圖中的曲線將會(huì)向上移動(dòng),對(duì)應(yīng)的可計(jì)算范圍縮小.

    在20<|z|<104范圍內(nèi)對(duì)多個(gè)不同的復(fù)變量z進(jìn)行了數(shù)值計(jì)算,并進(jìn)行曲線擬合,得到如下經(jīng)驗(yàn)公式:

    (16)

    a=ln(|z|-0.5232|z|0.3+20.49)-2.997

    (17)

    (18)

    4 復(fù)數(shù)的大數(shù)計(jì)算方法

    在復(fù)變量情況下,Bessel函數(shù)值及其Taylor級(jí)數(shù)展開式的求和項(xiàng)在一個(gè)很大的動(dòng)態(tài)范圍內(nèi)變化.為了保證計(jì)算的順利進(jìn)行,需要對(duì)中間計(jì)算量的數(shù)值范圍進(jìn)行擴(kuò)展.為此,所有變量均采用復(fù)數(shù)形式儲(chǔ)存并參與運(yùn)算.

    (19)

    5 數(shù)值計(jì)算結(jié)果和分析

    為了驗(yàn)證本文提出的算法,我們將所得結(jié)果與du

    表1 Jn(z)數(shù)據(jù)比較

    6 結(jié)論

    本文介紹了第一類整數(shù)階Bessel函數(shù)的Taylor級(jí)數(shù)展開算法.為了保證高速收斂并且具有較高的計(jì)算精度,Taylor級(jí)數(shù)展開應(yīng)在n≥|z|2/4-1條件下使用;在保證7位有效數(shù)計(jì)算精度的情況下,給出了Taylor級(jí)數(shù)展開式適用范圍的經(jīng)驗(yàn)公式;此外,在算法中采用了指數(shù)形式以擴(kuò)展計(jì)算數(shù)值范圍.計(jì)算結(jié)果與現(xiàn)有的軟件進(jìn)行了比較,表明該算法具有較好的效果.

    遞推算法與級(jí)數(shù)展開算法結(jié)合使用可有效提高計(jì)算效率.此外,本文中給出的算法也可以應(yīng)用于計(jì)算非整數(shù)階Bessel函數(shù)以及改進(jìn)的Bessel函數(shù).

    [1]G Gouesbet,G Gréhan.Generalized Lorenz-Mie Theories[M].Springer-Verlag Berlin Heidel-berg,2011.

    [2]JA Lock.Improved Gaussian beam-scattering algorithm[J].Applied Optics,1995,34 (3):599-570.

    [3]GN Watson.A Treatise on the Theory of Bessel Functions[M].Cambridge University Press,Cambridge,1958.

    [4]M Abramowitz,IA Stegun.Handbook of Mathematical Functions[M].10th ed.,National Bureau of Standards,Washington,DC,1972.

    [5]S Zhang,J Jin.Computation of Special Functions[M].Wiley,New York,1996.

    [6]M Goldstein,RM Thaler.Recurrence techniques for the calculation of Bessel functions[J].Mathematical Tables & Other Aids to Computation,1959,13(66):102-102.

    [7]W Gautschi.Computational aspects of three-term recurrence relations[J].Society for Industrial and Applied Mathematics,1967,9(1):24-28.

    [8]DE Amos.Algorithm 644:A portable package for Bessel functions of a complex argument and nonnegative order[J].Acm Transactions on Mathematical Software,1986,12:265-273.

    [9]CF du Toit.The numerical computation of Bessel functions of the first and second kind for integer orders and complex arguments[J],IEEE Antennas & Propagation Magazine,1990,38(9):1341-1349.

    [10]CF du Toit.Bessel functions Jn(z) and Yn(z) for integer order and complex argument[J].Computer Physics Communications,1993,78(1-2):181-189.

    [11]CF du Toit.Evaluation of some algorithms and programs for the computation of integer-order Bessel functions of the first and second kind with complex arguments[J].IEEE Antennas & Propagation Magazine,1993,35(3):19-25.

    [12]KG Valeev,OY Kostinskii.Calculation of Bessel functions by using continued fractions[J].Ukrainian Mathematical Journal,1997,47(12):1949-1950.

    [13]HA Yousif,R Melka.Bessel function of the first kind with complex argument[J].Computer Physics Communications,1997,106(3):199-206.

    [14]EJ Rothwell.Computation of the logarithm of Bessel functions of complex argument[J].Communications in Numerical Methods in Engineering,2005,21(10):597-605.

    [15]LW Cai.On the computation of spherical Bessel functions of complex arguments[J].Computer Physics Communications,2011,182(3):663-668.

    [16]張善杰,唐漢.任意實(shí)數(shù)階復(fù)宗量第一類和第二類Bessel函數(shù)的精確計(jì)算[J].電子學(xué)報(bào),1996,24(3):77-80.

    Zhang Shanjie,Tang Han.Accurate computation of Bessel functions of the first and second kinds with arbitrary real order and complex argument[J].Acta Electronica Sinica,1996,24(3):77-80.(in Chinese)

    馮佳計(jì) 男,碩士研究生,1990年5月生于江西省撫州市東鄉(xiāng)縣.主要從事光學(xué)測(cè)試方面的研究.

    E-mail:jiaji19900505@163.com

    沈建琪(通信作者) 男,1965年7月出生,浙江桐鄉(xiāng)人,博士生導(dǎo)師,分別于1987年和1990年在華東師范大學(xué)獲得理學(xué)學(xué)士和碩士學(xué)位,2003年在德國(guó)Cottbus工大獲得工學(xué)博士學(xué)位.主要從事光散射顆粒測(cè)試技術(shù)研究.

    E-mail:jqshenk@163.com

    Computation of the Integer Order Bessel Functions of First Kind with Small Arguments

    FENG Jia-ji,JIA Xiao-wei,SHEN Jian-qi

    (CollegeofScience,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China)

    Algorithm based on the backward recurrence for computing the integer order Bessel functions of the first kind is stable and efficient.However,the orders of the starting points should be high enough and the normalization is required.In this paper,we introduce an algorithm based on the Taylor series expansion (TSE),in which all the quantities involved are expressed in the exponential format so as to expand the numeric range of calculation.Comparison against du Toit’s algorithm as well as MATLAB and Mathematica shows that our algorithm is stable and reliable.

    Bessel function;Taylor series expansion;exponential scaling

    2015-11-24;

    2016-02-25;責(zé)任編輯:藍(lán)紅杰

    國(guó)家自然科學(xué)基金(No.NSFC 51476104)

    TP301.6

    A

    0372-2112 (2016)11-2720-06

    ??學(xué)報(bào)URL:http://www.ejournal.org.cn

    10.3969/j.issn.0372-2112.2016.11.022

    猜你喜歡
    展開式相位角階數(shù)
    例說二項(xiàng)展開式的通項(xiàng)公式應(yīng)用
    關(guān)于無窮小階數(shù)的幾點(diǎn)注記
    確定有限級(jí)數(shù)解的階數(shù)上界的一種n階展開方法
    泰勒展開式在函數(shù)中的應(yīng)用
    開封地區(qū)健康成人相位角及其影響因素
    相位角對(duì)容性耦合電非對(duì)稱放電特性的影響?
    共軸剛性旋翼直升機(jī)旋翼控制相位角問題分析
    精確檢測(cè)電網(wǎng)電壓矢量相位角方法的研究
    對(duì)一道冪級(jí)數(shù)展開式例題的思考
    一種新的多址信道有效階數(shù)估計(jì)算法*
    国产av精品麻豆| 秋霞在线观看毛片| 日韩中文字幕视频在线看片 | 女的被弄到高潮叫床怎么办| 少妇人妻 视频| 亚洲无线观看免费| 国产淫片久久久久久久久| 亚洲综合色惰| av天堂中文字幕网| 在线观看免费视频网站a站| 欧美3d第一页| 舔av片在线| www.色视频.com| 你懂的网址亚洲精品在线观看| 久久亚洲国产成人精品v| 亚洲图色成人| 在线免费观看不下载黄p国产| 久久 成人 亚洲| 国产一区亚洲一区在线观看| 国产午夜精品久久久久久一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久久久久久大尺度免费视频| 在线 av 中文字幕| 在线观看美女被高潮喷水网站| 搡女人真爽免费视频火全软件| 精品人妻偷拍中文字幕| 国产成人freesex在线| 亚洲第一av免费看| 欧美成人精品欧美一级黄| 国产老妇伦熟女老妇高清| 国产综合精华液| 欧美精品国产亚洲| 又大又黄又爽视频免费| 亚洲精品亚洲一区二区| 伦精品一区二区三区| 最近最新中文字幕大全电影3| 久久精品国产亚洲网站| 午夜老司机福利剧场| 国产精品女同一区二区软件| 日韩中文字幕视频在线看片 | 女性被躁到高潮视频| 嫩草影院新地址| 免费不卡的大黄色大毛片视频在线观看| 国产精品偷伦视频观看了| 精品熟女少妇av免费看| 免费看光身美女| 免费少妇av软件| 美女高潮的动态| 国模一区二区三区四区视频| 伦理电影免费视频| 爱豆传媒免费全集在线观看| 在线观看免费日韩欧美大片 | 久久 成人 亚洲| 我要看日韩黄色一级片| 国产成人91sexporn| 国产 一区精品| 日韩av不卡免费在线播放| 高清午夜精品一区二区三区| 国产一区二区三区综合在线观看 | 亚洲国产色片| 97热精品久久久久久| 免费观看性生交大片5| 免费久久久久久久精品成人欧美视频 | 啦啦啦视频在线资源免费观看| 亚洲精品乱久久久久久| 少妇丰满av| 一区二区三区精品91| 亚洲精品aⅴ在线观看| 国产精品久久久久久久电影| 免费大片黄手机在线观看| 国产在线免费精品| 老师上课跳d突然被开到最大视频| 精品一区二区三区视频在线| 91在线精品国自产拍蜜月| 色吧在线观看| 亚洲av国产av综合av卡| 欧美激情极品国产一区二区三区 | 能在线免费看毛片的网站| 中国三级夫妇交换| 欧美亚洲 丝袜 人妻 在线| 国产在线免费精品| 伊人久久国产一区二区| 在现免费观看毛片| 免费av中文字幕在线| 国产精品av视频在线免费观看| 国产精品一区二区性色av| 精品国产露脸久久av麻豆| 久久久久国产网址| 日韩电影二区| 夜夜爽夜夜爽视频| 国产久久久一区二区三区| 亚洲不卡免费看| 亚洲av在线观看美女高潮| 精品一品国产午夜福利视频| 成年女人在线观看亚洲视频| 极品少妇高潮喷水抽搐| 九九爱精品视频在线观看| 欧美3d第一页| 色视频www国产| 亚洲精品成人av观看孕妇| 人体艺术视频欧美日本| 岛国毛片在线播放| 麻豆成人av视频| 在线播放无遮挡| 小蜜桃在线观看免费完整版高清| 一级爰片在线观看| 亚洲第一av免费看| 欧美激情极品国产一区二区三区 | 麻豆精品久久久久久蜜桃| 免费在线观看成人毛片| 国产黄色视频一区二区在线观看| 下体分泌物呈黄色| 成年人午夜在线观看视频| 国语对白做爰xxxⅹ性视频网站| 国产成人免费无遮挡视频| 国产黄片视频在线免费观看| 搡老乐熟女国产| 美女xxoo啪啪120秒动态图| 国产精品免费大片| 夫妻午夜视频| 男女国产视频网站| 九色成人免费人妻av| 少妇人妻 视频| 亚洲美女视频黄频| 久久久久久伊人网av| 国产成人免费观看mmmm| 日韩伦理黄色片| 精华霜和精华液先用哪个| 国产永久视频网站| 欧美xxⅹ黑人| 亚洲欧美一区二区三区黑人 | 深爱激情五月婷婷| 丰满乱子伦码专区| 亚洲国产精品一区三区| 国产视频首页在线观看| 色哟哟·www| 激情 狠狠 欧美| 成人影院久久| 亚洲精品国产成人久久av| 亚洲精品国产色婷婷电影| 久久久久人妻精品一区果冻| 中国国产av一级| 亚洲av电影在线观看一区二区三区| 美女xxoo啪啪120秒动态图| 午夜免费鲁丝| 日本爱情动作片www.在线观看| 亚州av有码| 国产精品久久久久成人av| 亚洲国产欧美在线一区| 国产精品人妻久久久影院| 777米奇影视久久| 国产精品欧美亚洲77777| 99九九线精品视频在线观看视频| 久久人人爽av亚洲精品天堂 | 成人综合一区亚洲| 欧美变态另类bdsm刘玥| 91狼人影院| 亚洲色图av天堂| 在线天堂最新版资源| 在线看a的网站| 亚洲电影在线观看av| 日本免费在线观看一区| 久久人妻熟女aⅴ| 汤姆久久久久久久影院中文字幕| 亚洲精品日韩av片在线观看| 26uuu在线亚洲综合色| 亚洲美女视频黄频| 亚洲精品色激情综合| 永久免费av网站大全| 亚洲最大成人中文| 九九在线视频观看精品| 高清黄色对白视频在线免费看 | 国产精品.久久久| 国产精品久久久久久久电影| 久久国产精品大桥未久av | 一个人看视频在线观看www免费| h日本视频在线播放| 日本黄大片高清| 婷婷色麻豆天堂久久| 免费人妻精品一区二区三区视频| 亚洲成人中文字幕在线播放| 搡女人真爽免费视频火全软件| 亚洲怡红院男人天堂| 啦啦啦在线观看免费高清www| 久久ye,这里只有精品| 丰满迷人的少妇在线观看| 免费观看无遮挡的男女| 我要看黄色一级片免费的| 少妇人妻精品综合一区二区| 内地一区二区视频在线| 成人午夜精彩视频在线观看| 十分钟在线观看高清视频www | 女的被弄到高潮叫床怎么办| 欧美成人精品欧美一级黄| 干丝袜人妻中文字幕| 韩国av在线不卡| 三级国产精品片| 青春草国产在线视频| 亚洲三级黄色毛片| 在线看a的网站| 欧美日韩视频精品一区| 亚洲一区二区三区欧美精品| 亚洲精品视频女| 久久久亚洲精品成人影院| 中文字幕av成人在线电影| 午夜福利网站1000一区二区三区| 两个人的视频大全免费| 欧美日韩综合久久久久久| 亚洲欧美成人精品一区二区| 91狼人影院| 伦理电影大哥的女人| .国产精品久久| 欧美三级亚洲精品| 国产女主播在线喷水免费视频网站| 欧美精品人与动牲交sv欧美| 日韩欧美 国产精品| 国产乱人视频| 日日啪夜夜撸| 亚洲熟女精品中文字幕| 3wmmmm亚洲av在线观看| 99久久精品热视频| 一本—道久久a久久精品蜜桃钙片| 十分钟在线观看高清视频www | 九九在线视频观看精品| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 激情五月婷婷亚洲| 国语对白做爰xxxⅹ性视频网站| 国产中年淑女户外野战色| 中文字幕免费在线视频6| 黄色欧美视频在线观看| 成人高潮视频无遮挡免费网站| 亚洲av日韩在线播放| 一区在线观看完整版| 啦啦啦在线观看免费高清www| 在线观看美女被高潮喷水网站| 中文字幕精品免费在线观看视频 | 夜夜骑夜夜射夜夜干| 三级国产精品片| 久久国产精品男人的天堂亚洲 | 欧美一级a爱片免费观看看| 人妻制服诱惑在线中文字幕| 日韩一本色道免费dvd| 99久久精品热视频| 国产精品精品国产色婷婷| 精品久久久精品久久久| 麻豆乱淫一区二区| 久久久久久久亚洲中文字幕| 伦精品一区二区三区| 男女啪啪激烈高潮av片| 日韩中字成人| 免费久久久久久久精品成人欧美视频 | 国产成人午夜福利电影在线观看| 欧美精品人与动牲交sv欧美| 日日撸夜夜添| 天堂8中文在线网| 国产人妻一区二区三区在| av在线app专区| av专区在线播放| 久久国产精品大桥未久av | 免费av不卡在线播放| 久久久午夜欧美精品| 色吧在线观看| 久久97久久精品| 国产 精品1| 乱码一卡2卡4卡精品| 欧美人与善性xxx| 男的添女的下面高潮视频| 日韩制服骚丝袜av| 久久综合国产亚洲精品| 精品一区二区三区视频在线| 蜜桃久久精品国产亚洲av| 97超碰精品成人国产| 一级毛片黄色毛片免费观看视频| av免费观看日本| 国国产精品蜜臀av免费| 精品久久久精品久久久| 久热久热在线精品观看| 国产欧美日韩精品一区二区| .国产精品久久| 亚洲精品日韩av片在线观看| 亚洲欧洲日产国产| 欧美成人a在线观看| 性色av一级| 亚洲av中文字字幕乱码综合| 九九在线视频观看精品| 看非洲黑人一级黄片| 国产精品一区二区在线不卡| 国产精品嫩草影院av在线观看| 国产深夜福利视频在线观看| 精品少妇黑人巨大在线播放| 在线天堂最新版资源| 精品久久久精品久久久| 妹子高潮喷水视频| 亚洲人成网站在线观看播放| 日本黄大片高清| 99热网站在线观看| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲网站| 免费人成在线观看视频色| 国产av精品麻豆| 又大又黄又爽视频免费| 老司机影院成人| 秋霞伦理黄片| 最近中文字幕2019免费版| 久久久久精品久久久久真实原创| 欧美一区二区亚洲| 午夜激情福利司机影院| 亚洲三级黄色毛片| 免费观看在线日韩| 亚洲综合色惰| 六月丁香七月| 狠狠精品人妻久久久久久综合| 久久av网站| 蜜臀久久99精品久久宅男| 麻豆国产97在线/欧美| 国产精品熟女久久久久浪| 97在线人人人人妻| 亚洲国产精品999| 亚洲精品视频女| 黑人猛操日本美女一级片| 亚洲国产色片| 男女无遮挡免费网站观看| 成年美女黄网站色视频大全免费 | 亚洲精品456在线播放app| 欧美日本视频| 久久鲁丝午夜福利片| 一个人看的www免费观看视频| 观看av在线不卡| 黄色怎么调成土黄色| 国产视频内射| 尾随美女入室| 日本与韩国留学比较| 久久影院123| a级毛色黄片| 国产精品不卡视频一区二区| 久久久午夜欧美精品| 亚洲精品国产成人久久av| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 永久免费av网站大全| 高清av免费在线| av女优亚洲男人天堂| 精品一品国产午夜福利视频| 两个人的视频大全免费| 一级片'在线观看视频| 男人爽女人下面视频在线观看| 日韩视频在线欧美| 纵有疾风起免费观看全集完整版| 中文精品一卡2卡3卡4更新| 国产精品久久久久成人av| 少妇精品久久久久久久| 国产有黄有色有爽视频| 各种免费的搞黄视频| 中文字幕精品免费在线观看视频 | 在线观看人妻少妇| 免费观看av网站的网址| 又粗又硬又长又爽又黄的视频| 18禁裸乳无遮挡免费网站照片| 乱码一卡2卡4卡精品| 免费观看av网站的网址| 啦啦啦在线观看免费高清www| 国产精品久久久久久精品电影小说 | 国产一区二区三区综合在线观看 | 亚洲国产av新网站| 日韩欧美精品免费久久| tube8黄色片| 男女免费视频国产| 亚洲欧洲日产国产| 久久久久性生活片| 国产伦在线观看视频一区| 日本av免费视频播放| 久久青草综合色| 成年av动漫网址| 最近中文字幕2019免费版| 一本色道久久久久久精品综合| 国语对白做爰xxxⅹ性视频网站| 一个人免费看片子| 大话2 男鬼变身卡| 一级毛片电影观看| 色视频www国产| 黄色一级大片看看| 精品一区二区免费观看| 毛片女人毛片| 久久久精品免费免费高清| 国产黄色免费在线视频| 18禁在线播放成人免费| 嘟嘟电影网在线观看| 国产精品一区二区在线观看99| 国产精品一及| 五月天丁香电影| 久热这里只有精品99| 国产精品成人在线| 91精品国产国语对白视频| 国产色婷婷99| 汤姆久久久久久久影院中文字幕| 少妇裸体淫交视频免费看高清| 纯流量卡能插随身wifi吗| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看| 久久久久精品性色| 97超碰精品成人国产| 精华霜和精华液先用哪个| 国产精品久久久久久av不卡| 国产精品麻豆人妻色哟哟久久| 91久久精品电影网| 综合色丁香网| 尤物成人国产欧美一区二区三区| av在线蜜桃| 国产视频首页在线观看| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 大片电影免费在线观看免费| 三级国产精品片| 成人免费观看视频高清| 精品午夜福利在线看| 免费黄色在线免费观看| 久久鲁丝午夜福利片| 亚洲精品乱久久久久久| 亚洲精品第二区| 啦啦啦中文免费视频观看日本| 又爽又黄a免费视频| 欧美精品人与动牲交sv欧美| 超碰av人人做人人爽久久| 亚洲不卡免费看| 麻豆成人av视频| 五月伊人婷婷丁香| av国产精品久久久久影院| 午夜激情福利司机影院| www.av在线官网国产| 色吧在线观看| 亚洲国产成人一精品久久久| 91午夜精品亚洲一区二区三区| 久久韩国三级中文字幕| 国产深夜福利视频在线观看| 亚洲av日韩在线播放| 欧美日韩视频高清一区二区三区二| 国产爽快片一区二区三区| 身体一侧抽搐| 亚洲色图综合在线观看| 成人一区二区视频在线观看| 亚洲精品亚洲一区二区| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 永久网站在线| 成人影院久久| 国产欧美日韩一区二区三区在线 | 亚洲精品乱久久久久久| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 亚洲av国产av综合av卡| 国精品久久久久久国模美| 在线观看人妻少妇| 久久热精品热| 亚洲精品久久午夜乱码| 直男gayav资源| 亚洲av日韩在线播放| 激情五月婷婷亚洲| 亚洲av不卡在线观看| 99精国产麻豆久久婷婷| 国产亚洲午夜精品一区二区久久| 欧美3d第一页| 国产高清不卡午夜福利| 国产精品免费大片| 亚洲精品中文字幕在线视频 | 国产中年淑女户外野战色| 国产欧美亚洲国产| 亚洲精品色激情综合| 免费在线观看成人毛片| 国产精品一区二区三区四区免费观看| 国产精品国产三级国产av玫瑰| 久久精品夜色国产| 色综合色国产| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 国产69精品久久久久777片| 少妇高潮的动态图| 又粗又硬又长又爽又黄的视频| 色视频www国产| 最近中文字幕2019免费版| 一级毛片电影观看| 精品99又大又爽又粗少妇毛片| 身体一侧抽搐| 国产黄片视频在线免费观看| 一个人免费看片子| 国产精品福利在线免费观看| av国产久精品久网站免费入址| a级一级毛片免费在线观看| 少妇丰满av| 熟女人妻精品中文字幕| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 国产探花极品一区二区| 人人妻人人添人人爽欧美一区卜 | 免费观看av网站的网址| 日韩成人av中文字幕在线观看| 国产免费又黄又爽又色| 亚洲欧洲国产日韩| 免费在线观看成人毛片| 亚洲欧美精品自产自拍| 美女福利国产在线 | 国产免费福利视频在线观看| 久久6这里有精品| 人人妻人人看人人澡| 国产男人的电影天堂91| 亚洲精品乱码久久久v下载方式| 国产极品天堂在线| 五月天丁香电影| 免费av不卡在线播放| 色视频在线一区二区三区| 一级毛片aaaaaa免费看小| 国产成人91sexporn| 香蕉精品网在线| 午夜日本视频在线| 国产精品嫩草影院av在线观看| 黄色日韩在线| 黄色怎么调成土黄色| 少妇的逼水好多| xxx大片免费视频| 99久久人妻综合| 国产 一区精品| 99热6这里只有精品| 寂寞人妻少妇视频99o| 一区二区三区免费毛片| 日韩,欧美,国产一区二区三区| 狂野欧美激情性xxxx在线观看| 精品久久久噜噜| 七月丁香在线播放| 免费黄网站久久成人精品| 国产黄色免费在线视频| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99| 在线观看美女被高潮喷水网站| 另类亚洲欧美激情| 国产国拍精品亚洲av在线观看| 国产日韩欧美亚洲二区| 美女主播在线视频| 色视频www国产| 亚洲精品乱码久久久久久按摩| 日韩亚洲欧美综合| 黑丝袜美女国产一区| 天天躁日日操中文字幕| 中文资源天堂在线| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 国产亚洲最大av| 欧美老熟妇乱子伦牲交| 舔av片在线| 色5月婷婷丁香| 久久久久久久精品精品| 在线精品无人区一区二区三 | 国产午夜精品一二区理论片| 免费av中文字幕在线| 国产伦精品一区二区三区四那| 一边亲一边摸免费视频| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 久久久久久久久久久丰满| 亚洲av成人精品一二三区| 国产精品伦人一区二区| 欧美精品一区二区免费开放| 午夜激情福利司机影院| 男女边摸边吃奶| 丝袜脚勾引网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品456在线播放app| 欧美日韩在线观看h| 妹子高潮喷水视频| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| 久久久午夜欧美精品| 久久久久性生活片| 日产精品乱码卡一卡2卡三| 99热全是精品| 一区二区三区精品91| 欧美3d第一页| 亚洲天堂av无毛| 亚洲精品456在线播放app| 亚洲天堂av无毛| 纯流量卡能插随身wifi吗| 99热网站在线观看| 亚洲第一av免费看| 国国产精品蜜臀av免费| 日韩视频在线欧美| 免费久久久久久久精品成人欧美视频 | 亚洲熟女精品中文字幕| 久热久热在线精品观看| 婷婷色综合大香蕉| 99热这里只有精品一区| 日本wwww免费看| 中文在线观看免费www的网站| 色吧在线观看| 99热这里只有是精品50| 国产亚洲欧美精品永久| 国产伦精品一区二区三区四那| 欧美少妇被猛烈插入视频| 日韩在线高清观看一区二区三区| 在线观看国产h片| 国产精品99久久久久久久久| 国产黄片视频在线免费观看| 国产精品一区二区三区四区免费观看| 久久韩国三级中文字幕| 黑丝袜美女国产一区| 高清欧美精品videossex| 人妻夜夜爽99麻豆av| 国产精品99久久99久久久不卡 | 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 久久精品国产自在天天线| 国产深夜福利视频在线观看| 高清日韩中文字幕在线| 在线观看一区二区三区激情| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 91久久精品国产一区二区三区|