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

    基于慣性測(cè)量器件的無(wú)線(xiàn)步態(tài)分析平臺(tái)

    2016-10-12 02:17:50鳴,森,洋,
    關(guān)鍵詞:腳跟步態(tài)角速度

    姜   鳴, 仇   森, 趙   洋, 張   志

    ( 1.東莞理工學(xué)院 電子工程學(xué)院, 廣東 東莞 523808;2.大連理工大學(xué) 控制科學(xué)與工程學(xué)院, 遼寧 大連 116024 )

    ?

    基于慣性測(cè)量器件的無(wú)線(xiàn)步態(tài)分析平臺(tái)

    姜 鳴*1,仇 森2,趙 洋1,張 志1

    ( 1.東莞理工學(xué)院 電子工程學(xué)院, 廣東 東莞523808;2.大連理工大學(xué) 控制科學(xué)與工程學(xué)院, 遼寧 大連116024 )

    步態(tài)定量測(cè)量方法應(yīng)用于許多領(lǐng)域,如臨床醫(yī)學(xué)、雙足機(jī)器人控制等.采用慣性測(cè)量單元結(jié)合無(wú)線(xiàn)傳感器網(wǎng)絡(luò)建立了一個(gè)步態(tài)分析平臺(tái),將兩個(gè)無(wú)線(xiàn)慣性測(cè)量單元傳感器節(jié)點(diǎn)分別綁定在左右雙側(cè)腳踝,以同時(shí)采集雙腳運(yùn)動(dòng)過(guò)程中的加速度和角速度信號(hào),并將其通過(guò)無(wú)線(xiàn)方式發(fā)送到遠(yuǎn)程終端.通過(guò)模式識(shí)別、時(shí)間序列分析、閾值檢測(cè)和零速修正等多種數(shù)據(jù)融合方法計(jì)算步態(tài)參數(shù),并通過(guò)融合雙足傳感器數(shù)據(jù)得到雙支撐相、雙腳步行周期等重要的雙足步態(tài)參數(shù),其中雙支撐相參數(shù)對(duì)人體日常動(dòng)作的識(shí)別有重要意義.實(shí)驗(yàn)結(jié)果顯示該研究具有較高的計(jì)算精度.

    步態(tài)分析;傳感器網(wǎng)絡(luò);數(shù)據(jù)融合;步行周期劃分;雙腳運(yùn)動(dòng)參數(shù)

    0 引 言

    步態(tài)分析主要研究人步行的機(jī)制,通過(guò)計(jì)算得到步態(tài)參數(shù).步態(tài)分析對(duì)臨床醫(yī)學(xué)、雙足機(jī)器人步態(tài)規(guī)劃等都有重要的意義.目前步態(tài)分析可采用多種方法.傳統(tǒng)臨床測(cè)量借助秒表和米尺,使用目測(cè)以及足印法得到步態(tài)信息,這種方法難以準(zhǔn)確定量劃分步行周期,人為因素造成的誤差大.文獻(xiàn)[1-3]使用視頻處理方法分析步態(tài),通過(guò)記錄人走路過(guò)程中的視頻圖像計(jì)算步態(tài)參數(shù).這類(lèi)步態(tài)識(shí)別方法需要多臺(tái)攝像機(jī)從不同角度同時(shí)記錄數(shù)據(jù),圖像分析算法復(fù)雜,采集的視頻數(shù)據(jù)量龐大,同時(shí)環(huán)境光線(xiàn)容易影響分析結(jié)果.文獻(xiàn)[4-6]使用鞋底放置的壓力墊,通過(guò)采集步行過(guò)程中的壓力變化來(lái)計(jì)算步態(tài)周期.由于壓力信號(hào)僅能反映腳與地面接觸時(shí)間段的壓力變化規(guī)律,這類(lèi)方法不能很好地跟蹤連續(xù)運(yùn)動(dòng),計(jì)算的步態(tài)參數(shù)不全面.文獻(xiàn)[7]建立了一種大型模擬行走系統(tǒng),用來(lái)分析人步態(tài)的一般規(guī)律,這種大型步態(tài)分析平臺(tái)非常復(fù)雜且昂貴,使用者需要經(jīng)過(guò)專(zhuān)業(yè)培訓(xùn).

    隨著慣性系統(tǒng)向高精度、高集成度和低成本發(fā)展,更多研究開(kāi)始使用慣性系統(tǒng)進(jìn)行步態(tài)分析,通過(guò)加速度信號(hào)和角速度信號(hào)來(lái)分析步態(tài)規(guī)律以及計(jì)算步態(tài)參數(shù).文獻(xiàn)[8]將一個(gè)慣性傳感器綁定在人的右腳上,采集步行過(guò)程中右腳的角速度和加速度信號(hào),并計(jì)算步態(tài)參數(shù).文獻(xiàn)[9]結(jié)合壓力傳感器和陀螺儀對(duì)步態(tài)周期進(jìn)行劃分.目前采用慣性系統(tǒng)來(lái)測(cè)量步態(tài)運(yùn)動(dòng)存在一些不足.首先,慣性器件的固定方式與分析結(jié)果緊密相關(guān),綁定位置設(shè)計(jì)不合適會(huì)極大降低計(jì)算結(jié)果準(zhǔn)確性.另外,在人步行過(guò)程中,被測(cè)物體(小腿、腳背、腳踝等)的慣性坐標(biāo)系與世界坐標(biāo)系的相對(duì)位置不停發(fā)生變化[8],給加速度和角速度的數(shù)據(jù)融合方法帶來(lái)挑戰(zhàn).再有,目前多數(shù)基于慣性系統(tǒng)的步態(tài)分析方法都采用一個(gè)慣性器件綁定在一只腳上[8-12],由于步行過(guò)程是兩只腳配合的協(xié)調(diào)運(yùn)動(dòng),單腳運(yùn)動(dòng)規(guī)律不能完全反映人體步態(tài)信息.文獻(xiàn)[13]使用兩個(gè)加速度計(jì)分析左右兩個(gè)腳的步態(tài),但是由于沒(méi)有使用陀螺儀,該方法缺乏有效數(shù)據(jù)融合過(guò)程,步態(tài)信息不完整.總之,目前多數(shù)步態(tài)分析的研究都是采用單一慣性節(jié)點(diǎn)分析單腳運(yùn)動(dòng),缺乏分析雙腳協(xié)調(diào)運(yùn)動(dòng)規(guī)律.為了解決這個(gè)問(wèn)題,本文采用兩組慣性測(cè)量器件,分別綁定在左右兩側(cè)腳踝,采集步行過(guò)程中的加速度信號(hào)和角速度信號(hào),通過(guò)多傳感器數(shù)據(jù)融合算法計(jì)算出雙足步態(tài)數(shù)據(jù).

    1 測(cè)量裝置

    本文使用兩個(gè)慣性測(cè)量單元(inertial measurement unit,IMU),分別采集左右腳的慣性力.采用MSP430F149對(duì)IMU進(jìn)行數(shù)據(jù)采集,CPU時(shí)鐘頻率為1 MHz,使用無(wú)線(xiàn)射頻芯片CC2420將采集到的數(shù)據(jù)通過(guò)2.4 GHz無(wú)線(xiàn)信道與接收節(jié)點(diǎn)通信,將步行過(guò)程中雙腳的加速度和角速度數(shù)據(jù)傳回PC端,在PC端計(jì)算步態(tài)參數(shù).IMU傳感器節(jié)點(diǎn)選取一個(gè)三軸加速度計(jì)(LIS344ALH)和一個(gè)單軸陀螺儀(LY530AL),器件性能指標(biāo)如表1所示.

    表1 慣性器件的性能指標(biāo)

    加速度計(jì)和陀螺儀通過(guò)12位AD轉(zhuǎn)換得到數(shù)字量,采樣前使用低通濾波器濾去電源帶來(lái)的50 Hz頻率的干擾.慣性器件校正時(shí),將加速度計(jì)的每個(gè)軸與地面垂直,通過(guò)校正使其輸出+1g或-1g加速度量.使用陀螺儀專(zhuān)用校正控制轉(zhuǎn)臺(tái),將陀螺儀分別轉(zhuǎn)過(guò)90°、180°、270°和360°,記錄陀螺儀輸出曲線(xiàn),通過(guò)校正使其積分結(jié)果得到正確角度.

    實(shí)驗(yàn)中使用特制繃帶將兩個(gè)IMU傳感器節(jié)點(diǎn)綁定在志愿者左右兩個(gè)腳踝上,固定方式如圖1所示.在處理綁定誤差的問(wèn)題上使用了文獻(xiàn)[14]和[15]中所提到的傳感器初始對(duì)準(zhǔn)方法.

    在器件固定過(guò)程中,從PC端實(shí)時(shí)監(jiān)視器件輸出,使加速度計(jì)X軸與重力方向相反,輸出加速度量為-1g,并保證慣性坐標(biāo)系與矢狀平面平行.根據(jù)本文設(shè)計(jì)的慣性器件結(jié)構(gòu),這時(shí)陀螺儀測(cè)量軸與矢狀平面垂直.

    圖1 IMU傳感器節(jié)點(diǎn)的綁定方式

    2 數(shù)據(jù)分析

    本文使用圖2所示的流程計(jì)算單腳和雙腳的步態(tài)參數(shù).傳感器節(jié)點(diǎn)綁定在使用者的腳踝關(guān)節(jié)處,用于測(cè)試使用者行走過(guò)程中腳部的三維加速度和三維角速度.同時(shí)上位機(jī)控制使用者行走過(guò)程中傳感器數(shù)據(jù)采集的開(kāi)始和結(jié)束,采集完成后,通過(guò)數(shù)據(jù)濾波和校準(zhǔn),使用基于多傳感器數(shù)據(jù)融合的步態(tài)分析算法,其中包括融合加速度和角速度數(shù)據(jù)、融合雙足的傳感器步態(tài)數(shù)據(jù);然后分別計(jì)算出單足步態(tài)參數(shù)和雙足步態(tài)參數(shù),其中包括時(shí)間參數(shù)和空間參數(shù).

    圖2 步態(tài)數(shù)據(jù)處理流程

    2.1加速度與角速度數(shù)據(jù)融合

    (1)

    2.1.2坐標(biāo)系轉(zhuǎn)換步行過(guò)程中,慣性器件與地面參考位置一直發(fā)生著變化,如圖3所示.由于步態(tài)信息包含在世界坐標(biāo)系內(nèi)(圖3中的矢狀平面),需要將慣性系下的數(shù)值投影到世界坐標(biāo)系下,從而計(jì)算出步長(zhǎng)等步態(tài)參數(shù)[8].此處需要說(shuō)明的是,由于主要的空間步態(tài)參數(shù),比如步長(zhǎng)、步速和足角都體現(xiàn)在矢狀平面,本文的研究重點(diǎn)就是在人體矢狀平面內(nèi)的運(yùn)動(dòng)信息.

    圖3 走動(dòng)過(guò)程中的兩個(gè)坐標(biāo)系關(guān)系

    圖3中角度θ表示兩個(gè)坐標(biāo)系的偏差角度,由于陀螺儀測(cè)量軸與矢狀平面垂直,θ可由角速度計(jì)算得到,計(jì)算過(guò)程詳見(jiàn)2.1.4節(jié).根據(jù)θ(t)通過(guò)式(2)所示的線(xiàn)性轉(zhuǎn)換將加速度計(jì)輸出的x軸和y軸加速度信號(hào)ax(t)和ay(t)投影到世界坐標(biāo)系的Y軸內(nèi),得到步行方向上的加速度信號(hào)aY(t):

    aY(t)=ay(t)cosθ(t)+ax(t)sinθ(t)

    (2)

    2.1.3單足步行時(shí)相劃分一個(gè)步態(tài)周期分為擺動(dòng)相周期Tsw和站立相周期Tsp,擺動(dòng)相是腳與地面不接觸的時(shí)間,站立相是腳與地面接觸的時(shí)間.站立相又分為3個(gè)階段,分別是腳跟離地Tho、腳跟擊地Ths,以及完全站立相Tff.在2.1.1節(jié)得到步行周期劃分的基礎(chǔ)上,根據(jù)Y軸加速度信號(hào)進(jìn)一步對(duì)單足步行時(shí)相做劃分.本文采用滑動(dòng)方差方法分別使用式(3)和(4)計(jì)算y軸加速度在滑動(dòng)窗口內(nèi)的均值和方差:

    (3)

    (4)

    其中i是采樣序列號(hào);M是窗口大小,本文中M=30.圖4所示是實(shí)際采樣過(guò)程中的一個(gè)步行周期.

    圖4 一個(gè)步行周期內(nèi)步行時(shí)相的劃分

    每一步開(kāi)始時(shí),在滑動(dòng)方差大于0.01的情況下,認(rèn)為是腳跟離地時(shí)刻(tho);當(dāng)角速度第一次變正時(shí),認(rèn)為是擺動(dòng)相開(kāi)始的時(shí)間點(diǎn)(tsw);由于腳跟著地時(shí)y軸加速度會(huì)產(chǎn)生明顯脈沖,在發(fā)生脈沖時(shí)認(rèn)為是腳跟擊地的時(shí)刻(ths);當(dāng)滑動(dòng)方差小于0.01時(shí),認(rèn)為是完全站立開(kāi)始時(shí)刻(tff).一個(gè)步行周期內(nèi)的步行時(shí)相由式(5)~(8)來(lái)計(jì)算:

    Tho(i)=tsw(i)-tho(i)

    (5)

    Tsw(i)=ths(i)-tsw(i)

    (6)

    Ths(i)=tff(i)-ths(i)

    (7)

    (8)

    2.1.4角度的計(jì)算步行過(guò)程中單足轉(zhuǎn)過(guò)的角度θ(t)可以根據(jù)下式計(jì)算:

    θ(t)=∫ttho(i)θ.(τ)dτ+θinit(i)

    (9)

    初始角度θinit通過(guò)式(10)計(jì)算,本文中令N=10,方括號(hào)代表離散點(diǎn)的序列號(hào).

    (10)

    根據(jù)式(9),一個(gè)步行周期內(nèi)的足俯角θpitch可由下式計(jì)算:

    θpitch(i)=∫ths(i)tho(i)θ.(t)dt+θinit(i)

    (11)

    本文中規(guī)定順時(shí)針為負(fù)方向,逆時(shí)針為正方向,因此在腳跟離地期發(fā)生的負(fù)角度由式(12)計(jì)算,在擺動(dòng)期發(fā)生的正角度由式(13)計(jì)算:

    θpitch(-)=∫tsw(i)tho(i)θ.Z(t)dt+θinit(i)

    (12)

    (13)

    2.1.5步伐長(zhǎng)度計(jì)算在一個(gè)步行周期內(nèi),通過(guò)對(duì)步行方向上(世界坐標(biāo)系內(nèi))的加速度信號(hào)aY(t) 進(jìn)行積分可以得到步伐長(zhǎng)度,如式(14)所示,[ ]表示積分是數(shù)值積分,使用梯形法積分:

    (14)

    2.2雙腳傳感器數(shù)據(jù)融合

    人直立行走是身體協(xié)調(diào)運(yùn)動(dòng)的過(guò)程,需要左右腳協(xié)調(diào)配合.已有研究大多使用單個(gè)節(jié)點(diǎn)采集單腳步行數(shù)據(jù),僅能分析單腳步態(tài)信息.本文設(shè)計(jì)兩個(gè)節(jié)點(diǎn)同步采集左右腳的步行數(shù)據(jù),可以分析雙腳步行過(guò)程的相關(guān)性,反映兩腳協(xié)調(diào)運(yùn)動(dòng)的規(guī)律.兩個(gè)腳協(xié)調(diào)運(yùn)動(dòng)時(shí),步態(tài)周期分為單支撐相(single stand phase,SSP)、雙支撐相(dual stand phase,DSP),以及擺動(dòng)相(swing phase,SW).雙支撐相是指兩個(gè)腳同時(shí)與地面接觸,單支撐相是指只有一只腳與地面接觸.從左腳腳跟離地到開(kāi)始擺動(dòng)相的期間是左腳前雙支撐相;隨后左腳發(fā)生擺動(dòng)相,同時(shí)右腳處于單支撐相;然后左腳腳跟著地,并且身體重心向前轉(zhuǎn)變右腳腳跟離地,這個(gè)期間是第二次雙支撐相;最后左腳處于單支撐相,同時(shí)右腳離地發(fā)生擺動(dòng)相.根據(jù)這個(gè)規(guī)律可以得到雙腳的步行時(shí)相劃分,如圖5所示.通過(guò)計(jì)算雙支撐相可以有效區(qū)分“走路”和“跑步”兩種動(dòng)作,對(duì)人體運(yùn)動(dòng)識(shí)別具有重要作用.

    圖5 一個(gè)步態(tài)周期內(nèi)的雙足步行時(shí)相劃分

    3 實(shí)驗(yàn)過(guò)程和實(shí)驗(yàn)結(jié)果

    實(shí)驗(yàn)中采集了6名志愿者(3男3女)的步態(tài)數(shù)據(jù),并對(duì)他們的步態(tài)數(shù)據(jù)進(jìn)行分析.圖6所示為某志愿者正常行走時(shí)的左腳x軸加速度、y軸加速度,以及角速度的信號(hào)曲線(xiàn).從圖6中可以看到,x軸方向有一個(gè)初始1g的重力分量,隨后數(shù)值累計(jì)變化2g左右,這個(gè)變化反映出人腳在垂直方向上的運(yùn)動(dòng).y軸(慣性坐標(biāo)系)與步行方向直接相關(guān),y軸加速度從-2g變化到2g,變化范圍4g左右.從圖6中還可以看到,在每一步腳跟擊地時(shí),y軸加速度在短時(shí)間內(nèi)(不超過(guò)0.01 s)有一個(gè)劇烈的變化,可以從2g直接躍變到-2g.角速度在擺動(dòng)期達(dá)到峰值,大概在300°/s.對(duì)比加速度信號(hào),角速度信號(hào)變化比較緩和,在腳跟擊地時(shí)沒(méi)有過(guò)大的沖擊.通過(guò)實(shí)驗(yàn)數(shù)據(jù)可以看出本文設(shè)計(jì)的步行周期劃分方法的有效性.文獻(xiàn)[16]將慣性器件固定在腳背上,腳跟擊地的加速度變化明顯小于本文設(shè)計(jì)方法,本文方法更能準(zhǔn)確反映步行過(guò)程中腳跟的運(yùn)動(dòng)情況.

    (a) x軸加速度信號(hào)

    (b)y軸加速度信號(hào)

    (c) 角速度信號(hào)

    圖6某志愿者左腳步態(tài)信號(hào)

    Fig.6Gait signal of left foot from a volunteer

    圖7所示是某志愿者正常步行時(shí)采集的左右腳的角速度曲線(xiàn),從中可以看到兩個(gè)腳交替運(yùn)動(dòng)的協(xié)調(diào)過(guò)程.

    表2所示是使用本文所提方法計(jì)算得到的某志愿者的單腳步行時(shí)相,本文給出了從起步開(kāi)始到第10步的數(shù)據(jù),圖8所示為每個(gè)步行時(shí)相占整個(gè)步態(tài)周期的百分比R.通過(guò)表2和圖8可以發(fā)現(xiàn),左右腳時(shí)相基本對(duì)稱(chēng),右腳擺動(dòng)相比左腳稍大一點(diǎn),而腳跟離地稍小一點(diǎn).本文計(jì)算結(jié)果與文獻(xiàn)[17-19]相比基本一致,但是本文中Ths所占百分比偏小(13%),而Tho偏大(25%).之所以存在差異,首先是因?yàn)楸疚闹袘T性器件綁定的位置更靠近腳踝;其次不同志愿者的步態(tài)不完全一樣,因此產(chǎn)生這種差異也是合理的.

    圖7 某志愿者雙腳步態(tài)信號(hào)

    表2 某志愿者步行時(shí)相計(jì)算結(jié)果

    圖8 左腳步行時(shí)相占步行周期的百分比

    表3所示是使用本文所提方法計(jì)算得到的某志愿者的步長(zhǎng)和足俯角的計(jì)算結(jié)果.從表3中可以看到,志愿者平均步長(zhǎng)為1.204±0.219 m;足俯角平均負(fù)角度-38.474°±7.878°,這個(gè)角度反映了在矢狀平面內(nèi)腳跟向后抬起的角度;平均正角度為65.686°±8.711°,正角度反映了走路的過(guò)程中主要在擺動(dòng)相期間腳逆時(shí)針擺動(dòng)的角度.值得注意的是第一步的步態(tài)參數(shù)對(duì)比后面參數(shù)差異明顯,例如第一步的負(fù)角度和正角度只有-17.381° 和36.295°,而其他步的參數(shù)都在-35°和60°附近.產(chǎn)生這個(gè)現(xiàn)象是因?yàn)橹驹刚咴陂_(kāi)始采集數(shù)據(jù)時(shí)是靜止的,本文給出的平均值統(tǒng)計(jì)中沒(méi)有包括第一步起步的數(shù)據(jù).

    表3 某志愿者步長(zhǎng)和足俯角的計(jì)算結(jié)果

    4 結(jié) 語(yǔ)

    本文提出的雙足步態(tài)分析系統(tǒng)被用來(lái)測(cè)量健康人的正常步態(tài),通過(guò)無(wú)線(xiàn)傳感器節(jié)點(diǎn)采集到的數(shù)據(jù),再使用上面提到的算法分析計(jì)算得到正常人的步態(tài)參數(shù).為了對(duì)測(cè)量結(jié)果的準(zhǔn)確性進(jìn)行一定的評(píng)估,本文在實(shí)際測(cè)量中使用了其他方法來(lái)測(cè)量對(duì)比系統(tǒng)的計(jì)算結(jié)果.比如在實(shí)際測(cè)量時(shí)采用秒表對(duì)步行周期的測(cè)量做了簡(jiǎn)單的對(duì)比,采用視頻的方法來(lái)估算俯仰角的大小,以及用米尺測(cè)量人的步長(zhǎng)大小.對(duì)比分析表明本文系統(tǒng)分析計(jì)算出來(lái)的參數(shù)具有一定的精度,尤其是在步態(tài)周期的測(cè)量上有很大的優(yōu)越性.并且比起其他的方法,本文提出的系統(tǒng)測(cè)量方法簡(jiǎn)單,容易被應(yīng)用到實(shí)際測(cè)量中.

    [1]ZHANG Xin, Fan G. Dual gait generative models for human motion estimation from a single camera [J]. IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, 2010, 40(4):1034-1049.

    [2]LIU Zong-yi, Sarkar S. Improved gait recognition by gait dynamics normalization [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2006, 28(6):863-876.

    [3]ZHANG Rong, Vogler C, Metaxas D. Human gait recognition [C]// Conference on Computer Vision and Pattern Recognition Workshop, 2004. CVPRW ′04. Piscataway:IEEE, 2004.

    [4]Bae Joon-bum, Kong Kyoung-chul, Nancy B,etal. A mobile gait monitoring system for gait analysis [C] // 2009 IEEE International Conference on Rehabilitation Robotics, ICORR 2009. Kyoto:IEEE Computer Society, 2009.

    [5]Bamberg S J M, Carson R J, Stoddard G,etal. The lower extremity ambulation feedback system for analysis of gait asymmetries:preliminary design and validation results [J]. JPO Journal of Prosthetics and Orthotics, 2010, 22(1):31-36.

    [6]Cuadrado J, Rosa P-V, Urbano L,etal. A force-based approach for joint efforts estimation during the double support phase of gait [J]. Procedia IUTAM, 2011, 2:26-34.

    [7]Shiozawa N, Arima S, Makikawa M. Virtual walkway system and prediction of gait mode transition for the control of the gait simulator [C] // 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Piscataway:IEEE, 2004:2699-2702.

    [8]Sabatini A M, Martelloni C, Scapellato S,etal. Assessment of walking features from foot inertial sensing [J]. IEEE Transactions on Biomedical Engineering, 2005, 52(3):486-494.

    [9]Pappas I P, Popovic M R, Keller T,etal. A reliable gait phase detection system [J]. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2001, 9(2):113-125.

    [10]Veltink P H, Slycke P, Hemssems J,etal. Three dimensional inertial sensing of foot movements for automatic tuning of a two-channel implantable drop-foot stimulator [J]. Medical Engineering & Physics, 2003, 25(1):21-28.

    [11]O′donovan K J, Kamnik R, O′keeffe D T,etal. An inertial and magnetic sensor based technique for joint angle measurement [J]. Journal of Biomechanics, 2007, 40(12):2604-2611.

    [12]WEN Shi-guang, WANG Fei, WU Cheng-dong. Realtime gait kinematics classification using LDA and SVM [C] // Proceedings of the 2011 Chinese Control and Decision Conference (CCDC). Shenyang:IEEE, 2011:592-595.

    [13]Lee J A, Cho S H, Lee J W,etal. Wearable accelerometer system for measuring the temporal parameters of gait [C] // 2007 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vols 1-16. Piscataway:IEEE, 2007:483-486.

    [14]JIANG Ming, SHANG Hong, WANG Zhe-long,etal. A method to deal with installation errors of wearable accelerometers for human activity recognition [J]. Physiological Measurement, 2011, 32(3):347-358.

    [15]QIU Sen, WANG Zhe-long, ZHAO Hong-yu,etal. Using distributed wearable sensors to measure and evaluate human lower limb motions [J]. IEEE Transactions on Instrumentation and Measurement, 2016, 65(4):939-950.

    [16]Alvarez J C, González R C, Alvarez D,etal. Multisensor approach to walking distance estimation with foot inertial sensing [C] // 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Piscataway:IEEE, 2007.

    [17]QIU Sen, YANG Yu, HOU Ji-jian,etal. Ambulatory estimation of 3D walking trajectory and knee joint angle using MARG sensors [C] // 2014 Fourth International Conference on Innovative Computing Technology (Intech). Piscataway:IEEE, 2014:191-196.

    [18]WANG Zhe-long, QIU Shen, CAO Zhong-kai,etal. Quantitative assessment of dual gait analysis based on inertial sensors with body sensor network [J]. Sensor Review, 2013, 33(1):48-56.

    [19]Godha S, Lachapelle G. Foot mounted inertial system for pedestrian navigation [J]. Measurement Science and Technology, 2008, 19(7):1-7.

    Wireless gait analysis platform based on inertial measurement device

    JIANGMing*1,QIUSen2,ZHAOYang1,ZHANGZhi1

    ( 1.School of Electronic Engineering, Dongguan University of Technology, Dongguan 523808, China;2.School of Control Science and Engineering, Dalian University of Technology, Dalian 116024, China )

    The method of quantitative measurement of gait makes a profound influence on many fields, such as clinical medicine, biped robot control strategy and so on. A gait analysis platform is established based on inertial measurement unit (IMU) combining wireless sensor network. Two IMU sensor nodes are bound on the left and right ankles so that the acceleration and angular velocity can be obtained from both sides at the same time and sent to remote terminal. Based on pattern recognition, time series analysis, threshold detection, zero velocity correction and various data fusion methods, the gait parameters are calculated. Furthermore, dual gait parameters, such as the dual stand phase and dual gait period can be obtained by the proposed sensor data fusion method. Dual stand phase detection has an important significance to the human activity recognition. The experimental results show that this study achieves high calculation accuracy.

    gait analysis; sensor network; data fusion; gait phase division; dual gait parameters

    1000-8608(2016)05-0518-07

    2016-05-10;

    2016-07-28.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51407031);廣東省自然科學(xué)基金資助項(xiàng)目(2016A030313134);廣東省高等學(xué)?!皠?chuàng)新強(qiáng)校工程”創(chuàng)新項(xiàng)目(2014KQNCX221);東莞市社會(huì)科技發(fā)展項(xiàng)目(2013108101007).

    姜 鳴*(1982-),男,博士,E-mail:jiangm@dgut.edu.cn.

    TH82

    A

    10.7511/dllgxb201605012

    猜你喜歡
    腳跟步態(tài)角速度
    小螞蟻與“三角步態(tài)”
    羊的吻
    遼河(2019年7期)2019-09-03 05:13:09
    基于面部和步態(tài)識(shí)別的兒童走失尋回系統(tǒng)
    電子制作(2018年18期)2018-11-14 01:48:04
    基于Kinect的學(xué)步期幼兒自然步態(tài)提取
    “腳跟疼”小心足底筋膜炎
    拓展游戲
    圓周運(yùn)動(dòng)角速度測(cè)量方法賞析
    半捷聯(lián)雷達(dá)導(dǎo)引頭視線(xiàn)角速度提取
    基于構(gòu)架點(diǎn)頭角速度的軌道垂向長(zhǎng)波不平順在線(xiàn)檢測(cè)
    便宜假肢能模仿正常步態(tài)
    少妇熟女欧美另类| 久久精品国产鲁丝片午夜精品| 久久国产亚洲av麻豆专区| av电影中文网址| 韩国高清视频一区二区三区| 热99久久久久精品小说推荐| 亚洲精品视频女| 我的女老师完整版在线观看| 精品少妇内射三级| 熟妇人妻不卡中文字幕| 免费av不卡在线播放| 两性夫妻黄色片 | 欧美 日韩 精品 国产| 一区二区日韩欧美中文字幕 | 搡老乐熟女国产| 久久久久网色| 久久久精品免费免费高清| 亚洲欧美一区二区三区黑人 | 日本欧美国产在线视频| 国产精品久久久av美女十八| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 高清不卡的av网站| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 亚洲伊人久久精品综合| 国产在视频线精品| 久久精品熟女亚洲av麻豆精品| 欧美xxⅹ黑人| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产色片| 赤兔流量卡办理| 天天操日日干夜夜撸| 18禁裸乳无遮挡动漫免费视频| av国产久精品久网站免费入址| 午夜精品国产一区二区电影| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 最后的刺客免费高清国语| 99久久人妻综合| 国产白丝娇喘喷水9色精品| 中文字幕人妻熟女乱码| 久久 成人 亚洲| 欧美 日韩 精品 国产| 成年人免费黄色播放视频| 街头女战士在线观看网站| 黑人猛操日本美女一级片| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| 亚洲精品自拍成人| 亚洲情色 制服丝袜| 一边摸一边做爽爽视频免费| 十八禁网站网址无遮挡| 91aial.com中文字幕在线观看| 精品第一国产精品| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 哪个播放器可以免费观看大片| 欧美日韩视频高清一区二区三区二| 激情视频va一区二区三区| 欧美老熟妇乱子伦牲交| 欧美人与性动交α欧美精品济南到 | 亚洲精品日韩在线中文字幕| 夫妻午夜视频| 妹子高潮喷水视频| 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 日本欧美国产在线视频| 亚洲欧美成人精品一区二区| 欧美人与性动交α欧美精品济南到 | 日韩av免费高清视频| 亚洲精品日本国产第一区| 精品一区二区三卡| 国产精品偷伦视频观看了| 精品第一国产精品| 老熟女久久久| 边亲边吃奶的免费视频| 国产精品女同一区二区软件| 欧美成人午夜免费资源| 熟女电影av网| 色哟哟·www| 欧美xxxx性猛交bbbb| av视频免费观看在线观看| 一级,二级,三级黄色视频| 国产精品一区二区在线不卡| 中文字幕亚洲精品专区| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 亚洲国产日韩一区二区| 中国三级夫妇交换| 一级黄片播放器| 97超碰精品成人国产| 一区二区三区精品91| 久久综合国产亚洲精品| 国产深夜福利视频在线观看| 日本黄大片高清| 国产精品免费大片| 婷婷成人精品国产| 日本免费在线观看一区| av天堂久久9| 热re99久久精品国产66热6| 捣出白浆h1v1| 日本vs欧美在线观看视频| 我要看黄色一级片免费的| 日本欧美视频一区| 午夜av观看不卡| 免费观看性生交大片5| 久久人人97超碰香蕉20202| 夫妻性生交免费视频一级片| 少妇的逼好多水| 亚洲精品第二区| 女人精品久久久久毛片| 极品人妻少妇av视频| 亚洲精品久久午夜乱码| 亚洲欧美中文字幕日韩二区| 久久精品国产a三级三级三级| 亚洲国产精品999| 精品人妻偷拍中文字幕| 免费不卡的大黄色大毛片视频在线观看| 91在线精品国自产拍蜜月| 97精品久久久久久久久久精品| 精品亚洲成a人片在线观看| 欧美激情极品国产一区二区三区 | 97精品久久久久久久久久精品| 丝袜喷水一区| 色94色欧美一区二区| 日韩中字成人| 在线观看国产h片| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 九色成人免费人妻av| 久久人人爽人人片av| 欧美激情国产日韩精品一区| 亚洲av成人精品一二三区| 国内精品宾馆在线| 成年人免费黄色播放视频| 国产成人免费观看mmmm| 男人爽女人下面视频在线观看| 国产精品人妻久久久久久| av免费观看日本| 国产免费一区二区三区四区乱码| 赤兔流量卡办理| 免费观看a级毛片全部| 久久午夜福利片| 精品一区二区三卡| 免费观看a级毛片全部| 人妻人人澡人人爽人人| 深夜精品福利| 视频区图区小说| 春色校园在线视频观看| 日韩在线高清观看一区二区三区| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 蜜桃国产av成人99| 高清黄色对白视频在线免费看| 久久ye,这里只有精品| 三上悠亚av全集在线观看| 欧美日韩综合久久久久久| 久久人人爽人人爽人人片va| 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠久久av| 午夜免费鲁丝| 日韩视频在线欧美| 欧美日韩视频精品一区| 久久国内精品自在自线图片| 欧美+日韩+精品| 精品久久久久久电影网| 欧美成人午夜精品| 日本免费在线观看一区| 亚洲 欧美一区二区三区| 丝袜喷水一区| 日韩中字成人| 国产精品熟女久久久久浪| 国产av一区二区精品久久| 美女主播在线视频| 水蜜桃什么品种好| 五月开心婷婷网| 中文精品一卡2卡3卡4更新| 韩国av在线不卡| 女人久久www免费人成看片| 久久久国产一区二区| 丰满饥渴人妻一区二区三| 亚洲人成77777在线视频| 免费久久久久久久精品成人欧美视频 | 国产午夜精品一二区理论片| 美女国产视频在线观看| 国产日韩欧美视频二区| 久久久久久久久久久免费av| 青春草亚洲视频在线观看| 久久精品国产亚洲av天美| 午夜福利网站1000一区二区三区| 国产成人免费无遮挡视频| 美女中出高潮动态图| 日本色播在线视频| 国产精品99久久99久久久不卡 | 亚洲,欧美精品.| 99热网站在线观看| 蜜桃国产av成人99| 男女边摸边吃奶| 91精品伊人久久大香线蕉| 精品一区二区三卡| 亚洲av.av天堂| 少妇熟女欧美另类| 男女下面插进去视频免费观看 | 中文字幕人妻丝袜制服| 在线观看美女被高潮喷水网站| 国产亚洲av片在线观看秒播厂| 1024视频免费在线观看| 高清欧美精品videossex| 亚洲精品一二三| 少妇人妻久久综合中文| 女的被弄到高潮叫床怎么办| 国产精品.久久久| 国产极品天堂在线| 插逼视频在线观看| 国产精品 国内视频| 午夜激情久久久久久久| 久久精品夜色国产| 日韩一区二区视频免费看| 这个男人来自地球电影免费观看 | 男女无遮挡免费网站观看| 亚洲av日韩在线播放| 国产国拍精品亚洲av在线观看| 成人无遮挡网站| av黄色大香蕉| 亚洲人成77777在线视频| av片东京热男人的天堂| 99久久综合免费| 久久久久久久国产电影| 人妻一区二区av| 18在线观看网站| 在线观看国产h片| 久久久久精品性色| 青春草视频在线免费观看| 男女啪啪激烈高潮av片| 精品一区在线观看国产| 国产黄色视频一区二区在线观看| 日日啪夜夜爽| 欧美最新免费一区二区三区| 香蕉国产在线看| 多毛熟女@视频| 国语对白做爰xxxⅹ性视频网站| 国产亚洲av片在线观看秒播厂| 99视频精品全部免费 在线| 精品久久蜜臀av无| 九色亚洲精品在线播放| 国内精品宾馆在线| 国产精品一二三区在线看| 黄色配什么色好看| 一区二区av电影网| 人妻人人澡人人爽人人| 伊人久久国产一区二区| 在线看a的网站| 91久久精品国产一区二区三区| 亚洲欧美日韩另类电影网站| 亚洲精品一二三| 只有这里有精品99| 91久久精品国产一区二区三区| 如何舔出高潮| 欧美激情 高清一区二区三区| 最近最新中文字幕大全免费视频 | 成人影院久久| 天堂俺去俺来也www色官网| 九草在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩卡通动漫| 另类精品久久| 18+在线观看网站| 国产成人精品久久久久久| av免费在线看不卡| 国产亚洲精品久久久com| 一级爰片在线观看| 汤姆久久久久久久影院中文字幕| 免费观看性生交大片5| 欧美 日韩 精品 国产| 国产高清不卡午夜福利| 亚洲国产av影院在线观看| 午夜91福利影院| 久久精品国产亚洲av天美| 亚洲成人av在线免费| 日韩大片免费观看网站| 欧美日韩国产mv在线观看视频| 国产色婷婷99| 欧美日韩视频精品一区| 有码 亚洲区| 久久久久久久国产电影| 永久免费av网站大全| 欧美性感艳星| 九九爱精品视频在线观看| 99精国产麻豆久久婷婷| 各种免费的搞黄视频| 最后的刺客免费高清国语| 国产午夜精品一二区理论片| 欧美人与性动交α欧美精品济南到 | 免费高清在线观看视频在线观看| 精品第一国产精品| 亚洲精品国产av蜜桃| 在线看a的网站| 亚洲av电影在线观看一区二区三区| 中文字幕制服av| 你懂的网址亚洲精品在线观看| 日韩大片免费观看网站| 欧美人与性动交α欧美软件 | 国产一区亚洲一区在线观看| 在线观看www视频免费| 看十八女毛片水多多多| 九色成人免费人妻av| av不卡在线播放| 久久狼人影院| 黄片无遮挡物在线观看| 久久午夜福利片| 久久久久久久亚洲中文字幕| 99国产综合亚洲精品| 亚洲美女视频黄频| 亚洲国产精品一区三区| 人妻系列 视频| 晚上一个人看的免费电影| 久久婷婷青草| 成年人免费黄色播放视频| 边亲边吃奶的免费视频| 国产精品三级大全| 欧美日韩av久久| av视频免费观看在线观看| 这个男人来自地球电影免费观看 | 久久狼人影院| 精品福利永久在线观看| 国产欧美另类精品又又久久亚洲欧美| videossex国产| 国产深夜福利视频在线观看| 成人亚洲精品一区在线观看| 中文字幕人妻熟女乱码| 午夜福利视频精品| 高清毛片免费看| 亚洲人成77777在线视频| 久久这里只有精品19| 国产精品久久久久久av不卡| 久久久久久人人人人人| 成人免费观看视频高清| 国产精品秋霞免费鲁丝片| 国产免费视频播放在线视频| 久久久久久人妻| 日本黄色日本黄色录像| 草草在线视频免费看| 丰满少妇做爰视频| 在线观看免费视频网站a站| 亚洲激情五月婷婷啪啪| 精品第一国产精品| 中文字幕免费在线视频6| 日本-黄色视频高清免费观看| 夫妻午夜视频| 欧美成人精品欧美一级黄| 久久久久国产网址| 99热网站在线观看| 老司机亚洲免费影院| 宅男免费午夜| 视频在线观看一区二区三区| 中文字幕最新亚洲高清| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 亚洲美女视频黄频| 春色校园在线视频观看| 国产淫语在线视频| 这个男人来自地球电影免费观看 | 午夜免费观看性视频| 只有这里有精品99| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 香蕉丝袜av| 亚洲国产欧美在线一区| 成人无遮挡网站| 亚洲国产毛片av蜜桃av| 国产成人精品久久久久久| 97在线视频观看| 久久久精品区二区三区| 国产免费现黄频在线看| 久久毛片免费看一区二区三区| 一级毛片黄色毛片免费观看视频| 国产成人免费观看mmmm| 大香蕉97超碰在线| 天天影视国产精品| 高清黄色对白视频在线免费看| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 国产精品久久久久久久电影| 久久久久久伊人网av| 免费黄色在线免费观看| 亚洲人成网站在线观看播放| 18在线观看网站| 欧美日韩视频高清一区二区三区二| 亚洲国产欧美日韩在线播放| 亚洲国产日韩一区二区| 精品国产一区二区久久| 欧美人与性动交α欧美软件 | 国产精品无大码| 9191精品国产免费久久| tube8黄色片| 91午夜精品亚洲一区二区三区| 国产综合精华液| 伊人久久国产一区二区| 欧美xxⅹ黑人| 欧美精品av麻豆av| 国产成人欧美| 国产白丝娇喘喷水9色精品| 亚洲国产欧美在线一区| 久久精品久久久久久久性| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 51国产日韩欧美| 九草在线视频观看| 久久久国产一区二区| 妹子高潮喷水视频| 99视频精品全部免费 在线| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的| 国产午夜精品一二区理论片| 欧美日韩视频精品一区| 18禁在线无遮挡免费观看视频| 亚洲伊人色综图| 国产永久视频网站| 午夜福利视频精品| 日韩精品免费视频一区二区三区 | 男女国产视频网站| 香蕉精品网在线| 中文字幕另类日韩欧美亚洲嫩草| 欧美另类一区| 中文字幕制服av| 纵有疾风起免费观看全集完整版| 国产日韩一区二区三区精品不卡| 丝瓜视频免费看黄片| 婷婷色麻豆天堂久久| 亚洲综合精品二区| 捣出白浆h1v1| xxx大片免费视频| 青春草国产在线视频| 国产成人精品在线电影| 国产欧美另类精品又又久久亚洲欧美| 国产精品99久久99久久久不卡 | 中文精品一卡2卡3卡4更新| 99热6这里只有精品| 日韩,欧美,国产一区二区三区| 国产精品人妻久久久影院| 一边摸一边做爽爽视频免费| 女性生殖器流出的白浆| 午夜激情久久久久久久| 国产精品久久久久久久久免| 日本免费在线观看一区| 人人澡人人妻人| 午夜福利视频在线观看免费| 国产女主播在线喷水免费视频网站| 欧美另类一区| 狠狠精品人妻久久久久久综合| 人体艺术视频欧美日本| 精品一区在线观看国产| 老熟女久久久| 丝瓜视频免费看黄片| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 婷婷成人精品国产| 国产成人精品在线电影| 九色成人免费人妻av| 婷婷色av中文字幕| 精品国产露脸久久av麻豆| 久久久久人妻精品一区果冻| 丝瓜视频免费看黄片| 夫妻性生交免费视频一级片| av一本久久久久| 高清毛片免费看| 日韩一本色道免费dvd| av国产久精品久网站免费入址| 久久精品熟女亚洲av麻豆精品| 亚洲色图 男人天堂 中文字幕 | 青春草国产在线视频| 人体艺术视频欧美日本| 久久影院123| 日日啪夜夜爽| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲国产看品久久| 免费播放大片免费观看视频在线观看| 一级毛片我不卡| 国产毛片在线视频| 国产一区亚洲一区在线观看| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| 青青草视频在线视频观看| 国产高清三级在线| 超色免费av| 欧美人与性动交α欧美软件 | 91久久精品国产一区二区三区| 亚洲av在线观看美女高潮| 免费高清在线观看视频在线观看| 亚洲精品av麻豆狂野| 男人添女人高潮全过程视频| 视频区图区小说| 亚洲国产精品成人久久小说| 女人久久www免费人成看片| 91久久精品国产一区二区三区| 亚洲国产看品久久| 国产成人免费观看mmmm| 国产精品久久久久久久久免| 美女大奶头黄色视频| 一级毛片我不卡| 国产在线视频一区二区| 亚洲内射少妇av| 18禁观看日本| 久久久久国产网址| 伦理电影免费视频| 夜夜爽夜夜爽视频| 黄片无遮挡物在线观看| 九色亚洲精品在线播放| 天堂中文最新版在线下载| 在线观看免费高清a一片| 高清视频免费观看一区二区| 一级a做视频免费观看| 国产在线一区二区三区精| 亚洲熟女精品中文字幕| 人妻少妇偷人精品九色| 在线看a的网站| 免费大片18禁| 久久久久久久久久成人| 最近中文字幕高清免费大全6| 老女人水多毛片| 人体艺术视频欧美日本| 咕卡用的链子| 成人免费观看视频高清| 我要看黄色一级片免费的| 边亲边吃奶的免费视频| 国产成人aa在线观看| 久久99精品国语久久久| 欧美成人午夜精品| 亚洲色图 男人天堂 中文字幕 | 国产精品人妻久久久久久| 国产高清国产精品国产三级| 成年动漫av网址| 精品久久国产蜜桃| 日韩一区二区视频免费看| 男女边吃奶边做爰视频| 成年女人在线观看亚洲视频| 日本猛色少妇xxxxx猛交久久| 午夜福利乱码中文字幕| 伦理电影免费视频| 午夜视频国产福利| 777米奇影视久久| 国产午夜精品一二区理论片| 天堂8中文在线网| 国产在线视频一区二区| 精品少妇久久久久久888优播| 91精品三级在线观看| 国产av国产精品国产| videosex国产| 又黄又粗又硬又大视频| 最近最新中文字幕大全免费视频 | 久久精品国产a三级三级三级| av片东京热男人的天堂| 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 成年女人在线观看亚洲视频| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 成人影院久久| 久久精品国产亚洲av涩爱| 久久国产亚洲av麻豆专区| 久久久久久久久久人人人人人人| 午夜福利视频精品| 亚洲成人av在线免费| 99久国产av精品国产电影| 国产女主播在线喷水免费视频网站| 精品一区二区三区视频在线| 汤姆久久久久久久影院中文字幕| 啦啦啦啦在线视频资源| 大香蕉久久成人网| 丰满饥渴人妻一区二区三| 99热国产这里只有精品6| 亚洲少妇的诱惑av| 国产精品人妻久久久影院| 亚洲人成网站在线观看播放| 久久97久久精品| 永久网站在线| 侵犯人妻中文字幕一二三四区| 日本与韩国留学比较| 纯流量卡能插随身wifi吗| 日韩大片免费观看网站| 久久毛片免费看一区二区三区| 啦啦啦中文免费视频观看日本| 肉色欧美久久久久久久蜜桃| 国产黄色视频一区二区在线观看| 国产精品一国产av| 国产色爽女视频免费观看| 欧美精品av麻豆av| 五月开心婷婷网| 国产精品熟女久久久久浪| 一级,二级,三级黄色视频| 久久久精品区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美一区视频在线观看| 视频中文字幕在线观看| 香蕉丝袜av| 欧美成人午夜精品| 人人妻人人澡人人看| 国产成人精品一,二区| 黑人欧美特级aaaaaa片| 美女脱内裤让男人舔精品视频| 97人妻天天添夜夜摸| 最黄视频免费看| 成人综合一区亚洲| 青青草视频在线视频观看| 香蕉丝袜av| 久久久久国产精品人妻一区二区|