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

    智能手機內(nèi)置加速度傳感器監(jiān)測走跑運動能量消耗的研究

    2018-12-04 09:21:44陳慶果袁川范江江陳恩格
    首都體育學(xué)院學(xué)報 2018年5期
    關(guān)鍵詞:能量消耗智能手機測量

    陳慶果 袁川 范江江 陳恩格

    摘 要:為探討走跑運動中手機內(nèi)置加速度傳感器原始信號的處理方法,在此基礎(chǔ)上建構(gòu)走跑運動的能耗模型并對其測量準(zhǔn)確性進(jìn)行評估。方法:45名受試者分別同時佩戴3部智能手機、1個能耗監(jiān)測儀和1套氣體代謝分析儀,在400 m田徑場上分別以4 、6 、8 和10 km/h的速度完成走跑運動各5 min。結(jié)果:在VM輸出上,手機(腰部)、手機(上臂)和手機(荷包)與GT3X的相關(guān)系數(shù)分別為0.963、0.952和0.941。在所建構(gòu)模型中,除10 km/h跑速下手機(荷包)和GT3X預(yù)測值90%CI落在校標(biāo)(K4b2)的等效區(qū)間外,其他均在等效區(qū)間內(nèi);在RMSE上,GT3X>手機(荷包)>手機(上臂)>手機(腰部)。3部手機在不同跑速下能耗預(yù)測值的ICC為0.872~0.975,均高度相關(guān)。各手機與GT3X在B-A圖上落在一致性限度區(qū)間以外的點的百分比均不高于5%。結(jié)論:設(shè)計的加速度信號預(yù)處理方法具有應(yīng)用價值,所使用的智能手機能夠準(zhǔn)確地預(yù)測走跑運動的能量消耗,并且腰部、上臂和荷包均是有效佩戴點。

    關(guān)鍵詞:智能手機;加速度傳感器;走跑運動;能量消耗;測量

    中圖分類號:G 804.2 學(xué)科代碼:040302 文獻(xiàn)標(biāo)識碼:A

    Abstract: Objective: The purpose of this work was to explore the processing approach on original signal from built-in accelerometer of smartphone, and develop EE predictive model of walking and running and evaluate on its accuracy of measurement. Methods: 45 subjects perform 20 minutes of filed walking at 4 km / h and 6 km / h and filed running at 8 km / h and 10km/h on the 400 meters track field. Three smart phones, an activity monitor and calorimetry analyzer were placed as recommended. Result: The correlation coefficient (VM) between phone (waist), phone (upper arm) and mobile phones (purse) and GT3X is respectively 0.963, 0.952 and 0.941. Except for EE of phone (purse) and GT3X at 10km/h,90%CI for a mean of the estimates is respectively fell within a proposed equivalence zone of the measured MET from OM. Order of RMSE size is GT3X> phone (purse)> Phone (upper arm)> arm (waist).ICC among each phone under every speed is highly relevant(0.872-0.975). The percentage of points out of limits of agreement from each device was under 5%. Conclusion: The methods of preprocessing acceleration signal designed by the study had applicable value. Smartphones could accurately assess EE of walking and running and the position of waist, upper arm and purse are effectively wearing points.

    Keywords: smartphone; accelerometer; walking and running; energy expenditure; assessment

    近年來,快走、慢跑的參與者數(shù)量呈井噴式增長,健身走、跑中的能量消耗問題也備受關(guān)注,以加速度傳感器為基礎(chǔ)的運動能耗監(jiān)測儀受到了越來越多健身者的青睞,相關(guān)產(chǎn)品如健身手環(huán)、健身手表等層出不窮,品牌型號更是紛繁復(fù)雜,其中美國Actigraph公司生產(chǎn)的系列產(chǎn)品在研究中被廣泛應(yīng)用 [1]。但這些運動能耗監(jiān)測設(shè)備在健身群體中很難大面積地普及,一是因為需要專門佩戴,耗費使用者額外的精力,二是因為價格不菲,普通健身者難以承擔(dān);因此,測量的效益問題有待解決[2]。而智能手機快速的普及和性能的不斷完善為解決這一問題提供了新思路。

    目前,加速度傳感器被普遍內(nèi)置于智能手機中,手機已具備運動能耗監(jiān)測產(chǎn)品的硬件條件,且不會給佩戴者造成任何的使用負(fù)擔(dān)和經(jīng)濟壓力,是一種便于推廣普及的運動能耗監(jiān)測設(shè)備。目前已有咕咚運動、樂動力等多款手機軟件可測量運動能耗,但這些軟件都是基于GPS(全球衛(wèi)星定位系統(tǒng))來獲取空間數(shù)據(jù),進(jìn)而獲取速度數(shù)據(jù)來推算運動能耗。受信號限制,該方法只能在室外使用,不便監(jiān)控家務(wù)和工作體力活動,并且在走跑活動中也存在較大的測量誤差。明鑫等的研究表明國內(nèi)的4款運動軟件預(yù)測值遠(yuǎn)低于校標(biāo)值[3],國外也有類似的研究結(jié)果[4],影響其測量準(zhǔn)確性的因素主要是:1)不能反映對運動能耗較敏感的垂直軸情況;2)采樣頻率低。

    現(xiàn)階段,基于智能手機內(nèi)置加速度傳感器的體力活動研究還主要集中在人機智能化領(lǐng)域,以探究動作識別為主[5-6],在健身走、跑的研究中也僅是集中在計步上,還缺乏對能耗測量的系統(tǒng)研究。雖然目前對運動能耗監(jiān)測儀測量準(zhǔn)確性的評估眾多,各種能耗模型也常見報告,但是原始加速度信號的處理方法是各品牌的商業(yè)機密;因此,有必要以智能手機內(nèi)置加速度傳感器為基礎(chǔ),探討健身走、跑的原始加速度信號處理過程,在此基礎(chǔ)上建構(gòu)能耗模型,并對其進(jìn)行測量準(zhǔn)確性評估,為智能手機在健身實踐和科學(xué)研究中的應(yīng)用奠定基礎(chǔ)。

    1 研究對象與方法

    1.1 研究對象

    受試者為45名在校大學(xué)生,其中男性23名,女性22名,平均年齡21.7歲,見表1。所有受試者均肢體健康,無心血管疾病史,測試前24 h內(nèi)無大強度體力活動。所有受試者均簽署知情同意書。

    1.2 研究方法

    1.2.1 測量方法

    1)測量儀器。

    受試者身高和體重采用恒康佳業(yè)HK-600身高體重儀測量,韓國VIVENTE-GOLD體成分儀測量其體脂率。

    使用3部紅米note2手機作為測試手機,其內(nèi)置的三軸加速度傳感型號為ACCELEROMETER,量程±2 g。選擇該型號手機是鑒于其正向研究結(jié)果具有良好的普適性。本研究編寫專門軟件收集手機加速度傳感器3個軸的原始數(shù)據(jù),采樣頻率為75 Hz,所獲數(shù)據(jù)以CSV格式保存在手機中。測試后,將原始數(shù)據(jù)下載后歸類整理,利用SPSS軟件、minitab軟件和MATLAB軟件進(jìn)行后期處理。

    為便于測量準(zhǔn)確性的橫向比較,使用美國產(chǎn)運動能耗監(jiān)測儀Actigraph-GT3X(以下簡稱GT3X)同步監(jiān)測健身走、跑的能量消耗,采樣頻率設(shè)為30 Hz,。使用配套的Actilife6.0軟件對GT3X進(jìn)行重置及后期數(shù)據(jù)處理。

    使用便攜式氣體代謝分析儀CosmedK4b2(以下簡稱K4b2)同步監(jiān)測活動時的能量消耗,其測量結(jié)果在同類研究中被廣泛視為金標(biāo)準(zhǔn)[7]。為保障測量的精度,每日測試前均進(jìn)行定標(biāo)。測試時以某一分鐘的第1秒為起始時間按壓K4b2開始鍵以記錄數(shù)據(jù),保證其與GT3X和智能手機數(shù)據(jù)的一致性。

    2) 測試方案。

    測試時每個受試者同時佩戴K4b2和3部智能手機以及1個GT3X。智能手機分別置于受試者的左側(cè)髖部髂棘處(運動腰帶固定)、左側(cè)上臂外上1/3處(Armband固定)和左側(cè)腹部荷包處。GT3X使用彈性腰帶置于受試者右側(cè)髖部、肚臍水平高度。

    走跑測試分為4個階段,每個階段分別為5 min,速度依次為4 (慢走)、6 (快走)、8 (慢跑)和10 km/h(跑),間隔時間依據(jù)心率的恢復(fù)情況而定,1~5 min不等。

    測試在標(biāo)準(zhǔn)400 m田徑場上進(jìn)行,跑道上每隔5 m放置一個標(biāo)志桶。測試時播放相應(yīng)的音頻,4個階段的提示音時間間隔分別為4.50、3.00、2.25和1.80 s,測試音每響一次受試者走(跑)完一個標(biāo)志桶(5 m),受試者在正式測試前會在工作人員的示范下熟悉節(jié)拍和速度。正式測試時根據(jù)標(biāo)志物和節(jié)拍調(diào)整速度,有工作人員陪同騎行,對受試者進(jìn)行實時提醒。

    1.2.2 數(shù)據(jù)處理

    自編程序在MATLAB7.0上運行以對手機加速度原始信號進(jìn)行處理,過程如下:1)進(jìn)行快速傅里葉變化觀察其主頻,根據(jù)主頻結(jié)果設(shè)置二階巴特沃茲帶通濾波器對信號進(jìn)行濾波;2)進(jìn)行修正處理以去除趨勢;3)將每軸的AC合成3軸的VM(vector magnitude, 矢量計數(shù))。

    本研究將45名受試者隨機分到建模組(30人)和驗證組(15人),并單獨建立相應(yīng)的統(tǒng)計分析數(shù)據(jù)庫。在模型構(gòu)建中,使用SPSS軟件對各跑速下手機、GT3X和K4b2的數(shù)據(jù)進(jìn)行描述性統(tǒng)計,采用Pearson相關(guān)考察各加速度設(shè)備之間,以及與K4b2之間輸出結(jié)果的一致性。智能手機預(yù)測模型建立的程序和方法參考FREEDSON VM3 Combination(2011)聯(lián)合方程(以下簡稱VM3C)的研制。

    在模型驗證中,使用MiniTAP軟件進(jìn)行等效性檢驗,以判斷預(yù)測METs和實測METs測量結(jié)果的一致性,標(biāo)準(zhǔn)為:預(yù)測均值90%置信區(qū)間是否落入實測均值的等效區(qū)間(?滋±10%?滋)[8]。如果落入則接受備擇假設(shè):上限<檢驗值/校標(biāo)均值<下限,可認(rèn)定2種測量方法等效;利用RMSE指標(biāo)考察各加速度設(shè)備預(yù)準(zhǔn)確性;利用ICC判斷智能手機能耗預(yù)測值之間的一致性;用Bland-Altman圖法進(jìn)一步分析各加速度設(shè)備與效標(biāo)的一致性。

    2 研究結(jié)果

    2.1 智能手機加速度傳感器信號的預(yù)處理

    常見的人體能耗監(jiān)測儀的輸出數(shù)據(jù)為count值,多數(shù)能耗預(yù)測模型也是以該數(shù)值為基礎(chǔ)建構(gòu)而來;但該數(shù)值并不是加速度傳感器產(chǎn)生的原始數(shù)據(jù),而是經(jīng)過預(yù)處理后的數(shù)據(jù),而處理過程是各廠商的商業(yè)機密,且處理的方法和流程均有所區(qū)別,震動盤上的實驗表明相同的活動各品牌的count值差異在300~3 000[9]。本研究借鑒相關(guān)研究成果,結(jié)合具體的使用條件對信號進(jìn)行處理,將其分為濾波、修正和數(shù)據(jù)合成3個環(huán)節(jié)。

    2.1.1 濾波

    在體力活動監(jiān)測中,濾波(wave filtering)的主要作用是將運動信號保留,將其他頻段噪聲和干擾排除,是提高測量精度的一項重要措施。本研究借鑒信號處理的常用方法,采用傅里葉轉(zhuǎn)換(FFT)技術(shù)從原始數(shù)據(jù)中得出頻譜信號,根據(jù)波峰觀察和先驗經(jīng)驗確定主頻,從而確定濾波器帶通頻率。圖1顯示3部手機在不同速度走跑下,三軸上的主頻率基本都在5 Hz以內(nèi),同時為排除低頻信號(直流漂移、電傳感器老化和溫度變化等)的干擾,將帶通頻率設(shè)計為0.2~5.0 Hz,參考孫泊等[3]的濾波實踐,采用二階巴特沃茲帶通濾波器對三軸的原始加速度數(shù)據(jù)分別進(jìn)行濾波。

    2.1.2 數(shù)據(jù)的修正

    鑒于智能手機在實際使用中的佩戴條件,運動軸和運動方向不能完全匹配,重力會產(chǎn)生分量影響各軸,參考SEKINE等[9]的建議,同時對3個軸過濾后的信號進(jìn)行去趨勢處理以消除各軸上重力的影響,公式為:?字t'=?字t-■t-5~t ,處理前后對比如圖2所示。

    2.1.3 數(shù)據(jù)的合成

    加速度傳感器用于運動能耗監(jiān)測時,其采樣頻率往往都不低于30 Hz,也就是說一個三軸加速度傳感器每秒就要生成不少于90個原始數(shù)據(jù),長時間的監(jiān)測會產(chǎn)生海量的數(shù)據(jù),這就需要對數(shù)據(jù)預(yù)先進(jìn)行合成。目前使用的計算指標(biāo)有單位時間內(nèi)的信號峰值、閾值計數(shù)、積分值、四分位差、峰度系數(shù)和偏度系數(shù)等,其中積分值最常使用。鑒于此,同時為克服信號雙向性帶來的積分消抵,本研究取絕對值后再計算積分。

    在實際應(yīng)用中,能耗儀能較為穩(wěn)定地固定在身體的某一部分,避免產(chǎn)生相對運動,而手機的固定條件不如能耗儀,傳感器3個運動軸與運動方向不能完全匹配,會受到重力分量的影響,單一軸的AC數(shù)值表征活動會產(chǎn)生較大的誤差。鑒于此,本研究采用三軸的合成矢量計數(shù)(VM):

    VM=(AC2x+AC2y+AC2z)1/2

    2.2 建模組走跑運動能耗和VM特征分析

    從表2可知,不同速度走跑時GT3X和各手機的VM在數(shù)值上有較大差異,但變化趨勢基本相同,均與校標(biāo)保持一致。進(jìn)一步使用Pearson相關(guān)分析后發(fā)現(xiàn)GT3X與各手機之間相關(guān)系數(shù)均在0.9以上(見表3),這不但表明GT3X和手機原始信號提供高度相似的信息,GT3X中能耗方程建立的方法和條件適用于本研究所使用的手機,還說明了本研究加速度信號處理方法是合理的。

    2.3 智能手機能耗方程的建構(gòu)

    最新版的Actilife6.0軟件中,共提供4個kcal能耗預(yù)測模型,王軍利的研究表明VM3C方程系統(tǒng)誤差較小[10],最為適合我國年輕人的能耗算法方程。該方程是Sasaki方程[11]和Williams方程[12]的組合,分為2部分:

    E/(kcal/min)=0.001 064×VM+0.087 512×BM-5.500 229 (VM≥2453);

    kcal/min=VM×0.000 019 1×BM(VM<2 453)。

    圖3~5為各智能手機的VM、體質(zhì)量和能耗的三維空間散點圖,各圖的散點均一定程度上存在著拐點現(xiàn)象,這也表明智能手機的能耗預(yù)測模型適宜于使用這種分段式方程。該方程建立思路為:1)建立METs的VM預(yù)測方程;2)以中等強度閾值3MET回代入方程求得相應(yīng)的VM閾值;3)然后以該點分別構(gòu)建低強度活動和中高強度活動預(yù)測方程。

    2.3.1 METs預(yù)測方程與3MET-VM臨界點

    由表4可知,各部智能手機的能耗預(yù)測模型中R2均高于0.9,且SEE較低,提示3個方程的擬合程度都較高。

    根據(jù)3個方程可以推算出3METs對應(yīng)的VM臨界點:手機(腰部)為86 counts/min、手機(上臂)為26 counts/min、手機(荷包)為120 counts/min。以上述切點為依據(jù)劃分活動為低強度活動和中高強度活動,并分別構(gòu)建相應(yīng)方程。

    2.3.2 低強度體力活動能耗方程

    VM3C聯(lián)合方程低強度體力活動的能耗預(yù)測公式為E/(kcal/min)=VM×0.0000191×BM。該公式是基于能量守恒的原理演算而來,并不是通過線性回歸建構(gòu)的方程;因此,本研究依據(jù)GT3X與各智能手機VM的比值,分別建立各手機低強度區(qū)間的能耗模型,見表5。

    2.3.3 中高強度體力活動能耗方程

    參考VM3C聯(lián)合方程中高強度體力活動方程建立的流程,將BM(體重)變量引入到模型中,得到各智能手機的預(yù)測模型見表6,各方程決定系數(shù)R2=0.908~0.919,均高于0.9,且SEE較低,3個模型均具有較高的擬合度。

    2.4 智能手機測量走跑運動能耗準(zhǔn)確性的檢驗

    判斷預(yù)測值和校標(biāo)值的一致性,傳統(tǒng)的方法為配對樣本t檢驗,但差異無統(tǒng)計學(xué)意義的結(jié)果只能得到一個“沒有結(jié)論”的結(jié)論,不能冒然得出2組數(shù)據(jù)等效或無差異[13];因此,本研究采用等效性檢驗中的可信區(qū)間法確定智能手機和GT3X的預(yù)測值與校標(biāo)值(K4b2測量結(jié)果)之間的等效性,依據(jù)為:預(yù)測均值90%置信區(qū)間是否落入實測均值的等效區(qū)間(?滋±10%?滋)[8]。從表7可知,4、6 、8 km/h,各手機以及GT3X預(yù)測數(shù)值90%的置信區(qū)間均落入校標(biāo)實測值的等效區(qū)間里,僅在10 km/h下,手機(荷包)和GT3X的測量與K4b2的測量不等效,90%置信區(qū)間分別為10.16~10.74、10.27~10.78落在了10.40~12.71的等效區(qū)間以外。從總的能耗上看,各加速度設(shè)備的置信區(qū)間也均落入等效區(qū)間里面,因此,基本可以認(rèn)定各加速度設(shè)備與K4b2之間能耗測量具有等效性。

    圖6顯示各跑速下手機與GT3X的均方根誤差(RMSE),整體上看,場地走的RMSE低于場地跑的RMSE。各加速度設(shè)備的RMSE均在10 km/h跑步中出現(xiàn)峰值。在各個速度上,手機(腰部)測量的RMSE均最小,分別為0.26、0.24、1.15和1.61,整體為0.81。而手機(荷包)隨著跑速的增加,其RMSE迅速增加,相比上一級速度增加的幅度分別為22.2%、149.1%和116.06%。從總能耗上看,各設(shè)備RMSE的排序為: GT3X>手機(荷包)>手機(上臂)>手機(腰部)。

    組內(nèi)相關(guān)系數(shù)(ICC)分析的結(jié)果顯示,手機(腰部)、手機(上臂)和手機(荷包)在不同速度下的能耗預(yù)測數(shù)值均高度相關(guān)。在以4、6 、8 和10 km/h的速度運動時,各手機能耗預(yù)測值之間的科隆巴赫系數(shù)分別為0.975、0.946、0.894和0.872,說明各智能手機能耗結(jié)果具有很好的一致性。

    通過Bland-Altman法,本研究分別計算了3部手機及GT3X與氣體代謝分析儀測量值的一致性限度區(qū)間(limits of agreement, mean ±1.96 s),手機(腰部)為(-1.7,2.1),手機(上臂)為(-2.0,2.0),手機(荷包)為(-2.4,2.1)和GT3X為(-2.8,1.7),各設(shè)備分別有7個、8個、9個和9個點在一致性限度區(qū)間以外,占總點數(shù)的比例分別為3.89%、4.44%、5%和5%,均不高于5%的標(biāo)準(zhǔn),表明其與氣體代謝分析儀的測量結(jié)果一致性好,與等效檢驗的結(jié)果相匹配。手機(荷包)和GT3X的測量值偏倚小于0,顯示整體上兩者略低估計能耗,而手機(腰部)的偏倚為0.2,又略微高估能耗,見表8。

    3 分析與討論

    3.1 智能手機與GT3X預(yù)測場地走跑能耗準(zhǔn)確性的對比分析

    按照等效性判斷標(biāo)準(zhǔn),除了10 km/h跑速下手機(荷包)和GT3X的90%置信區(qū)間落在校標(biāo)等效區(qū)間外,其他各跑速下,各手機及GT3X與K4b2的測量均等效??偟哪芰肯纳?,各智能手機、GT3X與K4b2的測量均具有等效性,通過進(jìn)一步分析測量誤差后發(fā)現(xiàn),GT3X在各個階段均低估能量消耗,這與王歡等[14]和王軍利等[10]的研究結(jié)果是一致的,GT3X的濾波頻率為0.25~2.5 Hz。Uemura等[15]認(rèn)為過窄的濾波頻率,影響其對于中高強度活動的測量;但不能認(rèn)為頻率設(shè)置不合理,因為GT3X運動能耗儀用來監(jiān)測日常的體力活動,走/跑活動只是一部分,濾波頻率的設(shè)置應(yīng)兼顧各項活動,不能因某類活動能耗預(yù)測誤差大就調(diào)整濾波的頻率,應(yīng)該從活動識別和能耗模型上入手,針對不同活動類型建立針對性的能耗方程[16]。

    在10 km/h下,手機測量的RMSE出現(xiàn)峰值,低于實際能耗,也同GT3X一樣,存在著高強度活動測量時的“平臺”現(xiàn)象[17],在走跑能耗監(jiān)測中,未來可以借用Crouter混合模型的思路,建立曲線回歸模型,解決走跑能耗預(yù)測的“平臺”問題[18]。

    但也應(yīng)該認(rèn)識到,GT3X使用的VM3C雖然說也是以不同速度走跑為運動方案建立的能耗模型,但其實驗對象、實驗方案有所不同;因此,并不能簡單地就通過等效性檢驗和RMSE分析就得出智能手機能耗監(jiān)測準(zhǔn)確性高于GT3X的結(jié)論, 還需要更多的研究去充實以得出概化性結(jié)論。本研究的結(jié)果只表明智能手機能夠準(zhǔn)確測量走跑時的能量消耗,可以作為GT3X的補充應(yīng)用到走跑類健身活動的指導(dǎo)和干預(yù)中。

    3.2 不同部位智能手機測量結(jié)果的一致性分析

    本研究中3部佩戴在不同部位的智能手機雖然說在輸出的VM上存在著顯著性的差異,但高度相關(guān)(P=0.962~0.984 ),表明不同部位的紅米note2手機提供高度相似的信息,在能耗監(jiān)測上,3部手機也均可以準(zhǔn)確地測量。這表明3個部位均可以用來佩戴手機監(jiān)測能量消耗,這主要是因為:1)場地走跑運動中上下肢的擺動幅度與人體行進(jìn)速度密切相關(guān),有研究就對佩戴在大腿、腰部、上臂和踝關(guān)節(jié)的三軸加速度計進(jìn)行比較后發(fā)現(xiàn)很難認(rèn)為哪個部位佩戴更有優(yōu)勢[19];2)智能手機內(nèi)置的傳感器為三軸傳感器,Jones等[20]的研究就表明佩戴位置對于三軸加速度計的影響要遠(yuǎn)遠(yuǎn)小于對于單軸加速度計的影響;3)本研究使用三軸的矢量合成計數(shù),而不是其中某一軸的積分值,能夠部分解決因放置方式不同導(dǎo)致的敏感軸與運動方向不一致的問題。Chen等[21]就認(rèn)為放置位置不同會導(dǎo)致三軸加速度計X、Y、Z 3個軸的偏轉(zhuǎn),但并不會對合成的VM造成影響。本研究的結(jié)果表明智能手機可以根據(jù)使用者的習(xí)慣放置在上衣荷包里,而這更符合普通人的佩戴習(xí)慣,有助于長時間進(jìn)行活動數(shù)據(jù)采集。此外,研究結(jié)果也支持跑步愛好者將手機使用Armband置于上臂處。

    本研究的結(jié)果雖然顯示部位對測量準(zhǔn)確性的影響小,但僅局限在監(jiān)測走跑這一類活動,在監(jiān)測其他非周期性運動或者靜力性活動時,佩戴部位卻很關(guān)鍵[22],加速度計均推薦戴在離身體質(zhì)心較近的部位[23-24]。Rosenberger等的研究表明count為積分值的回歸模型(一元回歸模型或者是基于動作識別的混合模型)預(yù)測活動能耗時,手腕部分的預(yù)測誤差顯著性高于髖部[25],而Staudermayer等采用機器學(xué)習(xí)算法集合多個特征指標(biāo)預(yù)測能耗表明腕部加速度計也能準(zhǔn)確預(yù)測能耗[26],另外一項研究[27]的結(jié)果也類似,這提示不同部位的加速度計需要有不同的特征指標(biāo),并組合不同的算法模型,這將是本領(lǐng)域研究的重點問題。

    3.3 手機加速度傳感器原始信號的處理

    研究者耳熟能詳?shù)母鱾€運動能耗監(jiān)測儀輸出的count值,實際并不是每個軸的加速度具體數(shù)值,而是對原始信號過濾、修正和合成而來的,是原始信號時域特征的提取,并沒有實際的生物學(xué)意義。采樣頻率的設(shè)置是信號采集初始步驟,遵循Nyquist標(biāo)準(zhǔn)[28]:所設(shè)頻率至少為運動中人體各環(huán)節(jié)最大頻率的2倍。在運動中身體的加速度信號一般低于10 Hz,日常體力活動中身體質(zhì)心處的運動頻率在0.3~3.5 Hz [29],因此人體軀干的運動信息采集不應(yīng)低于30 Hz [30],GT3X的最低采樣頻率為30 Hz。目前各品牌智能手機內(nèi)置的加速度傳感器的采樣頻率基本都可達(dá)到50 Hz,可以滿足體力活動能耗監(jiān)測的需求。

    濾波過程中設(shè)置帶通頻率是一個關(guān)鍵環(huán)節(jié),有效設(shè)置可以減少高頻噪聲干擾的同時降低低頻噪音的影響,避免所獲信號失真,不能反映運動情況。在實踐中,帶通頻率設(shè)置太大時各種低頻或高頻噪音難以消除,太小時則會造成運動信號遺漏。各品牌能耗儀的帶通頻率的設(shè)置主要是以獲取的活動數(shù)據(jù)為基礎(chǔ)進(jìn)行FFT頻譜分析,因此,獲取數(shù)據(jù)的活動項目不同,帶通區(qū)別有所差別,普遍在0.25~7 Hz [31]。本研究根據(jù)FFT頻譜分析的結(jié)果,結(jié)合具體的走跑特征確定手機的帶通頻率為0.2~5 Hz。而本研究后期的數(shù)據(jù)表明,智能手機VM的輸出也出現(xiàn)“平臺”現(xiàn)象[32-33],提示帶通頻率可進(jìn)行適當(dāng)調(diào)整,但調(diào)整的幅度需要深入的探究,避免非此即彼納入更多的無關(guān)信號。

    原始信號的合成是信號處理的最后環(huán)節(jié),而不同品牌能耗監(jiān)測儀提取的特征也不同,主要分為3類:1)閾值計數(shù)[34];2)時閾峰值[35];3)加速度積分值[36]。目前,加速度積分值的應(yīng)用最普遍,但不同品牌對加速度積分的計算也有所區(qū)別,有的只利用單軸(垂直軸)數(shù)據(jù),有的利用雙軸(垂直軸和矢狀軸)數(shù)據(jù)。Midorikawa等[33]認(rèn)為三軸合成的VM可以全面地反映體力活動中身體的活動情況,此外,還可以克服因為設(shè)備放置方式不同而造成的信號差異問題[37]。在實際使用中,智能手機往往是根據(jù)使用者的習(xí)慣放置,放置條件具有較大的不確定性,因此,在信號合成的時候應(yīng)使用三軸的數(shù)據(jù)合成VM來表征運動狀況。

    3.4 智能手機和GT3X內(nèi)置加速度傳感器的硬件基礎(chǔ)

    本研究采用紅米note2智能手機,市面價格一般為800元,內(nèi)置的傳感器型號是ACCELEROMETER,量程±2g、精度3.90E-3、功率為0.25 W,近似參數(shù)的傳感器市面價格10~20元,GT3X內(nèi)置的傳感器量程±6g,但是從測量的準(zhǔn)確性來說,更大的測量量程并沒有提高其測量的精度,這與人體走跑時候的加速度區(qū)間相關(guān)。Matthew等[38]的研究顯示走跑速度在1~5 m/s時人體骨盆處的加速度數(shù)值變化范圍為0.05~0.5g。本研究中即使跑速達(dá)到10 m/s,敏感軸的加速度數(shù)值也不超過1.5g,±2g量程能夠滿足測量場地走/跑的需要。GT3X內(nèi)置傳感器的其他參數(shù)未見報告,但從總體上看,GT3X輸出的VM數(shù)值與各智能手機的高度相似(0.941~0.984),兩者原始信號提供高度一致的信息,并且也與K4b2的測量結(jié)果的一致程度相當(dāng),說明了智能手機內(nèi)置加速度傳感器的配置能夠滿足體力活動測量的需要,此外,它不需要使用者單獨購買,不增加使用者額外的負(fù)擔(dān),是一種低成本高準(zhǔn)確性的運動能耗監(jiān)測工具。

    3.5 基于加速度傳感器的能耗方程分析

    無論是智能手機還是運動能耗監(jiān)測儀對運動能耗估算的準(zhǔn)確性除了與內(nèi)置的加速度傳感器的性能和原始數(shù)據(jù)的預(yù)處理的方法有關(guān)外,還與能量預(yù)測的方程密不可分,目前最為常見的模型建構(gòu)方法是以校標(biāo)為因變量,以count數(shù)值為自變量建立回歸方程。不僅不同能耗儀使用的方程不同,相同能耗儀也有多個方程可供選擇。例如GT3X軟件中就有8個METs和4個kcal的預(yù)測公式。不同的方程建立實驗條件也不同,有的建立在現(xiàn)實生活情景中[39],有的以跑臺或場地走跑為實驗方案[40];并且在實驗對象上,有的為兒童[41],有的為成年人[42]。本研究采用VM3C方程的思路建構(gòu)模型,該模型為分段模型,不同活動強度采用不同方程。該方程在國內(nèi)外廣泛使用,具有相對較好的信效度,但該方程測量不同活動類型的準(zhǔn)確性存在較大的差異。相關(guān)學(xué)者一致認(rèn)為:沒有一個方程能夠有效地預(yù)測各種類型的體力活動,只有建立在特定活動基礎(chǔ)上的模型才可能有較高的效度[43];因此,基于活動識別的能耗預(yù)測成為研究的熱點,該方法是對數(shù)據(jù)進(jìn)行預(yù)處理提取出相應(yīng)特征,使用分類器確定活動類型,從而調(diào)用相應(yīng)的能耗預(yù)測模型來預(yù)測能耗。Mikkel等[44]和Marco等[45]認(rèn)為在多類型活動的測試中,活動識別方法能夠顯著提高運動能耗的預(yù)測效度。

    4 結(jié)論與建議

    4.1 結(jié)論

    本研究選取的智能手機測量場地走跑能耗采用濾波(0.5~5 Hz二階巴特沃茲帶通濾波)、數(shù)據(jù)修正(?字t'=?字t-t-5~t)、數(shù)據(jù)合成[VM=(AC2x+AC2y+AC2z)1/2 ]的加速度信號處理方法具有應(yīng)用價值。以VM3C方程為模板建構(gòu)的智能手機能耗模型能夠準(zhǔn)確地測量本次研究驗證樣本的走跑活動能耗。腰部、上臂和荷包均是智能手機監(jiān)測走跑能耗時的有效佩戴點。由于受測試樣本數(shù)量的限制,本研究成果的概化性還需后續(xù)研究進(jìn)一步明確。

    4.2 建議

    本研究僅用一款手機作為測試設(shè)備,沒有對手機型號和品牌對測量結(jié)果的影響進(jìn)行評估;其次,對手機在荷包中的固定條件進(jìn)行了適當(dāng)限定,要求手機放置在左側(cè)腹部的荷包中,并且能夠相對固定避免產(chǎn)生較大的相對運動,這與使用者的習(xí)慣有出入;最后,本研究的方程專門預(yù)測走跑類活動,預(yù)測其他類型體力活動能耗時效度勢必降低。

    在以后的研究中,應(yīng)該進(jìn)一步豐富手機品牌的評測以得到概化性的結(jié)論,手機的選擇應(yīng)該更多地集中在中低端配置的手機,應(yīng)重視記錄原始數(shù)據(jù)的手機App,最好是專門設(shè)計可以自由修改關(guān)鍵參數(shù)的數(shù)據(jù)采集App。在能耗模型的輸入指標(biāo)的選擇上,應(yīng)該拓寬思路,輸入模型的count值除了用積分值代表外,還應(yīng)探討其他時域或頻域指標(biāo)的應(yīng)用價值,特別是某些具有活動類型識別價值的指標(biāo),并且選擇的指標(biāo)應(yīng)能夠匹配所使用的模型;而在建立能耗評估模型時,應(yīng)以多種類型的體力活動為方案,建構(gòu)基于活動識別技術(shù)的能耗評估方程,特別是隨著機器算法的興起,應(yīng)將諸如神經(jīng)網(wǎng)絡(luò)模型、支持向量機模型等引入到該領(lǐng)域,進(jìn)一步豐富該領(lǐng)域研究的同時切實提高能耗預(yù)測的準(zhǔn)確性。

    參考文獻(xiàn):

    [1] FREEDSON P,BOWLES H R,TROIANO R, et al. Assessment of physical activity using wearable monitors: recommendations for monitor calibration and use in the field[J]. Medicine and science in sports and exercise, 2012, 44(S1):1.

    [2] HEKLER E B,KLASNIJA P,TRAVER V, et al. Realizing effective behavioral management of health: The metamorphosis of behavioral science methods[J]. Pulse, IEEE, 2013, 4(5):29.

    [3]明鑫, 王斌. 手機運動軟件與 K4b2 測量步行運動能耗的比較研究[J]. 安徽體育科技, 2015, 36(3):55.

    [4] ZANDBERGEN P A,BARBEAU S J. Positional accuracy of assisted gps data from high-sensitivity gps-enabled mobile phones[J]. Journal of Navigation, 2011, 64(3):381.

    [5] MORILLO L M S,GONZALEZ-ABRIL L,RAMIREZ J A O, et al. Low energy physical activity recognition system on smartphones[J].Sensors, 2015, 15(3):5163.

    [6] SHOAIB M,BOSCH S,LNCEL O D,et al. Fusion of smartphone motion sensors for physical activity recognition[J]. Sensors, 2014, 14(6):10146.

    [7]湯強, 王香生, 盛蕾. 體力活動測量方法研究進(jìn)展[J]. 體育與科學(xué), 2008, 29(6):79.

    [8] NOLAN M, MITCHELL J R, DOYLE-BAKER P K. Validity of the Apple iPhoneR/iPod TouchR as an accelerometer-based physical activity monitor: a proof-of-concept study[J]. Journal of physical activity & health, 2014, 11(4):759.

    [9] SEKINE M, TAMURA T, OGAWA M, et al. Classification of acceleration waveform in a continuous walking record[C]//Austin: International Conference of the IEEE,1998:1523-1526.

    [10]王軍利, 張冰, 賈麗雅, 等. Actigraph (GT3X)加速度計測量我國19-29歲人群身體活動能耗的效度研究[J]. 體育科學(xué), 2012, 32(12):71.

    [11] SASAKI J E, JOHN D, FREEDSON P S. Validation and comparison of ActiGraph activity monitors[J].Journal of Science and Medicine in Sport, 2011, 14(5):411.

    [12] WILLIAMS R. Kcal estimates from activity counts using the Potential Energy Method[EB/OL].[2017-04-28].http://actigraphcorp.com/research-database.

    [13] 袁秉祥, 朱芳. 生物等效性檢驗的計算及計算程序[J]. 西安醫(yī)科大學(xué)學(xué)報, 2001, 22(2):164.

    [14]王歡, 王馨塘, 佟海青, 等. 三種加速度計測量多種身體活動的效度比較[J].體育科學(xué), 2014, 34(5):45.

    [15] UEMURA H, KATSUURA-KAMANO S, Yamaguchi M, et al. Abundant daily non-sedentary activity is associated with reduced prevalence of metabolic syndrome and insulin resistance[J]. Journal of endocrinological investigation, 2013, 36(11): 1069.

    [16] ELLIS K, KERR J, GODBOLE S, et al. Hip and Wrist Accelerometer Algorithms for Free-Living Behavior Classification[J]. Medicine and science in sports and exercise, 2016, 48(5):933.

    [17] SWARTZ A M, STRATH S J, BASSETT D R, et al. Estimation of energy expenditure using CSA accelerometers at hip and wrist sites[J]. Medicine and Science in Sports and Exercise, 2000, 32(9):450.

    [18] CROUTER S E, KUFFEL E, HAAS J D, et al. A refined 2-regression model for the actigraph accelerometer[J]. Medicine and science in sports and exercise, 2010, 42(5): 1029.

    [19] BOUTEN C V, SAUREN A A, VERDUIN M, et al. Effects of placement and orientation of body-fixed accelerometers on the assessment of energy expenditure during walking[J]. Medical & Biological Engineering & Computing, 1997, 35(1):50.

    [20] JONES S L, WOOD K ,Thompson R,et al. Effect ofmonitor placement on output from three different accelerometers[J].Medicine & Science in Sports & Exercise, 1999(31):142.

    [21] CHEN K Y, SUN M. Improving energy expenditure estimation by using a triaxial accelerometer[J].Journal of Applied Physiology, 1997, 83(6):2112.

    [22] MATTHEWS C E, HAGSTR?MER M, POBER D M, et al. Best practices for using physical activity monitors in population-based research[J]. Medicine and science in sports and exercise, 2012, 44(1):68.

    [23] TROST S G, MCIVER K L, PATE R R. Conducting accelerometer-based activity assessments in field-based research[J]. Medicine and science in sports and exercise, 2005, 37(11):531.

    [24] YANG C C, HSU Y L. A review of accelerometry-based wearable motion detectors for physical activity monitoring[J]. Sensors, 2010, 10(8):7772.

    [25] ROSENBERGER M E, HASKELL W L, ALBINALI F, et al. Estimating activity and sedentary behavior from an accelerometer on the hip or wrist[J]. Medicine and science in sports and exercise, 2013, 45(5):964.

    [26] STAUDENMAYER J, HE S, HICKEY A, et al. Methods to estimate aspects of physical activity and sedentary behavior from high-frequency wrist accelerometer measurements[J]. Journal of Applied Physiology, 2015, 119(4):396.

    [27] MONTOYE A H K, PIVARNIK J M, MUDD L M, et al. Wrist-independent energy expenditure prediction models from raw accelerometer data[J].Physiological Measurement, 2016, 37(10):1770.

    [28] VELHO L, FRERY A, GOMES J. Signal Theory[M].London:Springer,2009:13-55.

    [29] WELK G J. Use of accelerometry-based activity monitors to assess physical activity[J]. Physical activity assessments for health-related research, 2002:125.

    [30] MATHIE M J,COSTER A C F,LOVELL N H, et al. Accelerometry: providing an integrated, practical method for long-term, ambulatory monitoring of human movement[J]. Physiological measurement, 2004, 25(2):1.

    [31] UITERWAAL M,GLERUM E B C,BUSSER H J, et al. Ambulatory monitoring of physical activity in working situations, a validation study[J]. Journal of medical engineering & technology, 1998, 22(4):168.

    [32] BRAGE S,WEDDERKOPP N,F(xiàn)RANKS P W, et al. Reexamination of validity and reliability of the CSA monitor in walking and running[J]. Medicine and Science in Sports and Exercise, 2003, 35(8):1447.

    [33] LEVINE J A,BAUKOL P A,WESTERP K P. Validation of the Tracmortriaxial accelerometer system for walking[J]. Medicine and Science in Sports and Exercise, 2001, 33(9):1593.

    [34] HEINZ E A,KUNZE K S,SULISTVO S, et al. Experimental evaluation of variations in primary features used for accelerometric context recognition[M].Heidelberg:Springer, 2003:252.

    [35] VAN L K,GELLERSEN H W. Spine versus porcupine: A study in distributed wearable activity recognition[C]//Hiroshima:ISWC 2004, 2004:142-149.

    [36] KIANI K,SNIJDERA C J,GELSEMA E S. Computerized analysis of daily life motor activity for ambulatory monitoring[J]. Technology and Health Care, 1997, 5(4):307.

    [37] RAVI N, DANDEKAR N, MYSORE P, et al. Activity recognition from accelerometer data[C]//Austin:AAAI,2005:1541.

    [38] MATTHEW C E. Calibration of accelerometer output for adults[J]. Medicine and science in sports and exercise, 2005, 37(11):512.

    [39] CROUTER S E, CLOWERS K G, BASSETT D R. A novel method for using accelerometer data to predict energy expenditure[J]. Journal of applied physiology, 2006, 100(4): 1324.

    [40] NICHOLS J F, MORGAN C G, CHABOT L E, et al. Assessment of physical activity with the Computer Science and Applications, Inc., accelerometer: laboratory versus field validation[J]. Research quarterly for exercise and sport, 2000, 71(1):36.

    [41] BJORNSON K F. Physical activity monitoring in children and youths[J]. Pediatric Physical Therapy, 2005, 17(1): 37.

    [42] WESTERTERP K R, PLASPUI G. Physical activity and human energy expenditure[J]. Current Opinion in Clinical Nutrition & Metabolic Care, 2004, 7(6): 607.

    [43] HOWE C A, STAUDENMAYER J W, FREEDSON P S. Accelerometer prediction of energy expenditure: vector magnitude versus vertical axis[J]. Med Sci Sports Exerc, 2009, 41(12): 2199.

    [44] SCHNELLER M B, PEDERSEN M T, GUPTA N, et al. Validation of five minimally obstructive methods to estimate physical activity energy expenditure in young adults in semi-standardized settings[J]. Sensors, 2015, 15(3): 6133.

    [45] ALTINI M, PENDERS J, VULLERS R, et al. Estimating Energy Expenditure Using Body-Worn Accelerometers: A Comparison of Methods, Sensors Number and Positioning[J]. Biomedical & Health Informatics IEEE Journal, 2015, 19(1):219.

    猜你喜歡
    能量消耗智能手機測量
    智能手機是座礦
    太極拳連續(xù)“云手”運動強度及其能量消耗探究
    中年女性間歇習(xí)練太極拳的強度、能量消耗與間歇恢復(fù)探究分析
    沒別的可吃
    智能手機臉
    英語文摘(2020年5期)2020-09-21 09:26:30
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計算
    假如我是一部智能手機
    趣味(語文)(2018年8期)2018-11-15 08:53:00
    滑動摩擦力的測量與計算
    測量
    香蕉久久夜色| 老熟妇仑乱视频hdxx| 嫩草影视91久久| 91麻豆精品激情在线观看国产| 搡老妇女老女人老熟妇| 日本撒尿小便嘘嘘汇集6| 国内精品久久久久精免费| 99久久九九国产精品国产免费| 人人妻人人澡欧美一区二区| 天天添夜夜摸| 很黄的视频免费| 脱女人内裤的视频| 在线天堂最新版资源| 国产高清激情床上av| 乱人视频在线观看| 美女高潮的动态| 人人妻人人看人人澡| 国产高潮美女av| 18禁国产床啪视频网站| 国产亚洲av嫩草精品影院| 亚洲久久久久久中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产精品sss在线观看| 九色成人免费人妻av| 久久九九热精品免费| 国产综合懂色| 国产极品精品免费视频能看的| 久久久久久久久大av| 老鸭窝网址在线观看| 老司机午夜福利在线观看视频| netflix在线观看网站| 精品电影一区二区在线| 在线免费观看不下载黄p国产 | 成人精品一区二区免费| 老汉色∧v一级毛片| 国产主播在线观看一区二区| 成人无遮挡网站| 亚洲avbb在线观看| 国产乱人伦免费视频| 精华霜和精华液先用哪个| 亚洲性夜色夜夜综合| 18禁黄网站禁片午夜丰满| 怎么达到女性高潮| 嫩草影视91久久| 90打野战视频偷拍视频| 首页视频小说图片口味搜索| 国产精品一区二区三区四区免费观看 | a级毛片a级免费在线| 国产成人aa在线观看| 啦啦啦韩国在线观看视频| 日本黄大片高清| 国产成年人精品一区二区| 日韩欧美免费精品| 男插女下体视频免费在线播放| 亚洲国产中文字幕在线视频| 9191精品国产免费久久| 亚洲中文字幕日韩| 身体一侧抽搐| 欧美av亚洲av综合av国产av| 搡老岳熟女国产| 国产伦在线观看视频一区| 熟妇人妻久久中文字幕3abv| 黄色成人免费大全| 波野结衣二区三区在线 | 精品国产亚洲在线| 久久久久久九九精品二区国产| 亚洲av免费高清在线观看| 嫁个100分男人电影在线观看| www.www免费av| 国产亚洲精品久久久com| 欧美zozozo另类| 欧美成狂野欧美在线观看| 最近在线观看免费完整版| 国产午夜福利久久久久久| 亚洲国产高清在线一区二区三| 亚洲中文字幕一区二区三区有码在线看| 国内精品美女久久久久久| 精品不卡国产一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 91字幕亚洲| 久久精品国产亚洲av香蕉五月| www日本黄色视频网| 波多野结衣高清作品| 在线观看美女被高潮喷水网站 | 18禁黄网站禁片免费观看直播| 黄片小视频在线播放| 九九热线精品视视频播放| 国产精品久久久久久人妻精品电影| 高清毛片免费观看视频网站| 国产欧美日韩精品一区二区| 美女大奶头视频| 久久久久久久精品吃奶| 国产精品一区二区免费欧美| 99久久久亚洲精品蜜臀av| 女生性感内裤真人,穿戴方法视频| 在线免费观看不下载黄p国产 | 国产激情欧美一区二区| 久久人人精品亚洲av| 久久久久九九精品影院| 国产真人三级小视频在线观看| h日本视频在线播放| 亚洲七黄色美女视频| 丰满的人妻完整版| 亚洲精品成人久久久久久| ponron亚洲| 毛片女人毛片| 亚洲av一区综合| 一区二区三区免费毛片| 真人一进一出gif抽搐免费| e午夜精品久久久久久久| 午夜亚洲福利在线播放| 国产成人系列免费观看| 精品国产三级普通话版| 亚洲精品美女久久久久99蜜臀| 国产精品三级大全| 免费看美女性在线毛片视频| 国产真实乱freesex| www.熟女人妻精品国产| 日本黄色片子视频| 中文字幕熟女人妻在线| 亚洲av成人av| 免费大片18禁| 成人av在线播放网站| 黄色片一级片一级黄色片| 欧美中文日本在线观看视频| 色尼玛亚洲综合影院| 床上黄色一级片| 真实男女啪啪啪动态图| 又黄又粗又硬又大视频| 欧美日韩黄片免| 国产不卡一卡二| 免费大片18禁| 一级黄片播放器| 国内精品久久久久精免费| 老汉色∧v一级毛片| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 男插女下体视频免费在线播放| 欧美中文综合在线视频| 国产成+人综合+亚洲专区| 伊人久久大香线蕉亚洲五| 亚洲精品色激情综合| 日韩欧美精品免费久久 | 三级毛片av免费| 夜夜躁狠狠躁天天躁| 在线a可以看的网站| 美女免费视频网站| 中文字幕熟女人妻在线| 欧美色视频一区免费| 高清在线国产一区| 岛国在线观看网站| 亚洲国产精品999在线| a级毛片a级免费在线| 在线观看日韩欧美| 一本精品99久久精品77| 国语自产精品视频在线第100页| 国产不卡一卡二| 日韩欧美免费精品| 老司机午夜十八禁免费视频| 成人欧美大片| 三级国产精品欧美在线观看| 国产精品一及| 亚洲五月天丁香| 欧美日韩亚洲国产一区二区在线观看| 韩国av一区二区三区四区| 久久精品人妻少妇| 哪里可以看免费的av片| 青草久久国产| 国产精品久久电影中文字幕| 少妇高潮的动态图| 啦啦啦免费观看视频1| 日韩欧美免费精品| 国产97色在线日韩免费| netflix在线观看网站| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 久久性视频一级片| 男女那种视频在线观看| 亚洲人成网站在线播放欧美日韩| 男人舔女人下体高潮全视频| 亚洲av第一区精品v没综合| 99久久精品国产亚洲精品| 啦啦啦韩国在线观看视频| 国产老妇女一区| 国产私拍福利视频在线观看| 午夜老司机福利剧场| 高清在线国产一区| 在线a可以看的网站| 久久久成人免费电影| 女人十人毛片免费观看3o分钟| 欧美中文日本在线观看视频| 国产伦一二天堂av在线观看| 日本与韩国留学比较| 亚洲性夜色夜夜综合| aaaaa片日本免费| 午夜免费男女啪啪视频观看 | 国产伦人伦偷精品视频| xxxwww97欧美| 在线国产一区二区在线| 嫩草影视91久久| 听说在线观看完整版免费高清| 国产精品嫩草影院av在线观看 | 精品电影一区二区在线| 国产av不卡久久| 日韩欧美精品免费久久 | 特大巨黑吊av在线直播| 色综合婷婷激情| 国产一级毛片七仙女欲春2| 午夜影院日韩av| 国产av不卡久久| 午夜福利视频1000在线观看| 熟女少妇亚洲综合色aaa.| 一个人看视频在线观看www免费 | 男女视频在线观看网站免费| 成年人黄色毛片网站| 久久久精品大字幕| 一区二区三区激情视频| 搡老岳熟女国产| 变态另类丝袜制服| 好男人在线观看高清免费视频| 国产精品久久电影中文字幕| 一本综合久久免费| av中文乱码字幕在线| 黄色女人牲交| 亚洲精华国产精华精| 级片在线观看| 又紧又爽又黄一区二区| 亚洲精品亚洲一区二区| 国产一区二区激情短视频| 在线观看免费午夜福利视频| 亚洲人成网站在线播| 又黄又爽又免费观看的视频| 丁香欧美五月| 亚洲黑人精品在线| 日本一本二区三区精品| 午夜两性在线视频| 久久久成人免费电影| 一夜夜www| 麻豆国产av国片精品| 法律面前人人平等表现在哪些方面| 国产黄片美女视频| 美女大奶头视频| 乱人视频在线观看| 特大巨黑吊av在线直播| 黄色成人免费大全| 精品日产1卡2卡| 天堂av国产一区二区熟女人妻| 91字幕亚洲| 国产伦人伦偷精品视频| 国产亚洲欧美98| 韩国av一区二区三区四区| 亚洲精品成人久久久久久| 欧美成人a在线观看| 亚洲av电影在线进入| 午夜激情福利司机影院| 欧美乱妇无乱码| 中文字幕人妻丝袜一区二区| 欧美成人a在线观看| 女人被狂操c到高潮| 欧美绝顶高潮抽搐喷水| av女优亚洲男人天堂| 国产精品 欧美亚洲| 欧美区成人在线视频| 岛国视频午夜一区免费看| 国产极品精品免费视频能看的| 午夜福利在线观看免费完整高清在 | 亚洲av日韩精品久久久久久密| 极品教师在线免费播放| 欧美av亚洲av综合av国产av| 精品国产美女av久久久久小说| 精品国产亚洲在线| 国产伦精品一区二区三区视频9 | 一卡2卡三卡四卡精品乱码亚洲| 天天一区二区日本电影三级| 国产一区二区在线观看日韩 | 桃色一区二区三区在线观看| 99久久99久久久精品蜜桃| 欧美av亚洲av综合av国产av| 亚洲国产欧美人成| 精品熟女少妇八av免费久了| 国产激情偷乱视频一区二区| 老司机在亚洲福利影院| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 看黄色毛片网站| 欧美在线黄色| 毛片女人毛片| 中文字幕高清在线视频| 亚洲成人精品中文字幕电影| 女人十人毛片免费观看3o分钟| 精品久久久久久久久久久久久| 亚洲成人久久性| 国产一区二区三区在线臀色熟女| 国产精品嫩草影院av在线观看 | 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 18禁国产床啪视频网站| 国产色婷婷99| 999久久久精品免费观看国产| 亚洲av美国av| 国语自产精品视频在线第100页| 最好的美女福利视频网| 淫秽高清视频在线观看| 成人性生交大片免费视频hd| 国产精品久久电影中文字幕| 18禁裸乳无遮挡免费网站照片| 国产精品,欧美在线| 国模一区二区三区四区视频| 久久久国产成人精品二区| 午夜久久久久精精品| 一边摸一边抽搐一进一小说| 欧美色欧美亚洲另类二区| 中文在线观看免费www的网站| 亚洲最大成人中文| 婷婷六月久久综合丁香| 久9热在线精品视频| 国产亚洲精品av在线| 无遮挡黄片免费观看| 丝袜美腿在线中文| 亚洲天堂国产精品一区在线| 91久久精品电影网| 极品教师在线免费播放| 熟女少妇亚洲综合色aaa.| 国产99白浆流出| 99在线人妻在线中文字幕| 两人在一起打扑克的视频| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 性欧美人与动物交配| 三级男女做爰猛烈吃奶摸视频| 丰满人妻一区二区三区视频av | 久久久色成人| 国产精品99久久久久久久久| 国产精品一及| 亚洲成a人片在线一区二区| h日本视频在线播放| 色综合站精品国产| 亚洲国产精品久久男人天堂| 久久久久久人人人人人| 丝袜美腿在线中文| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 毛片女人毛片| 一个人看视频在线观看www免费 | 免费人成在线观看视频色| 18禁黄网站禁片免费观看直播| 好男人在线观看高清免费视频| 少妇熟女aⅴ在线视频| 欧美黑人欧美精品刺激| 少妇的丰满在线观看| 亚洲一区二区三区色噜噜| 一夜夜www| 国产精品电影一区二区三区| 日韩 欧美 亚洲 中文字幕| 好男人在线观看高清免费视频| 国产探花极品一区二区| 一级a爱片免费观看的视频| 窝窝影院91人妻| 麻豆成人av在线观看| 哪里可以看免费的av片| 欧美一区二区亚洲| 18禁美女被吸乳视频| 尤物成人国产欧美一区二区三区| 1000部很黄的大片| 欧美成人性av电影在线观看| 欧美日韩福利视频一区二区| 男女视频在线观看网站免费| 国产色婷婷99| 欧美一区二区精品小视频在线| 日本免费a在线| 成年女人毛片免费观看观看9| 又紧又爽又黄一区二区| 深夜精品福利| 美女高潮的动态| 香蕉久久夜色| 真实男女啪啪啪动态图| 亚洲中文字幕日韩| 成人特级av手机在线观看| 老司机福利观看| 国产午夜福利久久久久久| 综合色av麻豆| 国产精品 欧美亚洲| 久久欧美精品欧美久久欧美| 欧美乱妇无乱码| 亚洲国产精品合色在线| 亚洲av日韩精品久久久久久密| 高清在线国产一区| 欧美日本亚洲视频在线播放| 青草久久国产| 在线十欧美十亚洲十日本专区| 国产av不卡久久| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 欧美色视频一区免费| 国产精品一区二区三区四区免费观看 | 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 国产亚洲精品av在线| 99国产综合亚洲精品| 国产探花极品一区二区| 最近最新免费中文字幕在线| 欧美三级亚洲精品| 51国产日韩欧美| 欧美又色又爽又黄视频| 亚洲第一欧美日韩一区二区三区| 国产精品香港三级国产av潘金莲| 国产在视频线在精品| 美女高潮喷水抽搐中文字幕| 精品久久久久久久人妻蜜臀av| 一个人免费在线观看的高清视频| 日本熟妇午夜| 午夜福利在线观看吧| 欧美三级亚洲精品| 51国产日韩欧美| 日本与韩国留学比较| 国产v大片淫在线免费观看| 少妇丰满av| 久久伊人香网站| 香蕉av资源在线| 欧美一级a爱片免费观看看| 国产精品日韩av在线免费观看| 免费av观看视频| 超碰av人人做人人爽久久 | 五月伊人婷婷丁香| 精品久久久久久久久久久久久| 99精品欧美一区二区三区四区| 成人一区二区视频在线观看| 国内毛片毛片毛片毛片毛片| 99久久九九国产精品国产免费| 一边摸一边抽搐一进一小说| 成人av一区二区三区在线看| 免费在线观看成人毛片| 少妇丰满av| 最近最新中文字幕大全电影3| 欧美日韩国产亚洲二区| 不卡一级毛片| 国产精品一区二区三区四区久久| 一级黄片播放器| 九九在线视频观看精品| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看| 欧美在线一区亚洲| 亚洲在线自拍视频| 国内精品美女久久久久久| 免费看光身美女| 老汉色av国产亚洲站长工具| 中文在线观看免费www的网站| 美女免费视频网站| 在线免费观看不下载黄p国产 | 中文字幕熟女人妻在线| 亚洲在线自拍视频| 国产麻豆成人av免费视频| 久久精品国产亚洲av涩爱 | 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 女人高潮潮喷娇喘18禁视频| 久久久久国内视频| 亚洲专区国产一区二区| 久久久久久久精品吃奶| 成人高潮视频无遮挡免费网站| 亚洲狠狠婷婷综合久久图片| 欧美高清成人免费视频www| 国产精品三级大全| 小蜜桃在线观看免费完整版高清| 中文字幕熟女人妻在线| 国产伦在线观看视频一区| 俄罗斯特黄特色一大片| 欧美乱色亚洲激情| 欧美激情在线99| av福利片在线观看| 最近视频中文字幕2019在线8| 12—13女人毛片做爰片一| 精品久久久久久久末码| 19禁男女啪啪无遮挡网站| 国内精品美女久久久久久| 成人三级黄色视频| 亚洲成人中文字幕在线播放| 国产精品女同一区二区软件 | 两人在一起打扑克的视频| 久久精品国产综合久久久| 亚洲最大成人手机在线| 国产中年淑女户外野战色| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自拍偷在线| 观看美女的网站| 亚洲国产精品久久男人天堂| 国产精品99久久99久久久不卡| 91麻豆精品激情在线观看国产| 18禁黄网站禁片午夜丰满| 精品久久久久久久人妻蜜臀av| 女人高潮潮喷娇喘18禁视频| 嫩草影院入口| 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡| 俄罗斯特黄特色一大片| 久久久久精品国产欧美久久久| 1000部很黄的大片| 免费大片18禁| 欧美大码av| 波多野结衣高清无吗| 日日干狠狠操夜夜爽| 国产探花在线观看一区二区| 日本与韩国留学比较| 亚洲五月天丁香| 国产精品野战在线观看| 18禁黄网站禁片免费观看直播| 中国美女看黄片| 午夜日韩欧美国产| 色综合亚洲欧美另类图片| 99久久精品国产亚洲精品| 不卡一级毛片| 天天躁日日操中文字幕| 韩国av一区二区三区四区| 国产精品,欧美在线| 日本撒尿小便嘘嘘汇集6| av在线蜜桃| 日韩欧美精品免费久久 | 天天一区二区日本电影三级| 国产精品一区二区三区四区久久| 法律面前人人平等表现在哪些方面| 亚洲成人精品中文字幕电影| 亚洲成人久久爱视频| 一区二区三区激情视频| 99久久无色码亚洲精品果冻| 好看av亚洲va欧美ⅴa在| 国产成人福利小说| 亚洲在线自拍视频| 久久久成人免费电影| 婷婷丁香在线五月| 亚洲精品一区av在线观看| 国产真实乱freesex| 熟女电影av网| 成人无遮挡网站| 亚洲av五月六月丁香网| 少妇裸体淫交视频免费看高清| 欧美最黄视频在线播放免费| 久久中文看片网| 国产久久久一区二区三区| 色噜噜av男人的天堂激情| 免费在线观看成人毛片| 伊人久久大香线蕉亚洲五| 国产淫片久久久久久久久 | 老司机深夜福利视频在线观看| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 国产91精品成人一区二区三区| 亚洲avbb在线观看| 一级a爱片免费观看的视频| xxxwww97欧美| 欧美+亚洲+日韩+国产| 亚洲人与动物交配视频| 亚洲成人久久性| 午夜精品在线福利| 日本黄大片高清| 日本一二三区视频观看| 变态另类成人亚洲欧美熟女| 69av精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 舔av片在线| 9191精品国产免费久久| 色综合站精品国产| 亚洲欧美日韩高清在线视频| 97超视频在线观看视频| 国产又黄又爽又无遮挡在线| 国产蜜桃级精品一区二区三区| 啦啦啦免费观看视频1| 亚洲av一区综合| 18+在线观看网站| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 长腿黑丝高跟| 男女那种视频在线观看| 免费电影在线观看免费观看| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 国产又黄又爽又无遮挡在线| 亚洲,欧美精品.| 高潮久久久久久久久久久不卡| 99在线人妻在线中文字幕| 看免费av毛片| 变态另类丝袜制服| 波多野结衣巨乳人妻| 丰满乱子伦码专区| 色综合婷婷激情| 舔av片在线| 国产免费男女视频| 久久6这里有精品| 国产精品一区二区三区四区免费观看 | 在线播放无遮挡| 国产伦人伦偷精品视频| 又黄又爽又免费观看的视频| 欧美一区二区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 男女之事视频高清在线观看| 国产精品一区二区三区四区免费观看 | 亚洲18禁久久av| 中文字幕久久专区| 九九热线精品视视频播放| 久久精品亚洲精品国产色婷小说| 午夜亚洲福利在线播放| 亚洲av五月六月丁香网| 午夜免费激情av| 国产极品精品免费视频能看的| 成人性生交大片免费视频hd| 成人国产一区最新在线观看| 国产97色在线日韩免费| 国产精品香港三级国产av潘金莲| 亚洲精品乱码久久久v下载方式 | 日本在线视频免费播放| 国产精品一及| 国产激情欧美一区二区| a级毛片a级免费在线| 好看av亚洲va欧美ⅴa在| 日韩av在线大香蕉| 在线观看免费视频日本深夜|