CN104014076A - Analgesic waveform detection device and method - Google Patents
Analgesic waveform detection device and method Download PDFInfo
- Publication number
- CN104014076A CN104014076A CN201410160077.5A CN201410160077A CN104014076A CN 104014076 A CN104014076 A CN 104014076A CN 201410160077 A CN201410160077 A CN 201410160077A CN 104014076 A CN104014076 A CN 104014076A
- Authority
- CN
- China
- Prior art keywords
- waveform
- analgesia
- wave
- low frequency
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention discloses an analgesic waveform detection device which is characterized by comprising a waveform collection unit, an output unit, a waveform simulation unit and a waveform analysis unit. The invention also provides an analgesic waveform detection method. According the analgesic waveform generating process, ideal analgesic waveforms are simulated according to waveform configuration parameters in a signal simulation technology, and parameters of the analgesic waveforms are analyzed and calculated from output waveforms in a signal processing method. Thus, a user can visually observe difference between the configuration and calculation parameters as well as difference between ideal and practical waveforms, and further to determine whether a waveform generator works normally. The waveform detection device can be used as a matched self-checking program of an analgesic device to provide convenience for users and maintaining staff and greatly improve the waveform output security and validity of the analgesic device.
Description
Technical field
The present invention relates to analgesia field, particularly a kind of analgesia Wave form detector and method thereof.
Background technology
Open day is the patent of invention 201210488575.3 of 2012-11-26, discloses a kind of biofeedback type childbirth physics analgesic apparatus, has set forth the generating mode of analgesia waveform formula.This waveform has good clinical effectiveness to inhibiting pain in parturition.
This analgesia waveform is to utilize electronic devices and components and supporting Software Create to have accordingly the wave group of certain rule, and in specific implementation process, due to the build-in attribute of electronic devices and components and software, actual waveform and the ideal waveform generating has different; And that electronic devices and components itself can not ensure is just the same, there will be some beyond thought problems, finally cause output waveform inconsistent with expection.And also there is no a kind of whether consistent with expection device and method of waveform output that detects this equipment at present, in the time breaking down, be difficult for discovering more difficult definite source of failure.When occurring that analgesic effect is not good, being difficult to determine is formula problem, or individual variation problem, or software and hardware problem.While running into this class problem, clinical worker is often only sought the plant maintenance engineer of producer and is offered help.And Facilities Engineer also only has the oscillograph of utilization to detect output waveform, and utilize experience to judge that whether waveform is correct.But because this analgesia waveform formula is comparatively complicated, the time is long, checks by oscillograph, is difficult to go to judge that whether waveform is correct comprehensively, and need the comparatively personage of specialty to detect, the suitability is poor.
So, whether consistently with expection carry out detection waveform output in the urgent need to a kind of special Wave form detector at present.Can be used as on the one hand factory testing, plant maintenance instrument; Also can be used as on the other hand a kind of daily self-inspection instrument of clinical worker, whether accurate for the output of checkout facility waveform, the safety and efficacy that raising equipment uses.
Summary of the invention
Object of the present invention is just to address the above problem, and a kind of analgesia Wave form detector and method thereof are provided, and assisting users judges the accuracy of waveform output.
In order to achieve the above object, the present invention adopts following technical scheme: a kind of analgesia Wave form detector, is characterized in that: comprise waveform acquisition unit, output unit and wave simulation unit.Described waveform acquisition unit, wave simulation unit are connected with output unit respectively.
Described waveform acquisition unit is for the waveform that eases pain from analgesic apparatus collection;
Described wave simulation unit is for going out desirable waveform according to the configuration parameter simulation calculation that generates described analgesia waveform.
Described output unit is for exporting the analgesia waveform that described waveform acquisition unit collects, the waveform also going out for exporting described wave simulation unit simulation calculation.Described output unit comprises display, printer etc.
Further, the present invention also provides a kind of analgesia waveforms detection method based on described analgesia Wave form detector, it is characterized in that, comprises the following steps:
The waveform of steps A 1, the output of collection analgesic apparatus;
Steps A 2, according to the contoured configuration parameter that generates waveform described in steps A 1, simulation calculation goes out desirable waveform;
The actual analgesia waveform that steps A 3, output steps A 1 gather and the ideal waveform of steps A 2 emulation.
Further, steps A 2, it is as follows that described simulation calculation goes out desirable waveform step:
Step M1: emulation low frequency modulations ripple and carrier wave;
Step M2: utilize low frequency modulations ripple and the carrier wave of step M1 emulation to do AND operation, generate desirable output waveform.
Further, in step M1, emulation low frequency modulations ripple and carrier wave step are as follows:
M1, the sample rate of low frequency modulations ripple, carrier wave is set.
M2, according to the frequency of sample rate described in step m1 and each combination foundation waveform, counting and the numerical value of each point of each combination foundation waveform one-period is set, then generate low frequency modulations ripple according to the number of each combination foundation waveform and sequential combination.
M3, according to sample rate f s described in step m1 with carry wave frequency, counting and the numerical value of each point of carrier wave one-period is set, then according to the cycle-index of the duration calculation carrier wave of low frequency modulations ripple, generate carrier wave.
Further, the analgesia Wave form detector based on described, is characterized in that, also comprises waveform processing unit.Described waveform acquisition unit is also connected with output unit through described waveform processing unit.
Described waveform processing unit is for going out waveform configuration parameter from the actual analgesia resolve of wave shape of waveform acquisition unit collection;
Described output unit is also for exporting the waveform parameter of the contoured configuration parameter of described waveform processing element analysis calculating and the actual disposition of the described analgesia waveform of generation.
Further, the present invention also provides a kind of analgesia waveforms detection method based on described analgesia Wave form detector, it is characterized in that, comprises the following steps:
Step B1: the waveform that gathers analgesic apparatus output;
Step B2: go out waveform configuration parameter from output waveform analytical calculation;
Step B3, according to the contoured configuration parameter that generates waveform described in step B1, simulation calculation goes out desirable waveform;
The actual analgesia waveform of step B4, output step B1 collection and the ideal analgesia waveform of step B3 emulation; Export the contoured configuration parameter of step B2 calculating and the waveform parameter of actual disposition simultaneously.
Further, in step B2, go out waveform configuration parameter from output waveform analytical calculation, step is as follows:
Step X1: utilize A/D conversion that the simulation analgesia waveform transformation of collection is become to digital waveform, be designated as analgesia waveform S;
Step X2: S carries out pretreatment to analgesia waveform;
Step X3: calculate carrier wave and low frequency modulations wave parameter.
Further, step X2 carries out pretreatment to analgesia waveform S, comprises the following steps:
Step X21: obtain analgesia waveform S peak-to-peak value ± V, consider actual waveform output error, it is k*V that low and high level cut off value is set, 0<k<1;
Step X22: based on waveform S, be greater than+k*V is updated to 1, and be less than or equal to+k*V is updated to 0;
Step X23: obtain waveform S rising time point, and calculate successively the interval of adjacent rising edge, be designated as time period ordered series of numbers Δ t[1,2], Δ t[2,3] ... Δ t[n-1, n]; Δ t[n-1, n] n-1 interval to n rising edge of expression.
Further, step X3 calculates carrier wave and low frequency modulations wave parameter, comprises the following steps:
Step X31: according to the peak frequency of the minimum frequency of described analgesic apparatus carrier wave and low frequency modulations ripple, choose carrier wave and low frequency modulations separation wave period, be designated as Δ T;
Step X32: calculate carrier frequency.From time period Δ t[i-1, i] (i=2 ... n) in, choose the time period that is less than Δ T, to its averaging Δ Tp, carrier frequency f=1/ Δ Tp;
Step X33: the high level time that calculates low frequency modulations ripple.The time period that is less than Δ T of choosing from step X32, cumulative adjacent time period postscript is the new time period.The new time period that the high level time of low frequency modulations ripple is after adding up is added Δ Tp one year wave period;
Step X34: the low level time that calculates low frequency modulations ripple.From Δ t[i-1, i] (i=2 ... n) in, choose the time period that is greater than Δ T, the low level time of low frequency modulations ripple be this time period deduct again one carry wave period Δ Tp;
Step X35: calculate low frequency modulations wave parameter.According to the high and low level time of each low frequency modulations ripple of step X33, X34 calculating, calculate the frequency of each basic waveform.Calculate after the frequency of all basic waveform that are linked in sequence, calculate respectively the number of the basic waveform of adjacent and same frequency, extrapolate frequency and the number of each component basic waveform of low frequency modulations ripple.
The present invention has following advantage and effect with respect to prior art:
1, the invention provides a kind of means of the waveforms detection of easing pain, according to the generative process of analgesia waveform, can demodulate analgesia contoured configuration parameter, and simulate desirable analgesia waveform, user can observe between configuration parameter and calculating parameter accordingly more intuitively, difference between ideal waveform and actual waveform, and then judges that whether waveform generating work is normal.
2, Wave form detector of the present invention can be used as the supporting self-check program of described analgesic apparatus, facilitates user and attendant to use, and greatly improves safety and the effectiveness of described analgesic apparatus.
Brief description of the drawings
Fig. 1 is the internal element connection diagram of the analgesia Wave form detector of embodiment 1;
Fig. 2 (method: for Fig. 1) is the detecting step flow chart of the analgesia Wave form detector of embodiment 1;
Fig. 3 is the simulation waveform schematic diagram of the low frequency modulations ripple of embodiment 1;
Fig. 4 is the simulation waveform schematic diagram of the carrier wave of embodiment 1;
Fig. 5 is the simulation waveform schematic diagram of the one-level wave group of embodiment 1;
Fig. 6 is the internal element connection diagram of the analgesia Wave form detector of embodiment 2;
Fig. 7 is the detecting step flow chart of the analgesia Wave form detector of embodiment 2;
Fig. 8 is the analgesia waveform analysis waveform schematic diagram of embodiment 2;
Fig. 9 is the analgesia waveform analysis flow chart of steps of embodiment 2;
Detailed description of the invention
Below in conjunction with embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited to this.
The present invention is based on the open day patent of invention for 2012-11-26 " biofeedback type childbirth physics analgesic apparatus ", and a kind of devices and methods therefor that detects this analgesic apparatus output waveform is provided.First quote the generative process of this patent of invention analgesia waveform: rise time, high level time, fall time and low level time that basic waveform (modulation trapezoidal wave) is first set, several basic waveform are combined into basic wave group, then the carrier frequency of basic wave group is set, after modulation, generate one-level wave group, then by synthetic several fist-order wave groups secondary wave group output.
Ultimate principle of the present invention: according to the generative process of analgesia waveform, the one, utilize signal simulation technology to go out desirable analgesia waveform according to contoured configuration parameters simulation, the 2nd, utilize signal processing method to calculate the parameter of analgesia waveform from output waveform analysis meter.Analysis meter calculates analgesia contoured configuration parameter and simulates after desirable analgesia waveform, the difference of user accordingly between the parameter and the parameter that calculates by analysis of direct visual comparison actual disposition, between ideal waveform and actual waveform, and then judge that whether waveform generating work is normal.
It is to be noted, the carrier wave of selecting in patent of invention " the biofeedback type childbirth physics analgesic apparatus " embodiment that open day is 2012-11-26 is for sinusoidal wave, and selected carrier wave is square wave in the embodiment of the present invention, thus in the embodiment of the present invention analogous diagram waveform and embodiment accompanying drawing waveform different.
Embodiment 1
As the 1st embodiment of the present invention, as shown in Figure 1, a kind of analgesia Wave form detector, is characterized in that: comprise waveform acquisition unit 101, output unit 102 and wave simulation unit 103.Described waveform acquisition unit 101, wave simulation unit 103 are connected with output unit 102 respectively.
Described waveform acquisition unit 101 is for the waveform that eases pain from analgesic apparatus collection;
Described wave simulation unit 103 is for going out desirable waveform according to the configuration parameter simulation calculation that generates described analgesia waveform.
The analgesia waveform that described output unit 102 collects for exporting described waveform acquisition unit 101, the waveform also going out for exporting described wave simulation unit 103 simulation calculation.
Further, as shown in Figure 2, the present invention also provides a kind of analgesia waveforms detection method based on the Wave form detector that eases pain described in embodiment 1, it is characterized in that, comprises the following steps:
The waveform of steps A 1, the output of collection analgesic apparatus;
Steps A 2, according to the contoured configuration parameter that generates waveform described in steps A 1, simulation calculation goes out desirable waveform;
The actual analgesia waveform that steps A 3, output steps A 1 gather and the ideal waveform of steps A 2 emulation.
It should be noted that, can find out the patent of invention for 2012-11-26 " biofeedback type childbirth physics analgesic apparatus " description from open day, finally the analgesia waveform (SPA wave group) of output is to be combined to form by several fist-order waves group, one-level wave group is to form after the carrier wave by low frequency modulations wave group modulation corresponding frequencies, low frequency modulations wave group is to be combined to form by several basic waveform, and basic waveform is that the trapezoidal wave being formed by different low and high levels combines.Can find out according to above-mentioned analgesia waveform forming process, one-level wave group is the core of this analgesia waveform, it is to be formed by low frequency modulations ripple modulated carrier, and this modulated process is the Focal point and difficult point that whole analgesia waveform generates, and is also the important step that is prone to difference with analgesia waveform.Whether the emphasis of the present embodiment detects the generation of one-level wave group by described " Wave form detector " exactly normal, and then assists to judge that whether analgesic apparatus is working properly.
In the present embodiment, according to analgesic apparatus waveform generative process, generate one section of one-level wave group (as low frequency modulations wave group: 2(20Hz)+5(100Hz)+3(50Hz)+2(25Hz), carrier wave: square wave, frequency=1000Hz, above-mentioned method for expressing referring to open day be patent of invention " biofeedback type give a birth physics analgesic apparatus " 0088 section of related description of description of 2012-11-26), whether by detecting described one-level wave group, to detect the output of analgesic apparatus correct.
In the present embodiment, with low frequency modulations ripple: 2(20Hz)+5(100Hz)+3(50Hz)+2(25Hz) carrier wave: square wave, frequency=1KHz is example, simulation process is as follows:
Step M1: emulation low frequency modulations ripple and carrier wave;
Step M2: utilize low frequency modulations ripple and the carrier wave of step M1 emulation to do AND operation, generate desirable output waveform.
In step M1, emulation low frequency modulations ripple and carrier wave step are as follows:
M1, the sample rate f s of low frequency modulations ripple, carrier wave is set.
Taking the present embodiment as example, fs=5000.
M2, according to the frequency of sample rate f s described in step m1 and each combination foundation waveform, counting and the numerical value of each point of each combination foundation waveform one-period is set, then generate low frequency modulations ripple according to the number of each combination foundation waveform and sequential combination.
Taking the present embodiment as example, low frequency modulations ripple is combined by (n1=2) basic waveform 1, (n2=5) basic waveform 2, (n3=3) basic waveform 3, (n4=2) basic waveform 4, wherein, basic waveform 1 frequency f 1=20Hz, low level time be 25ms, rise time are 0, high level time is that 25ms, fall time are 0; Basic waveform 2 frequency f 2=100Hz, low level time be 5ms, rise time are 0, high level time is that 5ms, fall time are 0; Basic waveform 3 frequency f 3=50Hz, low level time be 10ms, rise time are 0, high level time is that 10ms, fall time are 0; Basic waveform 4 frequency f 4=25Hz, low level time be 20ms, rise time are 0, high level time is that 20ms, fall time are 0.
: the sampling number of the one-period T of basic waveform 1 is fs/f1=5000/20=250, to count is 125 to its low level, rise that to count be 0, high level is counted is 125, decline that to count be 0; The sampling number of the one-period T of basic waveform 2 is fs/f1=5000/100=50, and it is 25 that its low level is counted, rising is counted is 0, high level is counted is 25, decline is counted is 0; The sampling number of the one-period T of basic waveform 3 is fs/f1=5000/50=100, and it is 50 that its low level is counted, rising is counted is 0, high level is counted is 50, decline is counted is 0; The sampling number of the one-period T of basic waveform 4 is fs/f1=5000/25=200, and it is 100 that its low level is counted, rising is counted is 0, high level is counted is 100, decline is counted is 0.High level numerical value is 1, and low level numerical value is 0.Low frequency modulations ripple is made up of (n1=2) basic waveform 1, (n2=5) basic waveform 2, (n3=3) basic waveform 3, (n4=2) basic waveform 4 low frequency modulations ripple, total duration L1=n1/f1+n2/f2+n3/f3+n4/f4=2/20+5/100+3/50+2/25=0.29s of low frequency modulations ripple.
By the low frequency modulations ripple of above-mentioned emulation, simulation waveform as shown in Figure 3.
M3, according to sample rate f s described in step m1 with carry wave frequency, counting and the numerical value of each point of carrier wave one-period is set, then according to the cycle-index of the duration calculation carrier wave of low frequency modulations ripple, generate carrier wave.
Taking the present embodiment as example, carrier frequency is F=1000Hz, square wave, and total duration L2 must equate with low frequency modulations ripple.: the sampling number of the one-period T of carrier wave is fs/f1=5000/1000=5, its low level count be 2, high level count be 3, high level numerical value is 1, low level numerical value is-1; Cycle-index k=L1/ (1/1000)=0.29*1000=290 time.
By the carrier wave of above-mentioned emulation, simulation waveform as shown in Figure 4.
Above-mentioned low frequency modulations ripple and carrier wave do AND operation and generate one-level wave group.
In the present embodiment, above-mentioned low frequency modulations ripple (as shown in Figure 3) and carrier wave (as shown in Figure 4) do AND operation, obtain one-level wave group (as shown in Figure 5).
It is pointed out that in the present embodiment one-level wave group, the form of basic waveform and number of combinations are just for the implementation procedure of the method is described, irrelevant with reality use and clinical effectiveness, and consequent waveform is also just in order better to express the feature of figure.
Embodiment 2
Based on the analgesia Wave form detector described in embodiment 1, the 2nd embodiment of the present invention, as shown in Figure 6, a kind of analgesia Wave form detector, is characterized in that, also comprises waveform processing unit 601.Described waveform acquisition unit 101 is also connected with output unit 102 through described waveform processing unit 601.
Described waveform processing unit 601 goes out waveform configuration parameter for the actual analgesia resolve of wave shape gathering from waveform acquisition unit 101;
Described output unit 102 is also for exporting the waveform parameter of the contoured configuration parameter of described waveform processing element analysis calculating and the actual disposition of the described analgesia waveform of generation.
Further, the present invention also provides a kind of analgesia waveforms detection method based on the Wave form detector that eases pain described in embodiment 2, as shown in Figure 7, it is characterized in that, comprises the following steps:
Step B1: the waveform that gathers analgesic apparatus output;
Step B2: go out waveform configuration parameter from output waveform analytical calculation;
Step B3, according to the contoured configuration parameter that generates waveform described in step B1, simulation calculation goes out desirable waveform;
The actual analgesia waveform of step B4, output step B1 collection and the ideal analgesia waveform of step B3 emulation; Export the contoured configuration parameter of step B2 calculating and the waveform parameter of actual disposition simultaneously.
It should be noted that, the present embodiment 2 is on the basis of embodiment 1, by the reverse analysis to output waveform data, calculates waveform parameter value.Identical with embodiment 1, the present embodiment is also whether correct by detecting described one-level wave group if detecting the output of analgesic apparatus.
In the present embodiment, utilize equally one-level wave group (as low frequency modulations wave group: 5(5Hz)+20(100Hz)+10(2Hz)+20(10Hz), carrier wave: square wave, frequency=1000Hz) and the method that calculates waveform parameter value from output waveform is described.For better this method of explanation, the present embodiment is supposed the above-mentioned desirable one-level wave group of described analgesic apparatus output, considers the difference of Real output waveform simultaneously.
In the present embodiment, choose Fig. 5 waveform (by one-level wave group 2(20Hz)+5(100Hz)+3(50Hz)+2(25Hz), the simulation waveform of carrier wave: frequency=1KHz) in one section of waveform be designated as S, as shown in Fig. 8～9, analytical calculation goes out waveform parameter value, and step is as follows:
Step X1: utilize A/D conversion that the simulation analgesia waveform transformation of collection is become to digital waveform, be designated as analgesia waveform S;
Step X2: S carries out pretreatment to analgesia waveform;
Step X3: calculate carrier wave and low frequency modulations wave parameter.
Further, step X2 carries out pretreatment to analgesia waveform S, comprises the following steps:
Step X21: obtain analgesia waveform S peak-to-peak value ± V, consider actual waveform output error, it is k*V that low and high level cut off value is set, 0<k<1.The present embodiment k=0.4;
Step X22: based on waveform S, be greater than+k*V is updated to 1, and be less than or equal to+k*V is updated to 0;
Step X23: obtain waveform S rising time point, and calculate successively the interval of adjacent rising edge, be designated as time period ordered series of numbers Δ t[1,2], Δ t[2,3] ... Δ t[n-1, n]; Δ t[n-1, n] n-1 interval to n rising edge of expression.
Taking the present embodiment as example, as shown in Figure 8, under ideal waveform, Δ t[1,2]=Δ t[2,3]=Δ t[4,5]=Δ t[5,6]=1ms, Δ t[3,4]=4ms, Δ t[6,7]=7ms.
Further, step X3 calculates carrier wave and low frequency modulations wave parameter, comprises the following steps:
Step X31: according to the peak frequency of the minimum frequency of described analgesic apparatus carrier wave and low frequency modulations ripple, choose carrier wave and low frequency modulations separation wave period, be designated as Δ T.
Taking the present embodiment as example, carrier wave and low frequency modulations wave frequency separation are 500Hz, Δ T=2ms.
Step X32: calculate carrier frequency.From time period Δ t[i-1, i] (i=2 ... n) in, choose the time period that is less than Δ T, to its averaging Δ Tp, carrier frequency f=1/ Δ Tp;
Taking the present embodiment as example, the time period that is less than Δ T is Δ t[1,2], Δ t[2,3], Δ t[4,5], Δ t[5,6].Δ Tp=(Δ t[1,2]+Δ t[2,3]+Δ t[4,5]+Δ t[5,6])/4=1ms, carrier frequency f=1/ Δ Tp=1000Hz;
Step X33: the high level time that calculates low frequency modulations ripple.The time period that is less than Δ T of choosing from step X32, cumulative adjacent time period postscript is the new time period.The new time period that the high level time of low frequency modulations ripple is after adding up is added Δ Tp one year wave period.
Such as, the time period that is less than Δ T of choosing from step X32, there is m section a certain adjacent time period, respectively: Δ t[n-m, n-m+1], Δ t[n-m+1, n-m+2] ... Δ t[n-1, n], (n>m, m>1), according to time period sequence, cumulative adjacent time period postscript is new time period Δ t[n-m, n]=m* Δ Tp.The high level time TH[n-m of low frequency modulations ripple, n]=[n-(n-m)] * Δ Tp+ Δ Tp=(m+1) * Δ Tp.
Taking the present embodiment as example, one of them is less than Δ T and adjacent time period: Δ t[1,2], Δ t[2,3], by (Δ t[1,2]+Δ t[2,3]) as new time period Δ t[1,3]=(3-1) * Δ Tp=2* Δ Tp.The high level time TH[1 of low frequency modulations ripple, 3]=Δ t[1,3]+Δ Tp=3* Δ Tp=3ms.
Another one is less than Δ T and adjacent time period: Δ t[4,5], Δ t[5,6], by (Δ t[4,5]+Δ t[5,6]) as new time period Δ t[4,6]=(6-4) * Δ Tp=2* Δ Tp.The high level time TH[4 of low frequency modulations ripple, 6]=Δ t[4,6]+Δ Tp=3* Δ Tp=3ms.
Step X34: the low level time that calculates low frequency modulations ripple.From Δ t[i-1, i] (i=2 ... n) in, choose the time period that is greater than Δ T, the low level time of low frequency modulations ripple be this time period deduct again one carry wave period Δ Tp.
Such as, time period Δ t[k, k+1], the low level time TL[k of low frequency modulations ripple, k+1]=Δ t[k, k+1]-Δ Tp;
Taking the present embodiment as example, the time period that is less than Δ T is Δ t[3,4] and, the low level time TL[3 of low frequency modulations ripple, 4] and=Δ t[3,4] – Δ Tp=3ms;
Step X35: calculate low frequency modulations wave parameter.According to the high and low level time of each low frequency modulations ripple of step X33, X34 calculating, calculate the frequency of each basic waveform.Calculate after the frequency of all basic waveform that are linked in sequence, calculate respectively the number of the basic waveform of adjacent and same frequency, extrapolate frequency and the number of each component basic waveform of low frequency modulations ripple.
Taking the present embodiment Fig. 8 as example, the high and low level time of the each low frequency modulations ripple calculating according to step X33, X34 can cover whole waveform according to high and low level circulation, such as TH[1, and 3] high level, TH[3,4] low level, TH[4,6] high level, TH[6,7] low level etc.Each adjacent high and low level is combined into a basic waveform, calculates accordingly the frequency of each basic waveform.Basic waveform is the ingredient of low frequency modulations ripple, when calculating after the frequency of all basic waveform that are linked in sequence, just can calculate respectively the number of the basic waveform of adjacent and same frequency, just can extrapolate frequency and the number of each component basic waveform of low frequency modulations ripple.
It is pointed out that although the wave analyzing device described in the present embodiment 2 is the desirable simulation waveform based on shown in Fig. 8, the waveform of actual acquisition may be with it difference to some extent, this also affects the enforcement of wave analyzing device described in the present embodiment.On the contrary, whether accurate the wave analyzing device based on desirable simulation waveform research just, can detect waveform output, improved the output accuracy of waveform.
Waveform S1 has represented the frequency information of waveform, and after step X2～X3 Threshold segmentation, waveform S2 has not only continued to retain the frequency information of Real output waveform S1, and is that the parameter of calculating modulating wave and carrier wave is provided convenience.On the other hand, in the present embodiment, the low frequency modulations wave parameter of extrapolating according to above-mentioned X1～X9 method ignored basic waveform rise time and fall time parameter, have some systematic errors, but this does not affect the calculating to low frequency modulations ripple each component basic waveform frequency and number.Frequency and number be the core place of described analgesia waveform just, and the present embodiment is exactly the accuracy that accuracy by judging low frequency modulations ripple each component basic waveform frequency and number carrys out detection waveform output.
Main purpose of the present invention is the generative process according to waveform, simulate desirable analgesia waveform from input waveform parameter on the one hand, go out waveform configuration parameter from outfan resolve of wave shape on the one hand, they judge whether the output of analgesia waveform is correct then two groups of correction datas to be offered to user aid.Certainly, contoured configuration parameter is after software and hardware output, waveform and the desirable simulation waveform of actual output have certain gap, what the parameter value calculating from Real output waveform reverse push and reality arranged also has some gaps, and these all need, and user is follow-up further to be judged desirable and actual goodness of fit.Output unit of the present invention comprises the conventional output unit such as display, printer at present.
Above-described embodiment is preferably embodiment of the present invention; but embodiments of the present invention are not restricted to the described embodiments; other any do not deviate from change, the modification done under spirit of the present invention and principle, substitutes, combination, simplify; all should be equivalent substitute mode, within being included in protection scope of the present invention.
Claims (9)
1. an analgesia Wave form detector, is characterized in that: comprise waveform acquisition unit, output unit and wave simulation unit.Described waveform acquisition unit, wave simulation unit are connected with output unit respectively.
Described waveform acquisition unit is for the waveform that eases pain from analgesic apparatus collection;
Described wave simulation unit is for going out desirable waveform according to the configuration parameter simulation calculation that generates described analgesia waveform.
Described output unit is for exporting the analgesia waveform that described waveform acquisition unit collects, the waveform also going out for exporting described wave simulation unit simulation calculation.
2. analgesia Wave form detector according to claim 1, is characterized in that, also comprises waveform processing unit.Described waveform acquisition unit is also connected with output unit through described waveform processing unit.
Described waveform processing unit is for going out waveform configuration parameter from the actual analgesia resolve of wave shape of waveform acquisition unit collection;
Described output unit is also for exporting the waveform parameter of the contoured configuration parameter of described waveform processing element analysis calculating and the actual disposition of the described analgesia waveform of generation.
3. the analgesia waveforms detection method based on the Wave form detector that eases pain described in claim 1, is characterized in that, comprises the following steps:
The waveform of steps A 1, the output of collection analgesic apparatus;
Steps A 2, according to the contoured configuration parameter that generates waveform described in steps A 1, simulation calculation goes out desirable waveform;
The actual analgesia waveform that steps A 3, output steps A 1 gather and the ideal waveform of steps A 2 emulation.
4. based on analgesia waveforms detection method claimed in claim 3, it is characterized in that, in described steps A 2, simulation calculation goes out desirable waveform, comprises the following steps:
Step M1: emulation low frequency modulations ripple and carrier wave;
Step M2: utilize low frequency modulations ripple and the carrier wave of step M1 emulation to do AND operation, generate desirable output waveform.
5. based on analgesia waveforms detection method claimed in claim 4, it is characterized in that, emulation low frequency modulations ripple and carrier wave in described step M1, comprise the following steps:
M1, the sample rate of low frequency modulations ripple, carrier wave is set.
M2, according to the frequency of sample rate described in step m1 and each combination foundation waveform, counting and the numerical value of each point of each combination foundation waveform one-period is set, then generate low frequency modulations ripple according to the number of each combination foundation waveform and sequential combination.
M3, according to sample rate f s described in step m1 with carry wave frequency, counting and the numerical value of each point of carrier wave one-period is set, then according to the cycle-index of the duration calculation carrier wave of low frequency modulations ripple, generate carrier wave.
6. the analgesia waveforms detection method based on the Wave form detector that eases pain described in claim 2, is characterized in that, comprises the following steps:
Step B1: the waveform that gathers analgesic apparatus output;
Step B2: go out waveform configuration parameter from output waveform analytical calculation;
Step B3, according to the contoured configuration parameter that generates waveform described in step B1, simulation calculation goes out desirable waveform;
The actual analgesia waveform of step B4, output step B1 collection and the ideal analgesia waveform of step B3 emulation; Export the contoured configuration parameter of step B2 calculating and the waveform parameter of actual disposition simultaneously.
7. based on analgesia waveforms detection method claimed in claim 6, it is characterized in that, in described step B2, go out waveform parameter value from output waveform data analytical calculation, comprise the following steps:
Step X1: utilize A/D conversion that the simulation analgesia waveform transformation of collection is become to digital waveform, be designated as analgesia waveform S;
Step X2: S carries out pretreatment to analgesia waveform;
Step X3: calculate carrier wave and low frequency modulations wave parameter.
8. based on analgesia waveforms detection method claimed in claim 7, it is characterized in that, in described step X2, analgesia waveform S carried out to pretreatment, comprise the following steps:
Step X21: obtain analgesia waveform S peak-to-peak value ± V, consider actual waveform output error, it is k*V that low and high level cut off value is set, 0<k<1;
Step X22: based on waveform S, be greater than+k*V is updated to 1, and be less than or equal to+k*V is updated to 0;
Step X23: obtain waveform S rising time point, and calculate successively the interval of adjacent rising edge, be designated as time period ordered series of numbers Δ t[1,2], Δ t[2,3] ... Δ t[n-1, n]; Δ t[n-1, n] n-1 interval to n rising edge of expression.
9. based on analgesia waveforms detection method claimed in claim 7, it is characterized in that, in described step X3, calculate carrier wave and low frequency modulations wave parameter, comprise the following steps:
Step X31: according to the peak frequency of the minimum frequency of described analgesic apparatus carrier wave and low frequency modulations ripple, choose carrier wave and low frequency modulations separation wave period, be designated as Δ T;
Step X32: calculate carrier frequency.From time period Δ t[i-1, i] (i=2 ... n) in, choose the time period that is less than Δ T, to its averaging Δ Tp, carrier frequency f=1/ Δ Tp;
Step X33: the high level time that calculates low frequency modulations ripple.The time period that is less than Δ T of choosing from step X32, cumulative adjacent time period postscript is the new time period.The new time period that the high level time of low frequency modulations ripple is after adding up is added Δ Tp one year wave period;
Step X34: the low level time that calculates low frequency modulations ripple.From Δ t[i-1, i] (i=2 ... n) in, choose the time period that is greater than Δ T, the low level time of low frequency modulations ripple be this time period deduct again one carry wave period Δ Tp;
Step X35: calculate low frequency modulations wave parameter.According to the high and low level time of each low frequency modulations ripple of step X33, X34 calculating, calculate the frequency of each basic waveform.Calculate after the frequency of all basic waveform that are linked in sequence, calculate respectively the number of the basic waveform of adjacent and same frequency, extrapolate frequency and the number of each component basic waveform of low frequency modulations ripple.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410160077.5A CN104014076B (en) | 2014-04-21 | 2014-04-21 | A kind of analgesia Wave form detector and method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410160077.5A CN104014076B (en) | 2014-04-21 | 2014-04-21 | A kind of analgesia Wave form detector and method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104014076A true CN104014076A (en) | 2014-09-03 |
| CN104014076B CN104014076B (en) | 2016-02-24 |
Family
ID=51431251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410160077.5A Active CN104014076B (en) | 2014-04-21 | 2014-04-21 | A kind of analgesia Wave form detector and method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104014076B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105288837A (en) * | 2014-07-28 | 2016-02-03 | 广州三瑞医疗器械有限公司 | System and method for detecting labor analgesia waveforms |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4858142A (en) * | 1987-08-05 | 1989-08-15 | Tektronix, Inc. | Digitizer effective resolution measurement system using sinewave parameter estimation |
| US20030095099A1 (en) * | 2001-11-21 | 2003-05-22 | Azinger Frederick A. | Serial data cursor |
| CN101060445A (en) * | 2007-05-09 | 2007-10-24 | 华为技术有限公司 | Network flow simulation method and equipment, network flow test method and equipment |
| CN101493488A (en) * | 2008-01-22 | 2009-07-29 | 杭州士兰微电子股份有限公司 | Remote controller test system and method |
| CN101629995A (en) * | 2009-08-10 | 2010-01-20 | 华中电网有限公司 | Detection method of phasor measurement unit PMU |
| CN101769998A (en) * | 2010-01-18 | 2010-07-07 | 南通大学 | Switched reluctance motor flux linkage characteristic measurement method |
| CN102940934A (en) * | 2012-11-26 | 2013-02-27 | 广州三瑞医疗器械有限公司 | Biological feedback type delivery physical analgesia device |
| CN103135013A (en) * | 2012-12-28 | 2013-06-05 | 辽宁省电力有限公司 | Detecting and evaluating system of data acquisition unit |
-
2014
- 2014-04-21 CN CN201410160077.5A patent/CN104014076B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4858142A (en) * | 1987-08-05 | 1989-08-15 | Tektronix, Inc. | Digitizer effective resolution measurement system using sinewave parameter estimation |
| US20030095099A1 (en) * | 2001-11-21 | 2003-05-22 | Azinger Frederick A. | Serial data cursor |
| CN101060445A (en) * | 2007-05-09 | 2007-10-24 | 华为技术有限公司 | Network flow simulation method and equipment, network flow test method and equipment |
| CN101493488A (en) * | 2008-01-22 | 2009-07-29 | 杭州士兰微电子股份有限公司 | Remote controller test system and method |
| CN101629995A (en) * | 2009-08-10 | 2010-01-20 | 华中电网有限公司 | Detection method of phasor measurement unit PMU |
| CN101769998A (en) * | 2010-01-18 | 2010-07-07 | 南通大学 | Switched reluctance motor flux linkage characteristic measurement method |
| CN102940934A (en) * | 2012-11-26 | 2013-02-27 | 广州三瑞医疗器械有限公司 | Biological feedback type delivery physical analgesia device |
| CN103135013A (en) * | 2012-12-28 | 2013-06-05 | 辽宁省电力有限公司 | Detecting and evaluating system of data acquisition unit |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105288837A (en) * | 2014-07-28 | 2016-02-03 | 广州三瑞医疗器械有限公司 | System and method for detecting labor analgesia waveforms |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104014076B (en) | 2016-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| GB2537113A (en) | Apparatus and method for monitoring partial discharge | |
| CN102890258B (en) | Method for testing capture rate of parallel-structure digital storage oscilloscope | |
| CN103605062B (en) | Partial discharge signal trigger phase synchronous clock source | |
| DK2718673T3 (en) | A method of analysis of the power consumption of a system having a plurality of electrical devices | |
| CN103543329B (en) | A kind of high energy consumption smelting system electric energy loss measuring method | |
| KR101738076B1 (en) | Partial discharge remote diagnosis system for displaying reliability of partial discharge diagnosis | |
| CN105548739A (en) | Processing method of running state signal of arrester | |
| CN104198811B (en) | Method and device for measuring frequency of low frequency signal | |
| CN105510859A (en) | System and method for evaluating electronic transformer | |
| CN107530014A (en) | For system, method and the computerized algorithm of the change for measuring, showing and accurately detect electro physiology Evoked ptential | |
| CN105403820A (en) | On-line detection method of partial discharging signal of generator stator winding | |
| CN111934760A (en) | Signal processing transmission delay characteristic detection device and method and terminal equipment | |
| Dongfang | Dongfang Com Quantum Equations for LIGO Signal | |
| CN104014076B (en) | A kind of analgesia Wave form detector and method thereof | |
| CN108107394A (en) | Multiple-channel digital oscilloscope bandwidth flatness and consistency detecting method and its system | |
| CN104502875B (en) | Signal detection apparatus calibration method based on time and frequency parameter standard signal source mode | |
| KR20220075322A (en) | Analysis systems and methods for gear testing | |
| CN104597424B (en) | A kind of optical transformer electrical energy measurement verification platform and verification method | |
| CN104391209A (en) | Device for measuring line status | |
| CN106019042B (en) | A method of high-tension cable built-in electrical insulation defect type is identified based on characteristic spectrum | |
| CN203929911U (en) | The multi-functional analysis instrument of portable power | |
| Arandia-Perez et al. | Modeling automatic meter reading water demands as nonhomogeneous point processes | |
| CN105891637A (en) | Intelligent substation secondary device testing method and intelligent substation secondary device testing system | |
| CN103580654B (en) | Reactor period signal generator and its implementation and using method | |
| CN104156339A (en) | Method for recognizing weak periodic pulse signals through quadratic permutation entropy |
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 |