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

    基于恒星矢量差分的星敏感器星跟蹤算法研究

    2017-04-28 01:25:26毛曉楠孫少勇胡雄超
    上海航天 2017年2期
    關(guān)鍵詞:星點(diǎn)卡爾曼濾波矢量

    金 荷,毛曉楠,孫少勇,劉 露,胡雄超

    (1.中國航天科技集團(tuán)公司 紅外探測技術(shù)研發(fā)中心,上海 201109; 2.上海航天技術(shù)控制研究所,上海 201109)

    ?

    基于恒星矢量差分的星敏感器星跟蹤算法研究

    金 荷1, 2,毛曉楠1, 2,孫少勇1, 2,劉 露1, 2,胡雄超1, 2

    (1.中國航天科技集團(tuán)公司 紅外探測技術(shù)研發(fā)中心,上海 201109; 2.上海航天技術(shù)控制研究所,上海 201109)

    針對(duì)現(xiàn)有星跟蹤算法依賴外部信息、模型復(fù)雜等缺點(diǎn),結(jié)合衛(wèi)星平臺(tái)實(shí)際,提出了一種新的星敏感器內(nèi)部星跟蹤算法,利用星矢量與角速率間關(guān)系進(jìn)行星點(diǎn)預(yù)測。先根據(jù)當(dāng)前幀和前一幀的星矢量,用最小二乘原理估算出角速率,再對(duì)角速率進(jìn)行卡爾曼濾波以減小估計(jì)誤差,之后用濾波后的角速率預(yù)測下一幀的星矢量,將星矢量轉(zhuǎn)為星點(diǎn)坐標(biāo),最后在以該坐標(biāo)為中心的波門中提取星點(diǎn)。討論了角速率估計(jì)誤差、截?cái)嗾`差、星點(diǎn)成像隨機(jī)誤差和定姿星數(shù)等因素對(duì)算法精度的影響。該方法不依賴外部信息進(jìn)行預(yù)測,模型簡單,易于實(shí)現(xiàn)。實(shí)驗(yàn)驗(yàn)證在星敏感器角速率為0.6 (°)/s時(shí),平均跟蹤精度0.55像素,證明了其有效性。

    星跟蹤; 星矢量差分; 角速率估計(jì); 星點(diǎn)預(yù)測; 最小二乘法; 卡爾曼濾波; 截?cái)嗾`差; 星點(diǎn)成像誤差

    0 引言

    星敏感器是當(dāng)前精度最高的姿態(tài)敏感器,在衛(wèi)星姿態(tài)確定中被廣泛應(yīng)用。在全天捕獲成功后,星敏感器主要工作于星跟蹤模式。星跟蹤算法中首先進(jìn)行星點(diǎn)預(yù)測,以預(yù)測的星點(diǎn)位置為中心形成一定大小的波門,圖像處理軟件在該波門中提取星點(diǎn)位置信息,顯著減少了圖像處理的時(shí)間,提高了算法效率。當(dāng)星點(diǎn)預(yù)測誤差大時(shí),波門必須相應(yīng)增大,一方面增加了計(jì)算量,另一方面由于波門內(nèi)包含虛假目標(biāo)的概率增大,星點(diǎn)提取算法誤識(shí)別的可能性也增大,因此提高星點(diǎn)預(yù)測精度十分重要。星點(diǎn)預(yù)測算法分為兩類。一類是利用外部信息的星點(diǎn)預(yù)測方法,包括姿態(tài)與星圖信息,此類方法的研究已較充分[1-9]。文獻(xiàn)[1]研究了飛行器大角速率機(jī)動(dòng)時(shí)的星跟蹤算法,利用當(dāng)前幀和前一幀的瞬時(shí)姿態(tài)預(yù)測下一幀的姿態(tài),根據(jù)預(yù)測的姿態(tài)和星表信息計(jì)算出下一幀的星點(diǎn)信息。文獻(xiàn)[2]利用當(dāng)前幀姿態(tài)和陀螺輸出的角加速度積分推算下一時(shí)刻的星敏姿態(tài),再用該姿態(tài)信息和星表信息計(jì)算出下一幀的星點(diǎn)信息。此算法除姿態(tài)信息外還依賴陀螺的角加速度信息。第二類算法不依賴外部信息進(jìn)行星點(diǎn)預(yù)測,關(guān)于該類算法的研究較少。文獻(xiàn)[10]對(duì)星點(diǎn)運(yùn)動(dòng)軌跡建模,在該模型的基礎(chǔ)上預(yù)測星點(diǎn)位置并通過卡爾曼濾波濾除星點(diǎn)成像瞬時(shí)噪聲。文獻(xiàn)[11]提出了基于singer模型的卡爾曼濾波星點(diǎn)提取方法,該方法與文獻(xiàn)[10]的類似。這類算法專注于對(duì)星點(diǎn)軌跡的建模,在自主星跟蹤的同時(shí)濾除星點(diǎn)成像瞬時(shí)噪聲,但有兩個(gè)不足。第一,在星點(diǎn)軌跡的建模中需有星點(diǎn)機(jī)動(dòng)頻率和隨機(jī)加速度方差等難以獲得的先驗(yàn)信息。事實(shí)上在探測器不同位置處,星軌跡的形態(tài)并不相同,星點(diǎn)加速度并不一樣,這在大視場星敏感器和多頭星敏感器中尤為明顯,因此難以對(duì)所有的星軌跡精確建模。第二,對(duì)每顆恒星在探測器上的軌跡需在x、y向分別建立三維運(yùn)動(dòng)模型,在該模型的基礎(chǔ)上進(jìn)行星點(diǎn)預(yù)測和濾波。實(shí)際應(yīng)用中平均每幀星圖要利用星點(diǎn)15顆進(jìn)行定姿,故需對(duì)15條星軌跡模型分別進(jìn)行遞推和濾波,計(jì)算量大。

    為能利用盡可能少的外部信息進(jìn)行星跟蹤,并改善第二類方法中星點(diǎn)軌跡跟蹤模型復(fù)雜的缺陷,根據(jù)衛(wèi)星平臺(tái)實(shí)際,本文提出了一種利用星敏感器三維星矢量的星敏感器星跟蹤方法。該方法先用當(dāng)前幀和上一幀的星矢量估算角速率,再對(duì)角速率進(jìn)行卡爾曼濾波以減小角速率估算噪聲,用濾波后的角速率及當(dāng)前幀的星矢量信息預(yù)測下一幀的星矢量信息,計(jì)算出下一幀的星點(diǎn)預(yù)測值,最后在以預(yù)測星點(diǎn)為中心的波門中提取星點(diǎn)。該法不依賴導(dǎo)航星表與姿態(tài)四元數(shù)先驗(yàn)信息,在星敏感器出現(xiàn)短暫的數(shù)據(jù)無效時(shí),仍可高精度地預(yù)測星點(diǎn)位置,穩(wěn)定地完成星跟蹤,模型簡單、易實(shí)現(xiàn)。

    1 原理分析

    本文算法步驟具體如下。

    a)第一步,角速率估算。根據(jù)星敏感器星矢量與角速率間的關(guān)系,先由當(dāng)前幀星點(diǎn)及星敏感器參數(shù)計(jì)算出當(dāng)前幀星矢量,用當(dāng)前幀和上一幀的星矢量估算角速率。該角速率實(shí)際是上一幀時(shí)間點(diǎn)至當(dāng)前幀時(shí)間點(diǎn)間的角速率,當(dāng)前幀與下一幀間的角速率可用該值近似。

    b)第二步,為減小角速率估算的噪聲,對(duì)角速率進(jìn)行卡爾曼濾波。

    c)第三步,利用濾波后的角速率及當(dāng)前幀的星矢量信息預(yù)測下一幀的星矢量信息,由星矢量信息和星敏感器內(nèi)方位元素計(jì)算出下一幀的星點(diǎn)位置,即為星點(diǎn)預(yù)測值。

    d)第四步,在以預(yù)測星點(diǎn)為中心的波門中提取星點(diǎn)。

    1.1 單星矢量差分計(jì)算角速率

    星敏感器測量模型可表示為

    (1)

    式中:si(k)為時(shí)刻k第i個(gè)星矢量在星敏系中的表示;A(k)為時(shí)刻k星敏感器姿態(tài)矩陣;ri為時(shí)刻k第i個(gè)星矢量在慣性系中的表示;vi為測量誤差。

    根據(jù)SHUSTER的研究,測量誤差是方差為σiI的白噪聲,ri在忽略章動(dòng)時(shí)保持不變[12]。此處:I為33維的單位矩陣。對(duì)式(1)作差分可得

    si(k+1)-si(k)=

    [A(k+1)-A(k)]ri+vi(k+1)-vi(k)

    (2)

    假設(shè)時(shí)刻k與時(shí)刻k+1間的角速率為恒定,并忽略高階項(xiàng),則式(2)可簡化為

    si(k+1)-si(k)=

    -Δt[ω(k)]×A(k)ri+vi(k+1)-vi(k)

    (3)

    式中:Δt為兩個(gè)采樣時(shí)刻的間隔;ω(k)為時(shí)刻k的角速率,且

    (4)

    這在采樣頻率足夠高的條件下是成立的,而且實(shí)際工程中空間飛行器的姿態(tài)不會(huì)突變,多數(shù)情況下很穩(wěn)定,故此種假設(shè)符合實(shí)際。

    為實(shí)現(xiàn)不利用姿態(tài)信息,將式(1)代入式(3)消去A(k),可得

    (1/Δt) [si(k+1)-si(k)]=

    [si(k)]×ω(k)+wi(k)

    (5)

    式中:wi(k)為vi的函數(shù),且

    (6)

    式(5)揭示了角速率與星矢量的關(guān)系,這不需要依賴姿態(tài)信息等外部信息作為中間變量。

    假設(shè)vi為固定的噪聲過程,則wi(k)的方差可表示為

    (7)

    式中:σi為星敏感器對(duì)第i個(gè)星矢量的測量誤差。

    如所有時(shí)刻中,姿態(tài)變化遠(yuǎn)小于Nyquist頻率,有‖ω(k)‖Δt?π/10,即采樣周期相對(duì)角速率來說足夠小,成立

    (8)

    因ω(k)的2范數(shù)與(ω(k))×的2范數(shù)相等,故式(7)右端第二項(xiàng)為主導(dǎo)項(xiàng)。

    1.2 多星矢量最小二乘角速率估計(jì)

    由最小二乘原理和式(5),可求出角速率的估值

    (9)

    (10)

    式中:(xi(k),yi(k))為時(shí)刻k第i顆星在星敏探測器上的成像點(diǎn)坐標(biāo);(x0,y0)為星敏感器主點(diǎn)坐標(biāo);f為星敏感器焦距。

    該角速率的估值為時(shí)刻k至?xí)r刻k+1(上一個(gè)時(shí)間點(diǎn)至當(dāng)前時(shí)間點(diǎn))的角速率的估值,星點(diǎn)預(yù)測中需預(yù)測下一幀的星矢量(時(shí)刻k+2的星矢量),因此理論上需求出時(shí)刻k+1至?xí)r刻k+2的估值。假設(shè)角速率相對(duì)采樣周期緩變,可用前一個(gè)采樣周期的角速率估值近似本幀的角速率估值。

    式(5)兩邊同時(shí)乘以Δt,將si(k+1)的結(jié)果代入式(9),可得

    (11)

    估計(jì)誤差方差

    (12)

    1.3 卡爾曼濾波角速率平滑及星點(diǎn)預(yù)測

    對(duì)模型和噪聲特性已知的隨機(jī)過程,卡爾曼濾波可對(duì)模型狀態(tài)進(jìn)行最小方差線性估計(jì)[13]。根據(jù)角速率的線性模型和上文對(duì)角速率估計(jì)噪聲的分析,可用卡爾曼濾波對(duì)角速率進(jìn)行最優(yōu)線性估計(jì),減小角速率的噪聲,尤其是大角速率下的角速率估計(jì)誤差。當(dāng)采樣頻率足夠大時(shí),假設(shè)前后幀角速率不變,卡爾曼濾波的狀態(tài)方程為

    (13)

    測量方程、狀態(tài)估計(jì)方差的一步預(yù)測公式、卡爾曼濾波增益公式、狀態(tài)最優(yōu)線性估計(jì)和狀態(tài)估計(jì)方差更新公式分別為

    (14)

    (15)

    (16)

    (17)

    (18)

    式中:Z(k+1)為測量值,由式(9)給出;P(k)為狀態(tài)估計(jì)方差;Q(k)為狀態(tài)轉(zhuǎn)移噪聲;Kg(k)為卡爾曼濾波增益;R(k)為測量噪聲;I為單位陣。理論上,u(k)是均值和方差為零的白噪聲,但u(k)的方差Q取0時(shí)濾波易發(fā)散,故本文取其為0.001I3×3,另取v(k)的方差0.5I3×3。

    最后,根據(jù)式(5),用估計(jì)的角速率及當(dāng)前幀的星矢量可估得下一幀的星矢量,由星矢量和星敏感器參數(shù)可預(yù)測下一幀的星點(diǎn)位置,在以該點(diǎn)為中心的波門中提取星點(diǎn)。即

    (19)

    式中:第i個(gè)星矢量在下一幀的估值si=[si(1)si(2)si(3)]T。

    2 影響算法精度的因素

    定義本文算法的精度為星跟蹤算法中的星點(diǎn)預(yù)測精度,星點(diǎn)預(yù)測精度越高,星點(diǎn)提取波門越小,星點(diǎn)提取效率就越高,星跟蹤算法性能越好。本算法的誤差流如圖1所示。

    表征星矢量與角速率關(guān)系的式(5)忽略了姿態(tài)變化‖ω(k)‖Δt的高階項(xiàng),由此產(chǎn)生截?cái)嗾`差,該誤差通過影響角速率估值間接影響算法精度,同時(shí)因星點(diǎn)預(yù)測使用式(5),該誤差也會(huì)影響算法精度。

    由式(12)可知:星點(diǎn)成像瞬時(shí)誤差影響角速率估計(jì)誤差,會(huì)間接影響算法精度,同時(shí)由于星點(diǎn)預(yù)測時(shí)使用了上一幀的星矢量,也會(huì)引入星點(diǎn)成像瞬時(shí)誤差,該誤差也會(huì)直接影響算法精度;定姿星數(shù)會(huì)影響角速率估計(jì)誤差,間接影響算法精度。

    2.1 角速率估計(jì)誤差對(duì)算法精度的影響

    根據(jù)某星敏感器外場觀星的實(shí)測數(shù)據(jù),在角速率分別為0.2,0.4,0.6 (°)/s條件下角速率的估計(jì)誤差見表1。由表1可知:角速率估計(jì)誤差隨角速率增大而增大,由式(12),這是因?yàn)樾屈c(diǎn)成像誤差變大,定姿星數(shù)變少。

    表1 不同角速率下角速率估計(jì)誤差Tab.1 Error of angle rate estimation under different angle rate

    分析角速率估計(jì)誤差引起的星點(diǎn)預(yù)測誤差,三軸角速率誤差對(duì)預(yù)測星點(diǎn)坐標(biāo)的影響如圖2所示。根據(jù)星敏感器的原理,x、y向是等價(jià)的。由圖2(a)可知:x軸角速率誤差引起的星點(diǎn)預(yù)測誤差

    (20)

    式中:E為星點(diǎn)預(yù)測誤差;Δω為角速率誤差;dt為采樣時(shí)間;A為星敏感器視場;N為焦平面在x向或y向的像素。

    z軸角速率誤差,即滾動(dòng)角誤差引起的星點(diǎn)預(yù)測誤差與星點(diǎn)至主點(diǎn)(焦平面與光軸交點(diǎn))的距離有關(guān)。文獻(xiàn)[14]給出的焦平面內(nèi)星點(diǎn)至主點(diǎn)的平均距離

    (21)

    由圖2(b)可知:z軸角速率誤差引起的星點(diǎn)預(yù)測誤差

    (22)

    設(shè)星敏感器的圖像分辨率1 024 pix×1 024 pix,數(shù)據(jù)更新頻率5 Hz。根據(jù)表1,該型星敏感器運(yùn)動(dòng)速率為0.6 (°)/s時(shí),x軸角速率估計(jì)誤差為0.099 4 (°)/s,由式(20)可得該誤差引起的星點(diǎn)預(yù)測誤差為1.183 4 pix;y軸角速率估計(jì)誤差為0.206 3 (°)/s,由式(20)可得該誤差引起的星點(diǎn)預(yù)測誤差為2.456 0 pix;z軸角速率估計(jì)誤差為0.343 8 (°)/s,由式(22)可得該誤差引起的星點(diǎn)預(yù)測誤差為0.47 pix。由此可知:為更準(zhǔn)確地預(yù)測星點(diǎn)位置,實(shí)現(xiàn)高效的星跟蹤,有必要對(duì)角速率估計(jì)誤差進(jìn)行濾波。

    2.2 截?cái)嗾`差對(duì)精度的影響

    此處,影響指截?cái)嗾`差對(duì)精度的綜合影響,是直接影響和間接影響的總和。仿真分析不同‖ω(k)‖Δt時(shí)星跟蹤算法的星點(diǎn)預(yù)測誤差,結(jié)果見表2。由表2可知:隨著‖ω(k)‖Δt增大,平均星點(diǎn)預(yù)測誤差亦增大,表明本文算法在星敏感器數(shù)據(jù)更新頻率足夠大(即衛(wèi)星姿態(tài)相對(duì)采樣周期為緩變量)時(shí)可保證精度。

    表2 截?cái)嗾`差對(duì)精度的影響Tab.2 Influence to accuracy by truncation error

    2.3 星點(diǎn)成像隨機(jī)誤差對(duì)精度的影響

    此處的影響是指星點(diǎn)成像隨機(jī)誤差誤差對(duì)精度的綜合影響,是直接影響與間接影響的總和,當(dāng)‖ω(k)‖Δt=0.12,平均定姿星數(shù)為15時(shí),仿真所得不同星點(diǎn)成像隨機(jī)誤差對(duì)精度的影響見表3。由表3可知:星點(diǎn)成像誤差越大,星跟蹤精度就越差。濾波對(duì)該誤差的抑制有一定的效果。

    表3 星點(diǎn)成像隨機(jī)誤差對(duì)精度的影響Tab.3 Influence to accuracy by temporal error of star point imaging

    2.4 定姿星數(shù)對(duì)精度的影響

    當(dāng)‖ω‖Δt=0.12,星點(diǎn)成像隨機(jī)誤差為0.5pix時(shí),仿真分析用不同的定姿星數(shù)估計(jì)角速率和預(yù)測星點(diǎn)時(shí)對(duì)精度的影響,結(jié)果見表4。由表4可知:平均星點(diǎn)預(yù)測誤差隨定姿星數(shù)的減少而增加,角速率濾波對(duì)因定姿星減小而增大的誤差有較好的抑制作用。當(dāng)定姿星僅為2顆時(shí),星敏不能定姿,但仍用該算法進(jìn)行星跟蹤和姿態(tài)遞推。

    表4 定姿星數(shù)對(duì)精度的影響Tab.4 Influence to accuracy by number of stars used for attitude determination

    3 實(shí)驗(yàn)結(jié)果及分析

    3.1 角速率估計(jì)

    用某星敏感器的外場觀星試驗(yàn)對(duì)本文星跟蹤算法進(jìn)行驗(yàn)證,采集星敏感器在0.6 (°)/s條件下運(yùn)動(dòng)時(shí)的實(shí)測數(shù)據(jù),用本文方法預(yù)測星點(diǎn),實(shí)現(xiàn)星跟蹤,同時(shí)驗(yàn)證該方法的可行性和精度。實(shí)驗(yàn)所得星點(diǎn)運(yùn)動(dòng)軌跡如圖3所示。由傳統(tǒng)的用四元數(shù)差分所得角速率估計(jì)和本文算法使用前后幀星矢量估計(jì)角速率如圖4所示。

    由圖4可知:實(shí)驗(yàn)證明雖然未使用姿態(tài)信息,星矢量估計(jì)角速率的精度與用四元數(shù)差分估計(jì)的精度相當(dāng),估計(jì)角速率均值為0.603 3 (°)/s。

    3.2 角速率濾波及星點(diǎn)預(yù)測

    對(duì)角速率估值進(jìn)行卡爾曼濾波抑制噪聲,用濾波后的角速率預(yù)測星點(diǎn)位置。所得濾波前后的角速率誤差見表5,角速率如圖5所示。由表5和圖5可知:濾波對(duì)角速率估計(jì)噪聲有較好的抑制作用。

    表5 濾波前后角速率誤差Tab.5 Angle rate error before and after filtering

    星敏感器某幀的星點(diǎn)預(yù)測結(jié)果見表6。由表6可知:該預(yù)測方法基本可準(zhǔn)確預(yù)測出星點(diǎn)位置。角速率濾波前后的星點(diǎn)誤差見表7,星點(diǎn)誤差預(yù)測值如圖6所示。實(shí)驗(yàn)中采樣點(diǎn)212個(gè),其中濾波效果優(yōu)于不濾波的采樣點(diǎn)160個(gè),占75%;濾波前星點(diǎn)預(yù)測誤差大于1的采樣點(diǎn)39個(gè),其中濾波后預(yù)測誤差減小的38個(gè),占97%,平均預(yù)測誤差減小0.56 pix;

    表6 星敏感器某幀的預(yù)測星點(diǎn)坐標(biāo)Tab.6 Predicted coordinates of star points

    表7 濾波前后星點(diǎn)誤差Tab.7 Star points prediction error before and after filtering

    濾波前星點(diǎn)預(yù)測誤差大于0.5的采樣點(diǎn)122個(gè),其中濾波后預(yù)測誤差減小的106個(gè),占87%,平均預(yù)測誤差減小0.33 pix;角速率濾波后星點(diǎn)預(yù)測誤差峰值減小1 pix。這表明:卡爾曼濾波在角速率計(jì)算誤差較大時(shí)有較好的修正作用,特別是對(duì)誤差尖峰有很好的平滑作用。

    4 結(jié)束語

    本文提出了一種僅基于星敏感器坐標(biāo)系中星矢量坐標(biāo)計(jì)算角速率并實(shí)現(xiàn)星跟蹤的算法。該算法基于姿態(tài)和角速率相對(duì)采樣周期緩變的假設(shè),不依賴于姿態(tài)信息、星表或陀螺的角速率信息等外部信息,僅利用星點(diǎn)信息及星敏感器內(nèi)部參數(shù)就能實(shí)現(xiàn)星跟蹤。與其它不依賴于姿態(tài)信息的星跟蹤算法相比,其優(yōu)點(diǎn)是模型簡單、實(shí)現(xiàn)易、計(jì)算量小。實(shí)驗(yàn)驗(yàn)證在星敏感器角速率為0.6 (°)/s時(shí),平均跟蹤精度0.55 pix,證明了其有效性。該算法可作為現(xiàn)有星跟蹤算法的有效補(bǔ)充,當(dāng)星敏感器姿態(tài)無效時(shí)仍可進(jìn)行星點(diǎn)預(yù)測,實(shí)現(xiàn)星跟蹤,并能利用最后一幀有效姿態(tài)與角速率估值進(jìn)行姿態(tài)遞推,為下一次姿態(tài)捕獲成功提供有利條件。由于不同角速率時(shí)的角速率估計(jì)誤差各異,在角速率卡爾曼濾波中可根據(jù)不同的角速率水平調(diào)整噪聲估計(jì)方差,進(jìn)行自適應(yīng)卡爾曼濾波,進(jìn)一步提高角速率估計(jì)精度,從而提高星跟蹤精度。

    [1] 李葆華, 陳希軍, 鄭循江, 等. 星敏感器高動(dòng)態(tài)下自主星跟蹤算法[J]. 紅外與激光工程, 2012, 41(1): 191-195.

    [2] 王常虹, 吳志華, 李葆華. 一種大角加速度下的星跟蹤算法[J]. 中國慣性技術(shù)學(xué)報(bào), 2010, 18(6): 707-709.[3] 蔣明, 來逢昌, 葉以正, 等. 適用于星敏感器的星跟蹤算法研究[J]. 電子器件, 2009, 32(2): 426-428.

    [4] 江潔, 張廣軍, 李霄, 等. 一種快速的星敏感器星跟蹤方法研究[J]. 航空學(xué)報(bào), 2006, 27(5): 914-916.

    [5] 管宇, 張迎春, 王常虹. 一種適用于星敏感器的大角速度下的星跟蹤算法[J]. 中國慣性技術(shù)學(xué)報(bào), 2012, 20(4): 441-448.

    [6] 胡雄超, 吳永康, 閆曉軍, 等. 一種星敏感器標(biāo)片誤差修正方法研究[J]. 上海航天, 2016, 33(6): 88-92.

    [7] 余路偉, 毛曉楠, 金尚, 等. 星敏感器抗雜光圖像處理方法[J]. 上海航天, 2016, 33(4): 26-31.

    [8] LIEBE C C. Accuracy performance of star trackers——a tutorial[J]. IEEE Transactions on Aerospace and Electronic System, 2002, 38(2): 587-599.

    [9] 何家雄. 高精度全天時(shí)星敏感器關(guān)鍵技術(shù)研究[D]. 北京: 中國科學(xué)院研究生院, 2013.

    [10] 金雁, 江潔, 張廣軍. 高動(dòng)態(tài)星跟蹤方法[J]. 紅外與激光工程, 2013, 42(1): 213-217.

    [11] 劉海波. 基于星敏感器的無陀螺角速率測量新方法研究[D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 2013.

    [12] SHUSTER M D. Maximum likelihood estimation of spacecraft attitude[J]. Journal of the Astronautical Sciences, 189, 37(1): 79-88.

    [13] 秦永元, 張洪鉞, 汪叔華. 卡爾曼濾波與組合導(dǎo)航原理[M]. 西安: 西北工業(yè)大學(xué)出版社, 2012.

    [14] LIEBE C C. Accuracy performance of star trackers[J]. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(2): 587-599.

    Research on Star Tracking Algorithm Based on Star Vector Difference

    JIN He1, 2, MAO Xiao-nan1, 2, SUN Shao-yong1, 2, LIU Lu1, 2, HU Xiong-chao1, 2

    (1. Infrared Detection Technology Research & Development Center, China Aerospace Science and Technology Corporation, Shanghai 201109, China;2. Shanghai Institute of Spaceflight Control Technology, Shanghai 201109, China)

    According to the disadvantages of conventional algorithm which were dependent of outside information of attitude and complex model, an inner star prediction algorithm based on the relationship between star vector and angular rate was put forward considering the working condition of satellite platform in this paper. Firstly, angular rate was estimated according to the difference between current and last star vector based on the principle of least squares. Secondly, Kalman filtering was used to reduce the error of angular rate estimation in next sample period. Then, the filtered angular rate was applied to predict next star vectors which would be transformed to star positions. At last, center of sub window was applied for star extraction. The influence of the error of angular velocity estimation, truncation error, random error of star point imaging and star numbers for attitude determination on the accuracy of the algorithm was discussed. This method would make prediction in independent of attitude with simple model and was easy to be implemented. The real sky test had validated that this method had average tracking accuracy of 0.55 pixel at angular rate of 0.6 (°)/s.

    star tracking; star vector difference; angular rate estimation; star point prediction; least squares; Kalman filter; truncation error; star point imaging error

    1006-1630(2017)02-0099-07

    2016-11-10;

    2017-01-22

    國家重大科技專項(xiàng)工程(GFZX0401300206)

    金 荷(1992—),女,碩士,主要從事星敏感器誤差抑制研究。

    V448

    A

    10.19328/j.cnki.1006-1630.2017.02.010

    猜你喜歡
    星點(diǎn)卡爾曼濾波矢量
    矢量三角形法的應(yīng)用
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    基于矢量最優(yōu)估計(jì)的穩(wěn)健測向方法
    星點(diǎn)設(shè)計(jì)-效應(yīng)面法優(yōu)化雄黃乳膏劑的處方組成
    中成藥(2016年4期)2016-05-17 06:08:05
    一種基于數(shù)學(xué)形態(tài)學(xué)的星點(diǎn)提取方法
    三角形法則在動(dòng)態(tài)平衡問題中的應(yīng)用
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    基于擴(kuò)展卡爾曼濾波的PMSM無位置傳感器控制
    星點(diǎn)設(shè)計(jì)-效應(yīng)面法優(yōu)選止鼾顆粒成型工藝
    中成藥(2014年9期)2014-02-28 22:28:51
    基于自適應(yīng)卡爾曼濾波的新船舶試航系統(tǒng)
    欧美日韩精品成人综合77777| 欧美变态另类bdsm刘玥| 18禁在线无遮挡免费观看视频| 亚洲国产av新网站| 日韩,欧美,国产一区二区三区| 午夜免费鲁丝| 欧美区成人在线视频| 熟妇人妻不卡中文字幕| 18禁在线播放成人免费| 亚洲熟女精品中文字幕| 国产精品一区二区三区四区免费观看| 女的被弄到高潮叫床怎么办| 嫩草影院新地址| 精品人妻视频免费看| 最近中文字幕2019免费版| 亚洲第一av免费看| 国产精品久久久久久精品电影小说 | 色视频www国产| 激情 狠狠 欧美| 久久99蜜桃精品久久| 少妇高潮的动态图| av黄色大香蕉| av不卡在线播放| 激情五月婷婷亚洲| 在线观看免费日韩欧美大片 | 日韩强制内射视频| 久久久久精品久久久久真实原创| 国产高清不卡午夜福利| 国产色爽女视频免费观看| 亚洲精品乱久久久久久| 性高湖久久久久久久久免费观看| 一区二区三区免费毛片| 纵有疾风起免费观看全集完整版| 精品一品国产午夜福利视频| 国产中年淑女户外野战色| 国产爱豆传媒在线观看| 久久久久久久久久久丰满| 最后的刺客免费高清国语| 精品久久国产蜜桃| 天堂中文最新版在线下载| 国产精品久久久久久精品古装| 亚洲精品久久午夜乱码| 久久热精品热| 国产毛片在线视频| 九九爱精品视频在线观看| 99久久精品国产国产毛片| 国产精品一区二区在线不卡| 免费不卡的大黄色大毛片视频在线观看| 在现免费观看毛片| av播播在线观看一区| 亚洲国产成人一精品久久久| 深爱激情五月婷婷| 在线看a的网站| 一级爰片在线观看| 在现免费观看毛片| 各种免费的搞黄视频| 边亲边吃奶的免费视频| 国产成人freesex在线| 极品少妇高潮喷水抽搐| 亚洲欧美一区二区三区国产| 亚洲,一卡二卡三卡| 国产亚洲91精品色在线| 久久99蜜桃精品久久| 国产精品久久久久久久电影| 日韩视频在线欧美| 一区二区三区乱码不卡18| 日韩,欧美,国产一区二区三区| 亚洲精品国产av成人精品| 欧美日韩在线观看h| 亚洲无线观看免费| 国产在线一区二区三区精| 黄色视频在线播放观看不卡| 亚洲aⅴ乱码一区二区在线播放| 久久久精品94久久精品| 九九久久精品国产亚洲av麻豆| av在线观看视频网站免费| 欧美人与善性xxx| 久久99热这里只频精品6学生| 自拍欧美九色日韩亚洲蝌蚪91 | 一级毛片我不卡| 一级a做视频免费观看| 中文字幕av成人在线电影| 三级国产精品欧美在线观看| 日韩成人伦理影院| 在线 av 中文字幕| 色哟哟·www| 国产精品熟女久久久久浪| 国产白丝娇喘喷水9色精品| 内地一区二区视频在线| 久久人人爽人人片av| 一区二区三区乱码不卡18| 天堂8中文在线网| 18禁裸乳无遮挡动漫免费视频| 观看av在线不卡| 纵有疾风起免费观看全集完整版| 免费av不卡在线播放| 国产精品国产av在线观看| 久久精品国产亚洲av天美| 99久国产av精品国产电影| 在线观看美女被高潮喷水网站| 黑人高潮一二区| 亚洲成人av在线免费| 国产精品99久久久久久久久| 亚洲av在线观看美女高潮| 亚洲精品成人av观看孕妇| 国内精品宾馆在线| 丰满少妇做爰视频| 亚洲精品乱久久久久久| 男女边摸边吃奶| 午夜福利网站1000一区二区三区| 又爽又黄a免费视频| 人妻系列 视频| 国内精品宾馆在线| 夜夜骑夜夜射夜夜干| 久久午夜福利片| 国产黄片视频在线免费观看| 日本黄色日本黄色录像| 少妇精品久久久久久久| 国产视频内射| 成年免费大片在线观看| 狂野欧美激情性bbbbbb| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品久久久久久婷婷小说| 亚洲,欧美,日韩| 免费人成在线观看视频色| 美女视频免费永久观看网站| 久久午夜福利片| 夫妻午夜视频| 精品久久久久久久久亚洲| 在线观看av片永久免费下载| .国产精品久久| 日韩国内少妇激情av| 国产淫语在线视频| 国产日韩欧美亚洲二区| 亚洲国产成人一精品久久久| 亚洲成人一二三区av| 久久久久国产精品人妻一区二区| 国内精品宾馆在线| 欧美成人一区二区免费高清观看| 久久精品久久精品一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品自拍成人| 熟女电影av网| 小蜜桃在线观看免费完整版高清| 内地一区二区视频在线| 一级av片app| 乱系列少妇在线播放| 简卡轻食公司| 亚洲欧美一区二区三区国产| 国产成人aa在线观看| 97在线视频观看| 亚州av有码| 精品国产三级普通话版| 亚洲av国产av综合av卡| 寂寞人妻少妇视频99o| 欧美最新免费一区二区三区| 免费人成在线观看视频色| 国产在线一区二区三区精| 中文字幕制服av| 中文字幕精品免费在线观看视频 | av在线蜜桃| 日本黄色日本黄色录像| 国产探花极品一区二区| 亚洲性久久影院| 在线播放无遮挡| 91狼人影院| 亚洲最大成人中文| 伦理电影大哥的女人| 国产男女超爽视频在线观看| av天堂中文字幕网| 熟女电影av网| 三级国产精品欧美在线观看| 欧美精品一区二区大全| 久久久久久九九精品二区国产| 亚洲激情五月婷婷啪啪| av在线观看视频网站免费| av在线观看视频网站免费| 91aial.com中文字幕在线观看| av国产免费在线观看| 美女主播在线视频| 最新中文字幕久久久久| 日本wwww免费看| 久久久久久久久大av| 精品一区在线观看国产| 午夜激情久久久久久久| 99热国产这里只有精品6| 日韩中字成人| 爱豆传媒免费全集在线观看| 亚洲天堂av无毛| 男人舔奶头视频| 国产午夜精品久久久久久一区二区三区| 内射极品少妇av片p| 国产深夜福利视频在线观看| 国产淫片久久久久久久久| 国内精品宾馆在线| 女人十人毛片免费观看3o分钟| 王馨瑶露胸无遮挡在线观看| 精品少妇久久久久久888优播| 亚洲美女黄色视频免费看| 欧美日韩综合久久久久久| 亚洲无线观看免费| 国产 一区 欧美 日韩| 亚洲国产最新在线播放| 亚洲精品国产成人久久av| 最新中文字幕久久久久| 国产精品一及| 最新中文字幕久久久久| 精品亚洲乱码少妇综合久久| av不卡在线播放| 91精品伊人久久大香线蕉| 亚洲精品国产色婷婷电影| 亚洲欧美日韩无卡精品| 色综合色国产| 91精品一卡2卡3卡4卡| 亚洲av成人精品一区久久| 国产v大片淫在线免费观看| 国内精品宾馆在线| 高清不卡的av网站| 亚洲天堂av无毛| 边亲边吃奶的免费视频| 欧美最新免费一区二区三区| 免费看光身美女| 亚洲人成网站在线观看播放| 国产精品麻豆人妻色哟哟久久| 欧美97在线视频| 少妇精品久久久久久久| 亚洲国产成人一精品久久久| 亚洲人与动物交配视频| 国产成人aa在线观看| 超碰av人人做人人爽久久| 亚洲国产欧美人成| 三级国产精品欧美在线观看| av网站免费在线观看视频| 成年女人在线观看亚洲视频| 精品亚洲乱码少妇综合久久| 特大巨黑吊av在线直播| 欧美精品亚洲一区二区| 嫩草影院入口| 国产精品久久久久久av不卡| 欧美少妇被猛烈插入视频| 99久久精品热视频| av播播在线观看一区| av.在线天堂| 少妇高潮的动态图| 少妇的逼水好多| 在线观看一区二区三区激情| 国产色爽女视频免费观看| 中文字幕人妻熟人妻熟丝袜美| 一级av片app| 一二三四中文在线观看免费高清| 一级毛片电影观看| 亚洲精品日韩在线中文字幕| 久久这里有精品视频免费| 国产在线一区二区三区精| 性色av一级| 日韩 亚洲 欧美在线| 久久久久网色| 高清欧美精品videossex| 精品久久久久久电影网| 美女主播在线视频| 中文字幕av成人在线电影| 日本黄大片高清| 日本黄大片高清| 亚洲欧美一区二区三区黑人 | 老师上课跳d突然被开到最大视频| 成年女人在线观看亚洲视频| 国产乱人偷精品视频| 国产在线男女| 中文字幕免费在线视频6| 熟女人妻精品中文字幕| 我要看黄色一级片免费的| 性色avwww在线观看| 日韩欧美精品免费久久| 王馨瑶露胸无遮挡在线观看| av国产精品久久久久影院| 一区二区三区免费毛片| 国产国拍精品亚洲av在线观看| 国产精品国产三级国产av玫瑰| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品爽爽va在线观看网站| 一级毛片久久久久久久久女| 91久久精品国产一区二区三区| 国产精品欧美亚洲77777| 国产 一区精品| 亚洲熟女精品中文字幕| 亚洲av中文av极速乱| 亚洲精品久久久久久婷婷小说| 亚洲激情五月婷婷啪啪| 午夜福利高清视频| 天天躁夜夜躁狠狠久久av| 亚洲激情五月婷婷啪啪| 日产精品乱码卡一卡2卡三| 亚洲精品456在线播放app| 六月丁香七月| 国产免费视频播放在线视频| 一级毛片久久久久久久久女| 我要看日韩黄色一级片| 黑人猛操日本美女一级片| 国产亚洲午夜精品一区二区久久| 欧美少妇被猛烈插入视频| 免费久久久久久久精品成人欧美视频 | 亚洲av中文av极速乱| 精品国产露脸久久av麻豆| 人人妻人人澡人人爽人人夜夜| 国产有黄有色有爽视频| 黄色日韩在线| 日韩视频在线欧美| 国产又色又爽无遮挡免| 成人美女网站在线观看视频| 欧美一级a爱片免费观看看| 国产免费视频播放在线视频| 91精品国产国语对白视频| 亚洲,一卡二卡三卡| 国产精品女同一区二区软件| 久久综合国产亚洲精品| 最近中文字幕高清免费大全6| 啦啦啦中文免费视频观看日本| 人人妻人人添人人爽欧美一区卜 | 欧美激情极品国产一区二区三区 | 夜夜爽夜夜爽视频| 国产精品麻豆人妻色哟哟久久| 中文在线观看免费www的网站| 九九爱精品视频在线观看| 免费人成在线观看视频色| 亚洲欧美成人精品一区二区| 日韩国内少妇激情av| 波野结衣二区三区在线| 国产女主播在线喷水免费视频网站| 1000部很黄的大片| 亚洲精品自拍成人| 国产真实伦视频高清在线观看| 国产成人精品婷婷| 国产精品爽爽va在线观看网站| 偷拍熟女少妇极品色| 国精品久久久久久国模美| 一级毛片我不卡| 久久久久精品性色| 成年av动漫网址| 一级毛片aaaaaa免费看小| 久久久欧美国产精品| 欧美区成人在线视频| 婷婷色综合大香蕉| 男女边摸边吃奶| 日韩av在线免费看完整版不卡| 一区二区av电影网| 亚洲伊人久久精品综合| 国产精品国产三级国产av玫瑰| 男女国产视频网站| 丰满少妇做爰视频| 亚洲国产精品国产精品| kizo精华| 午夜福利在线在线| 菩萨蛮人人尽说江南好唐韦庄| 黑人猛操日本美女一级片| 天天躁日日操中文字幕| 亚洲成人手机| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 久久这里有精品视频免费| 天堂中文最新版在线下载| 国产 一区精品| 十分钟在线观看高清视频www | 韩国av在线不卡| 五月天丁香电影| 99热全是精品| 亚洲人成网站在线观看播放| 国产成人精品福利久久| av在线蜜桃| 自拍欧美九色日韩亚洲蝌蚪91 | 黑丝袜美女国产一区| 国产 一区精品| 爱豆传媒免费全集在线观看| 亚洲av电影在线观看一区二区三区| 国产精品一区www在线观看| 免费黄色在线免费观看| 午夜福利视频精品| 秋霞在线观看毛片| 高清在线视频一区二区三区| 国产伦理片在线播放av一区| 极品少妇高潮喷水抽搐| 国产成人a区在线观看| 不卡视频在线观看欧美| 国产一区二区三区av在线| 97超视频在线观看视频| 18禁在线播放成人免费| 干丝袜人妻中文字幕| 一级毛片aaaaaa免费看小| 久久久亚洲精品成人影院| 午夜福利网站1000一区二区三区| 五月伊人婷婷丁香| 国产亚洲午夜精品一区二区久久| 国国产精品蜜臀av免费| 午夜激情福利司机影院| 亚洲精品国产av蜜桃| 日韩一区二区视频免费看| 男的添女的下面高潮视频| 亚洲怡红院男人天堂| 精品久久久久久久末码| 国产精品国产三级国产专区5o| 久久久a久久爽久久v久久| 肉色欧美久久久久久久蜜桃| 深夜a级毛片| av在线播放精品| 老司机影院成人| 免费黄网站久久成人精品| av天堂中文字幕网| 久久国产乱子免费精品| av国产精品久久久久影院| 97精品久久久久久久久久精品| 纯流量卡能插随身wifi吗| 女人久久www免费人成看片| 日韩欧美 国产精品| 国产精品免费大片| 99久久精品国产国产毛片| 秋霞在线观看毛片| 亚洲国产日韩一区二区| 永久免费av网站大全| 多毛熟女@视频| 综合色丁香网| 最近的中文字幕免费完整| 另类亚洲欧美激情| 亚洲av日韩在线播放| 国国产精品蜜臀av免费| 亚洲四区av| 在线 av 中文字幕| 国产av一区二区精品久久 | 久久久久久久大尺度免费视频| 网址你懂的国产日韩在线| 日本av免费视频播放| 91久久精品国产一区二区成人| 久久精品国产亚洲网站| 亚洲国产毛片av蜜桃av| 亚洲精品aⅴ在线观看| 久久久久久九九精品二区国产| 久久国产乱子免费精品| 精品一品国产午夜福利视频| 亚洲精品成人av观看孕妇| 在线精品无人区一区二区三 | 街头女战士在线观看网站| 亚洲国产高清在线一区二区三| 国产男女超爽视频在线观看| 一区二区av电影网| 国内精品宾馆在线| 嘟嘟电影网在线观看| 久久久久久久久久久免费av| 免费播放大片免费观看视频在线观看| 久久午夜福利片| 国产成人精品一,二区| 午夜福利网站1000一区二区三区| 天堂中文最新版在线下载| 熟妇人妻不卡中文字幕| 日韩成人伦理影院| 久久精品人妻少妇| 少妇猛男粗大的猛烈进出视频| 国产v大片淫在线免费观看| videos熟女内射| 中文天堂在线官网| 久久久国产一区二区| 免费看av在线观看网站| a级一级毛片免费在线观看| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 在线免费十八禁| 久久久a久久爽久久v久久| 全区人妻精品视频| 国产精品爽爽va在线观看网站| 精品人妻熟女av久视频| 99国产精品免费福利视频| 国产精品国产三级专区第一集| 免费播放大片免费观看视频在线观看| 亚洲四区av| 精品亚洲成国产av| 亚洲激情五月婷婷啪啪| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频 | 国语对白做爰xxxⅹ性视频网站| 日本欧美视频一区| 国产精品久久久久久久电影| 国产精品av视频在线免费观看| 在线 av 中文字幕| 久热这里只有精品99| 美女主播在线视频| 大片免费播放器 马上看| 成人午夜精彩视频在线观看| h日本视频在线播放| 亚洲欧美成人综合另类久久久| 国产精品久久久久久精品电影小说 | 亚洲av福利一区| 在线观看免费高清a一片| 亚洲av中文字字幕乱码综合| 22中文网久久字幕| 天堂8中文在线网| 国产视频内射| 日本-黄色视频高清免费观看| 一本—道久久a久久精品蜜桃钙片| 欧美一区二区亚洲| 欧美丝袜亚洲另类| 男女无遮挡免费网站观看| 人人妻人人澡人人爽人人夜夜| 有码 亚洲区| 人妻系列 视频| 九九在线视频观看精品| 免费久久久久久久精品成人欧美视频 | 啦啦啦中文免费视频观看日本| 青春草视频在线免费观看| 男女免费视频国产| 99久久中文字幕三级久久日本| 美女高潮的动态| 国产毛片在线视频| 精品久久国产蜜桃| 在线免费观看不下载黄p国产| 黄色日韩在线| 菩萨蛮人人尽说江南好唐韦庄| 夜夜骑夜夜射夜夜干| 一本久久精品| 最近最新中文字幕大全电影3| 成人毛片a级毛片在线播放| 卡戴珊不雅视频在线播放| 国产亚洲av片在线观看秒播厂| 亚洲美女黄色视频免费看| 久久韩国三级中文字幕| 亚洲无线观看免费| 精品国产乱码久久久久久小说| 欧美3d第一页| 三级经典国产精品| 久久人人爽av亚洲精品天堂 | 国产 精品1| 免费av不卡在线播放| 九色成人免费人妻av| 直男gayav资源| 精品熟女少妇av免费看| 91久久精品国产一区二区成人| videos熟女内射| 欧美bdsm另类| 国产在视频线精品| 亚洲在久久综合| 日韩制服骚丝袜av| 欧美xxxx黑人xx丫x性爽| 国产精品蜜桃在线观看| 亚洲av成人精品一区久久| 两个人的视频大全免费| 成人18禁高潮啪啪吃奶动态图 | 国产欧美亚洲国产| 最近最新中文字幕免费大全7| 欧美精品亚洲一区二区| av视频免费观看在线观看| 只有这里有精品99| 日韩成人av中文字幕在线观看| 久久精品国产亚洲av天美| 自拍偷自拍亚洲精品老妇| av卡一久久| 在线精品无人区一区二区三 | 男女下面进入的视频免费午夜| 亚洲精品第二区| 久久久久久久久久久免费av| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 亚洲无线观看免费| 五月伊人婷婷丁香| 我要看黄色一级片免费的| 亚洲av二区三区四区| 亚洲国产av新网站| 在线亚洲精品国产二区图片欧美 | 亚洲四区av| 国产欧美日韩精品一区二区| 久久99热这里只有精品18| 免费看日本二区| 最近中文字幕2019免费版| 婷婷色综合www| 精品少妇黑人巨大在线播放| 久久青草综合色| 99热全是精品| 日韩av在线免费看完整版不卡| 一级毛片电影观看| 国产高清三级在线| 欧美xxⅹ黑人| 亚洲高清免费不卡视频| 人妻 亚洲 视频| 爱豆传媒免费全集在线观看| 中文字幕精品免费在线观看视频 | 午夜日本视频在线| 青春草亚洲视频在线观看| 亚洲不卡免费看| 三级经典国产精品| 一区二区三区四区激情视频| videos熟女内射| 搡女人真爽免费视频火全软件| 国产精品久久久久成人av| 国内少妇人妻偷人精品xxx网站| 欧美少妇被猛烈插入视频| 直男gayav资源| 亚洲怡红院男人天堂| 久久久久国产精品人妻一区二区| 欧美精品人与动牲交sv欧美| 亚洲高清免费不卡视频| 久久久久久久久久久免费av| av一本久久久久| 午夜免费男女啪啪视频观看| 免费观看av网站的网址| 中文字幕精品免费在线观看视频 | 七月丁香在线播放| 国产片特级美女逼逼视频| 亚洲在久久综合| 各种免费的搞黄视频| 少妇精品久久久久久久| 女性生殖器流出的白浆| 男男h啪啪无遮挡| av视频免费观看在线观看| 国产永久视频网站| av卡一久久| 一个人看的www免费观看视频| 成人影院久久| a级一级毛片免费在线观看|