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

    A Smart Wearable Ring Device for Sensing Hand Tremor of Parkinson’s Patients

    2021-04-28 05:03:10HaixiaYangYixianShenWeiZhuangChunmingGaoDongDaiandWeigongZhang

    Haixia Yang,Yixian Shen,Wei Zhuang,Chunming Gao,Dong Dai and Weigong Zhang,*

    1School of Instrument Science and Engineering,Southeast University,Nanjing,210096,China

    2School of Computer and Software,Nanjing University of Information Science and Technology,Nanjing,210044,China

    3School of Engineering and Technology,University of Washington,Tacoma,98402,USA

    4School of Cyber Science and Engineering,Southeast University,Nanjing,210096,China

    ABSTRACT Parkinson’s disease(PD)is a very common neurodegenerative disease that occurs mostly in the elderly.There are many main clinical manifestations of PD,such as tremor,bradykinesia,muscle rigidity,etc.Based on the current research on PD,the accurate and convenient detection of early symptoms is the key to detect PD.With the development of microelectronic and sensor technology,it is much easier to measure the barely noticeable tremor in just one hand for the early detection of Parkinson’s disease.In this paper,we present a smart wearable device for detecting hand tremor,in which MPU6050(MIDI Processing Unit)consisting of a 3-axis gyroscope and a 3-axis accelerometer is used to collect acceleration and angular velocity of fingers.By analyzing the time of specific finger movements,we successfully recognized the tremor signals with high accuracy.Meanwhile,with Bluetooth 4.0(Bluetooth Low Energy,BLE)and networking terminal ability,tremor data can be transferred to a monitoring device in real time with extremely low energy consumption.The experimental results have shown that the proposed device(smart ring)is convenient for long-term tremor detection which is vital for early detection and treatment for Parkinson’s disease.

    KEYWORDS Parkinson’s disease;tremor detection;IMU sensor;MPU6050;BLE

    1 Introduction

    Parkinson’s disease(PD)is a progressive nervous system disorder that affects movements of human body.Symptoms start gradually,sometimes starting with a barely noticeable tremor in just one hand.Tremors are common,but the disorder usually causes stiffness or bradykinesia[1].The elderly is at a higher risk of Parkinson’s disease compared with young people which causes great difficulties to their life.

    Tremor,especially static tremor,is usually the first symptom of Parkinson’s disease.The“static” refers to a state of quietness or relaxation of muscles of the body.These tremors look like as if a person is rubbing clay balls or counting money.They appear only occasionally at first,and get worse gradually.

    The main symptoms of Parkinson’s disease are distinguishable,especially the tremors which serve as the prelude and can be used to foresee an onslaught of Parkinson’s disease[2–4].In recent decades,the detection of tremor signals of Parkinson’s disease has attracted enthusiastic attentions and many detection methods have been proposed in academia and industry[5–7].

    Devices with builtin accelerometers and gyroscopes have been used to assist detecting tremors in PD patients.The study of using accelerometer to detect tremor signals is evolved gradually with the development of electronic techniques[8,9].In 2010,in the United States,Robert LeMoyne used builtin accelerometer in iPhone to detect tremor signals.By collecting and analyzing the waveform of tremors from a Parkinson’s patient,it is demonstrated that the accelerometer in iPhone can quantify the characteristics of tremor[10].Oess et al.[11]designed a smart glove equipped with gyroscope which can trace the motion of hand and analyze the motion waveform to detect tremor.Paola Pierlenoni designed a bracelet which is used to detect tremor signal;it uses the FFT(Fast Fourier Transform)and autoregressive parameter model spectrum to plot the power spectrum and identifies the features of tremor[12].Butson et al.[13]proposed a standardized method to classify motions based on motion programming DBS(Deep Brain Stimulation)parameters.In their study,subjects were asked to wear motion sensors while conducting DBS coded standard procedures to collect kinematics data.Researchers could use this method to better understand the parameters of motions.

    So far,it is a common belief that tremor signals embrace reliable characteristics for early detection of Parkinson’s disease[14,15].However,it remains a challenge about how to make the detection accurate and how to make the device easy to use for the elderly users.Researches are working to develop smaller and more convenient devices for tremor detection[16].In the future,it is inevitable that there will be more demands on small and convenient devices for PD’s early detection.

    In general,wearable devices[17–19]for PD detection have the following features in common:

    (1)Small size and light weight,so that the device can be used as convenient as possible.

    (2)High sensitivity.The device should be able to detect subtle tremors which can prelude the onslaught of PD.

    (3)Low power consumption.The battery capacity of microelectronic device is extremely limited.As a result,the operating power should be as low as possible.

    In this paper,we present our design and implementation of a small ring device to detect hand tremors for Parkinson’s disease patients.Section 2 presents our selection of the system hardware components and Section 3 presents the hardware circuit design.Section 4 explains the tremor detection algorithm design.Section 5 presents our system implementation and the experimental studies.In Section 6,we conclude with future work to fine-tune the smart ring device.

    2 System Hardware Architecture

    2.1 Off-the-Shelf Hardware List

    To minimize the size and cost of the device,we selected off-the-shelf hardware as listed in Tab.1 and explained below.

    Table 1:List of hardware materials

    CC2541 is a Bluetooth low energy and proprietary wireless MCU(Microcontroller),used for data transmission.MPU-6050(MIDI Processing Unit)is a motion tracking device designed for low power consumption,low cost and high-performance,used to obtain data of acceleration and angular velocity.Decoupling capacitor is used for signal filtering.Precision resistor(56 kohm,±1%)is a bias resistor.Main crystal oscillator(32 MHz)is to supply clock signals.Auxiliary crystal oscillator(32.768 kHz)works for low power consumption mode.The button battery(5 V)supplies power.LED lamp is used to indicate the state of device.Reset button is used to reset the device.10 KΩresistor is used to protect the circuit.

    2.2 MPU-6050

    2.2.1 A Brief Introduction of MPU-6050

    MPU6050 is an integrated chip with embedded accelerometer and gyroscope on the same silicon die that can detect acceleration,angular velocity and temperature.We only use acceleration and angular velocity in this study.The data collected were converted from analog signals to digital signals using a 16-bit ADC(Analog-to-Digital Converter).The integrated gyro detection diagram is shown in Fig.1.The pin diagram of MPU-6050 is shown in Fig.2.

    2.2.2 Configuring Register Parameters

    GYRO_CONFIG:XG_ST(self-testing X-axis gyroscope),YG_ST(self-testing Y-axis gyroscope),G_ST(self-testing Z-axis gyroscope)and FS_SEL(gyroscope range selection)are included in this register configuration part.When FS_SEL is 0,the output range is ±250?/s.When FS_SEL is 1,the output range is ±500?/s.When FS_SEL is 2,the output range is ±1000?/s.When FS_SEL is 3,the output range is ±2000?/s.

    ACCEL_CONFIG:XA_ST(self-testing X-axis accelerometer),YA_ST(self-testing Y-axis accelerometer),ZA_ST(self-testing Z-axis accelerometer)and AFS_SEL(accelerometer range selection)are included in this register configuration part.The accelerometer configuration is shown in Tab.2.

    Figure 1:Gyro detection diagram

    Figure 2:Pin diagram of MPU6050

    Table 2:Accelerometer configuration

    The registers used to store the triaxial acceleration values are:ACCEL_XOUT_H,ACCEL_XOUT_L,ACCEL_YOUT_H,ACCEL_YOUT_L,ACCEL_ZOUT_H,ACCEL_ZOUT_L.

    The registers used to store the value of gyroscope:GYRO_XOUT_H,GYRO_XOUT_L,GYRO_YOUT_H,GYRO_YOUT_L,GYRO_ZOUT_H,GYRO_ZOUT_L.

    2.2.3 Interrupt Mode

    The interrupt status register can read interrupt status in MPU-6050.The main interrupt sources are Free Fall Interrupt,Motion Interrupt,and Zero Motion Interrupt.

    Free Fall Interrupt:The free fall is judged by comparing the magnitude of triaxial acceleration with the threshold value.If the acceleration is close to the threshold value,it triggers the free fall interrupt and change the flag bit.If the acceleration is less than the threshold,it is ignored.

    Motion Interrupt:If the output value is greater than the threshold value,the counter increases by 1.Otherwise,the counter decreases by 1.When the value of counter reaches the threshold value,the motion interrupt will be triggered.When the flag bit is read,both free fall and motion interrupt become 0.

    Zero Motion Interrupt:If the output value is smaller than the threshold value,the counter will increase by 1.When the value of counter reaches the threshold value,the static interrupt is triggered.

    Interrupts can be triggered by the following events:Switching clock source,completing the DMP(Digital Motion Processor),reading data from FIFO(First-In,First-Out)data register and interrupt accelerometer.

    2.2.4 CC2541 Communication

    In addition to the IIC interface connected to MPU6050,CC2541 also has two powerful asynchronous serial communication interfaces that support multiple serial protocols:USART0 and USART1.These two interfaces can not only run asynchronous UART mode but also run synchronous SPI mode.

    The operation characteristics of UART mode are as follows:1)Three calibration methods are applicable,namely odd parity,even parity and no parity.2)Receive interrupts autonomously.3)Access memory autonomously.4)Check the error state through the checksum frame.

    UART has two interface selection methods:2-wire interface(using RXD and TXD),4-wire interface(using RXD,TXD,RTS and CTS).In USART0,P0_2 corresponds to RX,P0_3 to TX,P1_4 to RX,and P1_5 to TX.In USART1,P0_5 corresponds to RX,P0_4 to TX,P1_7 to RX,and P1_6 to TX.

    The UART operations mainly include UxCSR and UxUCR(x is 0 or 1,which is the USART number).Now take USART0 for example.U0CSR.MODE is used to set UART MODE.When U0CSR.RE is 0,the receiver is disabled.When U0CSR.RE is 1,the receiver is enabled.When U0CSR.FE is 0,it means no frame error detection.

    When U0CSR.FE is 1,it means the byte has received incorrect stop bit level.When U0CSR.ERR is 0,there is no parity error.When U0CSR.ERR is 1,there is parity error.U0CSR.RX_BYTE is 0 when no bytes are received.When U0CSR.RX_BYTE is 1,it is ready to receive bytes.When U0CSR.TX_BYTE is 0,the byte is not sent.When U0CSR.TX_BYTE is 1,the byte is written to the cache and the byte will be sent at last.When U0CSR.ACTIVE is 0,USART is idle.When U0CSR.ACTIVE is 1,USART is busy when sending or receiving.The table USART 0 UART Control-U0UCR(0xC4)is in Tab.3.

    Table 3:USART 0 UART control-U0UCR(0xC4)

    3 Hardware Circuit Design

    3.1 Overall Design Scheme and Circuit Diagram Analysis

    In this study,CC2541 is the main chip and MPU6050 is the data collection module which are to be integrated as a smart ring.Fig.3 shows the overall scheme diagram.The smart ring includes the sensor module,power supply,and the Bluetooth transmission unit.APP is developed as a smartphone application to display the data.Fig.4 shows the circuit of CC2541 module.

    Figure 3:Overall scheme diagram

    3.2 Design of CC2541 Reset Circuit

    CC2541 RESET pin can reset CC2541.The reset function will initialize the system to the original state when all registers are cleared.In this design,the reset button is used for future development.Fig.5 is the reset button circuit.

    Under normal conditions,the key will always be in a disconnected state,and the RESET pin will always be in a low-level state.When the key is pressed,RESET is directly connected to the power supply.At this time,RESET is 1 and the starts to reset.Of course,C and R determine the duration of the high level.The advantage of button reset is that it can be reset at any time.But there are also disadvantages.For instance,the operation is more complex and welding is more difficult.

    Figure 4:Circuit of CC2541 bluetooth module

    Figure 5:Reset button circuit

    3.3 Crystal Oscillator Circuit Design of CC2541

    Considering the low power consumption,a primary and a secondary crystal oscillator are added according to the user manual of CC2541.The main crystal oscillator is connected to the circuit of 22 and 23 pins.Fig.6 shows the crystal oscillator circuit.

    Figure 6:Crystal oscillator circuit

    The reason why the RC internal oscillator of 16 MHz is not used is that it has low accuracy which cannot be used in RF transceivers.The main frequency of auxiliary crystal oscillator is 32.768 KHz according to the manual,which is added to achieve low power consumption.In general,the crystal oscillator is matched with capacitors.The main crystal oscillator and the auxiliary crystal oscillator have two capacitors respectively according to the manual,which can effectively assist the crystal to vibrate and improve the stability of the circuit.

    3.4 CC2541 RF Circuit and Antenna Design

    The 25 and 26 pins of CC2541 are RF pins,and a 50 ohms impedance is the best match for RF communication according to the manual.Therefore,a device named 2450BM15A0002 of balun is used in this design.Considering the ceramic antenna is smaller than traditional PCB antenna,ceramic antenna is selected.Fig.7 is the RF circuit and antenna.

    Figure 7:RF circuit and antenna

    CC2541 and MPU6050 have SCL pin,SDA pin and IIC interface,which can receive and send data.NFC is used for communication between CC2541 and MPU6050.MPU6050 transmits the collected data to CC2541 using the two pins,and then the data will be transmitted to the mobile terminal using the Bluetooth module.

    3.5 Power Circuit Design

    Since 5 V is the reference voltage,the voltage stabilizing circuit Lm1117 is used.The advantage of LM1117 is the high voltage stability,accuracy,small temperature drift at high and low temperature and simple peripheral circuit.The output current can reach 0.8 A with 0.2% voltage stability and 0.4% load stability.A light emitting diode is added to the far right to indicate whether it is working to prevent the battery from connecting backwards.Fig.8 shows the power supply circuit.

    Figure 8:Power supply circuit

    3.6 Hardware Overall Analysis

    This design uses a CC2541 chip,which has low power consumption.The current is only 0.5 mA in external interrupt mode.In order to meet the requirements of low power consumption,32.768 kHz auxiliary crystal oscillator is used.

    4 Tremor Detection Algorithm Design

    4.1 Main Flow of System

    After calling the main function,the system will enter the initial state,and then continuously send the current state of the device in the while loop.The program flow diagram is shown in Fig.9.

    4.2 Internal Drive Function Flow of MPU6050

    After calling the MPU6050 drive function,it then tests readability of the data collected by MPU6050 and the stability of connection.It then sends all data of the triaxial acceleration and angular velocity to CC2541 by NFC.The flow diagram is shown in Fig.10.

    4.3 Inertial Data Collection

    4.3.1 Initialize and Test the Connection Fragment

    The algorithm for collecting inertial data is shown in Fig.11.

    At first,IIC is initialized and bIfConnected is set with a flag ahead of time.If AD0 of MPU6050 is in low level,then the connection is successful.And then assign a value of 1 to bIfConnected and initialize MPU6050.

    Using a for loop,all the data in the IIC device are read and stored in the array for invocation.Where Slave Address is the device address,REG_Address is the device’s internal register address,readNum is the maximum number of internal registers,and I2CSlaveBuffer[PORT_USED][i]is the array to be sent.

    Figure 9:System flow design

    The variables ax,ay,az receive the unprocessed triaxial acceleration data from registers and gx,gy,gz receive the unprocessed triaxial angular velocity data from the correspondent registers.These data will be continuously passed to the DMP motion chip for further processing and then being transmitted to the mobile client by the Bluetooth module.

    The data packet of the communication protocol usually has 20 bytes,including header(1 byte),effective data length(1 byte),function code(1 byte),effective data(0–16 bytes),checksum(1 byte,the sum of data before the checksum in the package),and complement byte(0–16 bytes).Firstly,check the following packet head to see whether it is 0×55.If so,continue the following operation and analyze whether the packet head of the transmitted packet is 0×52 or 0×51 respectively.If it is 0×51,it indicates that the packet is an acceleration packet.Otherwise,it is an angular velocity packet.

    Figure 10:Internal flow diagram

    Figure 11:Algorithm for collecting inertial data

    5 System Experiments and Results

    5.1 APP Design Test

    An smartphone application is designed as shown in Fig.12.In this interface,x,y and z in the diagram are displayed for acceleration and angular velocity.In the test process,click the“Select Bluetooth”button to select the corresponding Bluetooth.At this time,each axial waveform will be displayed on the screen,and the data will be saved into two files in the system with file extension “.txt”,with “ble6050a.txt” storing acceleration data,and “ble6050b.txt” storing angular velocity data.

    Figure 12:APP page

    5.2 Prototype of the Ring Device

    The smart ring implemented in this study uses a plastic platform,which does not produce much interference compared with using metal material.The Printed Circuit Board(PCB)contains a switch,MPU6050,CC2541.When the switch is open,a LED lamp will light up to indicate that the device is running.Fig.13 shows the prototype and the PCB diagram.

    5.3 Experiments and Results Analysis

    5.3.1 Measuring the Intensity of Finger Tremor

    In our experiment,hand tremor is simulated by healthy people using protocols based on the study as described in the research background.Three tests are conducted in this study.All participants wear the ring device in their fingers(see Fig.14).

    Sample 1:The tested person stays static.

    Sample 2:The tested person walks normally.

    Sample 3:The tested person simulates the tremor.

    During the test,we transmitted the data to the mobile APP according to the operation instructions and selected sample screenshots as shown in the figures.Fig.15 shows the waveform of three axis acceleration and angular velocity of Sample 1.Fig.16 shows the waveform of Sample 2.Fig.17 shows the waveform of Sample 3.

    Figure 13:Prototype and PCB diagram

    Figure 14:Wearing device on one finger

    Tabs.4–9 present the raw data of the three samples of the experiments,and Tabs.10–12 present the processed data,where acceleration of three-axis is processed by numerical value/32768×16 g(where g is the acceleration of gravity and taken as 9.8 m/s2),and angular velocity of three-axis is processed by numerical value/32768×2000(?/s).

    Tab.4 shows the acceleration raw data of Sample 1.

    Tab.5 shows the angular velocity raw data of Sample 1.

    Tab.6 shows the acceleration raw data of Sample 2.

    Tab.7 shows the velocity raw data of Sample 2.

    Tab.8 shows the acceleration raw data of Sample 3.

    Tab.9 shows the angular velocity raw data of Sample 3.

    Figure 15:Waveform of three axis acceleration and angular velocity of Sample 1

    Figure 16:Waveform of three axis acceleration and angular velocity of Sample 2

    Figure 17:Waveform of three axis acceleration and angular velocity of Sample 3

    Table 4:Acceleration raw data of Sample 1

    Table 5:Angular velocity raw data of Sample 1

    Table 6:Acceleration raw data of Sample 2

    Table 7:Angular velocity raw data of Sample 2

    Table 8:Acceleration raw data of Sample 3

    Tab.10 shows strength table after Sample 1 treatment.

    Tab.11 shows strength table after Sample 2 treatment.

    Tab.12 shows strength table after Sample 3 treatment.

    Table 9:Angular velocity raw data of Sample 3

    Table 10:Strength table after Sample 1 treatment

    Table 11:Strength table after Sample 2 treatment

    Table 12:Strength table after Sample 3 treatment

    Figs.18 and 19 show the intensity comparison of three samples for movements of static state,walking state,and simulated tremor by “counting money.” Examining the waveform of figure and the line graph of the processed intensity values of the acceleration and angular velocity,it can be found that acceleration of Sample 1 is stable at 275 and angular velocity is stable at 2000.The acceleration of Sample 3 fluctuates and so does the angular velocity.

    Figure 18:Acceleration intensity comparison diagram

    Figure 19:Angular velocity intensity comparison diagram

    5.3.2 Determine the Degree of Bradykinesia in Parkinson’s Patients

    To assess the severity of PD,the Unified Parkinson’s Disease Rating Scale(UPDRS)is most commonly used.In this experiment,measuring the movement time between the touching of thumb and the index finger is the most important part in order to analyze the symptoms of bradykinesia between Parkinson’s patients and normal people.In the experiment,three participants conducted as healthy people(Participant 1,Participant 2 and Participant 3)and three participants(Participant 4,Participant 5 and Participant 6)simulated PD of UPDRS level 1–3.All these participantss touched the index finger with thumb,and the movement time was measured respectively.

    The movement timeTbetween the touching of thumb and the index finger can be calculated using the collected datawhich represent the mutually perpendicular accelerated speed.Eq.(1)is the relationship between resultant acceleration and the three components.

    The movement timeTcan be measured by evaluating the periodPof angleθorφusing the Discrete-Time Fourier Transform Eq.(3).The sequence with N samples of angleθcan be denoted asθ[N]when only considering angleθ.

    where,ωrepresents the angular velocity ofθ.We then calculate theωwhich makes the |Θ(ω)|to achieve the maximum value.Then,T=P=1/ωwill be the movement time between the touching of thumb and the index finger.

    Fig.20 is the maximum finger opening.

    Fig.21 is the moderate opening of the finger.

    Fig.22 is the two-finger touch diagram.

    Tab.13 shows the movement time of each sample and the doctor’s ratings according to UPDRS.

    Figure 20:Maximum finger opening

    Figure 21:Moderate opening of the finger

    Figure 22:Two-fingers touch

    Table 13:Time recording table

    It can be seen that the average finger movement time is 307 ms,and the standard deviation is 55 ms.The higher the level of bradykinesia is,the more serious the disease is.According to the experimental results,the severity of symptomSis approximately linear with the finger touch timeT.The following Eq.(4)is a simple formula of the linear relationship between S and T:

    wherea,bare the coefficients which can be determined by experimental results.

    5.4 Error Analysis

    As can be seen in the experiments,the resting state waveform of Sample 1 is not very smooth which may be caused by the following factors:

    (1)The fingers cannot stay completely still.

    (2)The interference of surrounding environment exists.

    6 Conclusion

    In this paper,we presented the design and implementation of a smart wearable device(smart ring)for detecting hand tremor for Parkinson’s disease patients.We used MPU6050(MIDI Processing Unit)consisting of a 3-axis gyroscope and a 3-axis accelerometer to collect acceleration and angular velocity of fingers.By analyzing the time of specific finger movements,we successfully recognized the tremor signals with high accuracy.Meanwhile,with Bluetooth 4.0(Bluetooth Low Energy,BLE)and networking terminal ability,tremor data are transferred to a monitoring smartphone application in real time.The experimental results have shown that the proposed device is convenient for long-term tremor detection which is vital for early detection and treatment for Parkinson’s disease.As future work,we will continue our study to address the following issues we found with the product:(1)The selected 5 V button battery caused a noticeable amount of power consumption in the process of transmitting data,making the standby time of the smart ring shorter;(2)Although the size of the smart ring is small,the physical switch is big which should be improved;(3)Open source Android APP is selected so that the conclusion can be obtained only by refreshing the interface and then further analysis of data can be conducted to best utilize this smart device.

    Funding Statement:This work is supported by the National Natural Science Foundation of China(Grant Nos.61972207 and 61802196),Jiangsu Provincial Government Scholarship for Studying Abroad and the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)fund.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    欧美精品国产亚洲| 直男gayav资源| 亚洲欧洲日产国产| av免费观看日本| 3wmmmm亚洲av在线观看| 国产私拍福利视频在线观看| 国产高潮美女av| 欧美精品一区二区大全| 国产高潮美女av| 色5月婷婷丁香| 18+在线观看网站| 色尼玛亚洲综合影院| 国产淫语在线视频| 国产v大片淫在线免费观看| 欧美3d第一页| 亚洲国产最新在线播放| 天天躁日日操中文字幕| 国产男人的电影天堂91| 在线观看美女被高潮喷水网站| 国产三级在线视频| 蜜臀久久99精品久久宅男| 日产精品乱码卡一卡2卡三| 精品久久久久久久久亚洲| 精品人妻视频免费看| 国产高潮美女av| 国产欧美日韩精品一区二区| 可以在线观看毛片的网站| 在线免费观看不下载黄p国产| 国产成人a∨麻豆精品| 一区二区三区四区激情视频| 男女啪啪激烈高潮av片| 国产探花极品一区二区| 亚洲性久久影院| av福利片在线观看| 精品国产露脸久久av麻豆 | 18禁在线播放成人免费| 我要搜黄色片| 免费观看人在逋| 国产精品三级大全| 国产精品日韩av在线免费观看| 少妇的逼好多水| 一级毛片电影观看 | 能在线免费观看的黄片| 91精品国产九色| 一级毛片我不卡| 午夜福利在线观看吧| 日本wwww免费看| 精品久久久久久久久亚洲| 人妻制服诱惑在线中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 色综合色国产| 午夜精品国产一区二区电影 | 国产老妇女一区| 欧美区成人在线视频| 熟女电影av网| 熟女人妻精品中文字幕| 在线免费观看不下载黄p国产| 中文字幕制服av| 国产欧美另类精品又又久久亚洲欧美| 中文天堂在线官网| 亚洲第一区二区三区不卡| 亚洲精品乱码久久久久久按摩| 秋霞伦理黄片| 岛国毛片在线播放| 69人妻影院| 亚洲国产色片| av免费在线看不卡| 99久久中文字幕三级久久日本| 国产亚洲5aaaaa淫片| 91久久精品电影网| 人妻夜夜爽99麻豆av| 青春草视频在线免费观看| 久久久久国产网址| 老女人水多毛片| 国产 一区 欧美 日韩| 成人二区视频| 欧美成人免费av一区二区三区| 中文在线观看免费www的网站| a级毛片免费高清观看在线播放| 国产亚洲91精品色在线| 免费观看性生交大片5| 麻豆成人午夜福利视频| 国国产精品蜜臀av免费| 亚洲国产最新在线播放| 男人狂女人下面高潮的视频| 一区二区三区免费毛片| 在线免费观看不下载黄p国产| 日韩欧美精品免费久久| 亚洲最大成人手机在线| 精品久久久久久久久亚洲| 91久久精品国产一区二区三区| 国模一区二区三区四区视频| 麻豆久久精品国产亚洲av| 午夜福利在线观看免费完整高清在| 欧美三级亚洲精品| 丝袜喷水一区| 丝袜喷水一区| 尾随美女入室| 久久精品国产亚洲av涩爱| 国产熟女欧美一区二区| 成人毛片60女人毛片免费| 国产精品av视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 中文欧美无线码| 国产精品三级大全| 亚洲最大成人手机在线| 男女视频在线观看网站免费| 三级国产精品片| 亚洲欧美清纯卡通| a级毛色黄片| 青春草视频在线免费观看| 欧美丝袜亚洲另类| 国产精品蜜桃在线观看| 男女那种视频在线观看| 男人狂女人下面高潮的视频| 高清在线视频一区二区三区 | 国产精品福利在线免费观看| 午夜福利在线在线| av免费观看日本| 黑人高潮一二区| 少妇人妻一区二区三区视频| 久久婷婷人人爽人人干人人爱| 天堂中文最新版在线下载 | 中文精品一卡2卡3卡4更新| 日本wwww免费看| 欧美性猛交╳xxx乱大交人| 国产精品三级大全| 免费电影在线观看免费观看| 国产成人a∨麻豆精品| 最近2019中文字幕mv第一页| 两个人视频免费观看高清| 精品人妻熟女av久视频| 国产v大片淫在线免费观看| 99久久成人亚洲精品观看| 尤物成人国产欧美一区二区三区| 伦精品一区二区三区| 啦啦啦观看免费观看视频高清| 免费一级毛片在线播放高清视频| 国产精品国产高清国产av| 精品酒店卫生间| 日日啪夜夜撸| 色吧在线观看| 99久久无色码亚洲精品果冻| 成人鲁丝片一二三区免费| 亚洲精品乱码久久久v下载方式| 婷婷色av中文字幕| 波多野结衣巨乳人妻| 精品少妇黑人巨大在线播放 | 亚洲精品影视一区二区三区av| 神马国产精品三级电影在线观看| 国产伦精品一区二区三区视频9| 国产男人的电影天堂91| 美女cb高潮喷水在线观看| av国产久精品久网站免费入址| 久久人人爽人人片av| 国产色爽女视频免费观看| 韩国av在线不卡| 大又大粗又爽又黄少妇毛片口| 99久久精品国产国产毛片| 久久精品国产99精品国产亚洲性色| 成人综合一区亚洲| 国产老妇女一区| 插逼视频在线观看| 久久久亚洲精品成人影院| 亚洲av中文av极速乱| 亚洲经典国产精华液单| 日本五十路高清| 色综合亚洲欧美另类图片| 欧美精品国产亚洲| 精品一区二区三区人妻视频| 日本黄大片高清| 少妇丰满av| 床上黄色一级片| 男女边吃奶边做爰视频| 2021少妇久久久久久久久久久| 成年女人永久免费观看视频| 国产精品久久久久久精品电影小说 | 欧美丝袜亚洲另类| 精品久久久久久久人妻蜜臀av| 久久精品国产鲁丝片午夜精品| 两性午夜刺激爽爽歪歪视频在线观看| 插逼视频在线观看| 免费av毛片视频| 国产真实伦视频高清在线观看| 色综合亚洲欧美另类图片| 91av网一区二区| 99九九线精品视频在线观看视频| 男女视频在线观看网站免费| 麻豆成人午夜福利视频| 99热全是精品| 国产免费又黄又爽又色| 中国美白少妇内射xxxbb| 亚洲真实伦在线观看| 人体艺术视频欧美日本| 黄片无遮挡物在线观看| 97在线视频观看| 欧美一区二区亚洲| 日本免费a在线| 欧美成人a在线观看| 六月丁香七月| 欧美性猛交黑人性爽| 韩国av在线不卡| 一级av片app| 日本三级黄在线观看| 亚洲美女搞黄在线观看| 欧美激情国产日韩精品一区| 国产免费福利视频在线观看| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 禁无遮挡网站| 国产一级毛片在线| 色5月婷婷丁香| 午夜福利在线在线| 中文字幕熟女人妻在线| 久久人人爽人人片av| 亚洲人成网站在线观看播放| 国产精品一区二区性色av| 91久久精品电影网| 九九热线精品视视频播放| 亚洲五月天丁香| 国内精品宾馆在线| 精品人妻一区二区三区麻豆| 可以在线观看毛片的网站| 欧美日韩综合久久久久久| 国产老妇女一区| 99久久中文字幕三级久久日本| 亚洲人与动物交配视频| 欧美97在线视频| 永久网站在线| 婷婷六月久久综合丁香| 亚洲精品影视一区二区三区av| 99久久人妻综合| 午夜精品一区二区三区免费看| 国内精品美女久久久久久| 国内揄拍国产精品人妻在线| 亚洲精品日韩av片在线观看| 久久99热这里只频精品6学生 | 色尼玛亚洲综合影院| 免费av不卡在线播放| 高清午夜精品一区二区三区| 七月丁香在线播放| 国产高清国产精品国产三级 | 三级毛片av免费| 国产成人精品一,二区| 久久久久久久久久成人| 欧美精品一区二区大全| 国产极品精品免费视频能看的| 国产精品一及| 色哟哟·www| 亚洲人与动物交配视频| 老女人水多毛片| 日韩欧美在线乱码| 日产精品乱码卡一卡2卡三| 国产成人福利小说| 亚洲精品国产av成人精品| 国产又黄又爽又无遮挡在线| 免费无遮挡裸体视频| 精品久久久久久电影网 | 啦啦啦观看免费观看视频高清| 人体艺术视频欧美日本| 国产高清视频在线观看网站| 色播亚洲综合网| 国产精品伦人一区二区| 国产精品一区二区三区四区免费观看| 亚洲精品亚洲一区二区| 少妇裸体淫交视频免费看高清| 国产一级毛片在线| 国产激情偷乱视频一区二区| 精品一区二区三区人妻视频| 久久国产乱子免费精品| 男女视频在线观看网站免费| 国产单亲对白刺激| 亚洲欧美精品综合久久99| 小蜜桃在线观看免费完整版高清| 1000部很黄的大片| 亚洲av成人av| 非洲黑人性xxxx精品又粗又长| 久久精品熟女亚洲av麻豆精品 | 三级毛片av免费| 国产伦一二天堂av在线观看| 2021天堂中文幕一二区在线观| 国产精品爽爽va在线观看网站| 欧美xxxx性猛交bbbb| 欧美3d第一页| 一个人看的www免费观看视频| 婷婷六月久久综合丁香| 一级毛片电影观看 | 中国美白少妇内射xxxbb| 三级男女做爰猛烈吃奶摸视频| 日韩三级伦理在线观看| 国产成人91sexporn| 在现免费观看毛片| 午夜福利视频1000在线观看| 久久久欧美国产精品| 身体一侧抽搐| 亚洲av免费在线观看| 毛片女人毛片| 国产一区有黄有色的免费视频 | 直男gayav资源| 国产精品乱码一区二三区的特点| 免费av观看视频| 成人三级黄色视频| 亚洲精品国产成人久久av| 乱系列少妇在线播放| 亚洲在久久综合| av播播在线观看一区| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 色视频www国产| 欧美不卡视频在线免费观看| 免费观看a级毛片全部| 99久久精品热视频| 亚洲欧美日韩高清专用| 中文资源天堂在线| 国产免费福利视频在线观看| 欧美性猛交╳xxx乱大交人| 国产黄色视频一区二区在线观看 | 在线观看一区二区三区| 亚洲国产欧洲综合997久久,| 久99久视频精品免费| 国产精品不卡视频一区二区| 日韩成人伦理影院| 亚洲怡红院男人天堂| 亚洲五月天丁香| 日日摸夜夜添夜夜爱| 一卡2卡三卡四卡精品乱码亚洲| 蜜桃久久精品国产亚洲av| 人妻夜夜爽99麻豆av| 听说在线观看完整版免费高清| av卡一久久| 秋霞伦理黄片| eeuss影院久久| 青青草视频在线视频观看| 欧美3d第一页| 欧美三级亚洲精品| 人人妻人人澡欧美一区二区| 免费搜索国产男女视频| 国产69精品久久久久777片| 日本三级黄在线观看| 永久网站在线| 国产精品一及| 在线免费十八禁| 国产精品野战在线观看| 婷婷色综合大香蕉| 两个人视频免费观看高清| 级片在线观看| 中文亚洲av片在线观看爽| or卡值多少钱| 91在线精品国自产拍蜜月| 中文字幕久久专区| 三级经典国产精品| 欧美一区二区国产精品久久精品| 极品教师在线视频| 日本黄色视频三级网站网址| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看| av专区在线播放| 欧美一区二区精品小视频在线| 69av精品久久久久久| 99热6这里只有精品| 老司机影院毛片| 国国产精品蜜臀av免费| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 国产精品日韩av在线免费观看| 成年av动漫网址| 国产精品一区二区性色av| 久久99热这里只有精品18| 最后的刺客免费高清国语| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲午夜精品一区二区久久 | 三级男女做爰猛烈吃奶摸视频| 国产亚洲最大av| 又黄又爽又刺激的免费视频.| 人人妻人人澡人人爽人人夜夜 | 3wmmmm亚洲av在线观看| 日韩制服骚丝袜av| 久久精品国产亚洲网站| 亚洲成色77777| 精华霜和精华液先用哪个| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| 中文字幕久久专区| 成人国产麻豆网| 两个人视频免费观看高清| 美女被艹到高潮喷水动态| 天天躁夜夜躁狠狠久久av| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| av免费在线看不卡| 性色avwww在线观看| 欧美日本视频| 国产高清视频在线观看网站| 成人鲁丝片一二三区免费| 亚洲成人av在线免费| 欧美性感艳星| 欧美激情在线99| 看黄色毛片网站| 成人国产麻豆网| 成年免费大片在线观看| 熟女电影av网| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 欧美色视频一区免费| 内射极品少妇av片p| 99在线人妻在线中文字幕| av女优亚洲男人天堂| 小说图片视频综合网站| 在线观看美女被高潮喷水网站| 能在线免费观看的黄片| 亚洲经典国产精华液单| 国产精品国产三级国产av玫瑰| 18禁动态无遮挡网站| 国产成人freesex在线| 久久久久久大精品| 欧美成人午夜免费资源| 91aial.com中文字幕在线观看| 国产高清国产精品国产三级 | 伊人久久精品亚洲午夜| 亚洲欧美精品自产自拍| 色播亚洲综合网| 亚洲无线观看免费| 成人毛片a级毛片在线播放| 久久久精品欧美日韩精品| 99热这里只有是精品50| 日本熟妇午夜| or卡值多少钱| 亚洲不卡免费看| 亚洲自拍偷在线| 日日撸夜夜添| 亚洲精品乱码久久久v下载方式| 国产亚洲5aaaaa淫片| 91久久精品国产一区二区三区| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 亚洲精品影视一区二区三区av| 成年版毛片免费区| 男女边吃奶边做爰视频| 欧美3d第一页| 亚洲av二区三区四区| 欧美激情在线99| 久久久精品大字幕| 99久久精品热视频| 51国产日韩欧美| 亚洲第一区二区三区不卡| 久久午夜福利片| 欧美潮喷喷水| 欧美精品国产亚洲| 国产色爽女视频免费观看| 日韩人妻高清精品专区| 亚洲综合精品二区| 一级黄色大片毛片| 69av精品久久久久久| 麻豆一二三区av精品| 亚洲熟妇中文字幕五十中出| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 一级毛片久久久久久久久女| 精品久久久久久久人妻蜜臀av| 亚洲av福利一区| 内地一区二区视频在线| 天堂中文最新版在线下载 | 在线免费十八禁| 午夜福利在线观看免费完整高清在| videossex国产| 伦理电影大哥的女人| 国产av不卡久久| 国产 一区 欧美 日韩| 国产高潮美女av| 日本av手机在线免费观看| 内地一区二区视频在线| 我的老师免费观看完整版| 午夜老司机福利剧场| 午夜久久久久精精品| 久久亚洲精品不卡| 少妇的逼好多水| 中文字幕久久专区| 日韩三级伦理在线观看| 亚洲不卡免费看| 搡女人真爽免费视频火全软件| 最近的中文字幕免费完整| 国产熟女欧美一区二区| 天美传媒精品一区二区| 亚洲欧美成人精品一区二区| 男女国产视频网站| 欧美一级a爱片免费观看看| 欧美区成人在线视频| 亚洲国产精品成人综合色| 夜夜爽夜夜爽视频| 我的女老师完整版在线观看| 亚洲精品456在线播放app| 欧美一级a爱片免费观看看| 免费不卡的大黄色大毛片视频在线观看 | 插逼视频在线观看| 一个人观看的视频www高清免费观看| 国产精品一区二区三区四区免费观看| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 国产精品一区二区性色av| 国产一区二区在线av高清观看| 别揉我奶头 嗯啊视频| 久久久久精品久久久久真实原创| 久久国内精品自在自线图片| 亚洲成人av在线免费| 青青草视频在线视频观看| 又爽又黄a免费视频| 国产成人91sexporn| 一级av片app| 国产 一区精品| 国产乱来视频区| 国产乱人视频| 午夜激情欧美在线| 国产色婷婷99| 国产国拍精品亚洲av在线观看| 国产免费又黄又爽又色| 一个人免费在线观看电影| 中文资源天堂在线| 国产一区亚洲一区在线观看| 18禁动态无遮挡网站| 日本色播在线视频| 如何舔出高潮| 少妇丰满av| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 51国产日韩欧美| 亚洲激情五月婷婷啪啪| 国产乱来视频区| 男的添女的下面高潮视频| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕| 热99re8久久精品国产| 国产极品精品免费视频能看的| 91精品伊人久久大香线蕉| 国语自产精品视频在线第100页| 我的女老师完整版在线观看| 成人美女网站在线观看视频| 国产精品,欧美在线| 久久草成人影院| 国产一区有黄有色的免费视频 | 亚洲久久久久久中文字幕| 亚洲欧美成人综合另类久久久 | 少妇丰满av| 欧美高清性xxxxhd video| 国产乱人视频| 色视频www国产| 日韩人妻高清精品专区| 大香蕉97超碰在线| 国产精品综合久久久久久久免费| 一区二区三区四区激情视频| 97超碰精品成人国产| 色综合站精品国产| 国产午夜精品论理片| 久久久成人免费电影| 最近的中文字幕免费完整| 人人妻人人澡欧美一区二区| 欧美成人免费av一区二区三区| 老司机影院毛片| 1024手机看黄色片| 久久久亚洲精品成人影院| 能在线免费观看的黄片| 欧美三级亚洲精品| 日日干狠狠操夜夜爽| 三级男女做爰猛烈吃奶摸视频| 午夜精品在线福利| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久av不卡| 伊人久久精品亚洲午夜| 好男人视频免费观看在线| 国产毛片a区久久久久| 嫩草影院新地址| 欧美日韩一区二区视频在线观看视频在线 | 日韩亚洲欧美综合| 亚洲最大成人中文| 狠狠狠狠99中文字幕| 波野结衣二区三区在线| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 美女黄网站色视频| 精品久久久久久久人妻蜜臀av| 一夜夜www| 欧美97在线视频| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 天堂√8在线中文| 国产探花在线观看一区二区| 一级二级三级毛片免费看| 亚洲经典国产精华液单| 日韩亚洲欧美综合| 国产成人午夜福利电影在线观看| 一夜夜www| 高清午夜精品一区二区三区| 国产精品熟女久久久久浪| 成年版毛片免费区| 1024手机看黄色片| 深夜a级毛片| 极品教师在线视频| 日韩三级伦理在线观看| 成年免费大片在线观看| www.色视频.com| 欧美潮喷喷水| 高清av免费在线| 成人午夜精彩视频在线观看| 一区二区三区四区激情视频| 黄色配什么色好看| 国产亚洲精品久久久com| av在线观看视频网站免费| 在线播放国产精品三级| 国产精品av视频在线免费观看| 日本一本二区三区精品| 欧美性感艳星|