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

    基于遞推最小二乘法與模糊自適應(yīng)擴(kuò)展卡爾曼濾波相結(jié)合的車(chē)輛狀態(tài)估計(jì)

    2017-04-07 09:08:01魏民祥趙萬(wàn)忠張鳳嬌嚴(yán)明月
    中國(guó)機(jī)械工程 2017年6期
    關(guān)鍵詞:估計(jì)值協(xié)方差卡爾曼濾波

    汪 魏民祥 趙萬(wàn)忠 張鳳嬌,2 嚴(yán)明月

    1.南京航空航天大學(xué)能源與動(dòng)力學(xué)院,南京,2100162.常州工學(xué)院機(jī)械與車(chē)輛工程學(xué)院,常州,213002

    基于遞推最小二乘法與模糊自適應(yīng)擴(kuò)展卡爾曼濾波相結(jié)合的車(chē)輛狀態(tài)估計(jì)

    1.南京航空航天大學(xué)能源與動(dòng)力學(xué)院,南京,2100162.常州工學(xué)院機(jī)械與車(chē)輛工程學(xué)院,常州,213002

    針對(duì)汽車(chē)狀態(tài)估計(jì)中模型參數(shù)的變化和觀測(cè)噪聲的時(shí)變特性,提出了遞推最小二乘法與模糊自適應(yīng)擴(kuò)展卡爾曼濾波相結(jié)合的汽車(chē)狀態(tài)估計(jì)算法。為實(shí)現(xiàn)模型參數(shù)與觀測(cè)噪聲的實(shí)時(shí)更新,建立了基于三自由度非線性車(chē)輛動(dòng)力學(xué)模型的算法,首先利用遞推最小二乘法對(duì)汽車(chē)的總質(zhì)量進(jìn)行估計(jì),其次建立了模糊控制器對(duì)擴(kuò)展卡爾曼濾波的觀測(cè)噪聲進(jìn)行實(shí)時(shí)跟蹤。在搭建的CarSim與MATLAB/Simulink聯(lián)合仿真平臺(tái)中驗(yàn)證了該算法的有效性,結(jié)果表明該算法估計(jì)精度高于傳統(tǒng)擴(kuò)展卡爾曼濾波算法,研究結(jié)果為汽車(chē)的主動(dòng)安全控制提供了理論支持。

    汽車(chē)總質(zhì)量估計(jì);狀態(tài)估計(jì);遞推最小二乘法;模糊自適應(yīng)擴(kuò)展卡爾曼濾波

    0 引言

    隨著汽車(chē)工業(yè)的飛速發(fā)展,汽車(chē)主動(dòng)安全技術(shù)越來(lái)越受到人們的關(guān)注。對(duì)汽車(chē)主動(dòng)安全的研究幾乎都需要獲知汽車(chē)的當(dāng)前行駛狀態(tài),例如汽車(chē)的當(dāng)前車(chē)速、側(cè)向加速度、橫擺角速度和質(zhì)心側(cè)偏角等重要參數(shù)。對(duì)于上述關(guān)鍵狀態(tài)參數(shù)的獲取,目前研究中用得較多的方法是軟觀測(cè),即基于算法的估計(jì)[1]。

    目前汽車(chē)狀態(tài)估計(jì)器設(shè)計(jì)所采用的方法主要有擴(kuò)展卡爾曼濾波(EKF)[2-3]、無(wú)跡卡爾曼濾波[4-5]、神經(jīng)網(wǎng)絡(luò)[6-7]、模糊邏輯[8]等,它們都是對(duì)汽車(chē)控制系統(tǒng)中的關(guān)鍵控制變量(包括質(zhì)心側(cè)偏角、側(cè)向速度、橫擺角速度等)進(jìn)行估計(jì)。

    但這些方法都是基于模型參數(shù)(汽車(chē)總質(zhì)量)固定或算法參數(shù)(觀測(cè)噪聲協(xié)方差)固定來(lái)進(jìn)行狀態(tài)估計(jì)的,若參數(shù)改變將會(huì)導(dǎo)致結(jié)果精度降低,甚至可能會(huì)導(dǎo)致濾波發(fā)散。

    為了提高基于非線性汽車(chē)動(dòng)力學(xué)模型狀態(tài)估計(jì)算法的魯棒性與估計(jì)精度,本文提出一種遞推最小二乘(recursive least squares,RLS)法與模糊自適應(yīng)卡爾曼濾波(fuzzy adaptive extended Kalman filter,FAEKF)相結(jié)合的算法。

    1 非線性汽車(chē)動(dòng)力學(xué)模型

    本文估計(jì)算法中將車(chē)輛視為多個(gè)剛體組合,基于理論力學(xué)建立車(chē)輛動(dòng)力學(xué)模型并將其作為估計(jì)算法研究的對(duì)象。考慮建模的復(fù)雜程度和實(shí)時(shí)計(jì)算的需要,需引用縱向、側(cè)向和橫擺三個(gè)自由度的狀態(tài)變量來(lái)進(jìn)行估算,因此,在二自由度車(chē)輛模型的基礎(chǔ)上引入一個(gè)縱向運(yùn)動(dòng)自由度,使該模型具有橫擺、側(cè)向和縱向三個(gè)自由度[9],模型簡(jiǎn)圖見(jiàn)圖1。

    圖1 三自由度車(chē)輛動(dòng)力學(xué)模型Fig.1 3-DOF vehicle dynamics model

    圖1中,假設(shè)車(chē)輛是對(duì)稱(chēng)的,oxy是固定在車(chē)輛質(zhì)心處的坐標(biāo)系,x軸在汽車(chē)的縱向?qū)ΨQ(chēng)軸上,規(guī)定駕駛員的正向?yàn)閤軸正方向;y軸在水平面內(nèi),與x軸垂直,通過(guò)車(chē)輛質(zhì)心o點(diǎn),規(guī)定駕駛員的左向?yàn)閥軸正方向。其中,轉(zhuǎn)向盤(pán)轉(zhuǎn)角等于前輪轉(zhuǎn)角與轉(zhuǎn)向盤(pán)轉(zhuǎn)角到前輪轉(zhuǎn)角的傳動(dòng)比的乘積。對(duì)三自由度車(chē)輛模型建立方程如下:

    (1)

    (2)

    (3)

    (4)

    式中,vx為縱向速度;r為橫擺角速度;k1為前輪總側(cè)偏剛度;k2為后輪總側(cè)偏剛度;m為汽車(chē)總質(zhì)量;a為前軸到質(zhì)心的距離;b為后軸到質(zhì)心的距離;δ為前輪轉(zhuǎn)角;ax為縱向加速度;ay為側(cè)向加速度;Iz為繞z軸的轉(zhuǎn)動(dòng)慣量;β為質(zhì)心側(cè)偏角。

    2 基于RLS法的汽車(chē)總質(zhì)量估計(jì)

    RLS法是一種參數(shù)辨識(shí)方法,是最小二乘法的遞推形式,適用于單輸入單輸出(single input single output, SISO)動(dòng)態(tài)模型的在線參數(shù)辨識(shí),采用該方法不需要已知任何先驗(yàn)統(tǒng)計(jì)特性,并且計(jì)算量小、收斂速度快[10]。為更好地估計(jì)汽車(chē)總質(zhì)量,選擇在汽車(chē)平穩(wěn)行駛狀態(tài)時(shí)進(jìn)行估計(jì),考慮實(shí)際情況,在底盤(pán)主動(dòng)安全系統(tǒng)工作時(shí),應(yīng)當(dāng)有較好的汽車(chē)狀態(tài)估計(jì)輸入,即總質(zhì)量應(yīng)最先辨識(shí)出來(lái),故選擇汽車(chē)起步工況進(jìn)行參數(shù)辨識(shí),此時(shí)可認(rèn)為輪胎的側(cè)偏特性處在線性范圍,因此采用線性二自由度車(chē)輛模型[11]。汽車(chē)總質(zhì)量辨識(shí)模型如下:

    (5)

    RLS輸入輸出遞推方程如下[10]:

    (6)

    y(k)=φ(k)Tθ+e(k)

    (7)

    將式(6)改寫(xiě)成式(7),其中u(k)為輸入;y(k)為輸出;n為待辨識(shí)參數(shù)的總個(gè)數(shù);e(k)為由觀測(cè)噪聲或建模引起的誤差;θ為待估參數(shù)組成的向量,θ=(a1,a2,…,an,b1,b2,…,bn)T;φ(k)=[-y(k-1) … -y(k-n)u(k-1) …

    u(k-n)]T。下面針對(duì)辨識(shí)模型(式(5)),令θ=(m),φ(k)=[ay],進(jìn)行汽車(chē)總質(zhì)量估計(jì)。

    參數(shù)辨識(shí)增益:

    K(k)=P(k-1)φ(k)[φ(k)TP(k-1)φ(k)+1]-1

    參數(shù)辨識(shí)更新:

    辨識(shí)誤差更新:

    P(k)=[I-K(k)φ(k)T]P(k-1)

    3 模糊自適應(yīng)擴(kuò)展卡爾曼濾波

    根據(jù)式(1)~式(4),建立狀態(tài)空間方程:

    (8)

    其中,三維狀態(tài)向量x(t)=(r,β,vx)T;一維輸出向量y(t)=(ay);二維控制向量u(t)=(δ,ax)T;w(t)與v(t)是均值為零的白噪聲,且互不相關(guān)。

    將非線性函數(shù)f(x(t),u(t),w(t))、h(x(t),v(t))分別對(duì)狀態(tài)x(t)求偏導(dǎo),得到雅可比矩陣F(t)與H(t):

    至此完成了模型線性化。

    由于EKF是基于離散狀態(tài)空間方程的算法,故采用前向歐拉法將車(chē)輛連續(xù)狀態(tài)空間方程進(jìn)行離散,得到相應(yīng)離散系統(tǒng)的狀態(tài)空間方程:

    vx(k)=(r(k-1)β(k-1)vx(k-1)+
    ax(k))Δt+vx(k-1)

    其中,Δt為采樣時(shí)間步長(zhǎng)。觀測(cè)矩陣如下:

    傳統(tǒng)的EKF算法在濾波計(jì)算中,假設(shè)觀測(cè)噪聲協(xié)方差矩陣已知且固定不變。在汽車(chē)行駛過(guò)程中,傳感器的噪聲水平受不同的工況和外界環(huán)境的影響,傳統(tǒng)估計(jì)算法設(shè)定觀測(cè)噪聲協(xié)方差矩陣不變,這會(huì)影響EKF的估計(jì)精度。為此,建立一個(gè)模糊控制器對(duì)觀測(cè)噪聲協(xié)方差矩陣進(jìn)行實(shí)時(shí)調(diào)整,使EKF算法具有自適應(yīng)性。

    離散狀態(tài)空間方程為

    xk+1=Fkxk+Bkuk+wk
    zk=Hkxk+vk

    下面給出離散的EKF估計(jì)過(guò)程。

    狀態(tài)一步預(yù)測(cè):

    一步預(yù)測(cè)誤差方差陣:

    增益矩陣:

    狀態(tài)更新:

    方差更新:

    Pk=(I-KkHk)Pk|k-1

    取實(shí)際方差:

    理論方差:

    令方差的差值

    e=Pa-Pt

    (25)

    e={NB,NM,NS,Z,PS,PM,PB}

    U={PVS,PS,PSM,PM,PSB,PMB,PB,PVB}

    其中,NB為負(fù)大,NM為負(fù)中,NS為負(fù)小,Z為零,PS為正小,PM為正中,PB為正大,PVS為正的高等偏小,PSM為正的中等偏小,PSB為正的初等偏大,PMB為正的中等偏大,PVB為正的高等偏大。其輸入輸出隸屬函數(shù)如圖2~圖4所示,模糊規(guī)則如表1所示。

    圖2 輸入e的隸屬函數(shù)Fig.2 Membership function of Input variable e

    圖3 輸入的隸屬函數(shù)

    圖4 輸出U的隸屬函數(shù)Fig.4 Membership function of Input variables U

    表1 模糊規(guī)則

    定義k時(shí)刻觀測(cè)噪聲協(xié)方差的估計(jì)值為Rk,Rk=UR。通過(guò)模糊控制器不斷改變Rk的值達(dá)到EKF算法的自適應(yīng)。

    4 基于虛擬實(shí)驗(yàn)的算法驗(yàn)證

    4.1 汽車(chē)總質(zhì)量的估計(jì)驗(yàn)證

    在CarSim與MATLAB/Simulink環(huán)境下建立聯(lián)合仿真平臺(tái)。選擇普通前驅(qū)D級(jí)Sedan車(chē)型進(jìn)行仿真,整車(chē)的參數(shù)如表2所示。設(shè)定車(chē)輛速度為60 km/h。

    表2 整車(chē)參數(shù)

    仿真設(shè)定在固定附著系數(shù)為0.85的良好路面進(jìn)行蛇行試驗(yàn),前輪轉(zhuǎn)角變化如圖5所示,在1~2 s內(nèi)轉(zhuǎn)角為零,在2~8 s內(nèi)轉(zhuǎn)角呈幅值為10、周期為2 s的正弦變化,在8~10 s內(nèi)轉(zhuǎn)角為零??v向加速度變化值如圖6所示,由CarSim軟件輸出作為算法的輸入。

    圖5 前輪轉(zhuǎn)角變化圖Fig.5 The change of front wheel angle

    圖6 縱向加速度變化圖Fig.6 The change of longitudinal acceleration

    汽車(chē)總質(zhì)量初值設(shè)為1000 kg,實(shí)際值為1500 kg,基于RLS算法的汽車(chē)總質(zhì)量估計(jì)仿真結(jié)果如圖7所示。

    圖7 基于RLS的汽車(chē)總質(zhì)量的估計(jì)值based on RLS與實(shí)際值的對(duì)比Fig.7 Comparison between estimates and actual value of total vehicle quality estimation

    由圖7可知,基于RLS算法估計(jì)的汽車(chē)總質(zhì)量在初值為1000 kg的條件下能較快跟蹤到實(shí)際值,在9 s內(nèi)能實(shí)現(xiàn)實(shí)際值的良好跟蹤,且在之后波動(dòng)很小,最大偏差為23 kg,相對(duì)誤差為1.53%,此精度完全滿(mǎn)足控制要求??紤]實(shí)際情況,當(dāng)整車(chē)總質(zhì)量發(fā)生變化時(shí),在正常駕車(chē)情況下,總是汽車(chē)先停車(chē),然后上下人員或者裝卸貨物;當(dāng)汽車(chē)重新由停止?fàn)顟B(tài)啟動(dòng)而進(jìn)入加速行駛的過(guò)程中,程序又會(huì)自動(dòng)重新估計(jì)汽車(chē)總質(zhì)量,作為狀態(tài)估計(jì)的輸入,因此,當(dāng)汽車(chē)總質(zhì)量發(fā)生變化(即模型參數(shù)變化)時(shí),算法的狀態(tài)輸入能很好地跟蹤汽車(chē)總質(zhì)量的變化。由圖7可知,RLS較好地估計(jì)了汽車(chē)的總質(zhì)量,這為后續(xù)汽車(chē)狀態(tài)估計(jì)奠定了基礎(chǔ)。

    4.2 汽車(chē)關(guān)鍵狀態(tài)估計(jì)驗(yàn)證

    將CarSim軟件的狀態(tài)輸出結(jié)果作為實(shí)際值,將RLS與FAEKF相結(jié)合的算法估計(jì)結(jié)果,與EKF算法估計(jì)結(jié)果和實(shí)際值進(jìn)行對(duì)比,以驗(yàn)證本文優(yōu)化算法的有效性。

    在狀態(tài)估計(jì)過(guò)程中,不論是EKF算法還是RLS與FAEKF結(jié)合的算法,它們的過(guò)程噪聲協(xié)方差矩陣的值都固定不變且相同,即Q=diag(0.1,0.1,0.1)。EKF算法的觀測(cè)協(xié)方差矩陣值也固定不變, EKF觀測(cè)協(xié)方差矩陣值R=[0.001]。EKF算法的汽車(chē)總質(zhì)量為整備質(zhì)量,即為1370 kg。RLS與FAEKF結(jié)合的算法觀測(cè)協(xié)方差初值Rk=[0.001],且隨著時(shí)間的推移,觀測(cè)協(xié)方差Rk發(fā)生變化,變化規(guī)律如表3所示。FAEKF+RLS算法的汽車(chē)總質(zhì)量為RLS的估計(jì)值,仿真時(shí)間為10 s。汽車(chē)關(guān)鍵狀態(tài)估計(jì)結(jié)果如圖8~圖10所示。

    表3 觀測(cè)噪聲協(xié)方差變化規(guī)律

    圖8 縱向車(chē)速估算對(duì)比圖Fig.8 Comparison of longitudinal vehicle speed estimation

    圖9 質(zhì)心側(cè)偏角估算對(duì)比圖Fig.9 Comparison of sideslip angle estimation

    圖10 橫擺角速度估算對(duì)比圖Fig.10 Comparison of yaw rate estimation

    圖8是EKF算法與FAEKF+RLS算法對(duì)縱向車(chē)速的估計(jì)值與實(shí)際值的對(duì)比圖??梢钥闯觯現(xiàn)AEKF+RLS算法的縱向車(chē)速估計(jì)值與實(shí)際值一直很接近,估計(jì)效果比較理想,估計(jì)值與實(shí)際值最大瞬態(tài)偏差只有0.08m/s,相對(duì)誤差為0.04%;而EKF算法的估計(jì)值與實(shí)際值的最大瞬態(tài)偏差為0.36m/s,相對(duì)誤差為2.1%。在整個(gè)仿真時(shí)間內(nèi)EKF算法估計(jì)值偏差較大的主要原因是:沒(méi)有實(shí)時(shí)估計(jì)汽車(chē)總質(zhì)量,以整備質(zhì)量代替汽車(chē)總質(zhì)量;觀測(cè)協(xié)方差已更新,而EKF仍然按固定的觀測(cè)方差進(jìn)行循環(huán)計(jì)算,導(dǎo)致偏差較大。在2~8s內(nèi)偏差越來(lái)越大,主要是由于給定觀測(cè)協(xié)方差越來(lái)越大。在0~2s與8~10s內(nèi)EKF的觀測(cè)方差值與動(dòng)態(tài)給定值相同,所不同的是,車(chē)速估計(jì)值在0~1s內(nèi)有較大下降趨勢(shì),而在8~10s內(nèi)估計(jì)值與實(shí)際值接近。這種差異是輸入量縱向加速度的變化導(dǎo)致的,在圖6中0~1s內(nèi)縱向加速度急劇下降且為負(fù)值,即減速行駛造成EKF估計(jì)值比實(shí)際值低,在8s以后加速度基本接近零,使得EKF估計(jì)值與實(shí)際值偏差較小。

    圖9中,EKF算法的估計(jì)值與實(shí)際值最大瞬態(tài)偏差為0.83°,相對(duì)誤差為63.8%;FAEKF+RLS算法的估計(jì)值與實(shí)際值最大瞬態(tài)偏差為0.03°,相對(duì)誤差為1.76%。

    圖10中,EKF算法估計(jì)值與實(shí)際值最大瞬態(tài)偏差為27.3(°)/s,相對(duì)誤差為108%;FAEKF+RLS算法的估計(jì)值與實(shí)際值最大瞬態(tài)偏差為1.4(°)/s,相對(duì)誤差為4.6%。從圖9、圖10中可以看出,質(zhì)心側(cè)偏角與橫擺角速度的EKF算法估計(jì)值、FAEKF+RLS算法估計(jì)值與實(shí)際值相比,在0~2s與8~10s內(nèi)偏差都很小,主要是觀測(cè)協(xié)方差相同。這與縱向速度在0~2s內(nèi)與實(shí)際值偏差較大的情況比較,區(qū)別在于縱向速度直接與縱向加速度有關(guān),而橫擺角速度與質(zhì)心側(cè)偏角間接相關(guān),具體函數(shù)關(guān)系見(jiàn)式(1)~式(3)。在2~8s內(nèi),由于以整備質(zhì)量代替汽車(chē)總質(zhì)量,以及觀測(cè)協(xié)方差發(fā)生改變且變動(dòng)越來(lái)越大,所以造成EKF算法估計(jì)值的偏差越來(lái)越大,F(xiàn)AEKF+RLS算法由于對(duì)汽車(chē)總質(zhì)量進(jìn)行實(shí)時(shí)估計(jì)且對(duì)觀測(cè)方差及時(shí)更新,使得估計(jì)值與實(shí)際值偏差很小。以上分析說(shuō)明本文算法估計(jì)效果良好。

    5 結(jié)論

    (1)本文提出將RLS算法與FAEKF算法相結(jié)合進(jìn)行車(chē)輛狀態(tài)估計(jì),使估計(jì)誤差大大減小,達(dá)到與汽車(chē)總質(zhì)量恒定且觀測(cè)噪聲為定值狀態(tài)下的估計(jì)精度相當(dāng)?shù)男Ч?/p>

    (2)虛擬實(shí)驗(yàn)結(jié)果表明,RLS算法和FAEKF算法相結(jié)合可以很好地跟蹤虛擬實(shí)驗(yàn)值,優(yōu)于單純的EKF算法,具有較好的精度和魯棒性。

    [1] 張鳳嬌,魏民祥,趙萬(wàn)忠. 基于蟻群優(yōu)化UKF算法的汽車(chē)狀態(tài)估計(jì)[J].中國(guó)機(jī)械工程,2015,26(22) :3046-3050.ZHANGFengjiao,WEIMinxiang,ZHAOWanzhong.VehicleStateEstimationBasedonAntColonyOptimizationAlgorithm[J].ChinaMechanicalEngineering, 2015,26(22):3046-3050.

    [2]WENZELTA,BURNHAMKJ.DualExtendedKalmanFilterforVehicleStateandParameterEstimation[J].VehicleSystemDynamics, 2007, 42(2):153-171.

    [3] 宗長(zhǎng)富,胡丹,楊肖,等. 基于擴(kuò)展Kalman濾波的汽車(chē)行駛狀態(tài)估計(jì)[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版),2009,39(1): 7-11.ZONGChangfu,HUDan,YANGXiao,etal.VehicleDrivingStateEstimationBasedonExtendedKalmanFilter[J].JournalofJilinUniversity(EngineeringandTechnologyEdition), 2009, 39(1): 7-11.

    [4]ANTONOVS,FCHNA,KUGIA.UnscentedKalmanFilterforVehicleStateEstimation[J].VehicleSystemDynamics, 2011,49(9):1497-1520.

    [5] 趙又群,林棻.基于UKF算法的汽車(chē)狀態(tài)估計(jì)[J].中國(guó)機(jī)械工程,2010,21(5):615-619.ZHAOYouqun,LINFen.VehicleStateEstimationBasedonUnscentedKalmanFilterAlgorithm[J].ChinaMechanicalEngineering, 2010, 21(5): 615-619.

    [6]MCLZIS,SABBIONIE.OntheVehicleSideslipAngleEstimationthroughNeuralNetworks:NumericalandExperimentalResults[J].MechanicalSystemsandSignalProccssing, 2011,25(6):2005-2019.

    [7] 張小龍,李亮,李紅志,等.基于改進(jìn)RBF網(wǎng)的汽車(chē)側(cè)偏角估計(jì)方法試驗(yàn)研究[J].機(jī)械工程學(xué)報(bào),2010,46(22) : 105-110.ZHANGXiaolong,LILiang,LIHongzhi,etal.ExperimentalResearchonVehicleSideslipAngleEstimationBasedonImprovedRBFNeuralNetworks[J].ChineseJournalofMechanicalEngineering, 2010, 46(22): 105-110.

    [8] 施樹(shù)明,LUPKCRH,BREMMERP, 等.基于模糊邏輯的車(chē)輛側(cè)偏角估計(jì)方法[J].汽車(chē)工程,2005,27(4):426-430.SHIShuming,LUPKERH,BREMMERP,etal.EstimationofVehicleSideSlipAngleBasedonFuzzyLogic[J].AutomotiveEngineering, 2005, 27(4): 426-430.

    [9] 武鐘財(cái). 基于擴(kuò)展卡爾曼濾波的路面附著系數(shù)估計(jì)算法研究[D]. 長(zhǎng)春:吉林大學(xué),2008.WUZhongcai.ResearchontheAlgorithmoftheRoadFrictionCoefficientEstimatiomBasedontheExtendedKalmanFilter[D].Changchun:JilinUniversity, 2008.

    [10] 蔡季冰. 系統(tǒng)辨識(shí)[M]. 北京: 北京理工大學(xué)出版社,1989: 25-34.CAIJibing.SystemIdentification[M].Beijing:BeijingInstituteofTechnologyPress,1989: 25-34.

    [11] 余志生.汽車(chē)?yán)碚揫M]. 5版. 北京:機(jī)械工業(yè)出版社,2012 :144-146.YUZhisheng.VehicleDynamics[M]. 5thed.Beijing:MechanicalIndustryPress, 2012:144-146.

    (編輯 王旻玥)

    Vehicle State Estimation Based on Combined RLS and FAEKF

    WANG Yan1WEI Minxiang1ZHAO Wanzhong1ZHANG Fengjiao1,2YAN Mingyue1

    1.College of Energy and Power Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing,210016 2.School of Mechanical &Vehicle Engineering,Changzhou Institute of Technology, Changzhou,Jiangsu,213002

    For the problems of observation noise time-varying characteristics and model parameter variations in vehicle state estimation, a new algorithm which consisted of RLS method and FAEKF was proposed. The new algorithm was proposed based on 3-DOF nonlinear vehicle dynamics model in order to realize real time update of model parameters and observation noises. Firstly, the total mass of the vehicle was estimated by RLS. Then, a fuzzy controller was established to track the observation noises of extended Kalman filters. Finally, the algorithm was verified using CarSim and MATLAB/Simulink. Results show that the estimation accuracy of the new algorithm is higher than that of the traditional extended Kalman filter. It may provide theoretical support for the development of automobile active control systems.

    automobile total quality estimation; state estimation; recursive least squares(RLS); fuzzy adaptive extended Kalman filter(FAEKF)

    2016-08-09

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51375007);江蘇省自然科學(xué)基金資助項(xiàng)目(SBK2015022352);常州市科技計(jì)劃應(yīng)用基礎(chǔ)研究項(xiàng)目(CJ20159011)

    U461.6DOI:10.3969/j.issn.1004-132X.2017.06.019

    汪,男,1992年生。南京航空航天大學(xué)能源與動(dòng)力學(xué)院碩士研究生。主要研究方向?yàn)槠?chē)電子安全控制技術(shù)。E-mail:1305905618@qq.com。魏民祥,男,1963年生。南京航空航天大學(xué)能源與動(dòng)力學(xué)院教授、博士研究生導(dǎo)師。趙萬(wàn)忠,男,1982年生。南京航空航天大學(xué)能源與動(dòng)力學(xué)院教授。張鳳嬌,女,1978年生。常州工學(xué)院機(jī)械與車(chē)輛工程學(xué)院副教授,南京航空航天大學(xué)能源與動(dòng)力學(xué)院博士研究生。嚴(yán)明月,女,1992年生。南京航空航天大學(xué)能源與動(dòng)力學(xué)院碩士研究生。

    猜你喜歡
    估計(jì)值協(xié)方差卡爾曼濾波
    一道樣本的數(shù)字特征與頻率分布直方圖的交匯問(wèn)題
    統(tǒng)計(jì)信息
    2018年4月世界粗鋼產(chǎn)量表(續(xù))萬(wàn)噸
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    基于擴(kuò)展卡爾曼濾波的PMSM無(wú)位置傳感器控制
    基于自適應(yīng)卡爾曼濾波的新船舶試航系統(tǒng)
    2014年5月世界粗鋼產(chǎn)量表萬(wàn)噸
    亚洲激情五月婷婷啪啪| 亚洲中文av在线| 中文精品一卡2卡3卡4更新| 欧美精品一区二区大全| 亚洲人成77777在线视频| 国产精品麻豆人妻色哟哟久久| 天天躁夜夜躁狠狠躁躁| 五月天丁香电影| 国产免费一区二区三区四区乱码| av天堂久久9| 在线观看国产h片| 女警被强在线播放| 看十八女毛片水多多多| 涩涩av久久男人的天堂| 99热全是精品| netflix在线观看网站| 男人操女人黄网站| 亚洲精品久久成人aⅴ小说| 亚洲色图综合在线观看| 亚洲av日韩在线播放| 男男h啪啪无遮挡| 亚洲国产av影院在线观看| 两性夫妻黄色片| 电影成人av| 美女脱内裤让男人舔精品视频| 国产在线免费精品| 成人午夜精彩视频在线观看| 国产一区有黄有色的免费视频| 男人舔女人的私密视频| 妹子高潮喷水视频| 9热在线视频观看99| 美国免费a级毛片| 亚洲,一卡二卡三卡| 午夜免费鲁丝| 久久热在线av| 国产精品一区二区在线不卡| av片东京热男人的天堂| 18在线观看网站| 电影成人av| 欧美激情高清一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲av日韩精品久久久久久密 | 亚洲国产精品一区三区| 亚洲精品一二三| 女人被躁到高潮嗷嗷叫费观| 精品福利观看| 亚洲欧美清纯卡通| 中文字幕高清在线视频| www.熟女人妻精品国产| av片东京热男人的天堂| 国产极品粉嫩免费观看在线| 亚洲一区二区三区欧美精品| 汤姆久久久久久久影院中文字幕| 69精品国产乱码久久久| 久久久亚洲精品成人影院| 亚洲精品美女久久久久99蜜臀 | 一本—道久久a久久精品蜜桃钙片| 亚洲国产av新网站| 欧美亚洲 丝袜 人妻 在线| 精品免费久久久久久久清纯 | 久久精品熟女亚洲av麻豆精品| 黑人巨大精品欧美一区二区蜜桃| 国产精品国产三级国产专区5o| 在线 av 中文字幕| 欧美性长视频在线观看| 亚洲成人手机| 亚洲精品国产一区二区精华液| 日韩中文字幕视频在线看片| 女人精品久久久久毛片| 亚洲色图 男人天堂 中文字幕| 成人手机av| 日韩一本色道免费dvd| 久久女婷五月综合色啪小说| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜一区二区| 欧美激情 高清一区二区三区| 婷婷成人精品国产| 蜜桃国产av成人99| 男人爽女人下面视频在线观看| 欧美性长视频在线观看| 午夜福利影视在线免费观看| 精品少妇内射三级| 丝袜美足系列| 亚洲,一卡二卡三卡| 丝袜美足系列| 国产人伦9x9x在线观看| 亚洲综合色网址| 亚洲综合色网址| 亚洲精品国产色婷婷电影| 国产麻豆69| 91国产中文字幕| 黄色一级大片看看| 久久久欧美国产精品| 日韩 欧美 亚洲 中文字幕| 欧美老熟妇乱子伦牲交| 亚洲av美国av| 中文字幕人妻丝袜一区二区| 老汉色av国产亚洲站长工具| 看十八女毛片水多多多| 岛国毛片在线播放| 天天添夜夜摸| 99精国产麻豆久久婷婷| 精品免费久久久久久久清纯 | 黄片小视频在线播放| 欧美黄色片欧美黄色片| 午夜福利,免费看| a级片在线免费高清观看视频| 亚洲国产看品久久| 18禁国产床啪视频网站| 在线精品无人区一区二区三| 亚洲av美国av| 国产淫语在线视频| 国产免费又黄又爽又色| 国产又爽黄色视频| 国产亚洲av片在线观看秒播厂| 国产爽快片一区二区三区| 日韩制服丝袜自拍偷拍| 少妇精品久久久久久久| 一个人免费看片子| 久久精品熟女亚洲av麻豆精品| 亚洲人成网站在线观看播放| 久久av网站| 女人久久www免费人成看片| 美女高潮到喷水免费观看| 九草在线视频观看| 欧美精品人与动牲交sv欧美| 免费观看a级毛片全部| 欧美精品人与动牲交sv欧美| 亚洲国产精品一区二区三区在线| 日韩av不卡免费在线播放| 国产精品一区二区在线观看99| 成人亚洲精品一区在线观看| 这个男人来自地球电影免费观看| 亚洲美女黄色视频免费看| 欧美日韩一级在线毛片| 女性被躁到高潮视频| 国产成人免费无遮挡视频| 校园人妻丝袜中文字幕| 亚洲,一卡二卡三卡| 欧美人与性动交α欧美软件| 亚洲专区中文字幕在线| 国产在线视频一区二区| 超色免费av| 男女边吃奶边做爰视频| 精品卡一卡二卡四卡免费| 久久人人爽人人片av| 国产精品一二三区在线看| 久久久久久免费高清国产稀缺| 欧美 日韩 精品 国产| 欧美精品一区二区免费开放| 97人妻天天添夜夜摸| 欧美精品啪啪一区二区三区 | 亚洲精品乱久久久久久| 久久精品成人免费网站| 国产成人精品久久久久久| 亚洲国产精品一区二区三区在线| 99精国产麻豆久久婷婷| 国产伦理片在线播放av一区| 青春草视频在线免费观看| 男人添女人高潮全过程视频| 精品人妻一区二区三区麻豆| 久久久久国产精品人妻一区二区| 国产成人av激情在线播放| 欧美日韩亚洲国产一区二区在线观看 | 日本五十路高清| 国产黄色视频一区二区在线观看| 亚洲中文日韩欧美视频| 丁香六月欧美| av天堂久久9| 久久热在线av| 亚洲 欧美一区二区三区| 亚洲国产欧美在线一区| 欧美久久黑人一区二区| 制服人妻中文乱码| www.熟女人妻精品国产| 色视频在线一区二区三区| 久久久久久久精品精品| 亚洲欧美一区二区三区久久| 精品一区在线观看国产| 久久久久久亚洲精品国产蜜桃av| 中国美女看黄片| 精品国产乱码久久久久久男人| 国产xxxxx性猛交| 亚洲精品国产av蜜桃| 欧美日韩亚洲高清精品| 老司机亚洲免费影院| 亚洲色图综合在线观看| 精品久久蜜臀av无| 一区二区av电影网| 免费观看a级毛片全部| 在线观看免费午夜福利视频| 一级毛片黄色毛片免费观看视频| 国产欧美亚洲国产| 欧美日本中文国产一区发布| 考比视频在线观看| 亚洲自偷自拍图片 自拍| 看免费av毛片| 久久久久精品国产欧美久久久 | 又粗又硬又长又爽又黄的视频| 最新在线观看一区二区三区 | 亚洲欧美成人综合另类久久久| 电影成人av| 国产免费又黄又爽又色| 久久久精品免费免费高清| av欧美777| 一区福利在线观看| 国产精品免费视频内射| 精品少妇内射三级| 久久99热这里只频精品6学生| 国产97色在线日韩免费| 亚洲欧美成人综合另类久久久| 永久免费av网站大全| 超色免费av| 51午夜福利影视在线观看| 99re6热这里在线精品视频| 女人被躁到高潮嗷嗷叫费观| 一边摸一边做爽爽视频免费| 丁香六月欧美| 欧美精品av麻豆av| 在线观看免费日韩欧美大片| 视频在线观看一区二区三区| 精品久久久久久久毛片微露脸 | 大话2 男鬼变身卡| 久久精品国产亚洲av高清一级| 欧美变态另类bdsm刘玥| 1024香蕉在线观看| 老司机影院毛片| 99精国产麻豆久久婷婷| 亚洲精品美女久久久久99蜜臀 | 国产男人的电影天堂91| 国产精品人妻久久久影院| 亚洲av电影在线进入| 久久国产精品大桥未久av| 交换朋友夫妻互换小说| 亚洲天堂av无毛| 美女主播在线视频| xxx大片免费视频| 午夜免费鲁丝| 欧美 亚洲 国产 日韩一| 韩国高清视频一区二区三区| 精品国产超薄肉色丝袜足j| 精品少妇一区二区三区视频日本电影| 飞空精品影院首页| 久久久精品免费免费高清| 岛国毛片在线播放| 满18在线观看网站| 亚洲人成77777在线视频| 欧美激情高清一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲国产中文字幕在线视频| 黄片播放在线免费| 最近手机中文字幕大全| 婷婷成人精品国产| 精品视频人人做人人爽| 免费看不卡的av| 欧美日韩亚洲综合一区二区三区_| 亚洲人成电影免费在线| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩高清在线视频 | 国产精品免费视频内射| 晚上一个人看的免费电影| avwww免费| 午夜福利乱码中文字幕| 精品免费久久久久久久清纯 | 亚洲,欧美精品.| 亚洲精品美女久久久久99蜜臀 | 搡老乐熟女国产| 在线观看国产h片| 亚洲精品久久午夜乱码| 久久99热这里只频精品6学生| 久久久久久久久免费视频了| 国语对白做爰xxxⅹ性视频网站| 欧美日韩黄片免| 亚洲成人免费av在线播放| 精品人妻在线不人妻| 亚洲精品美女久久av网站| 日本午夜av视频| 国产男人的电影天堂91| 黄色a级毛片大全视频| 国产av国产精品国产| a级毛片黄视频| 久久人人爽人人片av| 亚洲伊人色综图| 久久久欧美国产精品| 满18在线观看网站| 亚洲国产av新网站| 黑人欧美特级aaaaaa片| 真人做人爱边吃奶动态| 777久久人妻少妇嫩草av网站| 国产一级毛片在线| 国产又爽黄色视频| 亚洲欧美一区二区三区黑人| 中国美女看黄片| 国产一区亚洲一区在线观看| 91精品国产国语对白视频| 中文乱码字字幕精品一区二区三区| 另类精品久久| 天天影视国产精品| 欧美人与善性xxx| 欧美日韩综合久久久久久| 午夜老司机福利片| 青春草亚洲视频在线观看| 亚洲精品国产av成人精品| 自线自在国产av| 日韩 欧美 亚洲 中文字幕| 国产黄色免费在线视频| 亚洲精品乱久久久久久| 亚洲 欧美一区二区三区| 欧美日韩国产mv在线观看视频| 久久精品亚洲熟妇少妇任你| 纯流量卡能插随身wifi吗| 国产主播在线观看一区二区 | 久久精品久久精品一区二区三区| 人妻人人澡人人爽人人| 亚洲欧美激情在线| 人人妻人人爽人人添夜夜欢视频| 老汉色av国产亚洲站长工具| 欧美在线一区亚洲| e午夜精品久久久久久久| 日本av手机在线免费观看| 91九色精品人成在线观看| 我要看黄色一级片免费的| 天天添夜夜摸| 中文字幕另类日韩欧美亚洲嫩草| 老司机影院毛片| 永久免费av网站大全| 免费不卡黄色视频| 亚洲欧美日韩高清在线视频 | 50天的宝宝边吃奶边哭怎么回事| 青青草视频在线视频观看| 国产成人免费无遮挡视频| 无限看片的www在线观看| 午夜福利一区二区在线看| 久久鲁丝午夜福利片| 国产精品香港三级国产av潘金莲 | 亚洲国产精品一区二区三区在线| 亚洲av成人精品一二三区| 黄色a级毛片大全视频| 下体分泌物呈黄色| 久久精品熟女亚洲av麻豆精品| 大话2 男鬼变身卡| 国产女主播在线喷水免费视频网站| cao死你这个sao货| 19禁男女啪啪无遮挡网站| 少妇裸体淫交视频免费看高清 | 久久国产精品大桥未久av| 亚洲av欧美aⅴ国产| 亚洲国产av影院在线观看| 国产精品欧美亚洲77777| 超碰成人久久| 人人妻,人人澡人人爽秒播 | 欧美激情极品国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 七月丁香在线播放| 精品高清国产在线一区| 黄色 视频免费看| 国产精品一区二区在线观看99| 在线观看免费午夜福利视频| 脱女人内裤的视频| 天天躁夜夜躁狠狠躁躁| 久久99一区二区三区| 777米奇影视久久| 国产精品国产三级国产专区5o| 国产亚洲欧美在线一区二区| 免费观看a级毛片全部| 热99国产精品久久久久久7| 91成人精品电影| 色网站视频免费| 国产精品av久久久久免费| 久久久久国产一级毛片高清牌| 人人澡人人妻人| 国产亚洲一区二区精品| 国产亚洲av片在线观看秒播厂| 五月开心婷婷网| 亚洲美女黄色视频免费看| 激情五月婷婷亚洲| 青青草视频在线视频观看| 中文乱码字字幕精品一区二区三区| 久久久国产欧美日韩av| 久久女婷五月综合色啪小说| 免费不卡黄色视频| 50天的宝宝边吃奶边哭怎么回事| 国产黄色视频一区二区在线观看| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃| 色94色欧美一区二区| 18禁黄网站禁片午夜丰满| 天天躁日日躁夜夜躁夜夜| 丁香六月天网| 亚洲久久久国产精品| 精品少妇黑人巨大在线播放| 久久精品成人免费网站| 女人精品久久久久毛片| 国产精品99久久99久久久不卡| 80岁老熟妇乱子伦牲交| 日本一区二区免费在线视频| 国精品久久久久久国模美| 国产精品一二三区在线看| 欧美黑人欧美精品刺激| 视频在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 久久久久久久大尺度免费视频| 考比视频在线观看| 人成视频在线观看免费观看| www.精华液| 乱人伦中国视频| 精品国产一区二区三区久久久樱花| 国产精品秋霞免费鲁丝片| 亚洲成人手机| 国产精品99久久99久久久不卡| kizo精华| 成人手机av| 多毛熟女@视频| 日本五十路高清| 国产欧美亚洲国产| 日韩一卡2卡3卡4卡2021年| 免费在线观看日本一区| 国产欧美日韩一区二区三区在线| 国产免费现黄频在线看| 日韩伦理黄色片| 国产99久久九九免费精品| 国产一区有黄有色的免费视频| 亚洲七黄色美女视频| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 美女高潮到喷水免费观看| 观看av在线不卡| 国产三级黄色录像| 熟女少妇亚洲综合色aaa.| 亚洲熟女毛片儿| 午夜两性在线视频| 久久99精品国语久久久| 免费av中文字幕在线| 一级黄片播放器| 免费在线观看完整版高清| 各种免费的搞黄视频| 丝袜脚勾引网站| 人妻 亚洲 视频| 国产成人av教育| 色视频在线一区二区三区| 亚洲三区欧美一区| 亚洲中文av在线| 久久久国产成人精品二区| 欧美国产精品va在线观看不卡| 两个人免费观看高清视频| 久久香蕉激情| 午夜两性在线视频| 久久久久久久久久黄片| 操出白浆在线播放| 搞女人的毛片| 亚洲中文av在线| 特大巨黑吊av在线直播 | 老司机靠b影院| 狠狠狠狠99中文字幕| 老司机福利观看| 久久欧美精品欧美久久欧美| 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 亚洲人成77777在线视频| 欧美色欧美亚洲另类二区| 白带黄色成豆腐渣| av在线天堂中文字幕| 亚洲国产精品久久男人天堂| 国产精品乱码一区二三区的特点| 国产精品亚洲美女久久久| 少妇熟女aⅴ在线视频| 亚洲,欧美精品.| 免费在线观看成人毛片| 久久久久免费精品人妻一区二区 | 国产午夜福利久久久久久| 国产黄a三级三级三级人| 最近最新免费中文字幕在线| 一本一本综合久久| 免费女性裸体啪啪无遮挡网站| 天堂动漫精品| 国产不卡一卡二| 亚洲中文字幕日韩| 一级毛片精品| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产亚洲av香蕉五月| 一夜夜www| 亚洲第一欧美日韩一区二区三区| 亚洲欧美精品综合一区二区三区| 男女午夜视频在线观看| 51午夜福利影视在线观看| 中文字幕最新亚洲高清| 人人妻人人澡欧美一区二区| 午夜免费激情av| 怎么达到女性高潮| 欧美亚洲日本最大视频资源| 丝袜美腿诱惑在线| 成年女人毛片免费观看观看9| 成年版毛片免费区| а√天堂www在线а√下载| 亚洲天堂国产精品一区在线| 一区二区三区激情视频| 亚洲国产欧美一区二区综合| 欧美最黄视频在线播放免费| 一a级毛片在线观看| av片东京热男人的天堂| 黄色成人免费大全| 琪琪午夜伦伦电影理论片6080| 最新美女视频免费是黄的| 777久久人妻少妇嫩草av网站| 欧美激情极品国产一区二区三区| 校园春色视频在线观看| 久久人人精品亚洲av| 青草久久国产| 女人爽到高潮嗷嗷叫在线视频| 色综合站精品国产| 99热这里只有精品一区 | 免费搜索国产男女视频| 高潮久久久久久久久久久不卡| 亚洲国产毛片av蜜桃av| 可以在线观看的亚洲视频| 亚洲色图av天堂| 久久精品影院6| 法律面前人人平等表现在哪些方面| 日韩欧美一区视频在线观看| 一区二区三区精品91| 搡老熟女国产l中国老女人| 免费看美女性在线毛片视频| 50天的宝宝边吃奶边哭怎么回事| 国产99白浆流出| 色精品久久人妻99蜜桃| 少妇熟女aⅴ在线视频| 国产精品久久久久久亚洲av鲁大| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 国产主播在线观看一区二区| 国产精品av久久久久免费| 欧美激情久久久久久爽电影| 免费女性裸体啪啪无遮挡网站| 中文字幕最新亚洲高清| 日韩欧美 国产精品| 久久久精品欧美日韩精品| 国产aⅴ精品一区二区三区波| 天天躁狠狠躁夜夜躁狠狠躁| 1024手机看黄色片| 窝窝影院91人妻| 男人舔女人下体高潮全视频| 男人操女人黄网站| 中文字幕av电影在线播放| 男人舔女人的私密视频| 69av精品久久久久久| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 亚洲七黄色美女视频| 不卡av一区二区三区| 黄色视频不卡| 看片在线看免费视频| 丁香欧美五月| 麻豆av在线久日| 久久久久久久久免费视频了| 亚洲男人天堂网一区| 最新在线观看一区二区三区| 999精品在线视频| 男人操女人黄网站| 精品福利观看| 久久久精品欧美日韩精品| 怎么达到女性高潮| 搡老岳熟女国产| 桃红色精品国产亚洲av| 色综合婷婷激情| 久久久久久久午夜电影| 久久久久精品国产欧美久久久| 91字幕亚洲| 又大又爽又粗| 国产亚洲精品av在线| 国产精品爽爽va在线观看网站 | 久久久久九九精品影院| 超碰成人久久| 18禁国产床啪视频网站| e午夜精品久久久久久久| 国产高清视频在线播放一区| 人成视频在线观看免费观看| 国产黄片美女视频| 亚洲无线在线观看| 精品久久久久久久久久免费视频| 中文字幕人成人乱码亚洲影| 国产视频内射| 亚洲国产精品sss在线观看| 中亚洲国语对白在线视频| 欧美乱色亚洲激情| 十八禁网站免费在线| 国产精品久久电影中文字幕| 久久精品国产清高在天天线| 久久久国产成人精品二区| 国产精品99久久99久久久不卡| 制服丝袜大香蕉在线| 欧美久久黑人一区二区| 国产在线观看jvid| 很黄的视频免费| 亚洲中文字幕一区二区三区有码在线看 | 免费无遮挡裸体视频| 成人亚洲精品一区在线观看| 精品国内亚洲2022精品成人| 久久香蕉国产精品| 黄片播放在线免费| 他把我摸到了高潮在线观看| 老司机福利观看| 韩国精品一区二区三区| 国产欧美日韩一区二区精品| 97碰自拍视频| 别揉我奶头~嗯~啊~动态视频| 国产黄片美女视频| 男女之事视频高清在线观看| 亚洲av成人一区二区三| 十八禁网站免费在线| 韩国av一区二区三区四区| 国产三级黄色录像|