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

    天體測(cè)量法探測(cè)系外行星?

    2016-06-27 08:14:12許偉維廖新浩周永宏許雪晴
    天文學(xué)報(bào) 2016年4期
    關(guān)鍵詞:外行星測(cè)量法根數(shù)

    許偉維 廖新浩 周永宏 許雪晴

    (1中國(guó)科學(xué)院上海天文臺(tái) 上海 200030) (2上??萍即髮W(xué)物質(zhì)學(xué)院 上海 201210) (3中國(guó)科學(xué)院行星科學(xué)重點(diǎn)實(shí)驗(yàn)室 上海 200030) (4中國(guó)科學(xué)院大學(xué) 北京 100049)

    天體測(cè)量法探測(cè)系外行星?

    許偉維1,2,3,4?廖新浩1,3?周永宏1,3許雪晴1,3

    (1中國(guó)科學(xué)院上海天文臺(tái) 上海 200030) (2上海科技大學(xué)物質(zhì)學(xué)院 上海 201210) (3中國(guó)科學(xué)院行星科學(xué)重點(diǎn)實(shí)驗(yàn)室 上海 200030) (4中國(guó)科學(xué)院大學(xué) 北京 100049)

    在目前已發(fā)現(xiàn)的系外行星中,絕大多數(shù)是由視向速度法和凌星法探測(cè)得到的,天體測(cè)量法僅發(fā)現(xiàn)了1顆.gaia衛(wèi)星數(shù)據(jù)即將發(fā)布,天體測(cè)量法將逐步成為系外行星探測(cè)的重要方法之一.基于天體測(cè)量法給出的恒星位置參數(shù)序列,討論了在求解行星質(zhì)量和軌道參數(shù)時(shí)涉及的動(dòng)力學(xué)條件方程計(jì)算問題,給出了具體微分改正公式,同時(shí)也進(jìn)行了必要的仿真模擬計(jì)算.建立的方法可以較容易地推廣到多行星系統(tǒng).

    系外行星,天體測(cè)量學(xué),天體力學(xué):軌道計(jì)算與定軌,方法:數(shù)據(jù)分析

    1 引言

    系外行星探測(cè)和研究是當(dāng)今國(guó)際天文學(xué)研究的熱點(diǎn).1995年,M ayor和Queloz采用視向速度法在飛馬座51附近發(fā)現(xiàn)了第1顆圍繞著主序恒星運(yùn)動(dòng)的木星質(zhì)量量級(jí)的系外行星飛馬座51 b[1].近年來,隨著天文探測(cè)技術(shù)的進(jìn)步,特別是Kep ler空間望遠(yuǎn)鏡的成功發(fā)射和觀測(cè),新的系外行星不斷被發(fā)現(xiàn),其數(shù)量明顯增多.根據(jù)系外行星網(wǎng)站(exop lanets.eu)統(tǒng)計(jì),截止2015年12月1日,已確認(rèn)發(fā)現(xiàn)2004顆系外行星,隸屬于1269個(gè)行星系統(tǒng),包括498個(gè)多行星系統(tǒng);此外,400多顆系外行星候選體需進(jìn)一步認(rèn)證.在這些已發(fā)現(xiàn)的系外行星中,絕大多數(shù)為氣態(tài)巨行星和熱木星(約占70%),部分為超級(jí)地球(20%)和地球質(zhì)量大小的行星(2%),但尚未確認(rèn)存在生命的類地行星系統(tǒng).

    系外行星的探測(cè)方法分為間接探測(cè)和直接成像,間接探測(cè)有視向速度法、天體測(cè)量法、脈沖星計(jì)時(shí)法、凌星法、微引力透鏡法等[2],目前實(shí)際觀測(cè)中采用最多的是視向速度法和凌星法.天體測(cè)量法至今只發(fā)現(xiàn)了一顆系外行星—HD 176051 b[3].2009年以前,由于HARPS(High Accuracy Radial velocity Planet Searcher)的應(yīng)用[4],視向速度法為主要的探測(cè)方法;2009年美國(guó)國(guó)家航空航天局(National Aeronautics and Space Adm inistration,NASA)建造的開普勒太空望遠(yuǎn)鏡發(fā)射升空,凌星法逐漸成為主流方法,目前該方法已發(fā)現(xiàn)超過1200顆系外行星(exop lanets.eu).2013年12月19日,歐洲空間局(European Space Agency,ESA)研制的gaia空間望遠(yuǎn)鏡升空,其目的是以前所未有的精度對(duì)銀河系內(nèi)數(shù)以十億計(jì)的恒星進(jìn)行觀測(cè),測(cè)量他們的位置、距離和運(yùn)動(dòng),其觀測(cè)效率將比同樣由ESA發(fā)射的依巴谷衛(wèi)星高出數(shù)百萬倍[5].因此,人們對(duì)于利用gaia開展太陽系與系外行星的探測(cè)寄予厚望.

    天體測(cè)量法是通過測(cè)量恒星空間位置變化來反演其附近行星的質(zhì)量和軌道參數(shù)的,相比于目前通常采用的視向速度法和凌星法,天體測(cè)量法有其獨(dú)特的優(yōu)勢(shì):(1)直接確定行星質(zhì)量.視向速度法也是通過測(cè)量恒星空間位置變化來反演其附近行星的軌道參數(shù)的,但因其觀測(cè)方程中僅包含了行星質(zhì)量m與行星軌道傾角的正弦sin i的乘積項(xiàng),所以只能推算出m sin i,不能直接確定行星的質(zhì)量.天體測(cè)量法對(duì)應(yīng)的觀測(cè)方程不僅包含了m cos i項(xiàng),而且包含m項(xiàng),故可以通過觀測(cè)資料同時(shí)直接推算出行星的質(zhì)量與軌道傾角[6].(2)尋找長(zhǎng)周期行星.根據(jù)開普勒第三定律可知,行星繞恒星運(yùn)動(dòng)的軌道半長(zhǎng)徑3次冪與行星軌道周期的平方成正比;同樣質(zhì)量的行星,行星軌道半長(zhǎng)徑越長(zhǎng),對(duì)于恒星繞系統(tǒng)質(zhì)心運(yùn)動(dòng)的擾動(dòng)就越大,因此越容易被天體測(cè)量法分辨出來.所以天體測(cè)量法對(duì)長(zhǎng)周期行星探測(cè)較為敏感.(3)檢測(cè)多星系統(tǒng)是否共面.在已知的系外行星中,大多數(shù)具有較大的軌道偏心率,其原因可能是由于受到一顆或多顆大質(zhì)量天體的引力作用[7],該作用力使得行星的軌道面不與恒星的赤道面共面,若在多星系統(tǒng)中,該作用力將使得各行星軌道不共面[8].研究行星質(zhì)量和相互夾角對(duì)于多星系統(tǒng)的動(dòng)力學(xué)模型、動(dòng)力學(xué)穩(wěn)定以及軌道共振具有重要意義.(4)太陽系附近尋找類地行星.利用空間干涉的天體測(cè)量能達(dá)到亞微角秒級(jí)的精度,對(duì)于造成恒星擾動(dòng)的行星質(zhì)量檢測(cè)可達(dá)到m⊕(地球質(zhì)量)量級(jí),滿足在太陽系附近尋找類地行星的要求[9].

    對(duì)于由1顆恒星和1顆行星組成的二體系外行星系統(tǒng),W right和Howard于2009年建立了基于軌道根數(shù)的天體測(cè)量法對(duì)應(yīng)的觀測(cè)方程[10],其在赤道坐標(biāo)系下,包含了視差、自行以及行星質(zhì)量和軌道根數(shù)等參數(shù).通過分析二體運(yùn)動(dòng)的解和參數(shù)估計(jì)理論,可以給出這些參數(shù)的最優(yōu)解.稍加分析即可發(fā)現(xiàn),基于軌道根數(shù)的天體測(cè)量方法,當(dāng)行星軌道偏心率e=0時(shí),因狀態(tài)轉(zhuǎn)移矩陣奇異而失效,特別是對(duì)于多行星系統(tǒng),基于軌道根數(shù)的狀態(tài)方程和狀態(tài)轉(zhuǎn)移矩陣計(jì)算都非常復(fù)雜,不便于實(shí)際應(yīng)用,因此需要建立簡(jiǎn)潔可行的確定系外行星軌道的天體測(cè)量方法.本文針對(duì)多體系統(tǒng)動(dòng)力學(xué)特性,提出采用坐標(biāo)速度描述,不僅可以避免e=0對(duì)應(yīng)的狀態(tài)轉(zhuǎn)移矩陣奇異問題,而且可以顯著降低整個(gè)計(jì)算的復(fù)雜性.本文第2節(jié)對(duì)gaia測(cè)量數(shù)據(jù)用于系外行星的確定進(jìn)行了必要的精度分析;第3節(jié)較全面地總結(jié)了基于軌道根數(shù)的天體測(cè)量方法涉及的理論基礎(chǔ);第4節(jié)給出了基于坐標(biāo)速度的天體測(cè)量方法對(duì)應(yīng)的觀測(cè)方程、運(yùn)動(dòng)方程和狀態(tài)轉(zhuǎn)移方程;第5節(jié)列出了仿真模擬計(jì)算結(jié)果并作了一些分析;第6節(jié)推導(dǎo)了行星軌道確定時(shí)初始估計(jì)滿足的約束條件,由此可以顯著提高軌道確定微分改正的收斂效率;第7節(jié)推導(dǎo)了三體問題(1顆恒星2顆行星)的計(jì)算公式;最后1節(jié)給出本文的結(jié)論,同時(shí)也討論存在的問題及其解決途徑.

    2 天體測(cè)量觀測(cè)精度分析

    當(dāng)恒星附近存在行星時(shí),由于引力作用,使得恒星與行星繞著它們的質(zhì)心以相同的角速度作橢圓運(yùn)動(dòng).假設(shè)恒星和行星圍繞系統(tǒng)質(zhì)心運(yùn)動(dòng)軌道半長(zhǎng)徑分別為as和ap,則在地面上觀察恒星軌道半長(zhǎng)徑時(shí),可產(chǎn)生一個(gè)角度2θ[11](如圖1).

    其中,Ms和mp分別為恒星和行星的質(zhì)量,G為萬有引力常數(shù),D為恒星與觀測(cè)者之間的距離,p為二體運(yùn)動(dòng)的軌道周期.角半徑θ也稱作恒星的天體測(cè)量特征角,表征了行星對(duì)恒星的擾動(dòng)大小,除行星質(zhì)量外,還與恒星的質(zhì)量、軌道周期及觀測(cè)距離有關(guān).根據(jù)目前已發(fā)現(xiàn)的系外行星,表1列出了它們的質(zhì)量與軌道周期分布范圍,以及所繞恒星的質(zhì)量和與觀測(cè)距離的分布范圍(參見網(wǎng)站exop lanets.eu),其中行星質(zhì)量單位為MJ(木星質(zhì)量),恒星質(zhì)量單位為M⊙(太陽質(zhì)量).針對(duì)上述參數(shù)分布范圍,本文選取了5種模型,模型中恒星質(zhì)量皆取為1 M⊙,行星質(zhì)量、軌道周期以及觀測(cè)距離不等,計(jì)算5種模型中恒星的天體測(cè)量特征角,見表2.

    圖1 二體運(yùn)動(dòng)的軌道示意圖Fig.1 The orbita l sketch of the tw o-body m otion

    表1 已知系外行星的部分參數(shù)范圍Tab le 1 The partial param eters’range of the know n exop lanets

    表2 幾種模型下恒星天體測(cè)量特征角Tab le 2 The astrom etric signatu re of the star in d ifferen t m odels

    由于gaia對(duì)于不同星等的恒星觀測(cè)精度不同,所以還需分析gaia對(duì)于系外行星的可觀測(cè)性.表3給出了對(duì)于太陽質(zhì)量級(jí)的恒星,當(dāng)其觀測(cè)距離在1000 pc內(nèi)時(shí)的星等及gaia對(duì)應(yīng)的位置測(cè)量精度,由表可知:在1000 pc內(nèi),恒星的視星等mV=14.75mag,gaia的觀測(cè)精度高于10μas,其值均小于表2中各模型的天體測(cè)量特征角,因此gaia可進(jìn)行有效觀測(cè).

    表3 太陽質(zhì)量級(jí)恒星在不同距離的視星等及gaia對(duì)其的觀測(cè)精度Tab le 3 The ap paren t m agn itude of the sun like stars at d ifferen t d istances and the accu racies of the gaia

    3 基于軌道根數(shù)的觀測(cè)方程

    利用天體測(cè)量法間接探測(cè)系外行星時(shí),通過觀測(cè)恒星的周期性位置變化來判斷其是否存在.對(duì)于1顆孤立的恒星,可通過5個(gè)天體測(cè)量觀測(cè)值來描述其在天球上的角位置,分別是赤道坐標(biāo)系下的赤經(jīng)α0、赤緯δ0(在某一特定的坐標(biāo)系下某一歷元時(shí)的值,現(xiàn)在一般都取在國(guó)際天球參考系下J2000.0的值),恒星自行相應(yīng)的正交分量μα、μδ,以及視差?;繞著恒星作橢圓運(yùn)動(dòng)的行星軌道可由6個(gè)軌道根數(shù)描述,即軌道半長(zhǎng)徑a、傾角i、升交點(diǎn)經(jīng)度?、近星點(diǎn)經(jīng)度ω、偏心率e及平近點(diǎn)角M;行星質(zhì)量m可由p和a通過Kep ler第三定律確定.因此,對(duì)一個(gè)含有n顆行星的系統(tǒng)進(jìn)行軌道擬合時(shí),需要5+7n個(gè)參數(shù).

    建立如圖2所示觀測(cè)示意圖,恒星的運(yùn)動(dòng)軌道處于O-X YZ坐標(biāo)系內(nèi),O-xyz系為觀測(cè)面所在坐標(biāo)系,z軸表示觀測(cè)者視線方向,x軸指向北極,y方向由右手法則確定,i、?、ω為3個(gè)方位角.

    圖2 二體運(yùn)動(dòng)的觀測(cè)示意圖Fig.2 T he sketch of the observation of the tw o-body m otion

    恒星在空間位置的變化反映了其橢圓軌道運(yùn)動(dòng)在觀測(cè)面上的投影,由赤經(jīng)和赤緯兩個(gè)位置參數(shù)決定,分別用?η(t)、?ξ(t)表示,對(duì)應(yīng)的觀測(cè)方程為

    等式左邊為觀測(cè)值,右邊各項(xiàng)依次為初始值、視差改正、自行改正以及天體擾動(dòng)項(xiàng),其中,?α、?δ表示視差改正參數(shù);?xs、?ys表示恒星運(yùn)動(dòng)在觀測(cè)面上兩個(gè)方向的投影,單位為角度;為了計(jì)算和公式表達(dá)的方便,選取適當(dāng)?shù)挠?jì)算單位是必要的,本文分別采取太陽質(zhì)量M⊙、太陽赤道半徑A⊙作為基本質(zhì)量和長(zhǎng)度單位,取引力常數(shù)G=1,則相應(yīng)的時(shí)間單位為

    故可通過觀測(cè)距離D將?xs、?ys歸算到系統(tǒng)長(zhǎng)度單位下,分別用xs、ys表示

    考慮恒星附近僅存在1顆行星情況,即n=1,由天體力學(xué)二體問題分析解可將xs、ys表示為

    其中,

    A、B、F、G稱為Thiele-Innes常數(shù)[13],可由軌道根數(shù)描述,

    恒星視差與自行改正可通過一定量天體測(cè)量數(shù)據(jù)計(jì)算出來,故在后面的計(jì)算中將視差與自行視為已知量.至此,我們可將觀測(cè)方程改寫為

    這里,η′(t)和ξ′(t)包含了自行和視差改正項(xiàng),已歸算為系統(tǒng)長(zhǎng)度單位.假設(shè)初始時(shí)刻恒星軌道根數(shù)和行星質(zhì)量的估計(jì)值為a?、i?、??、ω?、e?、M?及m?,為方便起見,用矢量將其表示為σ=(σ1,σ2,σ3,σ4,σ5,σ6,σ7).對(duì)上式在恒星初始運(yùn)動(dòng)位置處展開,有:

    其中,Pi(t)、Qi(t)依次表示對(duì)7個(gè)參數(shù)的微分結(jié)果.如果有k次觀測(cè),那么2k個(gè)觀測(cè)方程可以寫成矩陣形式

    其中,

    上述方程(10)包含7個(gè)未知量,因此問題可解至少需要觀測(cè)4次.

    將由估計(jì)理論首次解得的?σ(0)疊加到初始時(shí)刻近似軌道上,進(jìn)行迭代計(jì)算,當(dāng)?shù)趌次迭代給出的|?σ(l)|小于設(shè)定的收斂標(biāo)準(zhǔn)ε時(shí),求解完成.本文統(tǒng)一設(shè)定誤差限ε= 10?6.

    由上述獲得的恒星軌道根數(shù)與行星質(zhì)量的估計(jì)值,可以求出恒星初始時(shí)刻的坐標(biāo)→rs與速度˙→rs,由此根據(jù)二體關(guān)系式,可以求出行星此時(shí)對(duì)應(yīng)的坐標(biāo)→rp與速度˙→rp,

    對(duì)于多行星系統(tǒng),即n≥2,對(duì)應(yīng)的運(yùn)動(dòng)沒有分析解,故上述方法不再適用.

    4 基于坐標(biāo)速度的觀測(cè)方程

    觀測(cè)方程與軌道根數(shù)法對(duì)應(yīng)的恒星觀測(cè)方程相同,即(2)式.假定t=t0時(shí)刻恒星位置與速度和行星質(zhì)量的估計(jì)值分別為(t0)、(t0)、mp,其近似值(即初始估計(jì))為、、m?,兩者相差為?、?、?m,將觀測(cè)方程在、、m?處展開有

    當(dāng)t=t0時(shí),狀態(tài)轉(zhuǎn)移方程涉及的初始條件為

    這里,H為3階單位矩陣.

    5 軌道特征分析與仿真計(jì)算

    目前,gaia的巡天數(shù)據(jù)仍未釋放,因此為了驗(yàn)證上述方法的有效性,只能進(jìn)行仿真模擬計(jì)算.假定t0時(shí)刻1組恒星的軌道根數(shù)和行星質(zhì)量σ(t0),由該組參數(shù)可模擬出1組恒星位置的觀測(cè)數(shù)據(jù),然后用其反解σ(t0),由此可以驗(yàn)證和分析評(píng)估軌道確定的方法和精度.

    5.1 軌道特征的影響

    5.1.1 小偏心率

    對(duì)于小偏心率軌道,軌道根數(shù)法不能進(jìn)行軌道確定,但坐標(biāo)速度法仍有效,其原因是e=0時(shí)近星點(diǎn)經(jīng)度不能確定所致.具體體現(xiàn)在當(dāng)e=0時(shí),xs、ys可表示為:

    由上式可以看出

    即狀態(tài)轉(zhuǎn)移矩陣中有關(guān)ω和M的兩列相關(guān),導(dǎo)致狀態(tài)轉(zhuǎn)移矩陣奇異,故對(duì)于小偏心率的情況,軌道根數(shù)法失效.

    5.1.2 軌道傾角

    對(duì)于傾角i=0°或180°的軌道,兩種方法都不適用.

    軌道根數(shù)法中,由于xs、ys中包含的i項(xiàng)都以cos i出現(xiàn),故當(dāng)i=0°或180°時(shí),xs、ys對(duì)i的偏導(dǎo)數(shù)都為0,具體體現(xiàn)在

    坐標(biāo)速度法中,當(dāng)i=0°時(shí),中不顯含z分量,即z?和分量.故有

    兩種情況都導(dǎo)致線性方程組的系數(shù)矩陣奇異,無唯一解.

    此外,天體測(cè)量法描述的是恒星X、Y方向的運(yùn)動(dòng)在觀測(cè)面上的投影,對(duì)于Z方向沒有直接約束,這就造成了當(dāng)i?π?i時(shí),σ(t0)一樣能滿足觀測(cè)方程.但此時(shí)除傾角項(xiàng)外,升交點(diǎn)的經(jīng)度?和近地點(diǎn)的經(jīng)度ω也存在相應(yīng)的變化,這一點(diǎn)我們可從Thiele-Innes常數(shù)來理解.我們將其對(duì)應(yīng)關(guān)系整理如下:

    5.2 觀測(cè)誤差的影響

    觀測(cè)即有誤差,因此不同的仿真模型應(yīng)增加相應(yīng)的觀測(cè)誤差.以第2節(jié)中模型A為例,mV=5.63mag<10mag,此時(shí)gaia的觀測(cè)精度?θ=4μas,若將其視為天體測(cè)量觀測(cè)誤差,對(duì)于觀測(cè)距離D=15 pc,該誤差s的大小為

    將其歸算到系統(tǒng)長(zhǎng)度單位A⊙,即得該模型的觀測(cè)誤差?s為

    通過M atlab軟件可計(jì)算出任意兩組滿足均值為0、方差為0.013的高斯白噪聲,將其分別增加到恒星的x和y方向的觀測(cè)數(shù)據(jù)上,則可獲得含有觀測(cè)誤差的模擬觀測(cè)數(shù)據(jù).

    由含有誤差觀測(cè)數(shù)據(jù)可確定初始根數(shù)σ(t0)的估計(jì)值ˉσ(t0),兩者之差為?σ.表4列出了幾種模型對(duì)于不同觀測(cè)誤差和不同觀測(cè)次數(shù)時(shí)確定的初始軌道估計(jì)值誤差情況.模型A表示一個(gè)軌道周期為2 yr的木星繞一個(gè)距離地球15 pc的太陽作橢圓運(yùn)動(dòng)(見表2),此太陽在天球上位置變化最大角度(特征角)為100μas,考慮不同觀測(cè)誤差與觀測(cè)次數(shù)的組合,形成了3組模擬觀測(cè)數(shù)據(jù):A-1、A-2和A-3.模型A的初始軌道根數(shù)取為:a=1.587 au,e=0.4000,i=40.00°,?=80.00°,ω=50.00°,M=50.00°, m=9.552×10?4Ms;軌道確定誤差統(tǒng)計(jì)分為x、y方向的殘差及軌道根數(shù)的相對(duì)誤差.以A模型為例,對(duì)于特征角為100μas的系外行星系統(tǒng),考慮觀測(cè)隨機(jī)誤差分別為2μas和4μas以及觀測(cè)點(diǎn)數(shù)分別為20和40的情況.模型A-1對(duì)應(yīng)觀測(cè)誤差2μas,觀測(cè)點(diǎn)數(shù)20,由坐標(biāo)速度法確定的軌道殘差的x和y分量分別為23.75μas和24.20μas,由軌道根數(shù)法確定的軌道半長(zhǎng)徑和軌道偏心率的相對(duì)誤差分別為1.65%和0.92%,其他根數(shù)的誤差均小于10%,質(zhì)量確定的精度較高,可好于千分之一;模型A-2對(duì)應(yīng)觀測(cè)誤差2μas,觀測(cè)點(diǎn)數(shù)40,相較于A-1各根數(shù)的相對(duì)誤差均顯著減小,即擬合結(jié)果更優(yōu);模型A-3對(duì)應(yīng)觀測(cè)誤差4μas,觀測(cè)點(diǎn)數(shù)40,相較于A-2大部分根數(shù)的相對(duì)誤差有所增大.通過比較不同模型對(duì)應(yīng)的計(jì)算結(jié)果表明,觀測(cè)數(shù)據(jù)越多,相對(duì)誤差越小,軌道確定的擬合結(jié)果越可靠.

    表4 模型A-1、A-2和A-3的計(jì)算結(jié)果Tab le 4 T he ca lcu lation s of the m od els A-1,A-2,an d A-3

    5.3 初值的選取

    由于觀測(cè)方程(2)中關(guān)于恒星初始軌道根數(shù)或恒星初始坐標(biāo)速度函數(shù)的表示均是非線性的,故在軌道確定中需要采用微分迭代方法求解軌道估計(jì)值,因此需要選定恰當(dāng)?shù)牡踔?以使求解過程收斂.對(duì)于軌道初值的選取,最簡(jiǎn)單的方法是在一定區(qū)間內(nèi)對(duì)每個(gè)待解參數(shù)進(jìn)行掃描計(jì)算,當(dāng)?shù)踔当贿x取在參數(shù)收斂區(qū)間時(shí),便可給出估計(jì)值.本文討論的問題涉及的參數(shù)收斂區(qū)間是7維的,因此要獲得收斂的估計(jì)值,需要較大的掃描計(jì)算量.為了方便計(jì)算,采用固定?、ω、M、m進(jìn)行a、i、e 3維檢索,再固定a、i、e、m進(jìn)行?、ω、M 3維檢索.計(jì)算結(jié)果列于表5,表中列出的是相對(duì)誤差百分比值,需要說明的是其中角度量的誤差相對(duì)量取180,偏心率的誤差相對(duì)量取為1,主要是為了避免小偏心率作為相對(duì)量導(dǎo)致相對(duì)誤差失真.模型A-1和A-3的軌道半長(zhǎng)徑收斂區(qū)間約為5%,其他根數(shù)的收斂區(qū)間相對(duì)較大,一般大于10%,有些甚至可達(dá)20%.從初步的掃描計(jì)算結(jié)果看,軌道半長(zhǎng)徑a的迭代初值選取是影響計(jì)算收斂與否的關(guān)鍵因素,需要盡量靠近估計(jì)值.

    表5 模型A-1和A-3的收斂區(qū)間Tab le 5 The convergen ce in terval of the m od els A-1 and A-3

    6 關(guān)于初始估計(jì)取值的約束

    在上節(jié)討論中,為減少計(jì)算時(shí)間,對(duì)7維參數(shù)空間僅進(jìn)行了兩個(gè)3維參數(shù)空間掃描,且行星質(zhì)量均是固定的,沒有參與掃描計(jì)算.事實(shí)上,可通過恒星空間位置觀測(cè)資料時(shí)間序列對(duì)上述7維參數(shù)空間進(jìn)行適當(dāng)?shù)募s束,降低參數(shù)空間的維數(shù).具體做法如下:已知由天體測(cè)量方法給出的某顆恒星的位置時(shí)間序列xs(tj),ys(tj),j=1,2,···,k,其包含了運(yùn)動(dòng)周期信息,通過對(duì)其傅里葉譜分析(FFT)可獲得軌道周期p,再根據(jù)開普勒第三定律可以給出恒星軌道半長(zhǎng)徑a和行星質(zhì)量mp之間的關(guān)系式,

    一般認(rèn)為Ms的值可由天體物理方法獲得.

    傅里葉譜分析通常針對(duì)等間距數(shù)據(jù)(均勻采樣)進(jìn)行,故需要先對(duì)恒星的觀測(cè)數(shù)據(jù)進(jìn)行插值,再利用插值數(shù)據(jù)進(jìn)行譜分析.此外,還采用了另一種方法——LSQ(最小二乘擬合),將插值后的數(shù)據(jù)進(jìn)行最小二乘,擬合成一條正弦曲線,該正弦曲線的周期即為恒星運(yùn)動(dòng)周期.

    已知gaia的觀測(cè)年限為5 yr,故進(jìn)行傅里葉譜分析所選取的模型周期不應(yīng)該大于5 yr,否則無法觀測(cè)到一個(gè)完整周期的數(shù)據(jù);又根據(jù)奈奎斯特取樣定理,所選模型周期應(yīng)限制在2.5 yr內(nèi).此處,我們選取5.2節(jié)中所建立的模型進(jìn)行計(jì)算,圖3給出了A-3模型(任意選取)x、y方向誤差數(shù)據(jù)的分布及進(jìn)行傅里葉分析的結(jié)果,并將計(jì)算結(jié)果列于表6中.

    表6 模型A-1、A-2、A-3進(jìn)行最小二乘擬合和傅里葉分析的結(jié)果Tab le 6 T he LSQ an d FFT resu lts of the m od els A-1,A-2,and A-3

    圖3 模型A-3的x、y方向數(shù)據(jù)進(jìn)行傅里葉分析的結(jié)果Fig.3 T he resu lts of x and y com p onen ts of the m odel A-3 with FFT

    根據(jù)統(tǒng)計(jì),LSQ和FFT均能實(shí)現(xiàn)觀測(cè)數(shù)據(jù)的周期分析.模型A的精確周期為730 d,3個(gè)子模型進(jìn)行FFT的結(jié)果分別為720 d、716.67 d和716.67 d,明顯要好于LSQ的結(jié)果.

    對(duì)序列{xs(tj),ys(tj),j=1,2,···,k}在時(shí)間域上作平均,可以獲得如下關(guān)系式,假定已知行星質(zhì)量mp、恒星的平近點(diǎn)角M、恒星的升交點(diǎn)經(jīng)度?和恒星的近星點(diǎn)經(jīng)度ω,可以由上述3個(gè)約束關(guān)系式(21)–(23)式,近似給出恒星的軌道半長(zhǎng)徑、偏心率和軌道傾角,這樣將對(duì)7維參數(shù)空間掃描問題降低到4維參數(shù)空間掃描,可以顯著提高計(jì)算效率.當(dāng)然,在由觀測(cè)時(shí)間序列計(jì)算上述平均值時(shí),會(huì)產(chǎn)生計(jì)算誤差,將會(huì)直接影響到上述約束條件的精度.

    7 多行星系統(tǒng)的計(jì)算

    坐標(biāo)速度法可以作為天體測(cè)量法探測(cè)系外行星的基本方法,不僅是因?yàn)樗梢杂糜谛∑穆受壍烙?jì)算,而且較容易將其推廣到多行星系統(tǒng),求解各個(gè)行星的參數(shù)(位置、速度及質(zhì)量).下面僅推導(dǎo)由1顆恒星和2顆行星組成的三體系統(tǒng)中恒星觀測(cè)方程計(jì)算具體關(guān)系式,n體問題的推導(dǎo)思路類似.

    建立如圖4所示慣性坐標(biāo)系,恒星與兩顆行星質(zhì)量分別為Ms、m1和m2,O為系統(tǒng)的質(zhì)量中心,空間位置矢量如圖所示.

    圖4 三體運(yùn)動(dòng)坐標(biāo)示意圖Fig.4 T he coord inate’s sketch of the th ree-body m otion

    3個(gè)天體對(duì)應(yīng)的運(yùn)動(dòng)方程為

    因系統(tǒng)質(zhì)心為坐標(biāo)原點(diǎn),由質(zhì)心運(yùn)動(dòng)定理知:

    方程(25)左邊分別用η′(t)、ξ′(t)表示,可將其改寫為

    由(24)–(25)式可求出行星運(yùn)動(dòng)方程為

    狀態(tài)轉(zhuǎn)移方程涉及到的初始條件為:當(dāng)t=t0時(shí),

    這里,H為3階單位矩陣.

    8 結(jié)論與討論

    本文基于天體測(cè)量法給出的恒星空間位置參數(shù)序列,討論了在求解行星質(zhì)量和軌道參數(shù)時(shí)涉及的動(dòng)力學(xué)條件方程計(jì)算問題;對(duì)于采用Kepler軌道根數(shù)和采用坐標(biāo)速度描述形式分別給出了具體微分改正公式,同時(shí)也進(jìn)行了初步必要的仿真模擬計(jì)算;通過天體力學(xué)二體問題運(yùn)動(dòng)關(guān)系式,對(duì)恒星空間位置時(shí)間序列做平均分析,推導(dǎo)出了待估參數(shù)間存在的3個(gè)約束條件關(guān)系式,可以顯著提升計(jì)算效率.

    本文沒有對(duì)恒星自行和視差等天體測(cè)量參數(shù)做解算研究,以后若有實(shí)測(cè)數(shù)據(jù),可根據(jù)實(shí)際情況作具體處理:或單獨(dú)求解,或與行星軌道和質(zhì)量一起求解.

    [1]M ayor M,Q ueloz D.Natu re,1995,368:355

    [2]M ichael P.T he Exop lanet Handbook.New York:Cam b ridge Un iversity P ress,2011:3-5

    [3]Luisa G J,Lu is A,Barbara P.M NRAS,2014,443:260

    [4]姜明達(dá),肖東,朱永田.天文學(xué)進(jìn)展,2012,30:214

    [5]Dou J P,Zhu Y T,Ren D Q.Ch inese Jou rna l of Natu re,2014,36:124

    [6]張牛,季江徽.天文學(xué)進(jìn)展,2009,27:14

    [7]Lin D N C,Ida S.A p J,1997,477:781

    [8]Papaloizou J C B,Terquem C.M NRAS,2001,325:221

    [9]金文敬.天文學(xué)進(jìn)展,2013,31:375

    [10]W right J T,How ard A W.Ap JS,2009,182:205

    [11]Sara S.Exop lanets.A rizona:Un iversity of A rizona P ress,2010:327-329

    [12]M ichael P,Joel H,G aspar B,et a l.A p J,2014,797:14

    [13]Th iele T N.A stronom ische Nach rich ten,1883,104:245

    Exop lanet Detection by A strom etric M ethod

    XUWei-wei1,2,3,4LIAO Xin-hao1,3ZHOU Yong-hong1,3XU Xue-qing1,3

    (1 Shanghai A stronom ica l O bserva to ry,Chinese A cadem y of Scien ces,Shanghai 200030) (2 Schoo l of Physica l Scien ce an d Techno logy,ShanghaiTech Un iversity,Shanghai 201210) (3 K ey Labo ra to ry of P laneta ry Scien ces,Chinese A cadem y of Scien ces,Shanghai 200030) (4 Un iversity of Chinese A cadem y of Scien ces,Beijing 100049)

    Aswe known,the exoplanetsaremostly detected by themethodsof radial velocity and transit,only one is found by the astrometric method.As the data of the gaia to be released,astrom etry w ill becom e one of themost im portantmethod for detecting exoplanets gradually.Based on the position sequence of stars,this paper discusses the calculation of the equations of dynam ics conditions involved in solving themass and the orbit parameters of the planet.Due to the deficiency of the available theory(orbital elem ent m ethod),we put forward a new m ethod(coordinate velocity method).The differential correction formulae of the twomethodsare presented,aswell as the necessary simulation.In addition,themethod established in this paper can be app lied to themulti-planet system easily.

    exoplanet,astrometry,celestialmechanics:orbit calculation and determ ination,m ethods:data analysis

    P135;

    A

    10.15940/j.cnki.0001-5245.2016.04.004

    2015-12-08收到原稿,2016-03-04收到修改稿

    ?國(guó)家自然科學(xué)基金項(xiàng)目(11133004)和中國(guó)科學(xué)院先導(dǎo)B項(xiàng)目(XDBO9000000)資助

    ?1143880832@qq.com

    ?xhliao@shao.ac.cn

    猜你喜歡
    外行星測(cè)量法根數(shù)
    外行星
    更正
    尋找規(guī)律巧算根數(shù)
    玉米的胡須
    系外行星那些事——“呼啦圈”法
    基于比較測(cè)量法的冷卻循環(huán)水系統(tǒng)電導(dǎo)率檢測(cè)儀研究
    系外行星探索與發(fā)現(xiàn)
    太空探索(2016年4期)2016-07-12 15:17:52
    系外行星的探測(cè)
    太空探索(2016年3期)2016-07-12 09:58:42
    基于改進(jìn)型號(hào)第二婁無廳點(diǎn)根數(shù)的北斗CEO衛(wèi)星廣播星歷擬合算法及實(shí)現(xiàn)
    垂直面內(nèi)建立基線的特殊點(diǎn)位高程測(cè)量法
    .国产精品久久| 国产伦在线观看视频一区| 一区二区三区精品91| 亚洲av成人精品一区久久| 亚洲精品日本国产第一区| 大码成人一级视频| 亚洲电影在线观看av| 高清毛片免费看| 久久 成人 亚洲| a级毛片免费高清观看在线播放| av免费在线看不卡| 极品人妻少妇av视频| 亚洲真实伦在线观看| 这个男人来自地球电影免费观看 | freevideosex欧美| 色网站视频免费| 亚洲丝袜综合中文字幕| 99热这里只有精品一区| 亚洲精品日本国产第一区| 天美传媒精品一区二区| 国产亚洲91精品色在线| 国产精品无大码| 午夜免费男女啪啪视频观看| 免费在线观看成人毛片| 亚洲内射少妇av| 欧美日韩在线观看h| 久久久久久久久久久免费av| 久久久精品免费免费高清| 97精品久久久久久久久久精品| 久久青草综合色| 丰满乱子伦码专区| 又大又黄又爽视频免费| 日日啪夜夜撸| 国产有黄有色有爽视频| 久久精品国产亚洲av涩爱| 香蕉精品网在线| 一区二区三区四区激情视频| 成人国产av品久久久| 热99国产精品久久久久久7| av又黄又爽大尺度在线免费看| 国产 一区精品| 亚洲欧美一区二区三区黑人 | 亚洲欧美成人综合另类久久久| 国产一区二区三区综合在线观看 | av线在线观看网站| 国产 精品1| av在线播放精品| 国产极品天堂在线| 国产日韩欧美在线精品| 国产亚洲午夜精品一区二区久久| 亚洲精品自拍成人| 三级经典国产精品| 日韩成人伦理影院| 黑人巨大精品欧美一区二区蜜桃 | 80岁老熟妇乱子伦牲交| 亚洲av不卡在线观看| 国产成人免费无遮挡视频| 欧美国产精品一级二级三级 | 亚洲av福利一区| 交换朋友夫妻互换小说| 成年人免费黄色播放视频 | 九九爱精品视频在线观看| 男人和女人高潮做爰伦理| 国产综合精华液| 久久久久久久久大av| 久久精品国产自在天天线| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区黑人 | 午夜福利视频精品| 多毛熟女@视频| 久久99热这里只频精品6学生| 毛片一级片免费看久久久久| 啦啦啦中文免费视频观看日本| 毛片一级片免费看久久久久| 亚洲婷婷狠狠爱综合网| 乱系列少妇在线播放| 青春草国产在线视频| 中文字幕制服av| 男人狂女人下面高潮的视频| 在线亚洲精品国产二区图片欧美 | 欧美3d第一页| 一级黄片播放器| 欧美丝袜亚洲另类| 国产深夜福利视频在线观看| 日日啪夜夜爽| 日韩精品免费视频一区二区三区 | 亚洲熟女精品中文字幕| 国产欧美日韩综合在线一区二区 | 9色porny在线观看| 中文字幕av电影在线播放| av卡一久久| 亚洲av免费高清在线观看| 日韩亚洲欧美综合| 99热6这里只有精品| 婷婷色av中文字幕| 极品教师在线视频| av在线app专区| 精品一区二区三卡| 国产伦精品一区二区三区四那| 中国美白少妇内射xxxbb| 精品人妻熟女毛片av久久网站| 免费大片黄手机在线观看| 亚洲经典国产精华液单| h日本视频在线播放| 亚洲欧美一区二区三区黑人 | 日日啪夜夜撸| 精品久久久久久久久亚洲| 亚洲av.av天堂| 国产精品三级大全| 久久久久久久精品精品| 久久国产亚洲av麻豆专区| 国产成人aa在线观看| 日韩亚洲欧美综合| 一级av片app| 精品亚洲乱码少妇综合久久| 国产精品国产av在线观看| 免费看光身美女| 肉色欧美久久久久久久蜜桃| 国产av精品麻豆| 三级国产精品欧美在线观看| 街头女战士在线观看网站| 国产 一区精品| 桃花免费在线播放| 在线观看国产h片| 久久午夜福利片| 国产精品国产三级国产专区5o| 高清在线视频一区二区三区| 男女无遮挡免费网站观看| 日本黄色日本黄色录像| 久久国产精品大桥未久av | 卡戴珊不雅视频在线播放| 男人和女人高潮做爰伦理| 美女视频免费永久观看网站| 日本91视频免费播放| 国产精品人妻久久久久久| 啦啦啦啦在线视频资源| 一级毛片aaaaaa免费看小| 九九爱精品视频在线观看| 亚洲国产色片| 国产高清三级在线| 国产免费又黄又爽又色| 久久毛片免费看一区二区三区| 亚洲在久久综合| 欧美日韩综合久久久久久| 高清黄色对白视频在线免费看 | 天美传媒精品一区二区| 美女内射精品一级片tv| 久久国内精品自在自线图片| 亚洲欧美中文字幕日韩二区| 下体分泌物呈黄色| 久热久热在线精品观看| 夜夜爽夜夜爽视频| 十八禁网站网址无遮挡 | 中文字幕av电影在线播放| 99久国产av精品国产电影| 天堂俺去俺来也www色官网| 亚洲国产精品国产精品| 少妇被粗大的猛进出69影院 | 69精品国产乱码久久久| 亚洲一级一片aⅴ在线观看| 成人影院久久| 99九九在线精品视频 | 三级国产精品欧美在线观看| 久久99精品国语久久久| 精品久久久久久久久亚洲| 中国三级夫妇交换| 国产精品女同一区二区软件| 五月天丁香电影| 亚洲精品中文字幕在线视频 | 国产男女超爽视频在线观看| 久久久亚洲精品成人影院| 亚洲欧洲国产日韩| 国产亚洲午夜精品一区二区久久| 亚洲高清免费不卡视频| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| 亚洲第一av免费看| 国产成人一区二区在线| 寂寞人妻少妇视频99o| 国产精品无大码| 精品午夜福利在线看| 日韩成人av中文字幕在线观看| 精品一区二区三区视频在线| 蜜桃久久精品国产亚洲av| 人妻 亚洲 视频| 亚洲经典国产精华液单| 蜜臀久久99精品久久宅男| 亚洲自偷自拍三级| 亚洲色图综合在线观看| 男人爽女人下面视频在线观看| 亚洲人成网站在线播| 午夜免费观看性视频| 色94色欧美一区二区| 亚洲真实伦在线观看| 亚洲国产欧美在线一区| 亚洲欧美日韩卡通动漫| 亚洲欧美中文字幕日韩二区| 精品人妻一区二区三区麻豆| 久久精品熟女亚洲av麻豆精品| 国产白丝娇喘喷水9色精品| 最新的欧美精品一区二区| 最近最新中文字幕免费大全7| 欧美国产精品一级二级三级 | 国产日韩欧美在线精品| 欧美日本中文国产一区发布| 中文天堂在线官网| 美女国产视频在线观看| 日日爽夜夜爽网站| 国产中年淑女户外野战色| 又黄又爽又刺激的免费视频.| 国产伦理片在线播放av一区| 精品熟女少妇av免费看| 搡女人真爽免费视频火全软件| 亚洲精品一二三| 亚洲国产成人一精品久久久| 男人添女人高潮全过程视频| 黑丝袜美女国产一区| 26uuu在线亚洲综合色| 亚洲精品国产成人久久av| 老司机亚洲免费影院| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| 欧美日本中文国产一区发布| 日韩av不卡免费在线播放| 一本—道久久a久久精品蜜桃钙片| 久久人人爽av亚洲精品天堂| 在线免费观看不下载黄p国产| 观看免费一级毛片| 不卡视频在线观看欧美| 久久久亚洲精品成人影院| 亚洲自偷自拍三级| 18禁在线无遮挡免费观看视频| 三级国产精品片| 亚洲av成人精品一区久久| 黄色日韩在线| 日本色播在线视频| 狂野欧美激情性xxxx在线观看| 成人无遮挡网站| 国产91av在线免费观看| 亚洲三级黄色毛片| 亚洲丝袜综合中文字幕| 国产av码专区亚洲av| 精品国产乱码久久久久久小说| 97在线视频观看| 春色校园在线视频观看| 中文字幕亚洲精品专区| 另类亚洲欧美激情| 免费观看av网站的网址| 亚洲欧洲国产日韩| 卡戴珊不雅视频在线播放| 亚洲精品中文字幕在线视频 | 国产成人freesex在线| 欧美亚洲 丝袜 人妻 在线| √禁漫天堂资源中文www| 大香蕉97超碰在线| 久久久国产欧美日韩av| 免费av不卡在线播放| 国产欧美日韩综合在线一区二区 | 曰老女人黄片| 久久久久久久亚洲中文字幕| av福利片在线观看| 婷婷色综合大香蕉| 一级毛片久久久久久久久女| 午夜激情福利司机影院| av视频免费观看在线观看| 亚洲av日韩在线播放| 又粗又硬又长又爽又黄的视频| 国产精品熟女久久久久浪| 国产色婷婷99| 交换朋友夫妻互换小说| 精品久久国产蜜桃| 国产 精品1| 两个人免费观看高清视频 | 大香蕉97超碰在线| 在线观看一区二区三区激情| 国产高清不卡午夜福利| 一区二区av电影网| 嫩草影院入口| 99热这里只有精品一区| .国产精品久久| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 国产69精品久久久久777片| 国产又色又爽无遮挡免| 午夜日本视频在线| 国产日韩一区二区三区精品不卡 | 一级毛片久久久久久久久女| 免费看av在线观看网站| 黄色一级大片看看| 又大又黄又爽视频免费| 91久久精品国产一区二区成人| 欧美成人午夜免费资源| 国产高清不卡午夜福利| 亚洲av.av天堂| 中文字幕精品免费在线观看视频 | 午夜久久久在线观看| 国产淫片久久久久久久久| 人妻 亚洲 视频| 自拍偷自拍亚洲精品老妇| 男女国产视频网站| 国产日韩欧美视频二区| 在线精品无人区一区二区三| 少妇被粗大猛烈的视频| 亚洲,一卡二卡三卡| 丝瓜视频免费看黄片| 肉色欧美久久久久久久蜜桃| 精品久久久久久久久av| 久久精品国产亚洲网站| 亚洲,欧美,日韩| 人人妻人人澡人人看| 国产在视频线精品| 中文字幕久久专区| 人妻夜夜爽99麻豆av| 国国产精品蜜臀av免费| 久久99热6这里只有精品| 欧美+日韩+精品| 汤姆久久久久久久影院中文字幕| 精品久久久久久久久av| 三级国产精品片| 亚洲欧美日韩东京热| 99热国产这里只有精品6| 热re99久久国产66热| 欧美一级a爱片免费观看看| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 成人二区视频| 女人久久www免费人成看片| 一本色道久久久久久精品综合| 中文乱码字字幕精品一区二区三区| 秋霞伦理黄片| 日产精品乱码卡一卡2卡三| 日本av手机在线免费观看| 精品99又大又爽又粗少妇毛片| 久久久国产精品麻豆| 亚洲精品色激情综合| 少妇 在线观看| 内射极品少妇av片p| 一边亲一边摸免费视频| 青青草视频在线视频观看| 狠狠精品人妻久久久久久综合| 欧美日韩一区二区视频在线观看视频在线| 国产国拍精品亚洲av在线观看| 秋霞在线观看毛片| 国产精品人妻久久久久久| 亚洲熟女精品中文字幕| 亚洲成人手机| 国产高清三级在线| 亚洲av中文av极速乱| 校园人妻丝袜中文字幕| 亚洲不卡免费看| 99视频精品全部免费 在线| 大片电影免费在线观看免费| 色视频www国产| 日本黄色日本黄色录像| 99久久中文字幕三级久久日本| 国产成人aa在线观看| 精品少妇黑人巨大在线播放| 又大又黄又爽视频免费| 老司机亚洲免费影院| 麻豆成人av视频| 国产精品秋霞免费鲁丝片| 在现免费观看毛片| 国产精品人妻久久久久久| 最黄视频免费看| 日本欧美国产在线视频| av视频免费观看在线观看| 日韩中字成人| 久久久久精品性色| 免费看光身美女| 大码成人一级视频| 日韩一区二区三区影片| 成人免费观看视频高清| 国产一区亚洲一区在线观看| 亚洲人与动物交配视频| 国产成人精品福利久久| 亚洲自偷自拍三级| 丝袜在线中文字幕| 你懂的网址亚洲精品在线观看| 97在线视频观看| 老熟女久久久| 亚洲精品乱码久久久久久按摩| 久久久久久久精品精品| 97在线人人人人妻| 欧美精品人与动牲交sv欧美| 观看美女的网站| 又黄又爽又刺激的免费视频.| 亚洲欧美精品专区久久| 一本久久精品| 国产精品99久久久久久久久| av又黄又爽大尺度在线免费看| 久久精品久久久久久久性| .国产精品久久| 久久ye,这里只有精品| 久久精品国产鲁丝片午夜精品| 国产免费福利视频在线观看| 久久人人爽av亚洲精品天堂| 自线自在国产av| 亚洲精品亚洲一区二区| 一区二区三区免费毛片| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 国产成人aa在线观看| 女人久久www免费人成看片| 一级a做视频免费观看| 中文字幕免费在线视频6| 欧美精品人与动牲交sv欧美| 一二三四中文在线观看免费高清| 黄色欧美视频在线观看| 久久精品国产自在天天线| 曰老女人黄片| 色94色欧美一区二区| 国产日韩一区二区三区精品不卡 | 精品久久久精品久久久| 中文资源天堂在线| 亚洲av成人精品一二三区| 国产 精品1| 亚洲成人av在线免费| 我要看日韩黄色一级片| 亚洲人成网站在线播| 亚洲精品乱码久久久久久按摩| 午夜91福利影院| 久久久久久久久久成人| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 制服丝袜香蕉在线| 美女福利国产在线| 在线播放无遮挡| 啦啦啦在线观看免费高清www| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 成人无遮挡网站| 久久av网站| 亚洲国产欧美在线一区| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| 丰满少妇做爰视频| 欧美变态另类bdsm刘玥| 日韩不卡一区二区三区视频在线| 亚洲欧美日韩另类电影网站| 色网站视频免费| 免费观看性生交大片5| av视频免费观看在线观看| 亚洲av.av天堂| 如何舔出高潮| 男女边摸边吃奶| 亚洲综合精品二区| 欧美日本中文国产一区发布| 欧美另类一区| 成人美女网站在线观看视频| 久久精品国产自在天天线| 麻豆乱淫一区二区| 夜夜看夜夜爽夜夜摸| 在线观看美女被高潮喷水网站| 毛片一级片免费看久久久久| 精品少妇久久久久久888优播| 大又大粗又爽又黄少妇毛片口| 日韩伦理黄色片| 美女视频免费永久观看网站| 这个男人来自地球电影免费观看 | 国产午夜精品久久久久久一区二区三区| 亚洲欧美清纯卡通| 少妇高潮的动态图| 国产综合精华液| 人妻 亚洲 视频| 国产精品国产三级国产av玫瑰| 精品99又大又爽又粗少妇毛片| 国产午夜精品久久久久久一区二区三区| 男人爽女人下面视频在线观看| 久久午夜综合久久蜜桃| 老熟女久久久| 久久99精品国语久久久| 在线 av 中文字幕| 全区人妻精品视频| 国精品久久久久久国模美| 久久亚洲国产成人精品v| 国产男女内射视频| 在线观看人妻少妇| 26uuu在线亚洲综合色| 国产精品国产三级国产专区5o| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 久热这里只有精品99| 爱豆传媒免费全集在线观看| 中文字幕制服av| 亚洲精品视频女| 我要看黄色一级片免费的| 国产精品一区二区三区四区免费观看| 亚洲欧洲精品一区二区精品久久久 | 曰老女人黄片| 99精国产麻豆久久婷婷| 免费观看无遮挡的男女| 亚洲国产精品999| 男女国产视频网站| 一级爰片在线观看| 日本vs欧美在线观看视频 | 日本av免费视频播放| 在线免费观看不下载黄p国产| 午夜福利网站1000一区二区三区| 午夜老司机福利剧场| 免费大片黄手机在线观看| 综合色丁香网| 蜜桃在线观看..| 亚洲欧洲日产国产| 男人添女人高潮全过程视频| 97超视频在线观看视频| 乱码一卡2卡4卡精品| 日韩不卡一区二区三区视频在线| 国产精品久久久久久精品古装| 亚洲av中文av极速乱| av女优亚洲男人天堂| 人妻夜夜爽99麻豆av| 99久国产av精品国产电影| 亚洲电影在线观看av| 成人无遮挡网站| 免费av不卡在线播放| 亚洲高清免费不卡视频| 国产精品一区二区三区四区免费观看| 在线观看美女被高潮喷水网站| 午夜免费男女啪啪视频观看| 男女边吃奶边做爰视频| 美女中出高潮动态图| 又粗又硬又长又爽又黄的视频| 免费黄频网站在线观看国产| 久久久久精品久久久久真实原创| 特大巨黑吊av在线直播| 亚洲熟女精品中文字幕| 99热这里只有是精品在线观看| 久久99精品国语久久久| 另类精品久久| 国产一区二区三区av在线| 欧美丝袜亚洲另类| 国产一区二区三区av在线| 一级av片app| 日本色播在线视频| 天堂俺去俺来也www色官网| 黄片无遮挡物在线观看| 性色av一级| 亚洲精品日韩av片在线观看| 日本黄色片子视频| 欧美精品一区二区免费开放| 建设人人有责人人尽责人人享有的| 久久免费观看电影| 最近的中文字幕免费完整| 欧美精品人与动牲交sv欧美| 丝袜脚勾引网站| 我的老师免费观看完整版| 老司机影院成人| 精品人妻偷拍中文字幕| 少妇熟女欧美另类| 自拍偷自拍亚洲精品老妇| 啦啦啦啦在线视频资源| 蜜臀久久99精品久久宅男| 韩国av在线不卡| 亚洲一区二区三区欧美精品| 免费黄色在线免费观看| 久久99精品国语久久久| 啦啦啦啦在线视频资源| 国产一区二区在线观看日韩| 纯流量卡能插随身wifi吗| 一区二区av电影网| 国产在线一区二区三区精| 亚洲综合色惰| 在线观看一区二区三区激情| 久久人人爽人人爽人人片va| 青春草视频在线免费观看| 国产亚洲午夜精品一区二区久久| 赤兔流量卡办理| av天堂久久9| 一个人看视频在线观看www免费| 日本av手机在线免费观看| 大香蕉久久网| 毛片一级片免费看久久久久| 成人亚洲精品一区在线观看| 亚洲成人手机| 欧美区成人在线视频| 国产亚洲精品久久久com| 成人毛片60女人毛片免费| 下体分泌物呈黄色| 国产黄色免费在线视频| 妹子高潮喷水视频| av福利片在线| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 日韩大片免费观看网站| 成人综合一区亚洲| 一级a做视频免费观看| 日韩成人伦理影院| 黄色怎么调成土黄色| 国产在线视频一区二区| 91精品国产国语对白视频| 色网站视频免费| 色婷婷久久久亚洲欧美| 亚洲精品亚洲一区二区| 亚洲av综合色区一区| 免费播放大片免费观看视频在线观看| 中文在线观看免费www的网站| 中文精品一卡2卡3卡4更新| 美女中出高潮动态图| 黑人猛操日本美女一级片| 午夜福利,免费看| 啦啦啦视频在线资源免费观看| 国产精品一区二区三区四区免费观看| 精品卡一卡二卡四卡免费| 中文字幕人妻丝袜制服| 黄色视频在线播放观看不卡| freevideosex欧美| 极品教师在线视频| 99久久精品热视频| 免费少妇av软件| 伦理电影免费视频| 久久鲁丝午夜福利片| 王馨瑶露胸无遮挡在线观看| 成年美女黄网站色视频大全免费 | 国产精品久久久久久精品电影小说|