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

    基于正交多項(xiàng)式擬合的氣象數(shù)據(jù)處理方法*

    2013-09-11 09:14:10孫寶京
    關(guān)鍵詞:高空數(shù)據(jù)處理分段

    孫寶京

    (沈陽炮兵學(xué)院電子偵察系,遼寧沈陽 110867)

    基于正交多項(xiàng)式擬合的氣象數(shù)據(jù)處理方法*

    孫寶京

    (沈陽炮兵學(xué)院電子偵察系,遼寧沈陽 110867)

    針對(duì)高空氣象探測(cè)數(shù)據(jù)變化規(guī)律復(fù)雜、突變情況不可預(yù)測(cè)、數(shù)據(jù)量大等特點(diǎn),提出了基于離散數(shù)據(jù)正交多項(xiàng)式最小二乘的分段擬合方法.首先,從計(jì)算穩(wěn)定性角度,說明了采用離散數(shù)據(jù)的正交多項(xiàng)式最小二乘法的原因;其次,從曲線擬合的保形性角度,詳細(xì)闡述了二次分段擬合原則;最后,以具有典型特征的氣溫探測(cè)數(shù)據(jù)為例,證明了采用基于離散數(shù)據(jù)正交多項(xiàng)式最小二乘的分段擬合方法能夠獲得明顯優(yōu)于傳統(tǒng)內(nèi)插方法的擬合精度,提高高空氣象探測(cè)數(shù)據(jù)處理的精度和自動(dòng)化處理程度.

    離散數(shù)據(jù);正交多項(xiàng)式;最小二乘;分段擬合;氣象探測(cè)

    當(dāng)前氣象探測(cè)站整理高空氣象探測(cè)數(shù)據(jù)時(shí),主要采用的方法是將整個(gè)探測(cè)空間分層,計(jì)算規(guī)定高度的中間數(shù)值,作為計(jì)算層數(shù)據(jù),其他數(shù)據(jù)則采用依托計(jì)算層數(shù)據(jù)內(nèi)插的方法計(jì)算.上述計(jì)算方法使用高度上的數(shù)據(jù)均由計(jì)算層數(shù)據(jù)插值得到,存在一定的截?cái)嗾`差和舍入誤差;另外在較少的計(jì)算層數(shù)據(jù)中,必然會(huì)丟棄絕大多數(shù)探測(cè)數(shù)據(jù),造成探測(cè)資源的極大浪費(fèi).隨著氣象探測(cè)及作業(yè)裝備自動(dòng)化水平和數(shù)據(jù)處理能力的不斷提高,對(duì)高空氣象探測(cè)數(shù)據(jù)處理精度和時(shí)效性提出了更高的要求,這也使得從根本上改變氣象參數(shù)的計(jì)算方法成為可能[1].

    由于氣象要素隨時(shí)間、空間在不斷地變化,而且探測(cè)儀器本身也受到許多復(fù)雜因素的影響,因此各種高空氣象探測(cè)數(shù)據(jù)普遍具有數(shù)據(jù)量大、變化規(guī)律復(fù)雜、突變情況較多且不可預(yù)測(cè)等特點(diǎn)[2].以彈道氣象諸元中最重要的高空實(shí)時(shí)氣溫/高度曲線為例,當(dāng)探測(cè)高度達(dá)到2萬m時(shí),采用電子探空儀探測(cè)的數(shù)據(jù)超過2 000組,雖然高空實(shí)時(shí)氣溫總體呈下降趨勢(shì),但是在局部經(jīng)常出現(xiàn)不可預(yù)測(cè)的無規(guī)則變化.

    針對(duì)高空氣象探測(cè)數(shù)據(jù)特點(diǎn),筆者基于離散數(shù)據(jù)正交多項(xiàng)式最小二乘擬合思想,提出根據(jù)數(shù)據(jù)變化率的變化情況,分段擬合各種高空氣象探測(cè)數(shù)據(jù)曲線的方法.借助編程軟件的自動(dòng)處理,經(jīng)過多次試驗(yàn)證明,此方法可直接得出氣象通報(bào)所需要的氣(虛)溫-高度、氣壓-高度、風(fēng)速-高度和風(fēng)向-高度函數(shù)方程(高度為自變量),并且對(duì)高空氣象探測(cè)數(shù)據(jù)處理的精確性、實(shí)效性和計(jì)算簡(jiǎn)捷性都有了很大的提高.因?yàn)楦呖諏?shí)時(shí)氣溫/高度曲線在彈道氣象諸元中最為重要,所以筆者將以此作為說明實(shí)例.

    1 離散數(shù)據(jù)分段最小二乘擬合方法

    用插值方法所找出的近似曲線,通過所有己知點(diǎn),但是由觀測(cè)或?qū)嶒?yàn)所獲得的數(shù)據(jù),不可避免地含有誤差,這就使插值所找出的近似曲線保存所有觀測(cè)誤差,導(dǎo)致所得結(jié)果可能偏離了實(shí)際情況,出現(xiàn)曲線擬合的過度現(xiàn)象.而根據(jù)實(shí)際標(biāo)準(zhǔn)盡可能地通過節(jié)點(diǎn)的函數(shù)值近旁,以部分地抵消原始數(shù)據(jù)所包含的觀測(cè)誤差,則更具有實(shí)用價(jià)值,最小二乘擬合就是其中的一個(gè)重要方法[3].

    1.1 離散數(shù)據(jù)的正交多項(xiàng)式最小二乘擬合

    在線性空間由Φ=span(φ0,φ1,...,φn)求解最小二乘曲線擬合問題,經(jīng)常會(huì)出現(xiàn)法方程組是病態(tài)方程組的情況,特別是當(dāng)樣本組數(shù)n偏大時(shí)更是如此,同時(shí)在求解法方程組時(shí)系數(shù)矩陣或右端項(xiàng)微小的擾動(dòng)都可能導(dǎo)致解函數(shù)有很大的誤差.為了確保計(jì)算過程的穩(wěn)定性和精確性,筆者利用基于高空氣象探測(cè)數(shù)據(jù)的正交多項(xiàng)式作為擬合函數(shù)系,使法方程組的系數(shù)矩陣變?yōu)閷?duì)角矩陣,從而解得正交多項(xiàng)式的系數(shù)[4-5].

    根據(jù)點(diǎn)集{x1,x2,...,xm},構(gòu)造出相應(yīng)的正交多項(xiàng)式系{p0(x),p1(x),...,pn(x)}.

    由正交多項(xiàng)式三項(xiàng)遞推關(guān)系,最高次項(xiàng)系數(shù)為1的正交多項(xiàng)式系{pk(x)}(k=0,1,2,...,n)有如下遞推關(guān)系:

    其中pk(x)為最高項(xiàng)系數(shù)為1的k次多項(xiàng)式.由{pk(x)}正交性可知,

    1.2 二次分段擬合原則

    因一次完整綜合氣象探測(cè)的數(shù)據(jù)量大,若對(duì)整組數(shù)據(jù)進(jìn)行一次性擬合,則擬合精度和計(jì)算效果都比較差.同時(shí),為克服傳統(tǒng)分段曲線擬合方法中對(duì)數(shù)據(jù)點(diǎn)分段時(shí)經(jīng)驗(yàn)成分較多的缺點(diǎn),筆者根據(jù)數(shù)據(jù)曲線斜率的變化情況對(duì)曲線進(jìn)行分段擬合,并且加入了使相鄰曲線連續(xù),即曲線邊界點(diǎn)必須連續(xù)的約束條件[6-7].

    考慮到高度是單調(diào)遞增,為了適應(yīng)高空氣象數(shù)據(jù)曲線在局部隨高度變化的無規(guī)則性,采取二次分段原則.

    一次分段:在氣溫和高度曲線中,設(shè)Δti=(ti+1-ti)/(hi+1-h(huán)i),其中ti和hi分別為探測(cè)數(shù)據(jù)中第i個(gè)探測(cè)點(diǎn)的氣溫和高度,ti+1和hi+1分別為探測(cè)數(shù)據(jù)中第i+1個(gè)探測(cè)點(diǎn)的氣溫和高度.若Δti-1為負(fù),Δti,Δti+1,Δti+2,...,Δti+a(a≥1)為正,Δti+a+1,...,Δti+a+b(b≥1)為負(fù),Δti+a+b+1為正,則{ti,...,ti+a+b}為第i個(gè)分段,第i+1個(gè)分段從點(diǎn)ti+a+b開始.

    二次分段:在氣溫和高度曲線中,第i個(gè)分段已經(jīng)由一次分段確定,設(shè)δi+c=Δti+c-Δti+c-1(1≤c≤a+b-1).若δi+1,δi+2,...,δi+d(d≥1)均大于標(biāo)準(zhǔn)值δ0(為負(fù)數(shù)),δi+d+1≤δ0,且δi+d+2,...,δi+d+e(e≥1)小于或等于(-δ0),δi+d+e+1>(-δ0),則ti+d+e+1為分段點(diǎn),即{ti,...,ti+d+e+1}為一個(gè)分段,下一個(gè)分段以ti+d+e+1為起點(diǎn),從δi+d+e+2開始判斷;若δi+1≤δ0,且δi+2,...,δi+e(e≥1)小于或等于(-δ0),δi+e+1>(-δ0),則ti+e+1為分段點(diǎn),即{ti,...,ti+e+1}為一個(gè)分段,下一個(gè)分段以ti+e+1為起點(diǎn),從δi+e+2開始判斷[8-9].

    為使在單個(gè)分段內(nèi)擬合曲線在誤差允許范圍內(nèi)最大程度地保持?jǐn)?shù)據(jù)曲線的形狀,并且減少數(shù)據(jù)處理程序的計(jì)算量,實(shí)際計(jì)算過程對(duì)原始探測(cè)數(shù)據(jù)中心化計(jì)算后,采用如下做法進(jìn)行擬合:

    若a+b≥5,則用四階正交多項(xiàng)式為擬合函數(shù)系,即當(dāng)一個(gè)分段包含4個(gè)以上的點(diǎn)時(shí)用四階正交多項(xiàng)式擬合,

    若a+b=4,則用三階正交多項(xiàng)式為擬合函數(shù)系,即當(dāng)一個(gè)分段包含4個(gè)點(diǎn)時(shí)用三階正交多項(xiàng)式擬合,

    若a+b=3,則用二階正交多項(xiàng)式為擬合函數(shù)系,即當(dāng)一個(gè)分段包含3個(gè)點(diǎn)時(shí)用二階正交多項(xiàng)式擬合,

    同時(shí),使參與第i個(gè)分段擬合的第i個(gè)分段中的最后一個(gè)數(shù)據(jù)點(diǎn),作為第i+1個(gè)分段的第1個(gè)數(shù)據(jù)點(diǎn)參與第i+1個(gè)分段的擬合,這樣第i個(gè)分段的擬合函數(shù)和第i+1個(gè)分段的擬合函數(shù)在分界點(diǎn)上的值相等,從而實(shí)現(xiàn)所有分段擬合函數(shù)在整體上保持連續(xù)[10],即

    2 計(jì)算實(shí)例

    以編號(hào)為087233的GZZ8型電子探空儀探測(cè)數(shù)據(jù)為例.

    圖1是該組探測(cè)數(shù)據(jù)曲線中一部分,溫度t∈[-61.31,-62.27],高度H∈[15 125.89,1 5 494.56].根據(jù)文中的分段規(guī)則,數(shù)據(jù)處理系統(tǒng)自動(dòng)將此段數(shù)據(jù)曲線分為5段進(jìn)行擬合.

    圖1 溫度和高度探測(cè)數(shù)據(jù)分布

    第1個(gè)分段為前13個(gè)點(diǎn),h∈[15 125.89,15 230.23],擬合函數(shù)為

    平方誤差為0.08.

    第2個(gè)分段為后續(xù)7個(gè)點(diǎn),h∈[15 230.23,15 289.47],擬合函數(shù)為

    平方誤差為0.03.

    第3個(gè)分段為后續(xù)10個(gè)點(diǎn),h∈[15 289.47,15 397.85],擬合函數(shù)為

    平方誤差為0.02.

    第4個(gè)分段為后續(xù)5個(gè)點(diǎn),h∈[15 397.85,15 439.32],擬合函數(shù)為

    平方誤差為0.01.

    第5個(gè)分段為最后7個(gè)點(diǎn),h∈[15 439.32,15 494.56],擬合函數(shù)為

    平方誤差為0.03.

    根據(jù)傳統(tǒng)高空氣象探測(cè)數(shù)據(jù)處理方法,14 000~16 000m高度范圍內(nèi)的相關(guān)氣象數(shù)據(jù)從標(biāo)準(zhǔn)高度層14 000m和16 000m處的實(shí)際探測(cè)數(shù)據(jù)進(jìn)行線性內(nèi)插得到.在此例中,14 000m和16 000m處的氣溫分別為-56.16℃和-62.37℃,從而內(nèi)插得到高度為15 142.00,15 276.4,15 378.67,15 417.06,1 5 479.1m處的氣溫.各高度對(duì)應(yīng)的實(shí)際探測(cè)值、內(nèi)插結(jié)果和擬合結(jié)果見表1.由表1可見,采用基于離散數(shù)據(jù)正交多項(xiàng)式最小二乘的分段擬合法,可以獲得比傳統(tǒng)內(nèi)插法更好的擬合精度.

    表1 不同高度實(shí)際探測(cè)值、內(nèi)插值、擬合值對(duì)應(yīng)關(guān)系

    3 結(jié)語

    針對(duì)高空氣象探測(cè)數(shù)據(jù)量大、變化規(guī)律復(fù)雜、突變情況較多且不可預(yù)測(cè)等特點(diǎn),提出了基于離散數(shù)據(jù)正交多項(xiàng)式最小二乘的分段擬合法.經(jīng)過實(shí)踐證明,該方法能夠根據(jù)算法自動(dòng)分段,擬合出符合精度要求的曲線,在處理程序自動(dòng)計(jì)算過程中確保計(jì)算的穩(wěn)定性和高效性,避免傳統(tǒng)高空氣象探測(cè)數(shù)據(jù)處理方法中不斷累計(jì)的舍入誤差對(duì)結(jié)果的影響,提高高空氣象探測(cè)數(shù)據(jù)處理結(jié)果的精度和自動(dòng)化處理程度.

    [1] 孫寶京.炮兵防空兵氣象信息系統(tǒng)研究[D].沈陽:東北大學(xué),2004.

    [2] 曲延祿.外彈道氣象學(xué)概論[M].北京:氣象出版社,1987:59-176.

    [3] FRANK R.Scattered Data Interpolation:Tests of Some Methods[J].Mathematical Computation,1982,38:181-200.

    [4] MICCHELLI C.Interpolation of Scattered Data:Distance Matrix and Conditionally Positive Definite Functions[J].Constructive Approximation,1986,2:11-22.

    [5] 藺小林,蔣耀林.現(xiàn)代數(shù)值分析[M].北京:國防工業(yè)出版社,2004:212-255.

    [6] WU Zong-min,SCHABACK R.Shape Preserving Properties and Convergence of Univariate Multi-Quadric Quasti-Interpolation[J].ACTA Mathematice Application Sinica,1994,10:441-446.

    [7] GREGORY J A.Shape Preserving Spline Interpolation[J].Computer-Aided Design,1986,18(1):53-57.

    [8] KAUFMANN E,KLASS R.Smoothing Surfaces Using Reflection Lines for Families of Splines[J].Computer-Aided Design,1988,20(6):312-316.

    [9] HOSCHEK J.Smoothing Curves and Surfaces[J].Computer-Aided Geometric Design,1985,2(1-3):97-105.

    [10] JORG PRTERS.Local Smooth Surface Interpolation:A Classification[J].Computer-Aided Geometric Design,1 990(7):191-195.

    (責(zé)任編輯 向陽潔)

    Piecewise Curve Fitting Method of Processing Meteorological Data Based on the Orthogonal Polynomial

    SUN Bao-jing

    (Department of Electronic Reconnaissance,Shenyang Artillery Academy,Shenyang 110867,China)

    In accordance with the characterics of the large number and complex and unpredicted trend on change of the meteorological detecting data,the piecewise curve fitting method based on the least squares of orthogonal polynomial for scattered data is presented.Firstly,the reason of using the least squares of orthogonal polynomial for scattered data is introduced from the respect of the stability of data processing.Secondly,the principle of twice piecewise curve fitting is set out in detail from the respect of the shape preserving feature of curve fitting.Finally,as the temperature detecting data for example,it is proved clearly that the twice piecewise curve fitting method can make much better fitting precision than the traditional interpolation method,and improve the level of precision and automation of the meteorological detecting data processing.

    scattered data;orthogonal polynomial;least squares;piecewise curve fitting;meteorological detecting

    P412.2;TP274

    A

    10.3969/j.issn.1007-2985.2013.05.012

    1007-2985(2013)05-0049-04

    2013-01-24

    孫寶京(1968-),男,山東五蓮人,沈陽炮兵學(xué)院教授,主要從事彈道氣象學(xué)研究.

    猜你喜歡
    高空數(shù)據(jù)處理分段
    認(rèn)知診斷缺失數(shù)據(jù)處理方法的比較:零替換、多重插補(bǔ)與極大似然估計(jì)法*
    高空走繩
    一類連續(xù)和不連續(xù)分段線性系統(tǒng)的周期解研究
    ILWT-EEMD數(shù)據(jù)處理的ELM滾動(dòng)軸承故障診斷
    高空纜車
    分段計(jì)算時(shí)間
    不要高空拋物!
    高空莫拋物
    3米2分段大力士“大”在哪兒?
    太空探索(2016年9期)2016-07-12 10:00:04
    基于希爾伯特- 黃變換的去噪法在外測(cè)數(shù)據(jù)處理中的應(yīng)用
    av播播在线观看一区| 人人澡人人妻人| 婷婷色综合www| 欧美成人精品欧美一级黄| h视频一区二区三区| 国产欧美亚洲国产| 免费在线观看黄色视频的| 久久久久人妻精品一区果冻| 在线 av 中文字幕| 国产色婷婷99| 国产福利在线免费观看视频| 一边亲一边摸免费视频| 国产精品久久久久成人av| 美女视频免费永久观看网站| 免费在线观看完整版高清| 亚洲第一区二区三区不卡| 国产麻豆69| 久久99热这里只频精品6学生| 51国产日韩欧美| 成人综合一区亚洲| 一级毛片黄色毛片免费观看视频| 亚洲美女视频黄频| 久久人人爽av亚洲精品天堂| 国产精品熟女久久久久浪| 成年av动漫网址| 精品少妇黑人巨大在线播放| 免费观看a级毛片全部| 成人影院久久| 性色av一级| 国产成人午夜福利电影在线观看| 22中文网久久字幕| 九色成人免费人妻av| 我要看黄色一级片免费的| 亚洲国产精品国产精品| 中文字幕另类日韩欧美亚洲嫩草| 一本色道久久久久久精品综合| 大片免费播放器 马上看| 中文字幕亚洲精品专区| 免费高清在线观看日韩| www.色视频.com| 一级毛片 在线播放| 精品一区二区免费观看| 国产熟女午夜一区二区三区| av在线观看视频网站免费| 五月伊人婷婷丁香| 啦啦啦中文免费视频观看日本| 成年人午夜在线观看视频| 日本91视频免费播放| 五月开心婷婷网| 日韩中字成人| 国产综合精华液| 天堂俺去俺来也www色官网| 午夜视频国产福利| 亚洲精品日韩在线中文字幕| 久久久国产一区二区| 亚洲综合精品二区| 国产色婷婷99| 超色免费av| 中文字幕最新亚洲高清| av网站免费在线观看视频| 成年人免费黄色播放视频| videosex国产| 亚洲欧美一区二区三区黑人 | 久久久久精品性色| 妹子高潮喷水视频| 日韩电影二区| 九草在线视频观看| 久久久久久久久久人人人人人人| 看免费成人av毛片| 国产一区二区三区综合在线观看 | 一级毛片我不卡| 一边摸一边做爽爽视频免费| 三级国产精品片| 免费黄网站久久成人精品| av黄色大香蕉| 18禁观看日本| 男女无遮挡免费网站观看| 国产男女超爽视频在线观看| 交换朋友夫妻互换小说| 日韩熟女老妇一区二区性免费视频| 亚洲精华国产精华液的使用体验| 国产男女超爽视频在线观看| 肉色欧美久久久久久久蜜桃| 22中文网久久字幕| 亚洲精品乱久久久久久| 午夜免费鲁丝| 久久久久视频综合| 成人影院久久| 美女福利国产在线| 人妻少妇偷人精品九色| 99久久人妻综合| 亚洲性久久影院| 高清黄色对白视频在线免费看| 男人操女人黄网站| 各种免费的搞黄视频| 成人18禁高潮啪啪吃奶动态图| 欧美人与性动交α欧美软件 | 亚洲欧美精品自产自拍| 亚洲av电影在线进入| 91在线精品国自产拍蜜月| 在线天堂中文资源库| 国产成人精品在线电影| 2022亚洲国产成人精品| 桃花免费在线播放| 欧美日韩视频高清一区二区三区二| 国产精品熟女久久久久浪| 纯流量卡能插随身wifi吗| 秋霞伦理黄片| 亚洲人成77777在线视频| 成人手机av| 大片电影免费在线观看免费| 91精品伊人久久大香线蕉| 一级黄片播放器| 中国美白少妇内射xxxbb| 如日韩欧美国产精品一区二区三区| 久久精品久久精品一区二区三区| 国产乱人偷精品视频| 久久人妻熟女aⅴ| 久久人人爽人人爽人人片va| 三级国产精品片| 国产视频首页在线观看| 在现免费观看毛片| 99re6热这里在线精品视频| 日韩不卡一区二区三区视频在线| av电影中文网址| 国产又爽黄色视频| 久久久久久久久久久久大奶| 观看美女的网站| 高清毛片免费看| 插逼视频在线观看| 大码成人一级视频| 欧美xxxx性猛交bbbb| 中国美白少妇内射xxxbb| 日本-黄色视频高清免费观看| 免费播放大片免费观看视频在线观看| 国产免费一区二区三区四区乱码| 在线亚洲精品国产二区图片欧美| 国产精品.久久久| 色吧在线观看| 一二三四在线观看免费中文在 | 亚洲熟女精品中文字幕| 又黄又粗又硬又大视频| 午夜福利,免费看| 亚洲欧洲国产日韩| 亚洲精品第二区| 99re6热这里在线精品视频| 日本wwww免费看| 高清欧美精品videossex| 国产淫语在线视频| av福利片在线| 欧美精品亚洲一区二区| av有码第一页| 精品久久国产蜜桃| 99久久人妻综合| 男女下面插进去视频免费观看 | 国产免费现黄频在线看| √禁漫天堂资源中文www| 建设人人有责人人尽责人人享有的| 看免费成人av毛片| 免费av中文字幕在线| 国产亚洲一区二区精品| 成年美女黄网站色视频大全免费| 中国国产av一级| 久久久久久久大尺度免费视频| 国产精品一区二区在线观看99| 国产一区亚洲一区在线观看| 国产一级毛片在线| 欧美精品人与动牲交sv欧美| 美女内射精品一级片tv| 国产成人午夜福利电影在线观看| 国产成人免费观看mmmm| 侵犯人妻中文字幕一二三四区| 热99久久久久精品小说推荐| 日日撸夜夜添| 妹子高潮喷水视频| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 高清毛片免费看| 欧美精品一区二区大全| 高清不卡的av网站| 午夜免费男女啪啪视频观看| 亚洲av男天堂| 丰满饥渴人妻一区二区三| 国产成人精品在线电影| 久久综合国产亚洲精品| 国产色婷婷99| 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 免费看光身美女| 十分钟在线观看高清视频www| 亚洲欧美日韩另类电影网站| 成人亚洲精品一区在线观看| 色哟哟·www| 香蕉精品网在线| 国产精品久久久久久久电影| 日本爱情动作片www.在线观看| 免费黄频网站在线观看国产| 免费观看a级毛片全部| av电影中文网址| 亚洲欧美成人精品一区二区| 桃花免费在线播放| 欧美3d第一页| 国产女主播在线喷水免费视频网站| 黄片播放在线免费| 精品视频人人做人人爽| 男人操女人黄网站| 国产精品久久久久成人av| 午夜免费观看性视频| 五月伊人婷婷丁香| 一边摸一边做爽爽视频免费| 亚洲av日韩在线播放| 国产综合精华液| 国产一区有黄有色的免费视频| 日韩av不卡免费在线播放| 高清不卡的av网站| 久久亚洲国产成人精品v| 一级毛片 在线播放| 丰满饥渴人妻一区二区三| 狠狠婷婷综合久久久久久88av| 黑人欧美特级aaaaaa片| 精品国产一区二区三区四区第35| 精品国产一区二区久久| 久久久精品区二区三区| av网站免费在线观看视频| 91精品国产国语对白视频| 日韩电影二区| 亚洲欧美成人综合另类久久久| 国产精品一区www在线观看| 欧美激情极品国产一区二区三区 | 99re6热这里在线精品视频| 亚洲av福利一区| 国产成人欧美| 免费人妻精品一区二区三区视频| 人妻少妇偷人精品九色| 久久人妻熟女aⅴ| 久久99热这里只频精品6学生| 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品电影小说| 国产成人精品一,二区| 精品人妻在线不人妻| 女人久久www免费人成看片| 免费观看a级毛片全部| 曰老女人黄片| 美女脱内裤让男人舔精品视频| 免费大片18禁| 成年女人在线观看亚洲视频| 女人精品久久久久毛片| 天天躁夜夜躁狠狠躁躁| 亚洲熟女精品中文字幕| 黄网站色视频无遮挡免费观看| 插逼视频在线观看| 精品少妇久久久久久888优播| 免费人成在线观看视频色| 亚洲熟女精品中文字幕| 欧美日韩av久久| 女人久久www免费人成看片| 欧美精品av麻豆av| 男人添女人高潮全过程视频| 日本wwww免费看| 国产乱来视频区| 伦理电影大哥的女人| 激情五月婷婷亚洲| 欧美97在线视频| 久久这里只有精品19| 国产成人aa在线观看| 久久久久久久久久人人人人人人| 91精品国产国语对白视频| 考比视频在线观看| 爱豆传媒免费全集在线观看| 99久久人妻综合| av在线观看视频网站免费| 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 久久ye,这里只有精品| 欧美国产精品va在线观看不卡| 高清视频免费观看一区二区| 亚洲综合色网址| 女的被弄到高潮叫床怎么办| 日本av免费视频播放| 久久久久精品性色| 午夜福利在线观看免费完整高清在| 亚洲 欧美一区二区三区| 成年女人在线观看亚洲视频| 亚洲av福利一区| 自线自在国产av| 两性夫妻黄色片 | 女性生殖器流出的白浆| 日韩成人av中文字幕在线观看| 亚洲欧美色中文字幕在线| 久久久久精品久久久久真实原创| 最近2019中文字幕mv第一页| 午夜福利网站1000一区二区三区| 国产精品久久久久成人av| 岛国毛片在线播放| 91在线精品国自产拍蜜月| 久久毛片免费看一区二区三区| 一本色道久久久久久精品综合| 国产欧美日韩一区二区三区在线| 日本av免费视频播放| 国产成人aa在线观看| 亚洲欧美精品自产自拍| 亚洲四区av| 女人久久www免费人成看片| 青青草视频在线视频观看| 午夜免费观看性视频| 国产精品人妻久久久影院| 精品99又大又爽又粗少妇毛片| 在线观看免费日韩欧美大片| 久久久久国产精品人妻一区二区| 男的添女的下面高潮视频| 一本久久精品| 精品国产露脸久久av麻豆| 亚洲色图综合在线观看| 国产一区二区三区综合在线观看 | 人体艺术视频欧美日本| 男男h啪啪无遮挡| 午夜福利,免费看| 精品一区在线观看国产| 国产av一区二区精品久久| 少妇的丰满在线观看| 狠狠精品人妻久久久久久综合| 中文欧美无线码| 9热在线视频观看99| 亚洲欧美精品自产自拍| 在线 av 中文字幕| 久久99一区二区三区| 一本大道久久a久久精品| 亚洲中文av在线| 人人澡人人妻人| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 日本欧美视频一区| 亚洲综合色网址| 精品亚洲成国产av| 日韩中文字幕视频在线看片| 久久国产精品男人的天堂亚洲 | 亚洲精华国产精华液的使用体验| 精品一区在线观看国产| 宅男免费午夜| 中文字幕最新亚洲高清| 国产毛片在线视频| 国产精品女同一区二区软件| 国产不卡av网站在线观看| 欧美国产精品一级二级三级| 两个人免费观看高清视频| 国产男女内射视频| 男人添女人高潮全过程视频| 国产亚洲午夜精品一区二区久久| 美女中出高潮动态图| 欧美精品国产亚洲| 国产淫语在线视频| 激情五月婷婷亚洲| 免费在线观看完整版高清| 一级毛片电影观看| 免费看不卡的av| 97在线视频观看| 成年动漫av网址| 桃花免费在线播放| 日本欧美视频一区| 九色成人免费人妻av| 久久狼人影院| 一级a做视频免费观看| 精品国产一区二区三区四区第35| 多毛熟女@视频| 久久婷婷青草| 国产高清国产精品国产三级| 久久亚洲国产成人精品v| 中国国产av一级| 亚洲精华国产精华液的使用体验| 天堂8中文在线网| 亚洲,欧美精品.| 国产精品一区二区在线不卡| 午夜福利,免费看| 成人漫画全彩无遮挡| 亚洲精华国产精华液的使用体验| 久久热在线av| 欧美激情极品国产一区二区三区 | 成人毛片60女人毛片免费| av不卡在线播放| 国产熟女欧美一区二区| 婷婷色综合www| 午夜福利在线观看免费完整高清在| 一本大道久久a久久精品| 欧美xxⅹ黑人| 欧美日韩av久久| 亚洲精品乱码久久久久久按摩| 熟女人妻精品中文字幕| 91国产中文字幕| 妹子高潮喷水视频| 亚洲欧洲国产日韩| 精品亚洲成国产av| 国产男人的电影天堂91| 午夜福利乱码中文字幕| 成人影院久久| 日韩,欧美,国产一区二区三区| 视频中文字幕在线观看| 性高湖久久久久久久久免费观看| av国产久精品久网站免费入址| 国产高清不卡午夜福利| 2021少妇久久久久久久久久久| 在线天堂最新版资源| 欧美国产精品va在线观看不卡| 边亲边吃奶的免费视频| 亚洲五月色婷婷综合| 日韩av在线免费看完整版不卡| 亚洲五月色婷婷综合| 欧美亚洲 丝袜 人妻 在线| 99热这里只有是精品在线观看| 亚洲精品久久久久久婷婷小说| 日韩一区二区三区影片| 街头女战士在线观看网站| 欧美日韩亚洲高清精品| 国产一区二区在线观看av| 精品人妻一区二区三区麻豆| 精品午夜福利在线看| 青春草视频在线免费观看| 免费久久久久久久精品成人欧美视频 | 国产欧美日韩一区二区三区在线| 26uuu在线亚洲综合色| 丝瓜视频免费看黄片| 美女国产视频在线观看| 国产色婷婷99| 久久久久视频综合| 国产无遮挡羞羞视频在线观看| 久久久国产精品麻豆| 亚洲综合色网址| a级毛色黄片| 日本黄色日本黄色录像| 国产成人av激情在线播放| 成年人免费黄色播放视频| 日日摸夜夜添夜夜爱| 国产精品一区二区在线观看99| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 97精品久久久久久久久久精品| 日产精品乱码卡一卡2卡三| 香蕉精品网在线| 亚洲欧洲国产日韩| 夜夜爽夜夜爽视频| 国产亚洲最大av| 欧美+日韩+精品| av卡一久久| 99视频精品全部免费 在线| 欧美国产精品va在线观看不卡| 男女啪啪激烈高潮av片| 国产精品久久久久久精品古装| 草草在线视频免费看| 午夜福利视频在线观看免费| 一级片'在线观看视频| 少妇被粗大的猛进出69影院 | 午夜免费观看性视频| 搡老乐熟女国产| 免费看光身美女| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 国产精品一区www在线观看| 成人手机av| 国产深夜福利视频在线观看| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 国产成人av激情在线播放| av有码第一页| 国内精品宾馆在线| 日韩av在线免费看完整版不卡| 国产69精品久久久久777片| 26uuu在线亚洲综合色| 国产亚洲av片在线观看秒播厂| 多毛熟女@视频| 免费看光身美女| 久久午夜福利片| 成人毛片60女人毛片免费| 一边摸一边做爽爽视频免费| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 夫妻午夜视频| 精品酒店卫生间| 亚洲国产最新在线播放| 伊人亚洲综合成人网| 丰满饥渴人妻一区二区三| 国产国语露脸激情在线看| 丝袜喷水一区| 男女高潮啪啪啪动态图| 日韩熟女老妇一区二区性免费视频| 国内精品宾馆在线| 免费高清在线观看视频在线观看| 男人舔女人的私密视频| 亚洲成国产人片在线观看| 18+在线观看网站| 日本爱情动作片www.在线观看| 香蕉丝袜av| 久久人人爽人人片av| 日本午夜av视频| 国产色爽女视频免费观看| 欧美另类一区| 久久久国产精品麻豆| 亚洲国产精品成人久久小说| 亚洲精品,欧美精品| 自线自在国产av| 自拍欧美九色日韩亚洲蝌蚪91| kizo精华| 99久久综合免费| 母亲3免费完整高清在线观看 | 老女人水多毛片| 免费久久久久久久精品成人欧美视频 | 中文欧美无线码| 亚洲av在线观看美女高潮| 高清av免费在线| 亚洲成av片中文字幕在线观看 | 97精品久久久久久久久久精品| 欧美人与性动交α欧美精品济南到 | 晚上一个人看的免费电影| 国产精品国产三级国产av玫瑰| a级片在线免费高清观看视频| 久久韩国三级中文字幕| 在线看a的网站| 另类精品久久| 亚洲一级一片aⅴ在线观看| 欧美激情 高清一区二区三区| av国产久精品久网站免费入址| av网站免费在线观看视频| 在线亚洲精品国产二区图片欧美| 在线观看三级黄色| 久久鲁丝午夜福利片| 欧美 亚洲 国产 日韩一| av在线观看视频网站免费| 精品人妻一区二区三区麻豆| 黑人猛操日本美女一级片| 国产免费现黄频在线看| 少妇 在线观看| 成人国产麻豆网| 成人午夜精彩视频在线观看| 亚洲 欧美一区二区三区| 国产福利在线免费观看视频| 中文字幕亚洲精品专区| 777米奇影视久久| 亚洲,欧美,日韩| 一本大道久久a久久精品| 亚洲美女黄色视频免费看| 久久久久久久亚洲中文字幕| 亚洲情色 制服丝袜| 国产精品久久久久久久电影| 国产精品一二三区在线看| 亚洲一级一片aⅴ在线观看| 亚洲,欧美,日韩| 我要看黄色一级片免费的| 国产精品一国产av| 国产日韩欧美在线精品| 精品少妇久久久久久888优播| 最近最新中文字幕大全免费视频 | 满18在线观看网站| 超碰97精品在线观看| 日本wwww免费看| 成年动漫av网址| 亚洲精品久久午夜乱码| 纵有疾风起免费观看全集完整版| 成年av动漫网址| 日韩av免费高清视频| 国产白丝娇喘喷水9色精品| 国产精品99久久99久久久不卡 | 国产黄频视频在线观看| 精品人妻一区二区三区麻豆| 中国国产av一级| 午夜福利在线观看免费完整高清在| 妹子高潮喷水视频| 精品午夜福利在线看| 国产成人一区二区在线| 亚洲四区av| 日韩伦理黄色片| 亚洲人成网站在线观看播放| 精品少妇内射三级| av在线观看视频网站免费| 午夜影院在线不卡| 啦啦啦中文免费视频观看日本| 22中文网久久字幕| 国产日韩欧美亚洲二区| 免费大片18禁| 国产女主播在线喷水免费视频网站| 内地一区二区视频在线| 国产视频首页在线观看| 免费观看在线日韩| 久久久久久久国产电影| www.熟女人妻精品国产 | 亚洲精品第二区| av片东京热男人的天堂| 精品少妇黑人巨大在线播放| 国产 精品1| 久久久久精品久久久久真实原创| kizo精华| 老司机亚洲免费影院| 哪个播放器可以免费观看大片| 日韩欧美一区视频在线观看| 好男人视频免费观看在线| 亚洲精品,欧美精品| 日韩视频在线欧美| 亚洲婷婷狠狠爱综合网| 久久99蜜桃精品久久| 成年动漫av网址| 午夜老司机福利剧场| 国产一区二区三区av在线| 久久这里有精品视频免费| 一本久久精品| 69精品国产乱码久久久| 两个人免费观看高清视频| 亚洲综合精品二区| 日本色播在线视频| 赤兔流量卡办理| 一本久久精品| 美女主播在线视频| 国产不卡av网站在线观看| 在现免费观看毛片| 亚洲欧洲国产日韩| 亚洲,欧美,日韩| 狂野欧美激情性bbbbbb| 日韩大片免费观看网站|