CN108245749A - Respiratory flow method of counting, device, computer equipment and medium - Google Patents
Respiratory flow method of counting, device, computer equipment and medium Download PDFInfo
- Publication number
- CN108245749A CN108245749A CN201711488861.9A CN201711488861A CN108245749A CN 108245749 A CN108245749 A CN 108245749A CN 201711488861 A CN201711488861 A CN 201711488861A CN 108245749 A CN108245749 A CN 108245749A
- Authority
- CN
- China
- Prior art keywords
- sliding window
- data
- arithmetic mean
- instantaneous value
- window array
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0003—Accessories therefor, e.g. sensors, vibrators, negative pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0003—Accessories therefor, e.g. sensors, vibrators, negative pressure
- A61M2016/003—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter
- A61M2016/0033—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical
- A61M2016/0036—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical in the breathing tube and used in both inspiratory and expiratory phase
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M2205/00—General characteristics of the apparatus
- A61M2205/33—Controlling, regulating or measuring
- A61M2205/3327—Measuring
Abstract
The application provides a kind of respiratory flow method of counting, device, computer equipment and medium, wherein, method includes：Obtain respiratory flow sampled data, obtain original data on flows, original data on flows is written in the sliding window of default first data length, update calculates the arithmetic mean of instantaneous value of single order sliding window array, arithmetic mean of instantaneous value is written in the sliding window of default second data length, according to second order sliding window array, lung ventilator flow count results are obtained.In whole process, without complicated filtering process, using simple mean value calculation, and by sliding window mode accurate counting respiratory flow, data support is provided to improve lung ventilator response speed.
Description
Technical field
The present invention relates to the field of medical instrument technology, more particularly to respiratory flow method of counting, device, computer equipment
And medium.
Background technology
The main function of lung ventilator is assisted respiartion, and when user's air-breathing, lung ventilator provides higher pressure of inspiration(Pi)
(IPAP) user's air-breathing is assisted, when user exhales, pressure is switched to relatively low expiratory pressure from IPAP by lung ventilator in time
Power (EPAP) so that user smooth can breathe out, while can guarantee that minimum pressure avoids the occurrence of Airway Collapse again.
In order to guarantee the respiratory movement of follow-up user in time, IPAP is just provided immediately when air-breathing just starts, is being exhaled
EPAP is provided immediately when gas just starts, lung ventilator needs to analyze and determine the respiratory air flow waveform of user in real time, but due to
It easily by external interference when flow sensor acquires, therefore will appear shake, and then influence accurate judgement.In order to reduce interference
Influence to flow, it is common practice to which using filtering, currently available technology is typically low-pass filtering, high-pass filtering, passband filter
The methods of wave.
Although above-mentioned filtering method can realize filter effect, at complex disposal process, mass data
Reason, data handling procedure occupies resource, and mostly time-consuming, causes lung ventilator can not quick response.
Invention content
Based on this, it is necessary to count filtering method complex disposal process for flow in general lung ventilator, cause lung ventilator
The problem of response speed is slow provides a kind of simple respiratory flow method of counting of data processing, device, computer equipment and Jie
Matter improves lung ventilator response speed.
A kind of respiratory flow method of counting, including step：
Respiratory flow sampled data is obtained, obtains original data on flows；
Original data on flows is written in the sliding window of default first data length, obtains single order sliding window array；
Update calculates the arithmetic mean of instantaneous value of single order sliding window array, and arithmetic mean of instantaneous value is written default second data length
In sliding window, second order sliding window array is obtained；
According to second order sliding window array, lung ventilator flow count results are obtained.
Arithmetic mean of instantaneous value is written in the sliding window of default second data length in one of the embodiments, obtains two
The step of rank sliding window array, includes：
When second order sliding window array has write full data, judge that current time corresponding arithmetic mean of instantaneous value is slided with second order
Whether absolute difference is more than preset increments threshold value between newest arithmetic mean of instantaneous value in window array；
If more than the filling of current time corresponding arithmetic mean of instantaneous value then is updated to second order sliding window array；
If being not more than, newest arithmetic mean of instantaneous value in second order sliding window array is refilled with being updated to second order sliding window
Array.
Original data on flows is written in the sliding window of default first data length, obtained in one of the embodiments,
The step of single order sliding window array, includes：
According to time sequencing, original flow data update is written in the sliding window of default first data length, obtains one
Rank sliding window array.
In one of the embodiments, according to second order sliding window array, the step of obtaining lung ventilator flow count results it
Afterwards, it further includes：
According to the second order sliding window array of real-time update, judge whether to reach lung ventilator assisted inhalation trigger condition.
In one of the embodiments, according to the second order sliding window array of real-time update, judge whether to reach lung ventilator auxiliary
The step of helping inspiration trigger condition includes：
Latest data in second order sliding window array is obtained, using latest data as slope zequin；
According to time sequencing, searched for forward from slope zequin；
When the arithmetic mean of instantaneous value searched arithmetic mean of instantaneous value corresponding less than starting point, record and mark the arithmetic searched
Average value；
Calculate matched curve slope between slope zequin and the arithmetic mean of instantaneous value searched；
The arithmetic mean of instantaneous value searched is re-used as slope zequin, is returned according to time sequencing, is calculated from slope
The step of starting point is searched for forward；
When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, judgement, which reaches, exhales
Suction machine assisted inhalation trigger condition.
In one of the embodiments, when the matched curve slope value when previous calculating is less than the matched curve of last computation
During slope, judge to include the step of reaching lung ventilator assisted inhalation trigger condition：
When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, a upper breathing is obtained
The breathing amplitude in period；
When current search to the breathing amplitude of arithmetic mean of instantaneous value and a upper respiratory cycle meet preset condition when, judgement reaches
To lung ventilator assisted inhalation trigger condition.
The formula of preset condition is in one of the embodiments,：
FF [m] ＞ X+k*A
Wherein, FF [m] is the arithmetic mean of instantaneous value that current search arrives, and flow minimum when X is exhales, A was upper one breathing week
The breathing amplitude of phase, the value range of k is 0.1~0.6.
A kind of respiratory flow counting device, including：
Acquisition module for obtaining respiratory flow sampled data, obtains original data on flows；
First processing module for original data on flows to be written in the sliding window of default first data length, obtains one
Rank sliding window array；
For updating the arithmetic mean of instantaneous value for calculating single order sliding window array, arithmetic mean of instantaneous value is written for Second processing module
In the sliding window of default second data length, second order sliding window array is obtained；
Counting module, for according to second order sliding window array, obtaining lung ventilator flow count results.
The application respiratory flow method of counting and device obtain respiratory flow sampled data, obtain original data on flows, will
In the sliding window of original default first data length of data on flows write-in, update calculates the arithmetic average of single order sliding window array
Value arithmetic mean of instantaneous value is written in the sliding window of default second data length, according to second order sliding window array, obtains lung ventilator stream
Measure count results.In whole process, without complicated filtering process, using simple mean value calculation, and pass through sliding window side
Formula accurate counting respiratory flow provides data support to improve lung ventilator response speed.
In addition, the application also provides a kind of computer readable storage medium, computer program is stored thereon with, which is located
It manages when device performs and realizes such as the step of the above method.
In addition, the application also provides a kind of computer equipment, including memory, processor and storage on a memory and can
The computer program run on a processor, processor are realized when performing program such as the step of the above method.
Above computer readable storage medium storing program for executing and computer equipment when performing above-mentioned respiratory flow method of counting, obtain
Respiratory flow sampled data obtains original data on flows, and original data on flows is written to the sliding window of default first data length
In, update calculates the arithmetic mean of instantaneous value of single order sliding window array, and arithmetic mean of instantaneous value is written to the slip of default second data length
In window, according to second order sliding window array, lung ventilator flow count results are obtained.In whole process, without complicated filtering process,
Using simple mean value calculation, and by sliding window mode accurate counting respiratory flow, carried to improve lung ventilator response speed
It is supported for data.
Description of the drawings
Fig. 1 is conventional ventilators structure diagram；
Fig. 2 is the data on flows waveform diagram that lung ventilator samples before not filtering；
Fig. 3 is the flow diagram of the application respiratory flow first embodiment of method of counting；
Fig. 4 is data on flows waveform diagram after the processing of the application respiratory flow method of counting；
Fig. 5 is the flow diagram of the application respiratory flow second embodiment of method of counting；
Fig. 6 is the structure diagram of the one of embodiment of the application respiratory flow counting device.
Specific embodiment
For be described in detail the application respiratory flow method of counting, below in conjunction with Fig. 1 for lung ventilator basic structure and
Function illustrates.Fig. 1 is lung ventilator basic structure schematic diagram, in Fig. 1 it can be seen that lung ventilator is by wind turbine, gas circuit, flow
The components such as sensor, pressure sensor, processor, display screen form.Wind turbine is air source, generates gas-pressurized, defeated via gas circuit
It is sent to outside lung ventilator, then the gas after pressurization is transported to by user terminal by external breathing pipeline.Flow sensor is adopted
Collecting the flow value at lung ventilator gas outlet, pressure sensor acquires the pressure value at lung ventilator gas outlet, and display screen shows pressure,
The various detection informations such as flow, respiratory state.Lung ventilator can also detect that user's (patient) respiratory cycle (completes one and exhales a suction
Action), the parameters such as maximum flow and flow minimum in the respiratory cycle.Since user's respiratory is the mistake of a fluctuation
Journey so directly collected data on flows is the data of a fluctuation, specifically may refer to Fig. 2, the flow quantity of fluctuation is not
Good assisted respiartion function is realized conducive to lung ventilator, and therefore, traditional mode is using complicated filtering process, filtering process
Complicated data is needed to calculate, consumption lung ventilator mass data process resource.
Based on traditional approach there are drawbacks described above, as shown in figure 3, the application provides a kind of respiratory flow method of counting, packet
Include step：
S200：Respiratory flow sampled data is obtained, obtains original data on flows.
Respiratory flow sampled data is obtained, which can be the acquisition of lung ventilator itself real-time sampling or periphery
Ancillary equipment sampling inputs later.The sampled data got is the original data on flows of a fluctuation.Specifically, it exhales
It is to be obtained based on default sampling period real-time sampling, such as conventional lung ventilator sample frequency is to inhale traffic sampling data
1000HZ, i.e. lung ventilator can acquire 1000 data points in 1 second, this 1000 data points are with noisy flow waveform.
S400：Original data on flows is written in the sliding window of default first data length, obtains single order sliding window array.
The original data on flows that step S200 is obtained is written in the sliding window of default first data length, obtains first
Sliding window array.Default first data length is preset numerical value, according to actual conditions can set,
Specifically, the first data length can be any one integer value in 2-5000.The data length of single order sliding window array is
Fixed (being equal to default first data length), and initial data is as time sequencing continuous sampling is sent to default first number
According in the sliding window of length, therefore, data are as time sequencing is newer in single order sliding window array.It is original when write-in
It when data on flows deficiency presets the first data length, is written successively according to time sequencing, when the original data on flows of write-in is more than
During default first data length, per newly one original data on flows of write-in, then the original being at first written in single order sliding window array
Beginning data are abandoned, other data move forward one successively, then will be newly written original data on flows and be put into last position, are followed successively
Ring keeps single order sliding window array length to fix, data sequential update at any time.It below will be using default first data length as 50
For, data update situation in single order sliding window array is described in detail.When initial, there is no appoint in the sliding window of 50 data lengths
What data, the collected original datas on flows of step S100 are written to the slip of 50 data lengths according to acquisition sequencing successively
In window, obtain { 1,2,3 ... ... 50 }, wherein 1 number at first sampling be written to data in sliding window, when have ground 51 originals
When beginning data on flows samples, the initial data being written at first is abandoned, that is, abandons No. 1 original data on flows, { 2,3 ... ... 50 }
One is moved forward, and single order sliding windows are written into newest 51 original datas on flows, as single order sliding window latest data,
It is also a newest data, is updated to the single order sliding window array of { 2,3 ... ... 50,51 }.
S600：Update calculates the arithmetic mean of instantaneous value of single order sliding window array, and arithmetic mean of instantaneous value is written default second data
In the sliding window of length, second order sliding window array is obtained.
Arithmetic mean of instantaneous value refers to again divided by entire data include the original flow data accumulation that entire array includes
The number of original data on flows.As above-mentioned, single order sliding window array is a newer array of dynamic, therefore, is calculating it
Arithmetic mean of instantaneous value is also the newer process of dynamic, i.e., with time sequencing, original data on flows in single order sliding window array
It is updating, corresponding arithmetic mean of instantaneous value is also in continuous calculating, continuous renewal.Default second data length is preset
Value, according to actual conditions can set, can specifically select the arbitrary integer in the range of 2-2000.It needs
, it is noted that second order sliding window array is also the array of a fixed data length, it is suitable based on the time when wherein data write full
Sequence updates arithmetic mean of instantaneous value wherein included, and specific renewal process is similar with above-mentioned single order sliding window array, no longer superfluous herein
It states.
S800：According to second order sliding window array, lung ventilator flow count results are obtained.
Based on the second order sliding window array that step S600 is obtained, most of interference is filtered out, has retained the base of original waveform
Eigen specifically may refer to Fig. 4, can obtain accurate lung ventilator flow count results.
The application respiratory flow method of counting obtains respiratory flow sampled data, original data on flows is obtained, by primary flow
In the sliding window for measuring default first data length of data write-in, update calculates the arithmetic mean of instantaneous value of single order sliding window array, will calculate
In the sliding window of default second data length of art average value write-in, according to second order sliding window array, obtain lung ventilator flow and count
As a result.In whole process, without complicated filtering process, using simple mean value calculation, and it is accurate by sliding window mode
Respiratory flow is counted, data support is provided to improve lung ventilator response speed.
Arithmetic mean of instantaneous value is written in the sliding window of default second data length in one of the embodiments, obtains two
The step of rank sliding window array, includes：
Step 1：When second order sliding window array has write full data, judge current time corresponding arithmetic mean of instantaneous value with
Whether absolute difference is more than preset increments threshold value between newest arithmetic mean of instantaneous value in second order sliding window array.
Step 2：If more than the filling of current time corresponding arithmetic mean of instantaneous value then is updated to second order sliding window array.
Step 3：If being not more than, newest arithmetic mean of instantaneous value in second order sliding window array is refilled with being updated to two
Rank sliding window array.
As above-mentioned, original data on flows and its arithmetic mean of instantaneous value are with time change in single order sliding window array
, the arithmetic mean of instantaneous value obtained over time can be more and more, the data length of second order sliding window array be it is fixed,
At a time, it when second order sliding window array has write full data, then needs to be updated the data wherein having been written into, update
Equally by the way of above-mentioned single order sliding window array, write time earliest data (arithmetic mean of instantaneous value) are abandoned, it will be new
Write-in data is needed to be written to second order sliding window array end position.In particular it is necessary to the data of write-in are corresponded to according to current time
Arithmetic mean of instantaneous value and second order sliding window array between newest arithmetic mean of instantaneous value absolute difference whether be more than preset increments
Threshold value judges, when between newest arithmetic mean of instantaneous value in current time corresponding arithmetic mean of instantaneous value and second order sliding window array
When difference is more than preset increments threshold value, show original data on flows that present sample arrives with before there are larger fluctuation needs
The data are paid close attention to, therefore, the corresponding arithmetic mean of instantaneous value of current time is filled in the manner described above and is updated to second order sliding window number
Group；When newest arithmetic mean of instantaneous value difference is less than default increase in current time corresponding arithmetic mean of instantaneous value and second order sliding window array
When measuring threshold value, show that original data on flows that present sample arrives and fluctuation before are little, active user's respiratory flow is maintained at
One metastable range (compared to preset increments threshold value), at this point, by newest arithmetic average in second order sliding window array
Value is refilled with being updated in second order sliding window array, i.e., will be newest including 2 in second order sliding window array in the updated
Arithmetic mean of instantaneous value.Specifically, preset increments threshold value is preset numerical value, can need to carry out according to actual conditions
Setting, for example, it can be set to being 2 liters/min.
Entire second order sliding window array renewal process will be described in detail so that default second data length is 10 as an example below.
First, update calculates the arithmetic mean of instantaneous value of single order sliding window array, is written to second order sliding window array successively sequentially in time,
{ 1,2 ... ..., 10 } are obtained when writing full, when the 11st arithmetic mean of instantaneous value is calculated in update, calculate the 11st arithmetic average
The absolute value of difference between value and the 10th arithmetic average, judges whether the absolute value of the difference is more than preset increments threshold value, when
More than when, abandon No. 1 arithmetic mean of instantaneous value in former { 1,2 ... ..., 10 }, other data move forward one, then by the 11st successively
Arithmetic mean of instantaneous value is filled to last position, obtains { 2,3 ... ..., 10,11 } newer second order sliding window array；When no more than when, throw
Abandon No. 1 arithmetic mean of instantaneous value in former { 1,2 ... ..., 10 }, other data move forward one successively, then by the 10th arithmetic mean of instantaneous value
It refills to last position, obtains { 2,3 ... ..., 10,10 } newer second order sliding window array.
As shown in figure 5, it is further included after step S800 in one of the embodiments,：
S900：According to the second order sliding window array of real-time update, judge whether to reach lung ventilator assisted inhalation trigger condition.
The second order sliding window array of real-time update has been able to accurate characterization active user's respiratory capacity, based on the accurate number
According to, it can precisely judge whether to reach lung ventilator assisted inhalation trigger condition, that is, judge whether to need to switch to IPAP from EPAP,
Pressure is improved, assists user's air-breathing.
In one of the embodiments, according to the second order sliding window array of real-time update, judge whether to reach lung ventilator auxiliary
The step of helping inspiration trigger condition includes：
Step 1：Latest data in second order sliding window array is obtained, using latest data as slope zequin.
Step 2：According to time sequencing, searched for forward from slope zequin.
Step 3：When the arithmetic mean of instantaneous value searched arithmetic mean of instantaneous value corresponding less than starting point, record and mark search
The arithmetic mean of instantaneous value arrived.
Step 4：Calculate matched curve slope between slope zequin and the arithmetic mean of instantaneous value searched.
Step 5：The arithmetic mean of instantaneous value searched is re-used as slope zequin, is returned according to time sequencing, from oblique
The step of rate zequin is searched for forward.
Step 6：When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, judgement
Reach lung ventilator assisted inhalation trigger condition.
Data can be arranged in order sequentially in time in second order sliding window array, and the more early data of write-in are more forward, then
Latest data is last digit according to (arithmetic mean of instantaneous value) in second order sliding window array, using last digit according to as slope zequin, one
As for, in second order interaction window array latest data be characterize present flow rate value data.According to time sequencing, from the data meter
It counts and a little searches for the arithmetic mean of instantaneous value that different moments are written in second order sliding window array one by one forward, when there are less than newest number
According to arithmetic mean of instantaneous value when, record and mark the arithmetic mean of instantaneous value searched, zequin and the arithmetic mean of instantaneous value that searches it
Between matched curve slope, matched curve slope can characterize flow value (arithmetic mean of instantaneous value) variation tendency, the arithmetic that will be searched
Average value is re-used as starting point, the step of scanning for again before repetition, when in the adjacent twice fitting slope of curve, when previous
When the matched curve slope value of calculating is less than the matched curve slope of last computation, shows that flow is in and accelerate the stage, then
Judgement reaches lung ventilator assisted inhalation trigger condition, needs adherence pressure, assists user's air-breathing.
In one of the embodiments, when the matched curve slope value when previous calculating is less than the matched curve of last computation
During slope, judge to include the step of reaching lung ventilator assisted inhalation trigger condition：
When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, a upper breathing is obtained
The breathing amplitude in period, when current search to the breathing amplitude of arithmetic mean of instantaneous value and a upper respiratory cycle meet preset condition
When, judgement reaches lung ventilator assisted inhalation trigger condition.
Other than needing to consider changes in flow rate trend, it is also necessary to consider that user breathes amplitude within the single respiratory cycle,
Further precisely determine a need for adherence pressure auxiliary user's air-breathing.Respiratory cycle refers to that user completes air-breathing and primary
It exhales corresponding complete cycle, which can accurately be captured by lung ventilator, and lung ventilator can also monitor user in the respiratory cycle
Breathing amplitude, specifically, breathing amplitude is within the respiratory cycle in second order sliding window array when maximum value and user's air-breathing
The difference of minimum value.When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, in acquisition
The breathing amplitude of one respiratory cycle, when current search to arithmetic mean of instantaneous value and a upper respiratory cycle breathing amplitude meet it is default
During condition, judgement reaches lung ventilator assisted inhalation trigger condition.Specifically, preset condition needs to meet following formula：FF[m]
＞ X+k*A, in formula, F [m] is the arithmetic mean of instantaneous value that current search arrives, and flow minimum when X is exhales, A was upper one breathing week
The breathing amplitude of phase, the value range of k is 0.1~0.6.More specifically, the value of k is preferably 0.3.
As shown in fig. 6, a kind of respiratory flow counting device, including：
Acquisition module 200 for obtaining respiratory flow sampled data, obtains original data on flows；
First processing module 400 for original data on flows to be written in the sliding window of default first data length, obtains
Single order sliding window array；
For updating the arithmetic mean of instantaneous value for calculating single order sliding window array, arithmetic mean of instantaneous value is write for Second processing module 600
In the sliding window for entering default second data length, second order sliding window array is obtained；
Counting module 800, for according to second order sliding window array, obtaining lung ventilator flow count results.
The application respiratory flow counting device, acquisition module 200 obtain respiratory flow sampled data, obtain original flow number
According to first processing module 400 original data on flows is written in the sliding window of default first data length, Second processing module
600 updates calculate the arithmetic mean of instantaneous value of single order sliding window array, and arithmetic mean of instantaneous value is written to the slip of default second data length
In window, counting module 800 obtains lung ventilator flow count results according to second order sliding window array.In whole process, without complexity
Filtering process, using simple mean value calculation, and by sliding window mode accurate counting respiratory flow, to improve lung ventilator
Response speed provides data and supports.
Second processing module 600 is additionally operable to write full data when second order sliding window array in one of the embodiments,
When, judge that difference is absolute between newest arithmetic mean of instantaneous value in current time corresponding arithmetic mean of instantaneous value and second order sliding window array
Whether value is more than preset increments threshold value；If more than the filling of current time corresponding arithmetic mean of instantaneous value then is updated to second order and is slided
Window array；If being not more than, newest arithmetic mean of instantaneous value in second order sliding window array is refilled with being updated to second order sliding window
Array.
First processing module 400 is additionally operable to according to time sequencing in one of the embodiments, by original data on flows more
In the sliding window of new default first data length of write-in, single order sliding window array is obtained.
Respiratory flow counting device further includes in one of the embodiments,：
Judgment module for the second order sliding window array according to real-time update, judges whether to reach lung ventilator assisted inhalation
Trigger condition.
Judgment module is additionally operable to obtain latest data in second order sliding window array in one of the embodiments, will be newest
Data are as slope zequin；According to time sequencing, searched for forward from slope zequin；When the arithmetic mean of instantaneous value searched
During arithmetic mean of instantaneous value corresponding less than starting point, record and mark the arithmetic mean of instantaneous value searched；Calculate slope zequin and
Matched curve slope between the arithmetic mean of instantaneous value searched；The arithmetic mean of instantaneous value searched is re-used as slope zequin,
It returns according to time sequencing, the operation searched for forward from slope zequin；When the matched curve slope value when previous calculating is small
When the matched curve slope of last computation, judgement reaches lung ventilator assisted inhalation trigger condition.
Specifically, the fitting that judgment module is less than last computation into current matched curve slope value calculated of the trade is bent
During line slope, the operation that judgement reaches lung ventilator assisted inhalation trigger condition includes：When the matched curve of previous calculating it is oblique
When rate value is less than the matched curve slope of last computation, the breathing amplitude of a upper respiratory cycle was obtained；When the calculation that current search arrives
When art average value and the breathing amplitude of a upper respiratory cycle meet preset condition, judgement reaches lung ventilator assisted inhalation triggering item
Part.
The formula of preset condition is in one of the embodiments,：
FF [m] ＞ X+k*A
Wherein, FF [m] is the arithmetic mean of instantaneous value that current search arrives, and flow minimum when X is exhales, A was upper one breathing week
The breathing amplitude of phase, the value range of k is 0.1~0.6.
For the effect for further explaining in detail the technical solution of the application respiratory flow method of counting and device and its realizing
Fruit will be illustrated below using a concrete application example.In concrete application example, following processing step is specifically included：
1st, in lung ventilator flow sensor according to the frequency collection patient of 1000HZ respiratory flow data.
2nd, the data on flows arrived according to 1000HZ frequency samplings is obtained, as original data on flows.
3rd, in order to quickly calculating as a result, first constructing the sliding window of length for 50 data.With array F [i]
It represents, wherein i=1,2,3 ... 50.When data deficiencies 50 in sliding window array, then a data on flows is often collected
Just it is put into array in sequence, when data have reached 50, as soon as a data on flows is then often collected, by the 1st
A data are abandoned, other data move forward one successively, then newest collected data on flows is put into the 50th position.Whenever
Enter after sliding window array there are one data on flows, the arithmetic mean of instantaneous value that lung ventilator will calculate entire array, calculation formula isWherein N takes 50 in this programme.
4th, after the average value for calculating sliding window array, average value is saved in the sliding window number that a length is 10
It organizes in FF [i], i=0,1,2 ... 10, the specific logic that preserves is：As the average value avgF newly calculated, than in array most
New data FF [new] is 2 liters/min big, then updates avgF into sliding window array FF, otherwise, update FF [new] to cunning
It (means in dynamic window array FF and inputs a FF [new] again to sliding window array).According to above step, then final FF is reflected
The flow waveform curve gone out is as shown in figure 4, filtered out most of interference, and remain the essential characteristic of original waveform.
5th, after the filtering, the judgement that lung ventilator will be triggered using filtered data FF, detailed process are：From FF
The last one data of array start to search for forward, when find a flow-rate ratio present flow rate it is small when, then it is smaller to record this
Flow value and its corresponding label.It is smaller value such as to find j-th of data, then records FF [j] and j, then with minimum two
Multiplication calculates the slope of the matched curve between FF [j] and FF [end] (end=10 herein), it is contemplated that filtered data
Present step-like, and real curve should be circular arc type, therefore is considered as when matched curve slope using initial data,
It is fitted using F arrays.The calculation formula of slope is after fitting
6th, again since j, numerical value more smaller than FF [j] is found out forward, such as FF [m], is repeated above step, is calculated m
To the slope k mj of the fitting formula of j, if kmj<Kjend then illustrates that flow is at this time and accelerates the stage, and by preceding
The rule of update FF mentioned in text has met triggering it is found that FF [end] at least rises 4 liters/min than FF [m] at this time
One of condition.Other than above procedure, lung ventilator also by the lasting amplitude for calculating prior respiration cycle, that is, is triggering judgement
Afterwards, lung ventilator starts to search the peak flow values of appearance in the waveform formed in FF, and records numerical value Fmax, meanwhile,
Lung ventilator searches the flow minimum occurred, and records numerical value Fmin.After Fmax and Fmin is obtained, you can calculate
Go out the amplitude A=Fmax-Fmin of a respiratory cycle, finally, if FF [m]>Minimum value+0.3* prior respiration cycles during expiration
Amplitude A, then it can be assumed that for triggering (while meet 2 conditions).
In addition, the application also provides a kind of computer readable storage medium, computer program is stored thereon with, which is located
It manages when device performs and realizes such as the step of the above method.
In addition, the application also provides a kind of computer equipment, including memory, processor and storage on a memory and can
The computer program run on a processor, processor are realized when performing program such as the step of the above method.
Above computer readable storage medium storing program for executing and computer equipment when performing above-mentioned respiratory flow method of counting, obtain
Respiratory flow sampled data obtains original data on flows, and original data on flows is written to the sliding window of default first data length
In, update calculates the arithmetic mean of instantaneous value of single order sliding window array, and arithmetic mean of instantaneous value is written to the slip of default second data length
In window, according to second order sliding window array, lung ventilator flow count results are obtained.In whole process, without complicated filtering process,
Using simple mean value calculation, and by sliding window mode accurate counting respiratory flow, carried to improve lung ventilator response speed
It is supported for data.
Above example only expresses the several embodiments of the present invention, and description is more specific and detailed, but can not
Therefore it is construed as limiting the scope of the patent.It should be pointed out that for those of ordinary skill in the art,
Under the premise of not departing from present inventive concept, various modifications and improvements can be made, these belong to protection scope of the present invention.
Therefore, the protection domain of patent of the present invention should be determined by the appended claims.
Claims (10)
1. a kind of respiratory flow method of counting, which is characterized in that including step：
Respiratory flow sampled data is obtained, obtains original data on flows；
By in the sliding window of original default first data length of data on flows write-in, single order sliding window array is obtained；
Update calculates the arithmetic mean of instantaneous value of the single order sliding window array, and it is long that the arithmetic mean of instantaneous value is written default second data
In the sliding window of degree, second order sliding window array is obtained；
According to the second order sliding window array, lung ventilator flow count results are obtained.
2. respiratory flow method of counting according to claim 1, which is characterized in that described that the arithmetic mean of instantaneous value is written
In the sliding window of default second data length, the step of obtaining second order sliding window array, includes：
When the second order sliding window array has write full data, current time corresponding arithmetic mean of instantaneous value and the second order are judged
Whether absolute difference is more than preset increments threshold value between newest arithmetic mean of instantaneous value in sliding window array；
If more than the filling of current time corresponding arithmetic mean of instantaneous value then is updated to the second order sliding window array；
If being not more than, newest arithmetic mean of instantaneous value in second order sliding window array is refilled with being updated to the second order sliding window
Array.
3. respiratory flow method of counting according to claim 1, which is characterized in that described to write the original data on flows
In the sliding window for entering default first data length, obtain single order sliding window array the step of include：
According to time sequencing, by the sliding window of original default first data length of flow data update write-in, one is obtained
Rank sliding window array.
4. respiratory flow method of counting according to claim 1, which is characterized in that described according to the second order sliding window number
After the step of group, acquisition lung ventilator flow count results, further include：
According to the second order sliding window array of real-time update, judge whether to reach lung ventilator assisted inhalation trigger condition.
5. respiratory flow method of counting according to claim 4, which is characterized in that described according to described the two of real-time update
Rank sliding window array judges whether to include the step of reaching lung ventilator assisted inhalation trigger condition：
Latest data in second order sliding window array is obtained, using the latest data as slope zequin；
According to time sequencing, searched for forward from the slope zequin；
When the arithmetic mean of instantaneous value searched arithmetic mean of instantaneous value corresponding less than starting point, record and mark the arithmetic average searched
Value；
Calculate matched curve slope between the slope zequin and the arithmetic mean of instantaneous value searched；
By described search to arithmetic mean of instantaneous value be re-used as slope zequin, return it is described according to time sequencing, from described
The step of slope zequin is searched for forward；
When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, judgement reaches lung ventilator
Assisted inhalation trigger condition.
6. respiratory flow method of counting according to claim 5, which is characterized in that described when the fitting when previous calculating is bent
When line slope value is less than the matched curve slope of last computation, judge to wrap the step of reaching lung ventilator assisted inhalation trigger condition
It includes：
When the matched curve slope value when previous calculating is less than the matched curve slope of last computation, a upper respiratory cycle was obtained
Breathing amplitude；
When current search to the breathing amplitude of arithmetic mean of instantaneous value and a upper respiratory cycle meet preset condition when, judgement reaches
To lung ventilator assisted inhalation trigger condition.
7. respiratory flow method of counting according to claim 6, which is characterized in that the formula of the preset condition is：
FF [m] ＞ X+k*A
Wherein, FF [m] is the arithmetic mean of instantaneous value that current search arrives, and flow minimum when X is exhales, A was a upper respiratory cycle
Amplitude is breathed, the value range of k is 0.1~0.6.
8. a kind of respiratory flow counting device, which is characterized in that including：
Acquisition module for obtaining respiratory flow sampled data, obtains original data on flows；
First processing module, for by the sliding window of original default first data length of data on flows write-in, obtaining one
Rank sliding window array；
Second processing module, for updating the arithmetic mean of instantaneous value for calculating the single order sliding window array, by the arithmetic mean of instantaneous value
In the sliding window of default second data length of write-in, second order sliding window array is obtained；
Counting module, for according to the second order sliding window array, obtaining lung ventilator flow count results.
9. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is held by processor
The step of any one the method in such as claim 1-7 is realized during row.
10. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage
Computer program, which is characterized in that the processor is realized when performing described program as described in any one in claim 1-7
The step of method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711488861.9A CN108245749A (en) | 2017-12-29 | 2017-12-29 | Respiratory flow method of counting, device, computer equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711488861.9A CN108245749A (en) | 2017-12-29 | 2017-12-29 | Respiratory flow method of counting, device, computer equipment and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108245749A true CN108245749A (en) | 2018-07-06 |
Family
ID=62725259
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711488861.9A Pending CN108245749A (en) | 2017-12-29 | 2017-12-29 | Respiratory flow method of counting, device, computer equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108245749A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109597574A (en) * | 2018-11-27 | 2019-04-09 | 深圳市酷开网络科技有限公司 | Distributed data storage method, server and readable storage medium storing program for executing |
| CN110840430A (en) * | 2018-08-21 | 2020-02-28 | 北京万生人和科技有限公司 | Intra-abdominal pressure data screening method, computer-readable storage medium, and intra-abdominal pressure data screening device |
| CN110947062A (en) * | 2019-03-31 | 2020-04-03 | 湖南明康中锦医疗科技发展有限公司 | Respiration withdrawal judgment method of respiration support equipment |
| CN113065979A (en) * | 2021-03-22 | 2021-07-02 | 贵州电网有限责任公司 | Load report improving and self-checking method for dispatching automation system |
| CN113546260A (en) * | 2021-06-07 | 2021-10-26 | 中国人民解放军总医院第一医学中心 | Method, apparatus and storage medium for monitoring transnasal high flow usage anomalies |
| CN113827224A (en) * | 2021-08-31 | 2021-12-24 | 真健康（北京）医疗科技有限公司 | Respiratory flow information processing method and respiratory detection device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5303700A (en) * | 1991-07-04 | 1994-04-19 | Dragerwerk Aktiengesellschaft | Method for detecting the respiratory phases of a patient during an assisted ventilating process |
| CN105854141A (en) * | 2016-05-05 | 2016-08-17 | 湖南明康中锦医疗科技发展有限公司 | Respiratory triggering conversion method and device |
| US20160250427A1 (en) * | 2011-12-07 | 2016-09-01 | Covidien Lp | Methods and systems for adaptive base flow |
-
2017
- 2017-12-29 CN CN201711488861.9A patent/CN108245749A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5303700A (en) * | 1991-07-04 | 1994-04-19 | Dragerwerk Aktiengesellschaft | Method for detecting the respiratory phases of a patient during an assisted ventilating process |
| US20160250427A1 (en) * | 2011-12-07 | 2016-09-01 | Covidien Lp | Methods and systems for adaptive base flow |
| CN105854141A (en) * | 2016-05-05 | 2016-08-17 | 湖南明康中锦医疗科技发展有限公司 | Respiratory triggering conversion method and device |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110840430A (en) * | 2018-08-21 | 2020-02-28 | 北京万生人和科技有限公司 | Intra-abdominal pressure data screening method, computer-readable storage medium, and intra-abdominal pressure data screening device |
| CN110840430B (en) * | 2018-08-21 | 2022-09-13 | 北京万生人和科技有限公司 | Intra-abdominal pressure data screening method, computer-readable storage medium, and intra-abdominal pressure data screening device |
| CN109597574A (en) * | 2018-11-27 | 2019-04-09 | 深圳市酷开网络科技有限公司 | Distributed data storage method, server and readable storage medium storing program for executing |
| CN109597574B (en) * | 2018-11-27 | 2021-09-24 | 深圳市酷开网络科技股份有限公司 | Distributed data storage method, server and readable storage medium |
| CN110947062A (en) * | 2019-03-31 | 2020-04-03 | 湖南明康中锦医疗科技发展有限公司 | Respiration withdrawal judgment method of respiration support equipment |
| CN113065979A (en) * | 2021-03-22 | 2021-07-02 | 贵州电网有限责任公司 | Load report improving and self-checking method for dispatching automation system |
| CN113546260A (en) * | 2021-06-07 | 2021-10-26 | 中国人民解放军总医院第一医学中心 | Method, apparatus and storage medium for monitoring transnasal high flow usage anomalies |
| CN113546260B (en) * | 2021-06-07 | 2023-08-22 | 中国人民解放军总医院第一医学中心 | Method, apparatus and storage medium for monitoring for use anomalies in transnasal high flow |
| CN113827224A (en) * | 2021-08-31 | 2021-12-24 | 真健康（北京）医疗科技有限公司 | Respiratory flow information processing method and respiratory detection device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108245749A (en) | Respiratory flow method of counting, device, computer equipment and medium | |
| JP2015096212A (en) | Systems and methods for determining patient effort and/or respiratory parameters in ventilation system | |
| CN103458784B (en) | Automatization's spirography map analysis and explanation | |
| CN105101873A (en) | Sampling and storage registry device for breath gas analysis | |
| JP2009528857A5 (en) | ||
| CN108125683A (en) | The gas of the exhalation of certain volume is collected and analyzed in the case where compensating respiration parameter frequency | |
| CN103974736B (en) | Automatic patient's synchronization control for invasive ventilation | |
| EP2473106B1 (en) | System and method for quantifying lung compliance in a self-ventilating subject | |
| JP6876057B2 (en) | Simultaneous estimation of respiratory mechanics and patient effort by optimizing parameters | |
| WO2016188343A1 (en) | Method and device for estimating air resistance and compliance | |
| CN105662417B (en) | A kind of control method and device based on the pressure signal feature recognition sound of snoring | |
| CN107427259B (en) | System for performing histogram analysis of time-based capnography signals and method of operation thereof | |
| WO2016185470A1 (en) | Pulmonary function test devices and methods | |
| CN103284724A (en) | Method, arrangement and computer program product for respiratory gas monitoring of ventilated patients | |
| CN105854141B (en) | A kind of respiratory triggering conversion method and device | |
| CN107029326A (en) | A kind of breath signal decision algorithm for positive airway pressure machine | |
| EP1792566A2 (en) | Method and apparatus for producing an average signal characteristic profile from cyclically recurring signals | |
| CN106840734B (en) | Method and device for evaluating following performance of respirator and noninvasive respirator | |
| CN111565780A (en) | Anesthesia machine, anesthesia output concentration monitoring method, system, equipment and storage medium | |
| CN110025864A (en) | Show CO2The air regenerating device for CPR of trend | |
| CN109876262B (en) | Automatic detection method for effusion in breathing machine pipeline based on wavelet | |
| CN104740735B (en) | Expiration tidal volume VTE computational methods | |
| WO2021135029A1 (en) | Method, system, and apparatus for measuring airway elasticity employing oct apparatus, and medium | |
| US9649051B2 (en) | Method and system for analyzing a patients respiratory activity and corresponding uses | |
| JP3527302B2 (en) | FV waveform synthesis method for FV waveform analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20210730 |
|
| AD01 | Patent right deemed abandoned |