CN108703748A - Health state evaluating method based on pulse wave harmonic wave - Google Patents

Health state evaluating method based on pulse wave harmonic wave Download PDF

Info

Publication number
CN108703748A
CN108703748A CN201810628488.0A CN201810628488A CN108703748A CN 108703748 A CN108703748 A CN 108703748A CN 201810628488 A CN201810628488 A CN 201810628488A CN 108703748 A CN108703748 A CN 108703748A
Authority
CN
China
Prior art keywords
pulse wave
wave
harmonic
probability density
harmonic wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810628488.0A
Other languages
Chinese (zh)
Other versions
CN108703748B (en
Inventor
不公告发明人
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Qianmai Health Technology Co ltd
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201810628488.0A priority Critical patent/CN108703748B/en
Publication of CN108703748A publication Critical patent/CN108703748A/en
Application granted granted Critical
Publication of CN108703748B publication Critical patent/CN108703748B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4854Diagnosis based on concepts of traditional oriental medicine
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Cardiology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The present invention measures human pulse wave using photoelectricity volume pulses measure technology and carries out simple venation wavelength-division of fighting and cuts, characteristics extraction and sequence are carried out using statistical method respectively to each harmonic wave of continuous single pulse wave, phase space is constructed on this basis and calculates probability density function, the health status of measured's viscera function can be assessed according to the logarithm normal distribution state of each harmonic wave probability density.

Description

Health state evaluating method based on pulse wave harmonic wave
Technical field
The present invention relates to a kind of health state evaluating methods, specially establish pulse wave harmonic wave and collaterals of human and dirty A kind of mathematical method of internal organs health status incidence relation.
Background technology
Due to life style, operating pressure etc., current social inferior health incidence is high, however to inferior health There are great difficulties for the actual inspection and treatment of state.Since its symptom is more complicated, assessment mode is more subjective and is difficult to It is unified.It can not check that any organic disease, function sexually revise or metabolic disorder on the person in sub-health state, therefore base It can not effectively be made a definite diagnosis in general medical means and method, thus can not treat or treat and be lack of pertinence.
In terms of human circulation angle, pulse wave is caused by the periodical bounce of heart, is passed through as a kind of quasi-periodicity signal Blood vessel is transmitted to the whole body.So being influenced by a variety of physiology or pathological factor in its transmittance process, it may appear that speed, strong The difference of degree, i.e., so-called strong and weak rapid point of pulse condition.By in pulse wave imply physiology, pathological information parses can Efficient diagnosis is carried out to individual psychosomatic state.
From the perspective of Traditional Chinese Medicine diagnosis by feeling the pulse, the finger that can reflect health state is effectively extracted from pulse wave Mark establishes pulse wave harmonic wave and will be helpful to push to human-body sub-health shape with channels and collaterals and the associated mathematical method of internal organs health status The diagnosis of state, the evaluation and test of viscera health status realization.
Invention content
For the complexity for overcoming human body sub-health status to detect, present invention combination Chinese medicine collaterals of human is theoretical, designs number Method establishes the incidence relation between pulse wave harmonic wave and channels and collaterals and internal organs health status, to reach the mesh of health evaluation and test 's.
The present invention obtains the original pulse signal of human body privileged site using photoelectricity volume pulse transducer.Design bandpass filtering Device eliminates various interference and the noises such as baseline, and single pulse wave identification is carried out using difference clustering method.Each pulse wave is carried out Fourier transform obtains the harmonic component of waveform, takes preceding 12 component of degree n n statistics each harmonic time delays, while analyzing frequency multiplication spy Sign.Phase space reconfiguration is carried out to the harmonic characteristic value of continuous pulse wave, probability density function is calculated and carries out normal distribution fitting. According to blood vessel resonance theory, it is believed that pulse wave each harmonic corresponds respectively to the acupuncture point of 12 channels and collaterals in theory of traditional Chinese medical science and corresponding dirty Internal organs organ, the organ corresponding to these harmonic waves is followed successively by liver,kidney,spleen, lung, stomach, courage etc., so can be according to the normal distribution of harmonic wave State carries out the health condition analysis of each organ of human body.
The technical scheme is that:
A kind of health state evaluating method based on pulse wave harmonic wave, its step are as follows:
1)Pulse wave at human body privileged site is obtained using photoplethysmographic detection method.
2)Due to there are the interference of the high-frequency noises such as power frequency, needing to be removed in signal using low-pass filter.Equally It needs to eliminate the interference because of low-frequency noises such as the baseline drifts that human motion generates.Amid all these factors, we select band logical to filter The completion of wave device is filtered.
3)Difference is asked to the consecutive points of pulse wave, clustering is carried out to gained difference:To removing mistake because of calculus of differences The respective point of difference takes amplitude, carries out cluster again to amplitude and peak point can be obtained.Single pulse wave can be carried out according to peak value point value Separation.
4)Fourier transform is carried out to single pulse wave after separation and obtains frequency signal, is extracted in pulse frequency spectrum first 12 times The logarithm of harmonic amplitude is characterized value, by each harmonic characteristic value makeup time array of continuous pulse.
5) above-mentioned time series is subjected to phase space reconfiguration and seeks the probability density function corresponding to each harmonic.
6)Logarithm normal distribution fitting is carried out to the probability density of pulse wave harmonic wave.Since pulse wave each harmonic is right respectively The acupuncture point of 12 channels and collaterals in theory of traditional Chinese medical science and corresponding internal organs organ are answered, they are followed successively by liver,kidney,spleen, lung, stomach, courage etc..Institute The health state evaluation of each organ of human body can be carried out according to the normal distribution state of harmonic wave probability density.
The beneficial effects of the present invention are:
The present invention eliminates an anisotropic error factor with statistical method, with Mathematical Method pulse wave each harmonic with A kind of incidence relation is established between the evaluation and test of human body corresponding organ health status.
Experiment shows that the present invention is feasible in terms of the inferior health of human body is evaluated and tested with each organ lesion prediction.The present invention Effective theoretical and experimental basis is provided for subsequent analysis and implementation system.
Description of the drawings
Fig. 1 be the present invention experiment in from human finger acquire pulse wave signal;
Fig. 2 is pulse wave schematic diagram.Wherein, mark 1 is the main wave in oscillogram, and mark 2 is tidal wave or dicrotic pulse prewave, mark 3 For dicrotic wave, 4 dicrotic waves are marked;
Fig. 3 is the result that the experimental data of the present invention is fitted through normal distyribution function.
Specific implementation mode
The present invention uses photoplethysmographic detection method, and specific implementation mode is:It is placed in finger or specific pulse wave The light source of acquisition position side sends out the light of specific wavelength, due to the light intensity of the non-beating tissue resorption such as skin, muscle, soft tissue It does not change over, and the volumetric blood in artery is fluctuated with the diastole of heart and contraction cycle so that its light intensity absorbed Variation synchronous with heart arteriopalmus, is also at this time a fluctuation signal by the luminous intensity that photoreceiver detects, as shown in Figure 1. Since the movement of human body can cause motion artifacts, so the fluctuation signal detected can have baseline drift, while with power frequency And High-frequency Interference, need design filter to be removed.Since the energy of pulse wave is mainly distributed on 0.5Hz~4.0Hz, herein The realization of 9 rank IIR Butterworth filters is selected to be filtered.It is flatter in passband amplitude frequency curve according to Butterworth filter The characteristics of falling to zero faster with suppressed frequency band, it is respectively 0.5~4hz and 0.2~6hz to choose passband and stop-band frequency, 3db and 40db is passband maximum attenuation and minimum attenuation in stop band, can filter out interference and the pulse wave of high-frequency noise Frequency range.
Since the wave character of the pulse wave signal time domain between different people and under same person different pathological situation has Prodigious difference, these differences are presented as that time domain pulse wave will present the different rhythm and pace of moving things, but the peak point difference in magnitude of same individual It is different smaller.So the method for using difference cluster herein, i.e., seek difference to the consecutive points of pulse wave, in combination with clustering, Extreme point is judged to obtain peak point.Based on difference clustering method the specific steps are:
1)Calculate signal forward difference calibration value.When difference is more than or equal to zero(The ascending branch of pulse)When, it defines constant;When difference is small In zero(The descending branch of pulse)When, it is defined as zero.
2)It is the extreme point being likely to occur to record at calibration value change.
3)Calculate the difference of extreme point.If the difference of current point and the difference of latter point are not belonging to same class, can determine whether This point is possible extreme value points, error dot of the preliminary removal due to signal burr.
4)Clustering is carried out to the amplitude of induction signal to possible extreme value points, data are divided into two classes with k-means methods, Middle peak point is the classification of negligible amounts, is thus removed due to signal dicrotic wave(See Fig. 2)Peak point caused by error dot.
5)All qualified peak points are stored in array.
Efficiently separating for single pulse wave can be realized according to the above-mentioned minimum sought and maximum point.
Since the wave crest of the single pulse wave signal frequency spectrum in multiple tools very close period will appear certain deviation each other.With The increase of pulse wave signal number of cycles, each secondary non-zero harmonic wave can be correspondingly formed specific probability distribution, so at this to list Pulse wave carries out Fourier transform, for statistical analysis to continuous simple venation wave frequency spectrum signature of fighting, every in wave frequency spectrum of being fought with simple venation Characteristic value of the amplitude as pulse wave spectrum signal at a Frequency point carries out phase space reconfiguration, by reconstructing to these characteristic values Sequence can acquire probability density function using correlation integral.Probability density function calculating process can be analyzed to:
1)Input the matrix reconstructed;
2)Distance is asked to all the points in phase space;
3)Heaviside functions are called to calculate all Heaviside functional values and sum;
4)It calculates the value of correlation integral and seeks its density.
To obtain the statistical significance of integral density function, logarithm normal distribution fitting is carried out to it at this.In order to make experiment Data are as accurate as possible, distinguish the more pulse waves of multi collect, the youth that acquisition crowd is 21~27 years old, experiment to every subject When be in tranquility.After pulse wave is pre-processed, the probability density for taking out the data of 20s carries out lognormal fitting.Fig. 3 is One of which fitting experimental data figure, from left to right overtone order be sequentially increased, first three subharmonic of the first behavior is as a result, the second row It it is the 4th time to the 6th harmonic wave accordingly result, it is seen that each harmonic has preferable degree of fitting with logarithm normal distribution.
According to the organ resonance principle that professor Wang Weigong proposes, each organ, that is, heart of human body, liver,kidney,spleen, lung, stomach, Courage, bladder, large intestine, three cokes, small intestine, pericardium channel are corresponding in turn to the 0-11 subharmonic arteries and veins in the frequency spectrum of pulse wave signal from low to high Frequently, that is to say, that fundamental wave be the first multiplied frequency harmonic, by fundamental wave to the 4th harmonic wave corresponding to human body the five internal organs, behind each high order it is humorous Wave mainly corresponds to channels and collaterals corresponding to six internal organs and pericardium.
Analysis result shows in the present invention:In the frequency spectrum of pulse wave signal, the i.e. corresponding average palmic rate of fundamental frequency, and phase After wave group amplitude wave crest correspond to each multiplied frequency harmonic, corresponding harmonic frequency is the integral multiple of fundamental frequency, that is, palmic rate, this with The conclusion proposed in Wang Wei work blood vessel resonance theories is consistent.
The fitting degree that the present invention carries out pulse wave harmonic wave normal distribution has carried out significance analysis.By repeatedly trying It tests, calculates separately the sum for being fitted all harmonious wave fitting error amounts every time, obtain error of fitting minimum, m&gt when m=5;Error of fitting when 6 Increase.Meanwhile m&gt can be obtained by the judgement situation of normal approach;When=5, each harmonic corresponds to phase space reconstruction probability density and meets pair Number normal distribution.Harmonic wave is mainly corresponding with internal organs and passages through which vital energy circulates, when qi-blood circulating state is in health status, phase space reconstruction Probability density function shows logarithm normal distribution, shows probability density function and can be used for harmonic wave distribution detection, and harmonic wave divides The exception of cloth generally means that the change of health status, therefore can be applied to the health degree evaluation and test to each internal organs of human body.

Claims (1)

1. a kind of health state evaluating method based on pulse wave harmonic wave, it is characterised in that:
1)Human body privileged site pulse wave beat signals are detected, are filtered to obtain effective letter after removal interference noise Number, signal is split using difference threshold algorithm and identifies single pulse wave;
2)Fourier transform is carried out to single pulse wave and obtains spectrum signal, according to the SHG properties of pulse wave harmonic wave to continuous single Pulse wave harmonic amplitude carries out eigenvalue and sequence;Phase space reconfiguration is carried out to harmonic characteristic value with statistical method, Its probability density function is sought, the normal distribution state of corresponding harmonic wave probability density is analyzed;
3)Human heart, liver, taste, intestines courage functions are corresponded to according to the normal distribution state analysis of 12 subharmonic before pulse wave Health status.
CN201810628488.0A 2018-06-19 2018-06-19 Human health state evaluation method based on pulse wave harmonic waves Active CN108703748B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810628488.0A CN108703748B (en) 2018-06-19 2018-06-19 Human health state evaluation method based on pulse wave harmonic waves

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810628488.0A CN108703748B (en) 2018-06-19 2018-06-19 Human health state evaluation method based on pulse wave harmonic waves

Publications (2)

Publication Number Publication Date
CN108703748A true CN108703748A (en) 2018-10-26
CN108703748B CN108703748B (en) 2023-06-16

Family

ID=63872865

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810628488.0A Active CN108703748B (en) 2018-06-19 2018-06-19 Human health state evaluation method based on pulse wave harmonic waves

Country Status (1)

Country Link
CN (1) CN108703748B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110974158A (en) * 2019-10-23 2020-04-10 重庆特斯联智慧科技股份有限公司 Community patient identification and distress calling method, system method and system based on deep learning
WO2020031106A3 (en) * 2018-08-07 2020-06-11 Goldtech Sino Ltd Health maps for navigating a health space
CN114847889A (en) * 2022-04-20 2022-08-05 河南省肿瘤医院 Pulse reproduction method, device, storage medium, terminal equipment and system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1543911A (en) * 2003-11-13 2004-11-10 李士春 Pulse-taking sensing system for traditional Chinese medical science
US20080103398A1 (en) * 2006-10-31 2008-05-01 Chin Ming Huang Pulse Taking and Spectral Analysis and Display Instrument
CA2848498A1 (en) * 2013-06-20 2014-12-20 Mary Ann Weitnauer Systems and methods for extracting physiological characteristics using frequency harmonics
CN104665803A (en) * 2014-12-10 2015-06-03 上海理工大学 Atrial fibrillation detecting system based on intelligent platform
CN106540374A (en) * 2015-09-16 2017-03-29 赵光正 System capable of generating qi-blood resonance corresponding to channels and collaterals of human body in traditional Chinese medicine

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1543911A (en) * 2003-11-13 2004-11-10 李士春 Pulse-taking sensing system for traditional Chinese medical science
US20080103398A1 (en) * 2006-10-31 2008-05-01 Chin Ming Huang Pulse Taking and Spectral Analysis and Display Instrument
CA2848498A1 (en) * 2013-06-20 2014-12-20 Mary Ann Weitnauer Systems and methods for extracting physiological characteristics using frequency harmonics
CN104665803A (en) * 2014-12-10 2015-06-03 上海理工大学 Atrial fibrillation detecting system based on intelligent platform
CN106540374A (en) * 2015-09-16 2017-03-29 赵光正 System capable of generating qi-blood resonance corresponding to channels and collaterals of human body in traditional Chinese medicine

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020031106A3 (en) * 2018-08-07 2020-06-11 Goldtech Sino Ltd Health maps for navigating a health space
CN110974158A (en) * 2019-10-23 2020-04-10 重庆特斯联智慧科技股份有限公司 Community patient identification and distress calling method, system method and system based on deep learning
CN114847889A (en) * 2022-04-20 2022-08-05 河南省肿瘤医院 Pulse reproduction method, device, storage medium, terminal equipment and system
CN114847889B (en) * 2022-04-20 2023-07-18 河南省肿瘤医院 Pulse reproduction method, device, storage medium, terminal equipment and system

Also Published As

Publication number Publication date
CN108703748B (en) 2023-06-16

Similar Documents

Publication Publication Date Title
EP1121051B1 (en) Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
US6519486B1 (en) Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
CN110384491A (en) A kind of heart rate detection method based on common camera
Rangaprakash et al. Study of wrist pulse signals using time domain spatial features
CN108703748A (en) Health state evaluating method based on pulse wave harmonic wave
Al–Fahoum et al. A multiple signal classification approach for photoplethysmography signals in healthy and athletic subjects
Li et al. The correlation study of Cun, Guan and Chi position based on wrist pulse characteristics
Rojano et al. Singular value decomposition of the time-frequency distribution of PPG signals for motion artifact reduction
Luke et al. Motion artifact removal and feature extraction from PPG signals using efficient signal processing algorithms
Joshi et al. Arterial pulse system: modern methods for traditional Indian medicine
Ab Hamid et al. Methods of extracting feature from photoplethysmogram waveform for non-invasive diagnostic applications
Bhoi et al. Pre-processing of ppg signal with performance based methods
Chen et al. A fast ECG diagnosis using frequency-based compressive neural network
Liu et al. The differential method of phase space matrix for AF/VF discrimination application
More et al. Development of non-invasive diagnostic tool for diseases using Photo Plethysmography
Kadarmandalgi et al. Pulse diagnosis System for Nadi Pariksha using parametric and statistical analysis
Zhang et al. Human pulse recognition based on wavelet transform and BP network
Afkhami et al. A dynamic model of brain hemodynamics in near-infrared spectroscopy
He et al. De-noising of 3D pulse images by channel-weighted robust principal component analysis
Deepakfranklin Survey on Methods of Obtaining Biomedical Parameters from PPG Signal
Gbaoui et al. Decomposition of photoplethysmographical arterial pulse waves by independent component analysis: possibilities and limitations
Ekhlasi et al. Investigating the differences between atrial fibrillation and normal ECG signals using the EMD method
Burrello et al. Improving PPG-based Heart-Rate Monitoring with Synthetically Generated Data
Priyadarshini et al. Nonlinear Processing of Wrist Pulse Signals to Distinguish Diabetic and Non-Diabetic Subjects
Gao Wearable Sensing of Cardio-Pulmonary Function: Non-invasive Sensor Design and Statistical approaches to Signal Compression and Analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20181106

Address after: 200000 Shanghai, China (Shanghai) free trade zone, 14 Lane 1502, Lo Shan Road.

Applicant after: SHANGHAI JIUJIANTANG HEALTH MANAGEMENT Co.,Ltd.

Address before: 200433 302, 38 Lane 179 lane, Shun Shun Road, Yangpu District, Shanghai.

Applicant before: Sun Ying

Applicant before: Ye Ying

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230524

Address after: Building C, 888 Huanhu West 2nd Road, Lingang New Area, Pudong New Area, Shanghai, November 2013

Applicant after: Zhongyu Life (Shanghai) Medical Technology Development Co.,Ltd.

Address before: 200000 Shanghai, China (Shanghai) free trade zone, 14 Lane 1502, Lo Shan Road.

Applicant before: SHANGHAI JIUJIANTANG HEALTH MANAGEMENT Co.,Ltd.

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240327

Address after: March 2014 Room 929, Building 28, No. 6055 Jinhai Road, Fengxian District, Shanghai

Patentee after: Shanghai Qianmai Health Technology Co.,Ltd.

Country or region after: China

Address before: Building C, 888 Huanhu West 2nd Road, Lingang New Area, Pudong New Area, Shanghai, November 2013

Patentee before: Zhongyu Life (Shanghai) Medical Technology Development Co.,Ltd.

Country or region before: China