CN103989463B - A kind of radial artery pulse wave detection system based on fingertip pulse wave and method - Google Patents
A kind of radial artery pulse wave detection system based on fingertip pulse wave and method Download PDFInfo
- Publication number
- CN103989463B CN103989463B CN201410209464.3A CN201410209464A CN103989463B CN 103989463 B CN103989463 B CN 103989463B CN 201410209464 A CN201410209464 A CN 201410209464A CN 103989463 B CN103989463 B CN 103989463B
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- pulse
- radial
- radial artery
- finger tip
- 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.)
- Active
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Based on radial artery pulse wave detection system and the method for fingertip pulse wave, this system comprises finger tip pulse transducer, radial pulse sensor, pulse signal acquisition module, central processing module and input/output module; Pulse signal acquisition module comprises finger tip pulse signal acquisition module and radial pulse signal acquisition module; Finger tip pulse transducer, radial pulse sensor are placed in finger tip place and wrist place respectively, the outfan of finger tip pulse transducer, the outfan of radial pulse sensor connect the input of finger tip pulse signal acquisition module, the input of radial pulse signal acquisition module respectively, and the outfan of finger tip pulse signal acquisition module, the outfan of radial pulse signal acquisition module connect the input of central processing module respectively.Method of the present invention rebuilds the optimum waveform matching degree of the radial pulse waveform of radial pulse waveform and the actual measurement obtained all in 70% to 92% scope, is applicable to the collection of radial pulse wave datum.
Description
Technical field
The present invention relates to technical field of medical equipment, be specifically related to a kind of radial artery pulse wave detection system based on fingertip pulse wave and method.
Background technology
Radial artery is one in the some shallow table tremulous pulsies of human body, be compared to the pulse wave at other positions, the detection of radial artery pulse wave is more convenient, no matter be therefore the feeling the pulse of the traditional Chinese medical science, or the human blood glucose of doctor trained in Western medicine central vessel parameter, often always at radial artery place obtaining information.The information that radial artery place detects can reflect some physiology, the pathological characters of human body really to a certain extent, this point already confirm by the pulse-taking of the traditional medicine such as the traditional Chinese medical science, Tibetan medicine.
About the detection of radial artery pulse wave, now conventional method wears a pulse transducer at radial artery place.Operable kind of sensor has a lot, as piezoelectric type pulse transducer, pressure resistance type pulse transducer and acceleration pulse transducer etc.But compare with the measurement of fingertip pulse wave, the pulse wave measurement at radial artery place is more complicated, needs the position of finding radial artery accurately; And compared with wearable sensors on pointing, wrist place wearable sensors is larger to human influence.Although fingertip pulse wave gathers more convenient, compare with radial artery pulse wave, due to the impact that it is subject to microcirculqtory system in palm, the characteristic information amount contained is less.
Summary of the invention
For the deficiency that prior art exists, the invention provides a kind of radial artery pulse wave detection system based on fingertip pulse wave and method.
Technical scheme of the present invention is:
Based on a radial artery pulse wave detection system for fingertip pulse wave, comprise finger tip pulse transducer, radial pulse sensor, pulse signal acquisition module, central processing module and input/output module;
Described pulse signal acquisition module comprises finger tip pulse signal acquisition module and radial pulse signal acquisition module;
Described finger tip pulse transducer, radial pulse sensor are placed in finger tip place and wrist place respectively, the outfan of finger tip pulse transducer, the outfan of radial pulse sensor connect the input of finger tip pulse signal acquisition module, the input of radial pulse signal acquisition module respectively, and the outfan of finger tip pulse signal acquisition module, the outfan of radial pulse signal acquisition module connect the input of central processing module respectively;
Input/output module comprises button/touch screen, LCD screen and SD card, and button/touch screen connects the input of central processing module, and SD card and LCD screen are connected the outfan of central processing module respectively.
Described finger tip pulse signal acquisition module and radial pulse signal acquisition module include pre-amplification circuit, pulse signal acquisition cuicuit and pressure signal acquisition cuicuit;
Described pulse signal acquisition cuicuit comprises high-pass filtering circuit, the first low-pass filter circuit, the first second amplifying circuit and level lifting circuit;
Described pressure signal acquisition cuicuit comprises the second low-pass filter circuit, the second second amplifying circuit and inverter circuit;
The outfan of finger tip pulse transducer, the outfan of radial pulse sensor connects the input of pre-amplification circuit respectively, the outfan of pre-amplification circuit connects the input of high pass filter and the input of the second low pass filter respectively, the outfan of high pass filter connects the input of the first low pass filter, the outfan of the first low pass filter connects the input of the first second amplifying circuit, the outfan of the first second amplifying circuit connects the input of level lifting circuit, the outfan of level lifting circuit connects central processing module, the outfan of the second low pass filter connects the input of the second second amplifying circuit, the outfan of the second second amplifying circuit connects the input of inverter circuit, the outfan of inverter circuit connects central processing module.
The radial artery pulse wave detection system based on fingertip pulse wave described in employing carries out the method for radial artery pulse wave detection, comprises the following steps:
Step 1: respectively pressure calibration is carried out to finger tip pulse transducer, radial pulse sensor, optimum contact pressure when pulse wave being detected to the optimum contact pressure of finger tip applying and radial pulse sensor during setting finger tip pulse transducer detection pulse wave, wrist radial artery place applied;
Step 1.1: finger tip pulse transducer, radial pulse sensor gather fingertip pulse wave and radial artery pulse wave respectively;
Step 1.2: central processing module control finger tip pulse transducer applies different contact pressures to finger tip and radial pulse wave sensor opponent wrist radial artery place applies different contact pressures, gathers the fingertip pulse wave under different contact pressure state and radial artery pulse wave signal;
Step 1.3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process, and the signal after process is sent to central processing module;
Step 1.4: according to the amplitude of the fingertip pulse wave signal under the different contact pressures collected and the amplitude of radial artery pulse wave signal, determine the contact pressure that the pulse wave of amplitude maximum in the contact pressure that the pulse wave of amplitude maximum in fingertip pulse wave signal is corresponding and radial artery pulse wave signal is corresponding, these two contact pressures are set as respectively finger tip pulse transducer detect pulse wave time to finger tip apply optimum contact pressure and radial pulse sensor detect pulse wave time to wrist radial artery place apply optimum contact pressure;
Step 2: the fingertip pulse wave using finger tip pulse transducer collection experimenter under the optimum contact pressure when the finger tip pulse transducer set detects pulse wave to finger tip applying, uses the radial artery pulse wave of radial pulse sensor acquisition experimenter when the radial pulse sensor detection pulse wave set to wrist radial artery place under the optimum contact pressure applied;
Step 3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process;
Step 4: set up the radial artery pulse wave reconstruction model based on transfer function, this model be input as fingertip pulse wave, export as radial artery pulse wave;
Step 4.1: solve the frequency domain transfer function of radial artery pulse wave to fingertip pulse wave;
Step 4.1.1: the fingertip pulse wave using Welch ' s average cycle method to estimate and the cross-spectral density of radial artery pulse wave;
Step 4.1.2: use Welch ' s average cycle method to estimate the autopower spectral density of fingertip pulse wave;
Step 4.1.3: utilize cross-spectral density divided by autopower spectral density, namely obtains the frequency domain transfer function of fingertip pulse wave and radial artery pulse wave;
Step 4.2: determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Step 4.2.1: the radial artery pulse wave reconstruction model based on transfer function setting up the different orders under different window width respectively;
Step 4.2.2: utilize root-mean-square percent error and Waveform Matching degree assessment radial artery pulse wave reconstruction model precision, obtain the best window width of radial artery pulse wave reconstruction model and best order, namely determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Step 4.3: set up the radial artery pulse wave reconstruction model based on transfer function according to the window width of radial artery pulse wave reconstruction model and order;
Step 5: utilize and rebuild radial artery pulse wave based on the radial artery pulse wave reconstruction model of transfer function and the fingertip pulse wave of Real-time Collection, obtain radial artery pulse wave testing result, this testing result is stored into SD card and is shown by LCD screen.
Beneficial effect:
The optimum waveform matching degree of the radial pulse waveform of the radial pulse waveform that obtains based on radial artery pulse wave detection system and the method for fingertip pulse wave of the present invention and actual measurement is all in 70% to 92% scope, this device is applicable to the collection of radial pulse wave datum, compared with radial pulse checkout gear of the prior art, make the detection of radial artery pulse wave convenient.
Accompanying drawing explanation
Fig. 1 is the radial artery pulse wave detection system structured flowchart based on fingertip pulse wave of the specific embodiment of the invention;
Fig. 2 be the specific embodiment of the invention finger tip radial pulse signal acquisition module signal move towards schematic diagram;
Fig. 3 is the pre-amplification circuit schematic diagram of the specific embodiment of the invention;
Fig. 4 is the high-pass filtering circuit schematic diagram of the specific embodiment of the invention;
Fig. 5 is the first low-pass filter circuit schematic diagram of the specific embodiment of the invention;
Fig. 6 is the first second amplifying circuit schematic diagram of the specific embodiment of the invention;
Fig. 7 is the level lifting circuit theory diagrams of the specific embodiment of the invention;
Fig. 8 is the second low-pass filter circuit schematic diagram of the specific embodiment of the invention;
Fig. 9 is the second second amplifying circuit schematic diagram of the specific embodiment of the invention;
Figure 10 is the inverter circuit schematic diagram of the specific embodiment of the invention;
Figure 11 is the method flow diagram that the radial artery pulse wave of the specific embodiment of the invention detects;
Figure 12 is the SD card circuit theory diagrams of the specific embodiment of the invention;
Figure 13 is the transfer function modeling training and testing procedural block diagram of the specific embodiment of the invention;
Figure 14 is PRD and the FIT scattergram of the best model estimated result under the different window widths of the specific embodiment of the invention;
PRD and the FIT scattergram of Figure 15 is subjects's window width of the specific embodiment of the invention when being 1024 each order model estimated result;
Figure 16 is the working-flow figure of the specific embodiment of the invention.
Detailed description of the invention
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is elaborated.
As shown in Figure 1, based on the radial artery pulse wave detection system of fingertip pulse wave, comprise finger tip pulse transducer, radial pulse sensor, pulse signal acquisition module, central processing module and input/output module;
Pulse signal acquisition module comprises finger tip pulse signal acquisition module and radial pulse signal acquisition module;
Finger tip pulse transducer, radial pulse sensor are placed in finger tip place and wrist place respectively, the outfan of finger tip pulse transducer, the outfan of radial pulse sensor connect the input of finger tip pulse signal acquisition module, the input of radial pulse signal acquisition module respectively, and the outfan of finger tip pulse signal acquisition module, the outfan of radial pulse signal acquisition module connect the input of central processing module respectively;
Input/output module comprises button/touch screen, LCD screen and SD card, and button/touch screen connects the input of central processing module, and SD card and LCD screen are connected the outfan of central processing module respectively.
Finger tip pulse signal acquisition module and radial pulse signal acquisition module include pre-amplification circuit, pulse signal acquisition cuicuit and pressure signal acquisition cuicuit;
Pulse signal acquisition cuicuit comprises high-pass filtering circuit, the first low-pass filter circuit, the first second amplifying circuit and level lifting circuit;
Pressure signal acquisition cuicuit comprises the second low-pass filter circuit, the second second amplifying circuit and inverter circuit;
The outfan of finger tip pulse transducer, the outfan of radial pulse sensor passes through connector (as aviation plug, DB plug etc.) connect the input of pre-amplification circuit respectively, as shown in Figure 2, the outfan of pre-amplification circuit connects the input of high pass filter and the input of the second low pass filter respectively, the outfan of high pass filter connects the input of the first low pass filter, the outfan of the first low pass filter connects the input of the first second amplifying circuit, the outfan of the first second amplifying circuit connects the input of level lifting circuit, the outfan of level lifting circuit connects central processing module, the outfan of the second low pass filter connects the input of the second second amplifying circuit, the outfan of the second second amplifying circuit connects the input of inverter circuit, the outfan of inverter circuit connects central processing module, wherein, level lifting circuit exports pulse wave signal, inverter circuit output pressure signal.
The main core of pre-amplification circuit adopts instrumentation amplifier AD620, and its circuit as shown in Figure 3.Present embodiment, in order to obtain pulse signal (AC signal) and pressure signal (direct current signal) respectively, introduces pulse signal acquisition cuicuit and pressure signal acquisition cuicuit, these two circuit with the OPA2277 of TI company for master chip.
Pulse signal acquisition cuicuit comprises high-pass filtering circuit, the first low-pass filter circuit, the first second amplifying circuit and level lifting circuit.Signal trend is: first pre-amp output signal connects the flip-flop in high-pass filtering circuit (as shown in Figure 4) filtered signal; Then the high-frequency noise in the second low-pass filter circuit (as shown in Figure 5) filtering pulse signal is connected; Then access the first second amplifying circuit (as shown in Figure 6) pulse signal is amplified; The pulse wave signal of AD sampling condition is met after finally accessing level lifting circuit (as shown in Figure 7).
Pulse transducer pressure signal acquisition cuicuit comprises: the second low-pass filter circuit, the second second amplifying circuit and inverter circuit.Signal trend is: first pre-amp output signal connects the alternating component in the second low-pass filter circuit (as shown in Figure 8) filtered signal; Then access the second second amplifying circuit (as shown in Figure 9) pressure signal is amplified, finally access the pressure signal that inverter circuit (as shown in Figure 10) is met AD sampling condition.
Central processing module adopts arm processor, the AD sampling module input end of analog signal of the output termination arm processor of finger tip pulse signal acquisition module and radial pulse signal acquisition module, completes AD sample by A/D module built-in in ARM; Because pulse signal frequency range concentrates within the scope of 0-40Hz, according to sampling thheorem, in present embodiment, AD sample frequency is set as 100Hz, sampling precision is 12.
As shown in figure 12, SD card is communicated by SPI Interface realization with between arm processor SD card circuit, and adopt Pop-up SD card slot, in circuit, " C/D " foot is for detecting in slot whether insert SD card.When carrying out SD card initialization, SPI output clock frequency configuration is after the success of 330Khz, SD card initialization, when carrying out read-write operation to SD card data, is 16.5MHz by SPI output clock frequency configuration.
When carrying out data and storing, first the establishment of the complete paired data file of FAT16 file system in SD card is called, in present embodiment, the form of capture file and destination file is " dat " form, naming method is respectively " fpwav000 " and " rpwav000 ", after filename three bit digital represent respectively subjects check numbering individual position, ten and hundred, be convenient to the record of data.
Button or touch screen is used to realize the switching between each functional module of system in present embodiment.Comprise 5 (Key1, Key2, Key3, Key4, Reset) buttons in system, by the mode of outer interruption, in each keys interrupt service subprogram, realize corresponding function;
Present embodiment central processing module employing Samsung model is the ARM plate of S3C6410, it is inner containing touch screen interface and LCD screen interface, use four traditional line touch screen modes, by the judgement of the sampled result of the magnitude of voltage to point touching screen action generation, obtain coordinate information, at touch screen interrupt service subroutine, according to this coordinate information, perform and operate accordingly.
The LCD screen of present embodiment selects Qun Chuan company model to be the LCD liquid crystal display screen of AT070TN92 (7 cun of screens), adopts 18 (6: 6: 6) RGB true color patterns.
Can be asked for leave by the collection of LCD screen video data and data analysis in present embodiment, the pulse data waveform of the in real time numbering of the current subjects of display and actual measurement, for the measured data of data analysis and Waveform Matching degree (FIT) value between the radial artery pulse wave of testing result data and reconstruction and the radial artery pulse wave of actual measurement and root-mean-square percent error (PRD) value.
The character of LCD screen display is the dot matrix of 16*16, and character delivery mode is the horizontal delivery of a left high position, and data array is from left to right, from top to bottom; The waveform display area background color of data acquisition interface is set to black, grid lines are red dotted line (between points every three pixels), abscissa value represents the time, get ready and be spaced apart 10ms, every 100 pixel length representative 1s, vertical coordinate represents sample magnitude, unit is " volt ", two path pulse waves (radial artery and finger tip pulse) sample result value is mapped as the ordinate value in region 1 and region 2 respectively, and mapping relations are as shown in formula (1) and (2):
Y
1＝51+(4095-V
1)/22 (1)
Y
2＝259+(4095-V
2)/22 (2)
Wherein, Y
1and Y
2represent the ordinate value of viewing area 1 and viewing area 2 respectively, V
1and V
2represent two path data sample result value respectively.
Adopt the radial artery pulse wave detection system based on fingertip pulse wave to carry out the method for radial artery pulse wave detection, its flow process as shown in figure 11, comprises the following steps:
Step 1: respectively pressure calibration is carried out to finger tip pulse transducer, radial pulse sensor, optimum contact pressure when pulse wave being detected to the optimum contact pressure of finger tip applying and radial pulse sensor during setting finger tip pulse transducer detection pulse wave, wrist radial artery place applied;
Step 1.1: finger tip pulse transducer, radial pulse sensor gather fingertip pulse wave and radial artery pulse wave respectively;
Step 1.2: central processing module control finger tip pulse transducer applies different contact pressures to finger tip and radial pulse wave sensor opponent wrist radial artery place applies different contact pressures, gathers the fingertip pulse wave under different contact pressure state and radial artery pulse wave signal;
Step 1.3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process, and the signal after process is sent to central processing module;
Step 1.4: according to the different contact pressure (50g collected, 60g, ... 140g, the amplitude of the fingertip pulse wave signal 150g) and the amplitude of radial artery pulse wave signal, determine the contact pressure that the pulse wave of amplitude maximum in the contact pressure that the pulse wave of amplitude maximum in fingertip pulse wave signal is corresponding and radial artery pulse wave signal is corresponding, these two contact pressures are set as respectively finger tip pulse transducer detect pulse wave time to finger tip apply optimum contact pressure and radial pulse sensor detect pulse wave time to wrist radial artery place apply optimum contact pressure,
The finger tip pulse transducer of present embodiment and radial pulse sensor can select motor adding pressure type/inflating pressure formula sensor.First, system carries out zeroing operation to two sensors, after zeroing successful operation, completes the calibration to finger tip pulse transducer and radial pulse cell pressure respectively, under making it be operated in optimum pressure state.
Step 2: the fingertip pulse wave using finger tip pulse transducer collection experimenter under the optimum contact pressure when the finger tip pulse transducer set detects pulse wave to finger tip applying, uses the radial artery pulse wave of radial pulse sensor acquisition experimenter when the radial pulse sensor detection pulse wave set to wrist radial artery place under the optimum contact pressure applied;
Step 3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process;
Step 4: set up the radial artery pulse wave reconstruction model based on transfer function, this model be input as fingertip pulse wave, export as radial artery pulse wave;
Solving fingertip pulse wave in present embodiment to the square ratio juris of radial artery pulse wave frequency domain transfer function Txy (f) is: the power spectral density of cross-spectral density divided by input signal using input/output signal, and its expression formula is:
Txy(f)＝(Pyx(f))/(Pxx(f)) (3)
Wherein, Pyx (f) is the cross-spectral density between x and y, the power spectral density that Pxx (f) is x, and the two is estimated by Welch ' s average cycle method.Using radial artery pulse wave Pr as input, fingertip pulse wave Pf, as output, can try to achieve the frequency domain transfer function Trf of radial artery pulse wave to fingertip pulse wave.The time domain FIR model of different order can be tried to achieve according to the amplitude of the Trf tried to achieve and phase information.
Step 4.1: solve the frequency domain transfer function of radial artery pulse wave to fingertip pulse wave;
Step 4.1.1: the fingertip pulse wave using Welch ' s average cycle method to estimate and the cross-spectral density of radial artery pulse wave;
Step 4.1.2: use Welch ' s average cycle method to estimate the autopower spectral density of fingertip pulse wave;
Step 4.1.3: utilize cross-spectral density divided by autopower spectral density, namely obtains the frequency domain transfer function of fingertip pulse wave and radial artery pulse wave;
Knownly from above-mentioned steps solve radial artery pulse wave is Welch ' s average cycle method use to the core of the frequency domain transfer function of fingertip pulse wave.Welch ' s average cycle method is that the random fluctuation problem in order to overcome periodic method proposes.Random fluctuation refers to when the power spectrum of life cycle method signal calculated, and along with got signal sequence length is different, the periodogram obtained is different.Average cycle method carries out segmentation windowing to signal, first obtains the periodogram of every segment signal, be then averaged the periodogram of gained, obtain the power spectrum of signal.Namely random fluctuation can be reduced to signal windowing process.Therefore, the selection of window width on power spectrum signal to solve impact very large, and then to the frequency domain transfer function of input/output signal solve and the precision of FIR model plays decisive role.
Step 4.2: determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Step 4.2.1: the radial artery pulse wave reconstruction model based on transfer function setting up the different orders under different window width respectively;
Step 4.2.2: utilize root-mean-square percent error and Waveform Matching degree assessment radial artery pulse wave reconstruction model precision, obtain the best window width of radial artery pulse wave reconstruction model and best order, namely determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Present embodiment, in order to the performance of estimation model algorithm objectively, adopts the evaluation index conventional in compress ecg data field: root-mean-square percent error (PRD) and Waveform Matching degree (FIT).The computing formula of the two is:
Wherein, u (i) is primary signal, and v (i) is the signal estimated,
it is the average of primary signal u (i).
Present embodiment is in order to determine the window width of fingertip pulse wave to radial artery pulse wave transfer function model, 10 different window widths (128 are set up respectively for large sample subjects, 256,384,512, ..., 1024,1152,1280) the 1-110 rank FIR model under, the radial artery pulse wave Waveform Matching degree (FIT) of the radial artery waveform estimated with each model and actual measurement is evaluation index, FIR model corresponding under determining best window width and best order.
Under setting up different window width (128,256,384,512 ..., 1024,1152,1280) 1 ~ 110 rank FIR model, the FIT value of the FIR model estimated signal of the different orders under more same window width, with maximum FIT value for standard, obtain the order of the best FIR model under this window width, with best FIT value for standard, compare the best FIR model under the different window widths obtained, obtain the best window width of radial artery pulse wave reconstruction model and best order.
Figure 14 is PRD and the FIT scattergram of the best model estimated result under the different window width of a subjects.The best window width of the model of ten subjectss in large sample subjects and order as shown in table 1.
The window width of the best model of table 10 subjectss and order
Known by table 1 data, the window width of FIR model the best of ten subjectss be greater than 896 have seven people.And pulse wave analysis product (the SphygmoCor pulse wave analysis instrument etc. as Australia) advanced on market, often choose 10s pulse data and carry out analyzing and processing.The system sampling frequency of present embodiment is 100Hz, therefore corresponding 1000 sampled points of the data length of 10s.According to above-mentioned two standards, the window width based on the radial artery pulse wave reconstruction model of transfer function is chosen as 896,1024,1152 and 1280.Then comparison finds by experiment, determines best window width.
Present embodiment, after the best window width of Confirming model, sets up 1 ~ 110 rank FIR model under this window width respectively to large sample subjects.According to PRD and the FIT scattergram (as shown in figure 15, PRD and the FIT scattergram for each order model estimated result of 1 subjects) of each order model estimated result of subjects, set forth the concrete steps of the best order of Confirming model.
Present embodiment finds when setting up the radial artery pulse wave reconstruction model based on transfer function: when use fingertip pulse wave is as input signal, radial artery pulse wave as output signal, the model accuracy directly set up to far be worse than the model set up as output signal using radial artery pulse wave as input signal fingertip pulse wave get inverse after model accuracy.This is mainly because fingertip pulse wave is owing to being subject to the impact of the inner microcirculqtory system of palm, its quantity of information contained will be less than radial artery pulse wave, according to the knowledge of system identification, the quantity of information of input signal must more than the quantity of information of output signal, therefore first present embodiment solves the time domain FIR model of radial artery pulse wave to fingertip pulse wave, then according to the character of convolution, get inverse to gained model, obtain the FIR model of fingertip pulse wave to radial artery pulse wave.
Step 4.3: set up the radial artery pulse wave reconstruction model based on transfer function according to the window width of radial artery pulse wave reconstruction model and order;
Step 5: utilize and rebuild radial artery pulse wave based on the general radial artery pulse wave reconstruction model of transfer function and the fingertip pulse wave of Real-time Collection, obtain radial artery pulse wave testing result, this testing result is stored into SD card and is shown by LCD screen.
Figure 16 is the working-flow figure of present embodiment.Upon power-up of the system, first the relevant initialization operations such as I/O mouth, intervalometer, A/D module are performed, before system performs initialization operation to SPI interface and SD again, can first detect SD card whether to exist, if detect without SD card in system, buzzer rings please insert SD card to point out user, time in system without SD card three times, the main interface state of system will be in always, data acquisition and data analysis operation cannot be performed; SD card is there is in system, and after completing relevant initialization operation, data acquisition operations can be performed, now LCD liquid crystal display screen can show the data waveform collected in real time, after user judges that having collected 10s stablizes data, can stop gathering, now the last 10s data collected can be rebuild stored in the radial artery pulse wave performed in the capture file in SD card based on transfer function and fingertip pulse wave simultaneously by system, complete the reconstruction to radial artery pulse wave, reconstructed results is stored in SD card destination file and also shows on the lcd screen, complete one-time detection operation.
Claims (1)
1. a method for radial artery pulse wave detection, adopts a kind of radial artery pulse wave detection system based on fingertip pulse wave, comprises finger tip pulse transducer, radial pulse sensor, pulse signal acquisition module, central processing module and input/output module;
Described pulse signal acquisition module comprises finger tip pulse signal acquisition module and radial pulse signal acquisition module;
Described finger tip pulse transducer, radial pulse sensor are placed in finger tip place and wrist place respectively, the outfan of finger tip pulse transducer, the outfan of radial pulse sensor connect the input of finger tip pulse signal acquisition module, the input of radial pulse signal acquisition module respectively, and the outfan of finger tip pulse signal acquisition module, the outfan of radial pulse signal acquisition module connect the input of central processing module respectively;
Input/output module comprises button/touch screen, LCD screen and SD card, and button/touch screen connects the input of central processing module, and SD card and LCD screen are connected the outfan of central processing module respectively;
It is characterized in that: the method comprises the following steps:
Step 1: respectively pressure calibration is carried out to finger tip pulse transducer, radial pulse sensor, optimum contact pressure when pulse wave being detected to the optimum contact pressure of finger tip applying and radial pulse sensor during setting finger tip pulse transducer detection pulse wave, wrist radial artery place applied;
Step 1.1: finger tip pulse transducer, radial pulse sensor gather fingertip pulse wave and radial artery pulse wave respectively;
Step 1.2: central processing module control finger tip pulse transducer applies different contact pressures to finger tip and radial pulse wave sensor opponent wrist radial artery place applies different contact pressures, gathers the fingertip pulse wave under different contact pressure state and radial artery pulse wave signal;
Step 1.3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process, and the signal after process is sent to central processing module;
Step 1.4: according to the amplitude of the fingertip pulse wave signal under the different contact pressures collected and the amplitude of radial artery pulse wave signal, determine the contact pressure that the pulse wave of amplitude maximum in the contact pressure that the pulse wave of amplitude maximum in fingertip pulse wave signal is corresponding and radial artery pulse wave signal is corresponding, these two contact pressures are set as respectively finger tip pulse transducer detect pulse wave time to finger tip apply optimum contact pressure and radial pulse sensor detect pulse wave time to wrist radial artery place apply optimum contact pressure;
Step 2: the fingertip pulse wave using finger tip pulse transducer collection experimenter under the optimum contact pressure when the finger tip pulse transducer set detects pulse wave to finger tip applying, uses the radial artery pulse wave of radial pulse sensor acquisition experimenter when the radial pulse sensor detection pulse wave set to wrist radial artery place under the optimum contact pressure applied;
Step 3: finger tip pulse signal acquisition module and radial pulse signal acquisition module are amplified the fingertip pulse wave collected and radial artery pulse wave respectively, filtering and level lifting process;
Step 4: set up the radial artery pulse wave reconstruction model based on transfer function, this model be input as fingertip pulse wave, export as radial artery pulse wave;
Step 4.1: solve the frequency domain transfer function of radial artery pulse wave to fingertip pulse wave;
Step 4.1.1: the fingertip pulse wave using Welch ' s average cycle method to estimate and the cross-spectral density of radial artery pulse wave;
Step 4.1.2: use Welch ' s average cycle method to estimate the autopower spectral density of fingertip pulse wave;
Step 4.1.3: utilize cross-spectral density divided by autopower spectral density, namely obtains the frequency domain transfer function of fingertip pulse wave and radial artery pulse wave;
Step 4.2: determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Step 4.2.1: the radial artery pulse wave reconstruction model based on transfer function setting up the different orders under different window width respectively;
Step 4.2.2: utilize root-mean-square percent error and Waveform Matching degree assessment radial artery pulse wave reconstruction model precision, obtain the best window width of radial artery pulse wave reconstruction model and best order, namely determine the window width based on the radial artery pulse wave reconstruction model of transfer function and order;
Step 4.3: set up the radial artery pulse wave reconstruction model based on transfer function according to the window width of radial artery pulse wave reconstruction model and order;
Step 5: utilize and rebuild radial artery pulse wave based on the radial artery pulse wave reconstruction model of transfer function and the fingertip pulse wave of Real-time Collection, obtain radial artery pulse wave testing result, this testing result is stored into SD card and is shown by LCD screen.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410209464.3A CN103989463B (en) | 2014-05-16 | 2014-05-16 | A kind of radial artery pulse wave detection system based on fingertip pulse wave and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410209464.3A CN103989463B (en) | 2014-05-16 | 2014-05-16 | A kind of radial artery pulse wave detection system based on fingertip pulse wave and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103989463A CN103989463A (en) | 2014-08-20 |
| CN103989463B true CN103989463B (en) | 2015-09-30 |
Family
ID=51303985
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410209464.3A Active CN103989463B (en) | 2014-05-16 | 2014-05-16 | A kind of radial artery pulse wave detection system based on fingertip pulse wave and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103989463B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017051276A (en) * | 2015-09-07 | 2017-03-16 | オムロンヘルスケア株式会社 | Pulse wave detection device |
| CN112263230B (en) * | 2016-09-22 | 2023-05-09 | 上海潓美医疗科技有限公司 | Dynamic blood pressure monitoring system and method based on radial artery biosensor technology |
| CN106580273A (en) * | 2016-11-28 | 2017-04-26 | 中国科学院微电子研究所 | Pulse wave acquisition device and pulse wave acquisition calibration method |
| CN106725363B (en) * | 2016-11-28 | 2020-09-25 | 中国科学院微电子研究所 | Pulse wave acquisition device and pulse wave acquisition calibration method |
| CN106667455A (en) * | 2016-12-27 | 2017-05-17 | 合肥讯创信息科技有限公司 | Simulation pulse instrument |
| CN112826459B (en) * | 2021-01-08 | 2022-11-29 | 北京工业大学 | Pulse wave waveform reconstruction method and system based on convolution self-encoder |
| CN113545760B (en) * | 2021-07-29 | 2023-08-18 | 西安仙峒科技有限责任公司 | Traditional Chinese medicine pulse diagnosis tester based on modal analysis technology and application method thereof |
| CN113724879B (en) * | 2021-09-08 | 2024-02-09 | 上海中医药大学 | Method for establishing cardiovascular disease identification model by using mucoid optimization algorithm |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2505462B2 (en) * | 1987-06-10 | 1996-06-12 | コ−リン電子株式会社 | Pulse wave continuous measurement device |
| EP1433417B1 (en) * | 1995-05-12 | 2007-12-05 | Seiko Epson Corporation | Device for controlling a physiological state |
| JP3988674B2 (en) * | 2003-04-21 | 2007-10-10 | オムロンヘルスケア株式会社 | Pulse wave velocity information measuring device |
| CN1631319A (en) * | 2003-12-23 | 2005-06-29 | 罗志昌 | Blood flow parameter non-invasive detecting method based on finger tip volume pulsation blood flow wave |
| JP3940150B2 (en) * | 2005-03-03 | 2007-07-04 | シチズンホールディングス株式会社 | Caffres electronic blood pressure monitor |
-
2014
- 2014-05-16 CN CN201410209464.3A patent/CN103989463B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN103989463A (en) | 2014-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103989463B (en) | A kind of radial artery pulse wave detection system based on fingertip pulse wave and method | |
| Luo et al. | Possibility of quantifying TCM finger-reading sensations: I. Bi-sensing pulse diagnosis instrument | |
| CN101815466B (en) | A non-invasive device NADI TARANGINI useful for quantitative detection of arterial NADI pulse waveform | |
| CN101732040B (en) | Non-invasive multipath pulse wave detection device, system and analytical system | |
| CN101357062B (en) | Vital signs parameter evaluation device based on volume pulsation wave detection | |
| CN105433915B (en) | A kind of traditional Chinese pulse-diagnosis instrument and its pulse condition detection method | |
| CN107529646A (en) | A kind of non-contact type method for measuring heart rate and device based on Euler's image zoom | |
| CN107961001A (en) | Appraisal procedure, device and the atherosclerosis detector of Degree of arteriosclerosis | |
| CN110840405A (en) | Manufacturing method of flexible array sensor and pulse condition detection system based on flexible array sensor | |
| CN104757959B (en) | Pulse wave transmission velocity detecting method and system based on image foldover | |
| CN204121004U (en) | A kind of pulse diagnostic information acquisition device in conjunction with multiple pressure-type sensor | |
| CN202960481U (en) | Traditional Chinese medicine pulse condition acquisition device | |
| CN104546007B (en) | Anti-interference processing method and device for fetal movement detection | |
| EP3364856B1 (en) | Method and device for heart rate detection with multi-use capacitive touch sensors | |
| CN103393402A (en) | Portable burn and trauma surface binding bandage pressure monitor | |
| CN105078439B (en) | The method of display device, heart rate monitoring system and rhythm of the heart | |
| CN110327023A (en) | A kind of array of pressure sensors and diagnosis by feeling the pulse device | |
| CN110200642A (en) | A kind of measurement method and terminal of cognitive load and psychological pressure | |
| CN201244022Y (en) | Instrument for measuring pulse wave and analyzing physiological characteristic parameter | |
| CN201894645U (en) | Novel intelligent electrocardiogram test healthcare apparatus | |
| CN116369888B (en) | Non-contact heart rate variability data acquisition method and device | |
| US20150313489A1 (en) | Wearable Device and Method for Assessing Cardiac Wellness by Analyzing High Definition EKG Signals | |
| CN112089406A (en) | Detection hardware system for traditional Chinese medicine pulse-taking information | |
| CN104027097A (en) | Vascular function noninvasive detecting method and device | |
| CN104188651A (en) | Electrocardiogram monitoring device and control method of electrocardiogram monitoring device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |