• <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在线观看| 啦啦啦啦在线视频资源| 国产成人a∨麻豆精品| 男女国产视频网站| 午夜av观看不卡| 九九在线视频观看精品| 中文资源天堂在线| 免费少妇av软件| 日日摸夜夜添夜夜添av毛片| 91在线精品国自产拍蜜月| 狂野欧美激情性xxxx在线观看| 亚洲内射少妇av| 男女免费视频国产| 国产成人精品无人区| 日本欧美视频一区| 热re99久久精品国产66热6| 亚洲国产色片| 最后的刺客免费高清国语| 人人妻人人添人人爽欧美一区卜| 成人漫画全彩无遮挡| 九色成人免费人妻av| 女人精品久久久久毛片| 又大又黄又爽视频免费| 亚洲欧洲国产日韩| 如日韩欧美国产精品一区二区三区 | 男女无遮挡免费网站观看| 丝袜在线中文字幕| a级毛色黄片| 亚洲欧美中文字幕日韩二区| 国语对白做爰xxxⅹ性视频网站| 亚洲成人手机| 久久国产亚洲av麻豆专区| 亚洲国产精品一区三区| 一本—道久久a久久精品蜜桃钙片| 少妇丰满av| 国产一区二区在线观看日韩| 这个男人来自地球电影免费观看 | 国产69精品久久久久777片| 五月天丁香电影| 久久国产亚洲av麻豆专区| 卡戴珊不雅视频在线播放| 午夜免费鲁丝| 男女国产视频网站| 亚洲精品国产成人久久av| 日本猛色少妇xxxxx猛交久久| 亚州av有码| www.色视频.com| 色婷婷av一区二区三区视频| 久久久久精品久久久久真实原创| 又大又黄又爽视频免费| 日本色播在线视频| 精品人妻偷拍中文字幕| 亚洲电影在线观看av| 在线天堂最新版资源| 肉色欧美久久久久久久蜜桃| www.av在线官网国产| 亚洲人成网站在线播| 麻豆成人午夜福利视频| 亚洲四区av| 高清不卡的av网站| 91在线精品国自产拍蜜月| 久久午夜综合久久蜜桃| 久久ye,这里只有精品| 亚洲综合色惰| 欧美日韩在线观看h| 国产精品偷伦视频观看了| 国产午夜精品久久久久久一区二区三区| 免费高清在线观看视频在线观看| 久久久午夜欧美精品| 亚洲激情五月婷婷啪啪| 美女大奶头黄色视频| 日韩一本色道免费dvd| 欧美bdsm另类| a级毛片在线看网站| 天堂中文最新版在线下载| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩东京热| 亚洲av成人精品一区久久| 日韩电影二区| 91精品一卡2卡3卡4卡| 久久久久久久久久成人| 一个人看视频在线观看www免费| 国产精品不卡视频一区二区| 国产欧美日韩精品一区二区| 成人综合一区亚洲| 久久久欧美国产精品| 人妻 亚洲 视频| 久久久欧美国产精品| 成年人免费黄色播放视频 | 男女啪啪激烈高潮av片| 亚洲精品日韩在线中文字幕| 如日韩欧美国产精品一区二区三区 | 天天操日日干夜夜撸| av.在线天堂| 国产极品天堂在线| 亚洲精品视频女| 久久99精品国语久久久| 日本色播在线视频| 草草在线视频免费看| 亚洲欧美成人综合另类久久久| 少妇人妻 视频| 欧美日韩精品成人综合77777| 国产亚洲5aaaaa淫片| 国产午夜精品久久久久久一区二区三区| 免费人妻精品一区二区三区视频| 久久99蜜桃精品久久| 日本色播在线视频| 18禁裸乳无遮挡动漫免费视频| 精华霜和精华液先用哪个| 久久精品国产自在天天线| 久久精品国产自在天天线| 色吧在线观看| 啦啦啦在线观看免费高清www| 国产精品三级大全| 99九九线精品视频在线观看视频| 一级二级三级毛片免费看| 欧美国产精品一级二级三级 | 久久这里有精品视频免费| 亚洲人成网站在线观看播放| 午夜免费观看性视频| 日韩电影二区| 久久久欧美国产精品| 美女中出高潮动态图| 亚洲国产毛片av蜜桃av| 亚洲国产毛片av蜜桃av| 久久婷婷青草| 少妇的逼水好多| 国产探花极品一区二区| www.av在线官网国产| av福利片在线观看| 成人免费观看视频高清| av天堂中文字幕网| 肉色欧美久久久久久久蜜桃| 美女大奶头黄色视频| 一区在线观看完整版| av免费在线看不卡| 大码成人一级视频| 亚洲欧洲精品一区二区精品久久久 | 日本猛色少妇xxxxx猛交久久| 自拍欧美九色日韩亚洲蝌蚪91 | av在线app专区| av在线老鸭窝| 精品久久久噜噜| kizo精华| 欧美激情极品国产一区二区三区 | 国产精品一区二区在线不卡| 亚洲色图综合在线观看| videossex国产| 五月天丁香电影| 韩国高清视频一区二区三区| 国产爽快片一区二区三区| 午夜视频国产福利| 性高湖久久久久久久久免费观看| 久久久久久久久久久免费av| 国产av码专区亚洲av| 亚洲av国产av综合av卡| 最黄视频免费看| 少妇丰满av| 丰满迷人的少妇在线观看| 卡戴珊不雅视频在线播放| 亚洲国产精品成人久久小说| 国产精品久久久久久久久免| 免费少妇av软件| a级片在线免费高清观看视频| 欧美人与善性xxx| 久久国产精品大桥未久av | 99九九在线精品视频 | 国产免费一级a男人的天堂| 汤姆久久久久久久影院中文字幕| 欧美bdsm另类| 一本色道久久久久久精品综合| 九色成人免费人妻av| 偷拍熟女少妇极品色| 91精品国产国语对白视频| 一区二区三区免费毛片| 下体分泌物呈黄色| 国产精品久久久久久久久免| 久久久国产欧美日韩av| 亚洲精华国产精华液的使用体验| av福利片在线观看| 久久热精品热| av免费观看日本| 国产成人精品婷婷| av在线播放精品| a 毛片基地| 久久久久久久久久久丰满| 少妇被粗大的猛进出69影院 | 各种免费的搞黄视频| www.av在线官网国产| 99久久精品一区二区三区| 人人妻人人看人人澡| 女人久久www免费人成看片| 精品国产露脸久久av麻豆| 男人舔奶头视频| 亚洲av成人精品一二三区| 国产国拍精品亚洲av在线观看| 午夜免费鲁丝| 国产成人aa在线观看| 内地一区二区视频在线| 人妻 亚洲 视频| av福利片在线观看| 自拍偷自拍亚洲精品老妇| 日韩三级伦理在线观看| 天堂8中文在线网| 欧美日韩av久久| 久热久热在线精品观看| 一个人免费看片子| 国产成人午夜福利电影在线观看| 国产有黄有色有爽视频| 亚洲美女黄色视频免费看| 精品久久久久久电影网| 亚洲欧洲日产国产| 亚洲成人手机| 99热这里只有是精品在线观看| 一个人免费看片子| 五月开心婷婷网| 男的添女的下面高潮视频| 成人国产麻豆网| 日本猛色少妇xxxxx猛交久久| 国产免费视频播放在线视频| 99久国产av精品国产电影| 热re99久久精品国产66热6| 美女中出高潮动态图| 性高湖久久久久久久久免费观看| 欧美亚洲 丝袜 人妻 在线| 桃花免费在线播放| 久久国产亚洲av麻豆专区| a级毛色黄片| 夜夜爽夜夜爽视频| 黑丝袜美女国产一区| 亚洲熟女精品中文字幕| 多毛熟女@视频| 91精品国产国语对白视频| 麻豆成人午夜福利视频| 乱码一卡2卡4卡精品| 亚洲国产日韩一区二区| 我的老师免费观看完整版| 桃花免费在线播放| av国产精品久久久久影院| 国产日韩欧美视频二区| 日韩视频在线欧美| 波野结衣二区三区在线| 国产精品成人在线| 精品午夜福利在线看| av卡一久久| 极品教师在线视频| 国产精品人妻久久久久久| 日韩熟女老妇一区二区性免费视频| 一区在线观看完整版| 国产日韩一区二区三区精品不卡 | www.色视频.com| 久久韩国三级中文字幕| 性色av一级| kizo精华| 男女边摸边吃奶| 午夜激情久久久久久久| 嫩草影院入口| 欧美丝袜亚洲另类| 国产精品人妻久久久久久| 男的添女的下面高潮视频| xxx大片免费视频| 日韩免费高清中文字幕av| 亚洲av在线观看美女高潮| 99视频精品全部免费 在线| 亚洲怡红院男人天堂| 黄色毛片三级朝国网站 | 国产女主播在线喷水免费视频网站| 18禁裸乳无遮挡动漫免费视频| 国产成人免费无遮挡视频| 久久精品国产亚洲网站| 日韩精品免费视频一区二区三区 | 亚洲精品色激情综合| 一区二区三区免费毛片| 内地一区二区视频在线| 国产成人精品福利久久| 国产午夜精品久久久久久一区二区三区| 国产高清国产精品国产三级| 国产日韩欧美视频二区| 久久婷婷青草| 欧美区成人在线视频| av福利片在线| 啦啦啦视频在线资源免费观看| 特大巨黑吊av在线直播| 国产精品99久久久久久久久| 美女主播在线视频| 欧美日韩视频高清一区二区三区二| 热re99久久精品国产66热6| 午夜福利,免费看| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久| 日本欧美国产在线视频| 又爽又黄a免费视频| 亚洲成色77777| 久久久久久久久大av| 国产综合精华液| 国产精品国产三级国产专区5o| 成年女人在线观看亚洲视频| 91久久精品电影网| 精品午夜福利在线看| 亚洲精品日本国产第一区| 国产高清三级在线| h日本视频在线播放| 26uuu在线亚洲综合色| 观看美女的网站| a级毛色黄片| 免费看不卡的av| 欧美日韩一区二区视频在线观看视频在线| 国产有黄有色有爽视频| 亚洲一级一片aⅴ在线观看| 国语对白做爰xxxⅹ性视频网站| 久久精品国产鲁丝片午夜精品| 久久久久久久久久久免费av| 99热国产这里只有精品6| 亚洲色图综合在线观看| 女的被弄到高潮叫床怎么办| 人人妻人人爽人人添夜夜欢视频 | 精品久久久久久久久av| 曰老女人黄片| 99久久精品一区二区三区| 日韩av不卡免费在线播放| 五月玫瑰六月丁香| 美女内射精品一级片tv| av在线老鸭窝| av网站免费在线观看视频| 国产亚洲最大av| 成人二区视频| 在线精品无人区一区二区三| 只有这里有精品99| 亚洲国产精品一区三区| 香蕉精品网在线| 人妻 亚洲 视频| 99精国产麻豆久久婷婷| 99久久精品热视频| 亚洲精品国产成人久久av| 亚洲欧美日韩另类电影网站| 久久久精品免费免费高清| 2022亚洲国产成人精品| 91久久精品国产一区二区三区| 如何舔出高潮| 国产av码专区亚洲av| 在线观看三级黄色| 亚洲国产色片| 伦精品一区二区三区| 五月伊人婷婷丁香| 午夜老司机福利剧场| 亚洲国产毛片av蜜桃av| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 亚洲精品中文字幕在线视频 | av在线观看视频网站免费| 亚洲欧洲精品一区二区精品久久久 | 我要看黄色一级片免费的| 老女人水多毛片| 男女免费视频国产| 精品久久国产蜜桃| 少妇人妻精品综合一区二区| 国产伦在线观看视频一区| 日韩中字成人| 一区二区三区乱码不卡18| 一级爰片在线观看| 欧美精品一区二区大全| 精品一区在线观看国产| 国产无遮挡羞羞视频在线观看| 三级经典国产精品| 国产成人免费观看mmmm| 最后的刺客免费高清国语| 免费少妇av软件| 免费人妻精品一区二区三区视频| 免费观看性生交大片5| 丝袜在线中文字幕| 大香蕉久久网| 爱豆传媒免费全集在线观看| 国产淫片久久久久久久久| 亚洲久久久国产精品| 国产成人91sexporn| 插阴视频在线观看视频| 亚洲第一av免费看| 在线观看av片永久免费下载| 有码 亚洲区| 欧美丝袜亚洲另类| 久久女婷五月综合色啪小说| 美女国产视频在线观看| 国内精品宾馆在线| 99久久精品一区二区三区| 亚洲图色成人| 人体艺术视频欧美日本| 又粗又硬又长又爽又黄的视频| 少妇猛男粗大的猛烈进出视频| av专区在线播放| 亚洲第一区二区三区不卡| 午夜激情久久久久久久| 国产精品久久久久久精品电影小说| 在线观看一区二区三区激情| 大话2 男鬼变身卡| 一级毛片 在线播放| 91精品伊人久久大香线蕉| 午夜影院在线不卡| 欧美三级亚洲精品| 精品久久久久久久久av| 男男h啪啪无遮挡| 国产成人精品一,二区| 亚洲av二区三区四区| 一级a做视频免费观看| 黄色视频在线播放观看不卡| 日韩成人伦理影院| 国产极品天堂在线| 欧美少妇被猛烈插入视频| 中文字幕制服av| 久久久久久久久久久免费av| 亚洲国产欧美在线一区| 国产精品嫩草影院av在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产成人午夜福利电影在线观看| 国产精品福利在线免费观看| 校园人妻丝袜中文字幕| 久久久久久久久久久丰满| 国产精品久久久久久久电影| 亚洲熟女精品中文字幕| 高清欧美精品videossex| 精品久久久噜噜| 26uuu在线亚洲综合色| 精品国产国语对白av| 两个人免费观看高清视频 | 大片免费播放器 马上看| 久久精品国产a三级三级三级| 你懂的网址亚洲精品在线观看| 中国三级夫妇交换| 午夜老司机福利剧场| av在线老鸭窝| 2022亚洲国产成人精品| 美女xxoo啪啪120秒动态图| 亚洲精品第二区| 五月开心婷婷网| 日日撸夜夜添| 亚洲精品成人av观看孕妇| 欧美日韩精品成人综合77777| 久久综合国产亚洲精品| 天堂俺去俺来也www色官网| 午夜福利,免费看| 插逼视频在线观看| 视频区图区小说| 国产 精品1| 国产午夜精品一二区理论片| 日韩精品有码人妻一区| 丰满迷人的少妇在线观看| 91成人精品电影| 亚洲精品日韩在线中文字幕| 噜噜噜噜噜久久久久久91| 亚洲av成人精品一区久久| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频| 乱码一卡2卡4卡精品| 久久国产乱子免费精品| 一区二区三区四区激情视频| 99热这里只有是精品50| 我的老师免费观看完整版| 日日爽夜夜爽网站| 久久久a久久爽久久v久久| 精品人妻熟女毛片av久久网站| 日韩精品免费视频一区二区三区 | 只有这里有精品99| 王馨瑶露胸无遮挡在线观看| 精品亚洲成国产av| 九草在线视频观看| 成人国产麻豆网| 国内精品宾馆在线| 精品酒店卫生间| 各种免费的搞黄视频| 亚洲国产成人一精品久久久| 高清午夜精品一区二区三区| 欧美国产精品一级二级三级 | 色视频www国产| 免费人成在线观看视频色| 日韩,欧美,国产一区二区三区| 精品酒店卫生间| 亚洲性久久影院| 99精国产麻豆久久婷婷| 亚洲精品视频女| 亚洲欧洲精品一区二区精品久久久 | 有码 亚洲区| 亚洲人与动物交配视频| 成年美女黄网站色视频大全免费 | 色婷婷av一区二区三区视频| 最近2019中文字幕mv第一页| 人妻一区二区av| 精品国产一区二区久久| 性色avwww在线观看| 男女无遮挡免费网站观看| 永久免费av网站大全| 成人无遮挡网站| 一个人看视频在线观看www免费| 18禁裸乳无遮挡动漫免费视频| 成年av动漫网址| 国内少妇人妻偷人精品xxx网站| 亚洲欧美清纯卡通| 国产欧美日韩精品一区二区| 国产成人91sexporn| av女优亚洲男人天堂| 在线播放无遮挡| 国产精品秋霞免费鲁丝片| 久久午夜福利片| www.av在线官网国产| 啦啦啦在线观看免费高清www| 中文在线观看免费www的网站| 男人舔奶头视频| 久久久久精品性色| 国产在线视频一区二区| 亚洲av电影在线观看一区二区三区| 亚洲av在线观看美女高潮| 观看av在线不卡| 中文字幕精品免费在线观看视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av欧美aⅴ国产| 亚洲国产精品一区三区| 日韩一区二区视频免费看| 欧美日韩视频精品一区| 少妇熟女欧美另类| 伦理电影大哥的女人| 久久久久久久精品精品| 尾随美女入室| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 免费大片黄手机在线观看| 女人精品久久久久毛片| 黑人巨大精品欧美一区二区蜜桃 | 大片电影免费在线观看免费| 男人狂女人下面高潮的视频| 国产亚洲5aaaaa淫片| 国模一区二区三区四区视频| 久久久久国产网址| 国产亚洲午夜精品一区二区久久| 成人综合一区亚洲| 亚洲精品自拍成人| 国产精品99久久99久久久不卡 | 一本色道久久久久久精品综合| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 99九九在线精品视频 | 亚洲美女视频黄频| 久久人人爽人人片av| 国产有黄有色有爽视频| 亚洲国产最新在线播放| 99re6热这里在线精品视频| 在线观看三级黄色| 久久久久久久亚洲中文字幕| 国产精品国产三级专区第一集| 国产日韩欧美视频二区| 欧美xxⅹ黑人| 亚洲精品日本国产第一区| 成人亚洲精品一区在线观看| 男男h啪啪无遮挡| 高清视频免费观看一区二区| 人人澡人人妻人| 精品亚洲成国产av| 亚洲精品乱久久久久久| 免费看av在线观看网站| 永久免费av网站大全| 国产黄色视频一区二区在线观看| 免费av不卡在线播放| 国产成人a∨麻豆精品| 国产伦理片在线播放av一区| 久久久久人妻精品一区果冻| 在线观看一区二区三区激情| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产欧美日韩在线播放 | 亚洲精品国产av成人精品| 亚洲电影在线观看av| av在线播放精品| 免费av不卡在线播放| 综合色丁香网| 免费大片黄手机在线观看| 日本猛色少妇xxxxx猛交久久| 极品少妇高潮喷水抽搐| 美女主播在线视频| 亚洲欧美中文字幕日韩二区| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 午夜视频国产福利| 欧美日韩视频高清一区二区三区二| 日韩,欧美,国产一区二区三区| 99九九线精品视频在线观看视频| 内射极品少妇av片p| 国产亚洲欧美精品永久| 人人妻人人看人人澡| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 黄色怎么调成土黄色| 日韩成人伦理影院| 国产av国产精品国产| 久久亚洲国产成人精品v| 少妇裸体淫交视频免费看高清| 午夜视频国产福利| 午夜免费鲁丝| 亚洲av二区三区四区| 日日摸夜夜添夜夜爱| 九九在线视频观看精品| 成年人午夜在线观看视频| 欧美三级亚洲精品| 精品人妻熟女av久视频| 国产精品一区二区性色av| a 毛片基地| 女性生殖器流出的白浆| 能在线免费看毛片的网站| 在线看a的网站| 亚洲综合精品二区| 亚洲av欧美aⅴ国产| 99久国产av精品国产电影| 国产极品天堂在线| 18禁动态无遮挡网站| 少妇人妻一区二区三区视频|