• <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
    滑動摩擦力的測量與計算
    測量
    村上凉子中文字幕在线| 亚洲av日韩精品久久久久久密| 淫秽高清视频在线观看| 国产成人影院久久av| 色综合欧美亚洲国产小说| 少妇的逼水好多| 在线a可以看的网站| 日韩成人在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 青草久久国产| 又爽又黄a免费视频| 特级一级黄色大片| 午夜影院日韩av| 精品久久国产蜜桃| 精品熟女少妇八av免费久了| 日本在线视频免费播放| 又黄又爽又免费观看的视频| 18禁黄网站禁片免费观看直播| 午夜福利免费观看在线| 麻豆成人午夜福利视频| 国产色爽女视频免费观看| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩| 少妇裸体淫交视频免费看高清| 久久久国产成人免费| 国内久久婷婷六月综合欲色啪| 亚洲无线观看免费| 国内少妇人妻偷人精品xxx网站| 热99在线观看视频| 99久久无色码亚洲精品果冻| 国产精品久久电影中文字幕| 午夜免费男女啪啪视频观看 | 18+在线观看网站| 精品一区二区三区人妻视频| 日韩欧美三级三区| 给我免费播放毛片高清在线观看| 精品久久久久久久久久久久久| 日韩国内少妇激情av| 亚洲 国产 在线| 天堂动漫精品| 一本久久中文字幕| 最近最新中文字幕大全电影3| 婷婷亚洲欧美| av欧美777| 两个人视频免费观看高清| 最新在线观看一区二区三区| 欧美精品啪啪一区二区三区| 欧美bdsm另类| 欧美zozozo另类| 又爽又黄无遮挡网站| 中文字幕av在线有码专区| 三级国产精品欧美在线观看| 女生性感内裤真人,穿戴方法视频| 成年人黄色毛片网站| 欧美国产日韩亚洲一区| 变态另类成人亚洲欧美熟女| 亚洲,欧美精品.| 久久天躁狠狠躁夜夜2o2o| 欧美极品一区二区三区四区| 黄色日韩在线| 国产91精品成人一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲人成网站高清观看| 1000部很黄的大片| 男女国产视频网站| a级毛色黄片| 免费看不卡的av| 亚洲精品aⅴ在线观看| 国产成人精品婷婷| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区四那| 成人毛片60女人毛片免费| 啦啦啦在线观看免费高清www| 久久久久九九精品影院| 免费人成在线观看视频色| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 一区二区三区免费毛片| 直男gayav资源| 亚洲国产日韩一区二区| 亚洲美女视频黄频| 91久久精品国产一区二区成人| 国产精品一区www在线观看| 嘟嘟电影网在线观看| 涩涩av久久男人的天堂| 国产日韩欧美在线精品| 亚洲国产日韩一区二区| 日韩av不卡免费在线播放| 亚洲av日韩在线播放| 亚洲精品第二区| 国产精品久久久久久精品电影小说 | 熟女电影av网| 午夜日本视频在线| 精品一区二区免费观看| 舔av片在线| 一区二区三区免费毛片| a级毛色黄片| 国产一区有黄有色的免费视频| 亚洲精品一二三| 色网站视频免费| 亚洲激情五月婷婷啪啪| 亚洲国产av新网站| 嫩草影院精品99| 2021少妇久久久久久久久久久| 国产高清三级在线| 麻豆久久精品国产亚洲av| 天堂网av新在线| 免费观看av网站的网址| 男女国产视频网站| 干丝袜人妻中文字幕| 人人妻人人爽人人添夜夜欢视频 | 欧美激情国产日韩精品一区| 一区二区av电影网| 亚洲国产精品成人综合色| 欧美丝袜亚洲另类| 国产91av在线免费观看| 69人妻影院| 亚洲四区av| 国产亚洲一区二区精品| 精品一区二区三卡| 亚洲欧美精品自产自拍| 日产精品乱码卡一卡2卡三| 七月丁香在线播放| 久久韩国三级中文字幕| 婷婷色av中文字幕| 亚洲无线观看免费| 亚洲av中文av极速乱| 亚洲av免费在线观看| 大片免费播放器 马上看| 狂野欧美白嫩少妇大欣赏| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩东京热| 精品人妻偷拍中文字幕| 精品久久久久久电影网| 观看美女的网站| 成年人午夜在线观看视频| 欧美 日韩 精品 国产| 久热这里只有精品99| 最近中文字幕2019免费版| 欧美日韩在线观看h| 男女边吃奶边做爰视频| 国精品久久久久久国模美| 丝袜脚勾引网站| 欧美成人精品欧美一级黄| 蜜桃亚洲精品一区二区三区| 国产日韩欧美亚洲二区| 伦精品一区二区三区| 边亲边吃奶的免费视频| 日韩精品有码人妻一区| 美女xxoo啪啪120秒动态图| 亚洲国产精品成人久久小说| 精品亚洲乱码少妇综合久久| 91狼人影院| 久久久久久九九精品二区国产| 久久久成人免费电影| 99九九线精品视频在线观看视频| 在现免费观看毛片| 18禁动态无遮挡网站| 又黄又爽又刺激的免费视频.| 麻豆成人av视频| 欧美日韩在线观看h| 国产男人的电影天堂91| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看 | 成人高潮视频无遮挡免费网站| 少妇人妻 视频| 日本欧美国产在线视频| 在线精品无人区一区二区三 | 嫩草影院入口| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| 在线看a的网站| 亚洲天堂国产精品一区在线| 色5月婷婷丁香| 精品国产露脸久久av麻豆| 日韩欧美 国产精品| 一本色道久久久久久精品综合| 亚洲国产精品成人综合色| 国产精品熟女久久久久浪| 久久精品国产亚洲av天美| 国产探花极品一区二区| 午夜福利在线在线| 亚洲自偷自拍三级| 97热精品久久久久久| 亚洲精品国产av蜜桃| 亚洲国产精品999| 午夜福利高清视频| 一边亲一边摸免费视频| 18禁在线无遮挡免费观看视频| 91精品国产九色| 91精品一卡2卡3卡4卡| 男人舔奶头视频| 有码 亚洲区| 国产一级毛片在线| 纵有疾风起免费观看全集完整版| 色播亚洲综合网| 晚上一个人看的免费电影| 亚洲性久久影院| 伊人久久国产一区二区| 男插女下体视频免费在线播放| 久久久久性生活片| 大又大粗又爽又黄少妇毛片口| 最近中文字幕2019免费版| 国产精品秋霞免费鲁丝片| 一二三四中文在线观看免费高清| 蜜桃亚洲精品一区二区三区| 久久亚洲国产成人精品v| 成人毛片a级毛片在线播放| 美女主播在线视频| 久久久精品欧美日韩精品| 婷婷色综合www| 国产高清三级在线| 午夜福利在线观看免费完整高清在| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 乱系列少妇在线播放| 亚洲国产欧美人成| 国产一区亚洲一区在线观看| 观看美女的网站| 国产人妻一区二区三区在| 伊人久久国产一区二区| 中国美白少妇内射xxxbb| 日本熟妇午夜| 久久精品国产亚洲av涩爱| 免费观看在线日韩| 国产乱人视频| 国产精品久久久久久精品电影| 777米奇影视久久| 亚洲av一区综合| 2021少妇久久久久久久久久久| 久久久久久久久久久丰满| 日韩视频在线欧美| 成人综合一区亚洲| 中文资源天堂在线| 国产精品一区二区性色av| 一级片'在线观看视频| 天天躁日日操中文字幕| 亚洲av在线观看美女高潮| 国产毛片a区久久久久| 国产一区二区三区av在线| 激情五月婷婷亚洲| 国产成人a∨麻豆精品| 久久精品久久久久久久性| 国产又色又爽无遮挡免| av网站免费在线观看视频| 日韩制服骚丝袜av| 一二三四中文在线观看免费高清| 成年免费大片在线观看| 熟女av电影| 中文欧美无线码| 观看美女的网站| 如何舔出高潮| 国产欧美亚洲国产| 夜夜看夜夜爽夜夜摸| 国产黄频视频在线观看| 夫妻午夜视频| 亚洲人与动物交配视频| 亚洲在久久综合| 精品一区二区免费观看| 春色校园在线视频观看| 男女无遮挡免费网站观看| 男人舔奶头视频| 久久久久久久精品精品| 欧美极品一区二区三区四区| 久久精品熟女亚洲av麻豆精品| 成年女人看的毛片在线观看| 亚洲一级一片aⅴ在线观看| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 亚洲成色77777| 亚洲熟女精品中文字幕| 99久久精品国产国产毛片| 亚洲性久久影院| 爱豆传媒免费全集在线观看| 午夜精品国产一区二区电影 | 亚洲av.av天堂| 久久人人爽人人片av| 国产av码专区亚洲av| 五月玫瑰六月丁香| 韩国av在线不卡| 精品酒店卫生间| 国产在线一区二区三区精| 51国产日韩欧美| 大香蕉久久网| 各种免费的搞黄视频| 18禁在线无遮挡免费观看视频| 精品视频人人做人人爽| 日日啪夜夜撸| 久久人人爽av亚洲精品天堂 | 性插视频无遮挡在线免费观看| videos熟女内射| 三级国产精品欧美在线观看| 精品一区在线观看国产| 成人二区视频| 九九在线视频观看精品| 国产亚洲最大av| 女的被弄到高潮叫床怎么办| 身体一侧抽搐| 在线观看国产h片| 亚洲精品乱码久久久v下载方式| 久久女婷五月综合色啪小说 | 在线免费观看不下载黄p国产| 久久久久久久久久成人| 久久国产乱子免费精品| 亚洲va在线va天堂va国产| 天堂网av新在线| 在线精品无人区一区二区三 | 99久久精品热视频| 大片免费播放器 马上看| 久久久久性生活片| 免费看不卡的av| 97人妻精品一区二区三区麻豆| 免费观看性生交大片5| 观看免费一级毛片| 中文天堂在线官网| 亚洲国产欧美在线一区| 免费黄色在线免费观看| 久久这里有精品视频免费| 三级男女做爰猛烈吃奶摸视频| 一级片'在线观看视频| 99久国产av精品国产电影| 精品99又大又爽又粗少妇毛片| 人妻制服诱惑在线中文字幕| 精品久久久久久电影网| 亚洲精品国产成人久久av| 婷婷色综合www| 秋霞伦理黄片| 黄片wwwwww| 99久国产av精品国产电影| 国产精品国产av在线观看| 精华霜和精华液先用哪个| 国产精品嫩草影院av在线观看| 久热这里只有精品99| 久久久久久久久大av| 狂野欧美白嫩少妇大欣赏| 日韩,欧美,国产一区二区三区| 国产男女内射视频| 国产精品国产三级国产专区5o| 99精国产麻豆久久婷婷| 日本色播在线视频| 国产在视频线精品| 99re6热这里在线精品视频| 亚洲av欧美aⅴ国产| 久久99热这里只频精品6学生| 日本黄色片子视频| 九草在线视频观看| 国产成人福利小说| 在线免费观看不下载黄p国产| 午夜免费观看性视频| 男人舔奶头视频| 18禁裸乳无遮挡免费网站照片| 久久久精品免费免费高清| 午夜日本视频在线| 国语对白做爰xxxⅹ性视频网站| 99久久精品一区二区三区| 美女主播在线视频| av女优亚洲男人天堂| 成人高潮视频无遮挡免费网站| 亚洲精品国产色婷婷电影| 国产亚洲精品久久久com| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 卡戴珊不雅视频在线播放| 亚洲精华国产精华液的使用体验| 亚洲国产精品专区欧美| 久久久久国产精品人妻一区二区| 国产男女超爽视频在线观看| av在线天堂中文字幕| 久久久a久久爽久久v久久| 欧美日韩精品成人综合77777| 狂野欧美白嫩少妇大欣赏| 欧美区成人在线视频| 18禁在线播放成人免费| 精品久久久久久电影网| 午夜激情福利司机影院| 亚洲精品aⅴ在线观看| 久久久久久久久久成人| 日产精品乱码卡一卡2卡三| 国产精品秋霞免费鲁丝片| 久久久精品欧美日韩精品| 简卡轻食公司| 国产一级毛片在线| 欧美xxxx性猛交bbbb| av线在线观看网站| 一级毛片aaaaaa免费看小| 肉色欧美久久久久久久蜜桃 | 男女那种视频在线观看| 中文字幕久久专区| 热re99久久精品国产66热6| 晚上一个人看的免费电影| 性色av一级| 国产视频内射| 婷婷色av中文字幕| 成人美女网站在线观看视频| 日韩国内少妇激情av| 91午夜精品亚洲一区二区三区| 一区二区av电影网| 国产精品女同一区二区软件| 性色avwww在线观看| 国产精品久久久久久久久免| 国产探花在线观看一区二区| 国产午夜福利久久久久久| 一区二区av电影网| 日韩成人av中文字幕在线观看| 成人二区视频| 国产成人午夜福利电影在线观看| 欧美区成人在线视频| 免费av毛片视频| 午夜福利视频精品| 一级黄片播放器| 成人综合一区亚洲| 精品熟女少妇av免费看| 精品国产三级普通话版| 精品久久久久久久末码| 少妇熟女欧美另类| 日韩欧美一区视频在线观看 | 在线观看三级黄色| 国产精品久久久久久精品电影小说 | 日韩在线高清观看一区二区三区| 国产男女超爽视频在线观看| 哪个播放器可以免费观看大片| 日日撸夜夜添| 午夜免费鲁丝| 女的被弄到高潮叫床怎么办| 国产欧美日韩精品一区二区| 色网站视频免费| 少妇人妻久久综合中文| 黄色怎么调成土黄色| 久久精品久久久久久久性| 中文资源天堂在线| 日韩人妻高清精品专区| 有码 亚洲区| 亚洲精品成人av观看孕妇| 欧美变态另类bdsm刘玥| 久久精品人妻少妇| 精品国产三级普通话版| 91在线精品国自产拍蜜月| 狂野欧美白嫩少妇大欣赏| 国产 一区 欧美 日韩| av福利片在线观看| 亚洲伊人久久精品综合| 最近中文字幕高清免费大全6| 国产精品久久久久久久久免| 亚洲av二区三区四区| 欧美日韩在线观看h| 美女内射精品一级片tv| 狠狠精品人妻久久久久久综合| 亚洲av二区三区四区| 一级毛片 在线播放| 好男人视频免费观看在线| 九色成人免费人妻av| 亚洲高清免费不卡视频| 赤兔流量卡办理| 国产永久视频网站| 偷拍熟女少妇极品色| 高清在线视频一区二区三区| 亚洲国产精品成人久久小说| 亚洲久久久久久中文字幕| 特级一级黄色大片| 亚洲真实伦在线观看| 综合色丁香网| 建设人人有责人人尽责人人享有的 | 纵有疾风起免费观看全集完整版| 免费av毛片视频| av女优亚洲男人天堂| 尾随美女入室| 日韩av在线免费看完整版不卡| 亚洲欧美日韩卡通动漫| 欧美 日韩 精品 国产| 黄色日韩在线| 最近最新中文字幕免费大全7| 国产美女午夜福利| 欧美潮喷喷水| 免费播放大片免费观看视频在线观看| videossex国产| 99热这里只有是精品50| 制服丝袜香蕉在线| 九色成人免费人妻av| 性色av一级| 国产欧美亚洲国产| 自拍偷自拍亚洲精品老妇| 国产精品国产三级国产av玫瑰| 亚洲国产精品成人综合色| 国产精品国产三级国产专区5o| 国产伦理片在线播放av一区| 免费看a级黄色片| 成人无遮挡网站| 国产在线一区二区三区精| 大陆偷拍与自拍| 免费大片18禁| 中文字幕免费在线视频6| 亚洲经典国产精华液单| 有码 亚洲区| 免费看日本二区| 国产乱人视频| 色视频在线一区二区三区| 天堂网av新在线| 91aial.com中文字幕在线观看| 久久午夜福利片| 一个人看视频在线观看www免费| www.色视频.com| 少妇人妻精品综合一区二区| 亚洲成人av在线免费| 国产精品人妻久久久影院| 国产欧美日韩一区二区三区在线 | 97超碰精品成人国产| 搡女人真爽免费视频火全软件| 男人爽女人下面视频在线观看| 夫妻性生交免费视频一级片| 亚洲国产欧美人成| 亚洲天堂国产精品一区在线| 久久精品人妻少妇| 色播亚洲综合网| 亚洲在久久综合| 美女国产视频在线观看| 国产精品久久久久久精品古装| 亚洲av日韩在线播放| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 亚洲av成人精品一二三区| 国产亚洲一区二区精品| 亚洲第一区二区三区不卡| 久久人人爽av亚洲精品天堂 | av线在线观看网站| freevideosex欧美| 国产免费一区二区三区四区乱码| 欧美人与善性xxx| 亚洲美女视频黄频| 欧美丝袜亚洲另类| 国产男女超爽视频在线观看| 国产精品女同一区二区软件| 国产日韩欧美在线精品| 国产精品秋霞免费鲁丝片| 久久人人爽人人片av| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 国产又色又爽无遮挡免| 美女国产视频在线观看| 热re99久久精品国产66热6| 男人舔奶头视频| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区三区四区免费观看| 2018国产大陆天天弄谢| 亚洲精品乱码久久久v下载方式| 久久99蜜桃精品久久| 男插女下体视频免费在线播放| 久久久久久久久久人人人人人人| 成人国产av品久久久| 天美传媒精品一区二区| 男人添女人高潮全过程视频| 国产精品一及| 男人狂女人下面高潮的视频| av播播在线观看一区| 热99国产精品久久久久久7| 国产黄片美女视频| 日韩在线高清观看一区二区三区| 高清av免费在线| 国产成人精品一,二区| 一本色道久久久久久精品综合| 热re99久久精品国产66热6| 亚洲av国产av综合av卡| 久久久色成人| 久久久久久久久久人人人人人人| 一个人观看的视频www高清免费观看| 国产又色又爽无遮挡免| 噜噜噜噜噜久久久久久91| 亚洲人成网站高清观看| 亚洲av.av天堂| 干丝袜人妻中文字幕| 亚洲av中文av极速乱| 最近的中文字幕免费完整| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 偷拍熟女少妇极品色| 尾随美女入室| 麻豆国产97在线/欧美| 国产伦精品一区二区三区四那| 免费电影在线观看免费观看| 久久精品国产亚洲av天美| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 国产高潮美女av| 一级毛片我不卡| 精品国产一区二区三区久久久樱花 | 中文字幕免费在线视频6| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站 | av黄色大香蕉| 九九久久精品国产亚洲av麻豆| 免费看日本二区| 少妇被粗大猛烈的视频| 国产精品国产三级国产专区5o| 久久久亚洲精品成人影院| 丝袜脚勾引网站| av在线蜜桃| 黄色视频在线播放观看不卡| 免费观看av网站的网址| 日韩欧美精品免费久久| 精品酒店卫生间| 久久午夜福利片| 伦理电影大哥的女人| videossex国产| 国产精品.久久久| 一区二区三区免费毛片| 国产一区二区三区综合在线观看 | 亚洲丝袜综合中文字幕| 精品午夜福利在线看| 在线观看一区二区三区激情| www.色视频.com| 高清在线视频一区二区三区| 最近最新中文字幕免费大全7| 亚洲精品一区蜜桃| 亚洲综合色惰| 亚洲国产精品国产精品|