CN108042141A

CN108042141A - A kind of signal processing method and device

Info

Publication number: CN108042141A
Application number: CN201711146687.XA
Authority: CN
Inventors: 梁杰; 刘洪涛; 孟亚斌; 范欣薇; 龚梅军
Original assignee: Shenzhen H&T Intelligent Control Co Ltd
Current assignee: Shenzhen H&T Intelligent Control Co Ltd
Priority date: 2017-11-17
Filing date: 2017-11-17
Publication date: 2018-05-18
Anticipated expiration: 2037-11-17
Also published as: CN108042141B

Abstract

The embodiment of the invention discloses a kind of signal processing method and device, wherein method includes：Obtain target discrete signal, calculate the value of the fisrt feature parameter of the target discrete signal, the value of the fisrt feature parameter is used to represent the wave character of the target discrete signal, according to the value of the fisrt feature parameter, determine credit rating information of the target discrete signal as ballistocardiography signal, the credit rating information is used to represent the content of ballistocardiography signal described in the target discrete signal.Using the embodiment of the present invention, the credit rating information of Human Physiology information centre impact figure signal can be obtained, so as to which user be aided in judge the accuracy of heart activity information.

Description

A kind of signal processing method and device

Technical field

The present invention relates to field of computer technology more particularly to a kind of signal processing methods and device.

Background technology

With the raising of people's health consciousness, people are no longer satisfied with regular body for the demand in terms of medical treatment ＆ health health care Inspection is more desirable to obtain monitoring in real time, inadvertently can just get the health and fitness information of oneself.At present, noninductive hidden Covering property cardiac life sign vitro detection field, by ballistocardiography BCG (Ballistocardiogram) this signal Detection and research, it has already been proven that the physiologic information of cardiomotility can be got from BCG signals.For example, people can lead to Various sensor detection BCG signals are crossed, so as to obtain the physiologic information of human heart activity, common sensor has piezoelectric electro Cable, piezoresistive transducer, acceleration transducer, gyro sensor and piezoelectric film sensor etc., and to BCG signals at Reason and analysis can obtain the physiologic information of human heart activity.However, it is necessary to obtain the accurate heart activity information of human body Based on BCG signals best in quality, if containing much noise in BCG signals, live by the heart for handling with being obtained after analysis Dynamic information is exactly the inaccuracy of mistake.

All it is to think the signal got by various kinds of sensors at present for the monitoring of human heart activity health status In containing BCG signals obtain heart activity information to do data processing.However, in practical application, other east in addition to human body West is linked into sensor-based system, still can be considered as that human body accesses and obtains the heart activity information of corresponding mistake, and due to The reasons such as human motion cause to contain much noise in the BCG got signals, and the heart activity information that can make is inaccurate, And in addition to the medical staff of specialty, user oneself cannot judge the quality of BCG signals, so as to the heart activity information judged It is whether accurate, if to have reference value etc., this just loses the meaning of real-time body's cardiomotility health monitoring.

The content of the invention

The embodiment of the present invention provides a kind of signal processing method, can obtain Human Physiology information centre impact figure signal Credit rating information, so as to which user be aided in judge the accuracy of heart activity information.

In a first aspect, an embodiment of the present invention provides a kind of signal processing method, this method includes：

Obtain target discrete signal；

The value of the fisrt feature parameter of the target discrete signal is calculated, the value of the fisrt feature parameter is used to represent institute State the wave character of target discrete signal；

According to the value of the fisrt feature parameter, the quality of the target discrete signal as ballistocardiography signal etc. is determined Grade information, the credit rating information are used to represent the content of ballistocardiography signal described in the target discrete signal.

In a kind of possible design, it is close that the fisrt feature parameter includes approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate At least one of degree, signal amplitude variation coefficient, the approximate entropy are used to represent the information quantitative change of the target discrete signal Change degree, the frequency spectrum proportionality coefficient are used to represent that the spectrum energy of the target frequency bands of the target discrete signal accounts for the target The ratio of the frequency spectrum gross energy of discrete signal, the zero-crossing rate density are used to represent the intensive shape of fluctuation of the target discrete signal State, the signal amplitude variation coefficient are used to represent the average of peak value and the standard deviation of the peak value in the target discrete signal Ratio.

In a kind of possible design, the credit rating information includes multiple grades, the difference in the multiple grade Grade corresponds to the different content of ballistocardiography signal described in the target discrete signal, and the more high corresponding content of grade is higher；

The fisrt feature parameter includes approximate entropy, and the value according to the fisrt feature parameter determines the target Credit rating information of the discrete signal as ballistocardiography signal, including：

If the value of the approximate entropy in the first scope, determines quality of the target discrete signal as ballistocardiography signal Class information is the highest ranking in the multiple grade, and first scope is the impact of the heart described in the target discrete signal The fluctuation range of the approximate entropy of the target discrete signal when content of figure signal is more than first threshold；

If the value of the approximate entropy in the second scope, determines quality of the target discrete signal as ballistocardiography signal Class information is the lowest class in the multiple grade, and second scope is the impact of the heart described in the target discrete signal The fluctuation range of the approximate entropy of the target discrete signal when content of figure signal is less than second threshold；

Wherein, the second threshold is less than the first threshold.

In a kind of possible design, the fisrt feature parameter further includes frequency spectrum proportionality coefficient, described according to described The value of one characteristic parameter determines credit rating information of the target discrete signal as ballistocardiography signal, further includes：

If the value of the approximate entropy is in the 3rd scope, the initial grade of acquisition credit rating information, the 3rd scope The content of ballistocardiography signal described in the target discrete signal be more than the second threshold, and less than the first threshold when The fluctuation range of the approximate entropy of the target discrete signal；

According to the value of the frequency spectrum proportionality coefficient and the initial grade, the first alternative grade is determined, wherein, if the frequency The value for composing proportionality coefficient is more than the 3rd threshold value, and the first alternative grade is higher than the initial grade, if the frequency spectrum ratio system Several values is less than or equal to the 3rd threshold value, and the first alternative grade is identical with the initial grade；

According to the described first alternative grade, determine that the target discrete signal is believed as the credit rating of ballistocardiography signal Breath.

In a kind of possible design, the fisrt feature parameter further includes zero-crossing rate density, described according to described first Alternative grade determines credit rating information of the target discrete signal as ballistocardiography signal, including：

According to the value of the zero-crossing rate density and the first alternative grade, the second alternative grade is determined, wherein, if described The value of zero-crossing rate density is less than or equal to the 4th threshold value, and the second alternative grade is higher than the described first alternative grade, if institute The value for stating zero-crossing rate density is more than the 4th threshold value, and the second alternative grade is identical with the described first alternative grade, alternatively, The second alternative grade is less than the described first alternative grade；

According to the described second alternative grade, determine that the target discrete signal is believed as the credit rating of ballistocardiography signal Breath.

In a kind of possible design, the fisrt feature parameter further includes signal amplitude variation coefficient, described according to institute The second alternative grade is stated, determines credit rating information of the target discrete signal as ballistocardiography signal, including：

According to the value of the signal amplitude variation coefficient and the second alternative grade, the 3rd alternative grade is determined, wherein, If the value of the signal amplitude variation coefficient is less than or equal to the 5th threshold value, the 3rd alternative grade is standby higher than described second Grade is selected, if the value of the signal amplitude variation coefficient is more than the 5th threshold value, the 3rd alternative grade and described second Alternative grade is identical, alternatively, the 3rd alternative grade is less than the described second alternative grade；

The 3rd alternative grade is determined as the target discrete signal as the credit rating of ballistocardiography signal to believe Breath.

Second aspect, an embodiment of the present invention provides a kind of signal processing apparatus, which includes：

Acquisition module, for obtaining the target discrete signal containing wave characteristic；

Computing module, for calculating the value of the fisrt feature parameter of the target discrete signal, the fisrt feature parameter Value be used to represent the wave character of the target discrete signal；

Determining module for the value according to the fisrt feature parameter, determines that the target discrete signal is impacted as the heart The credit rating information of figure signal, the credit rating information are used to represent that ballistocardiography described in the target discrete signal is believed Number content.

The fisrt feature parameter includes approximate entropy, and the determining module is specifically used for：

If the value of the approximate entropy in the second scope, determines quality of the target discrete signal as ballistocardiography signal Class information is the lowest class in the multiple grade, and second scope is ballistocardiography described in the target discrete signal The fluctuation range of the approximate entropy of the target discrete signal when content of signal is less than second threshold；

Wherein, the second threshold is less than the first threshold.

In a kind of possible design, the fisrt feature parameter further includes frequency spectrum proportionality coefficient, and the determining module is also It is specifically used for：

In a kind of possible design, the fisrt feature parameter further includes zero-crossing rate density, and the determining module is specific For：

In a kind of possible design, the fisrt feature parameter further includes signal amplitude variation coefficient, the definite mould Block is specifically used for：

The third aspect, an embodiment of the present invention provides a kind of electronic equipment, which includes processor, input is set Standby, output equipment and memory, the processor, input equipment, output equipment and memory are connected with each other, wherein, it is described to deposit For reservoir for storing the computer program that terminal is supported to perform the above method, the computer program includes program instruction, described Processor is arranged to call described program instruction, the method for performing above-mentioned first aspect.

Fourth aspect, an embodiment of the present invention provides a kind of computer readable storage medium, the computer storage media Computer program is stored with, the computer program includes program instruction, and described program instruction makes institute when being executed by a processor State the method that processor performs above-mentioned first aspect.

The embodiment of the present invention calculates the fisrt feature parameter of the target discrete signal by obtaining target discrete signal Value, the value of the fisrt feature parameter are used to represent the wave character of the target discrete signal, according to the value of the fisrt feature parameter, Determine credit rating information of the target discrete signal as ballistocardiography signal, which is used to represent the target The content of the ballistocardiography signal in discrete signal can obtain the credit rating letter of Human Physiology information centre impact figure signal Breath, so as to which user be aided in judge the accuracy of heart activity information.

Description of the drawings

Technical solution in order to illustrate the embodiments of the present invention more clearly, below will be to needed in embodiment description Attached drawing is briefly described, it should be apparent that, the accompanying drawings in the following description is some embodiments of the present invention, general for this field For logical technical staff, without creative efforts, other attached drawings are can also be obtained according to these attached drawings.

Fig. 1 is a kind of schematic flow diagram of signal processing method provided in an embodiment of the present invention；

Fig. 2 is a kind of level of signal schematic diagram corresponding with signal quality；

Fig. 3 is the schematic flow diagram of another signal processing method provided in an embodiment of the present invention；

Fig. 4 is a kind of schematic block diagram of signal processing apparatus provided in an embodiment of the present invention；

Fig. 5 is the schematic block diagram of a kind of electronic equipment provided in an embodiment of the present invention.

Specific embodiment

Below in conjunction with the attached drawing in the embodiment of the present invention, the technical solution in the embodiment of the present invention is carried out clear, complete Site preparation describes, it is clear that described embodiment is part of the embodiment of the present invention, instead of all the embodiments.Based on this hair Embodiment in bright, the every other implementation that those of ordinary skill in the art are obtained without making creative work Example, belongs to the scope of protection of the invention.

It should be appreciated that the term " first ", " second " " in description and claims of this specification and the attached drawing Three " etc. be for distinguishing different objects rather than for describing particular order.In addition, term " comprising " and " having " and it Any deformation, it is intended that cover non-exclusive include.Such as it contains the process of series of steps or unit, method, be The step of system, product or equipment are not limited to list or unit, but optionally further include the step of not listing or list Member is optionally further included for the intrinsic other steps of these processes, method, product or equipment or unit.

It is also understood that referenced herein " embodiment " it is meant that describe in conjunction with the embodiments special characteristic, structure or Characteristic may be embodied at least one embodiment of the present invention.Each position in the description shows that the phrase might not Identical embodiment is each meant, nor the independent or alternative embodiment with other embodiments mutual exclusion.People in the art Member explicitly and implicitly understands that embodiment described herein can be combined with other embodiments.

It will be further appreciated that the term "and/or" used in description of the invention and the appended claims is Refer to any combinations and all possible combinations of one or more of the associated item listed, and including these combinations.

Below in conjunction with attached drawing 1- attached drawings 5, signal processing method provided in an embodiment of the present invention and device are carried out detailed It introduces.

It is that the embodiment of the present invention provides a kind of schematic flow diagram of signal processing method, in the embodiment of the present invention referring to Fig. 1 Method can realize including mobile phone, tablet computer, desktop computer, IPAD etc. have signal processing function terminal on, such as Shown in figure, which may include：

S101 obtains target discrete signal；

In embodiments of the present invention, terminal can obtain the target discrete signal containing wave characteristic, target discrete letter Number can include at least two target discrete data, the wave characteristic can be used to indicate that the target discrete signal waveform exist It fluctuates amplitude of the waveform of the either target discrete signal for straight line or the target discrete signal and is not all zero etc..It is optional , which can be filtered treated target discrete signal, the target after the filtering process from Scattered signal can include at least two target discrete data, such as filtering process can be that cutoff frequency is having for 1 hertz (Hz) Limit for length's unit impulse response (Finite Impulse Response, FIR) high-pass filter and/or cutoff frequency are 40Hz's FIR low pass filter, the FIR high-pass filters of 1Hz are mainly used for filtering out the baseline drift of below 1Hz in the target discrete signal It moves, the FIR low pass filter of 40Hz is mainly used for filtering out in the target discrete signal 40Hz with noise of super band etc..

Optionally, terminal can gather the object time by piezoresistive transducer, piezoelectric cable, piezoelectric film sensor etc. In the range of the initial data for including at least two discrete datas, which may range from default time range, such as 1 Minute, and the undulating value of the initial data is detected, which is used to represent the fluctuation of the initial data, if the undulating value is big In or equal to targets threshold, the initial data is determined there are fluctuation, and the initial data is filtered, and is contained There is the target discrete signal of wave characteristic；If the undulating value is less than targets threshold, the initial data is determined without fluctuation, and is determined No signal or signal waveform are straight line in the initial data.For example, the targets threshold can be 0.1, if the original number detected According to undulating value be more than or equal to 0.1, determine the initial data there are fluctuation, and the initial data be filtered Obtain the target discrete signal containing wave characteristic；If the undulating value of the initial data detected is less than 0.1, the original number is determined According to no fluctuation, and it can represent that no signal or signal waveform are straight line in the initial data with forms such as number or letters. Wherein, the undulating value for detecting the initial data can be by detecting the difference of two neighboring discrete data in the initial data come real Existing, for example, the common n data of the initial data, are represented with x (n), n is the natural number more than or equal to 2, makes x₁(n)=x (n)-x (n-1), x is made₂(n)=x₁(n+1)-x₁(n), undulating value Fluctuationrate：

Wherein, numberof (x₁(n)=0 x) is represented₁(n)=0 number, n-1 represent x₁(n) total number, numberof(x₂(n)=0 x) is represented₂(n)=0 number, n-2 represent x₂(n) total number.

S102, calculates the value of the fisrt feature parameter of the target discrete signal, and the value of the fisrt feature parameter is used for Represent the wave character of the target discrete signal；

In embodiments of the present invention, terminal can calculate the value of the fisrt feature parameter of above-mentioned target discrete signal, this The value of one characteristic parameter is used to represent the wave character of the target discrete signal.Wherein, which can include near Like at least one of entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, signal amplitude variation coefficient.The approximate entropy can be used for table Show the information content variation degree of the target discrete signal；The frequency spectrum proportionality coefficient can be used to indicate that the mesh of the target discrete signal Mark frequency range spectrum energy account for the target discrete signal frequency spectrum gross energy ratio, for example, target frequency bands can be 1Hz~ 20Hz (spectrum energy of ballistocardiography signal is concentrated mainly on 1Hz~20Hz), the frequency spectrum proportionality coefficient can be expressed as the mesh Mark the ratio of the frequency spectrum gross energy of the spectrum energy and the target discrete signal between 1Hz~20Hz of discrete signal；The zero-crossing rate Density is used to represent the fluctuation high density state of the target discrete signal；The signal amplitude variation coefficient is used to represent the target discrete The ratio of the standard deviation of the average and peak value of peak value in signal.For example, if fisrt feature parameter is approximate entropy, terminal can be counted Count stating the value of the approximate entropy of target discrete signal in；If fisrt feature parameter is approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate are close Degree, it is close that terminal can calculate the value of the approximate entropy of above-mentioned target discrete signal, the value of frequency spectrum proportionality coefficient and zero-crossing rate respectively The value of degree.

Optionally, if fisrt feature parameter is frequency spectrum proportionality coefficient, terminal can carry out above-mentioned target discrete signal fast Fast Fourier transformation obtains the frequency spectrum of the target discrete signal, calculates the spectrum energy of the target frequency bands of the target discrete signal With the ratio of the frequency spectrum gross energy of the target discrete signal, which is determined as to the value of frequency spectrum proportionality coefficient.For example, terminal-pair Above-mentioned target discrete signal carries out Fast Fourier Transform (FFT), obtains the frequency spectrum of the target discrete signal, frequency spectrum proportionality coefficient be Spectrum energy between 1Hz~20Hz accounts for the ratio of frequency spectrum gross energy, can be represented with PSDr：

Wherein, PSD (f>1Hz&f<20Hz) represent the spectrum energy between 1Hz~20Hz of the target discrete signal, PSD Represent the frequency spectrum gross energy of the target discrete signal.

Optionally, if fisrt feature parameter is zero-crossing rate density, terminal can calculate the zero passage of above-mentioned target discrete signal The value of rate density.Such as zero-crossing rate density can be represented with ZCRD：

Wherein, x_n(m) above-mentioned target discrete signal is represented, N represents target discrete data in above-mentioned target discrete signal Number, N are the natural number more than or equal to 2,

Optionally, if fisrt feature parameter is signal amplitude variation coefficient, terminal can calculate above-mentioned target discrete signal Signal amplitude variation coefficient value, which can be also used for the peak value for reflecting above-mentioned target discrete signal Variation degree.For example, signal amplitude variation coefficient can be represented with following formula：

Wherein, CV represents signal amplitude variation coefficient, and STD (Peaks) represents the average of above-mentioned target discrete signal peak, Mean (Peaks) represents the standard deviation of above-mentioned target discrete signal peak.

S103 according to the value of the fisrt feature parameter, determines the target discrete signal as ballistocardiography signal Credit rating information, the credit rating information are used to represent containing for ballistocardiography signal described in the target discrete signal Amount.

In embodiments of the present invention, it is close can to include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate for above-mentioned fisrt feature parameter Degree, at least one of signal amplitude variation coefficient, terminal can be according to the values of the above-mentioned fisrt feature parameter being calculated, really Fixed credit rating information of the above-mentioned target discrete signal as ballistocardiography signal, passes through the credit rating of the ballistocardiography signal Information, it can be determined that the accuracy of heart activity information.For example, if above-mentioned fisrt feature parameter is approximate entropy, terminal is according to upper The value for the approximate entropy being calculated is stated, determines credit rating information of the above-mentioned target discrete signal as ballistocardiography signal；If Above-mentioned fisrt feature parameter is approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, and terminal is according to the above-mentioned approximate entropy being calculated Value, the value of the value of frequency spectrum proportionality coefficient and zero-crossing rate density, determine above-mentioned target discrete signal as ballistocardiography signal Credit rating information.

Wherein, which can be used to indicate that the content of the ballistocardiography signal in the target discrete signal. The credit rating information can include multiple grades, and the plurality of grade can use number, letter, symbol or number and the group of letter The expressions such as conjunction, such as grade 1, grade 2 etc. or grade A, grade B etc. or grade A1, grade B2 etc., the present invention Embodiment does not limit.Different brackets in the plurality of grade corresponds to the difference of the ballistocardiography signal in the target discrete signal Content, and the more high corresponding content of grade is higher, the more high corresponding quality of content is better, therefore the more high corresponding target discrete of grade Signal is also better as the quality of ballistocardiography signal.As shown in Fig. 2, Fig. 2 is that a kind of level of signal is corresponding with signal quality Schematic diagram, grade 1 can represent that no signal or signal are straight line, and credit rating information can be divided into 5 grades, use grade respectively 2,3,4,5,6 represent, and the bigger expression higher grade of number, and the more high corresponding target discrete signal of grade is believed as ballistocardiography Number quality it is also better, grade 2 can represent that the target discrete signal for periodic vibration signal, has strong industrial frequency noise, i.e. the mesh Mark poor quality of the discrete signal as ballistocardiography signal；Grade 3 can represent the target discrete signal as ballistocardiography Signal it is of poor quality；Class 4 can represent the target discrete signal as the second-rate of ballistocardiography signal；Class 5 can be with Represent that the target discrete signal is general as the quality of ballistocardiography signal；Class 6 can represent the target discrete signal conduct The quality of ballistocardiography signal is fine.

The embodiment of the present invention calculates the of the target discrete signal by obtaining the target discrete signal containing wave characteristic The value of one characteristic parameter, the value of the fisrt feature parameter are used to represent the wave character of the target discrete signal, according to this first The value of characteristic parameter determines credit rating information of the target discrete signal as ballistocardiography signal, the credit rating information For representing the content of the ballistocardiography signal in the target discrete signal, Human Physiology information centre impact figure letter can be obtained Number credit rating information, so as to which user be aided in judge the accuracy of heart activity information.

Fig. 3 is refer to, is the schematic flow diagram of another signal processing method provided in an embodiment of the present invention, the present invention is real Applying the method in example can realize in the terminal including mobile phone, tablet computer, desktop computer, IPAD etc. with signal processing function On, as shown in the figure, the signal processing method may include：

S301 obtains target discrete signal；

S302, calculates the value of the fisrt feature parameter of the target discrete signal, and the value of the fisrt feature parameter is used for Represent the wave character of the target discrete signal；

Step S301- steps of embodiment of the present invention S302 refer to the embodiment step S101- step S102 of Fig. 1, herein It repeats no more.

S303, if the value of the approximate entropy in the range of first, determines that the target discrete signal is believed as ballistocardiography Number credit rating information be the highest ranking in the multiple grade, first scope is institute in the target discrete signal The fluctuation range of the approximate entropy of target discrete signal when stating the content of ballistocardiography signal more than first threshold；

In embodiments of the present invention, it is close can to include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate for above-mentioned fisrt feature parameter At least one of degree, signal amplitude variation coefficient.If above-mentioned fisrt feature parameter includes approximate entropy, terminal can calculate the mesh The value of the approximate entropy of discrete signal is marked, if the value of the approximate entropy in the first scope, determines the target discrete signal as heart impact The credit rating information of figure signal is the highest ranking in the plurality of grade, which corresponds to the target discrete signal conduct Ballistocardiography signal it is best in quality.This first may range from the ballistocardiography signal in the target discrete signal content it is big The fluctuation range of the approximate entropy of target discrete signal when first threshold.For example, first threshold is 0.8, which can Think the fluctuation of the approximate entropy of the target discrete signal when content of the ballistocardiography signal in the target discrete signal is more than 0.8 Scope, the fluctuation range if more than the approximate entropy of target discrete signal when 0.8 are 0.01~0.15, then the first scope is 0.01 ~0.15.

Wherein, the credit rating information can include multiple grades, the plurality of grade can use number, letter, symbol or The expressions such as the combination of number and letter, such as grade 1, grade 2 etc. or grade A, grade B etc. or grade A1, etc. Grade B2 etc., the embodiment of the present invention does not limit.Different brackets in the plurality of grade corresponds to the heart in the target discrete signal The different content of figure signal is impacted, and the more high corresponding content of grade is higher, the more high corresponding quality of content is better, therefore grade is got over High corresponding target discrete signal is also better as the quality of ballistocardiography signal.It is closed so that the grade in Fig. 2 is corresponding with quality Exemplified by system, credit rating information can be divided into 5 grades, be represented respectively with grade 2,3,4,5,6, and the bigger expression grade of number Higher, the more high corresponding target discrete signal of grade is also better as the quality of ballistocardiography signal, the first scope for 0.01~ 0.15, terminal calculates the value of the approximate entropy of the target discrete signal as 0.1, the value 0.1 of the approximate entropy the first scope 0.01~ Between 0.15, it is the highest in the plurality of grade to determine the target discrete signal as the credit rating information of ballistocardiography signal Class 6, class 6 can represent that the target discrete signal is fine as the quality of ballistocardiography signal.

S304, if the value of the approximate entropy in the second scope, determines the target discrete signal as ballistocardiography signal Credit rating information be the lowest class in the multiple grade, second scope is described in the target discrete signal The fluctuation range of the approximate entropy of the target discrete signal when content of ballistocardiography signal is less than second threshold；

In embodiments of the present invention, it is close can to include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate for above-mentioned fisrt feature parameter At least one of degree, signal amplitude variation coefficient.If above-mentioned fisrt feature parameter includes approximate entropy, terminal can calculate the mesh The value of the approximate entropy of discrete signal is marked, if the value of the approximate entropy in the second scope, determines the target discrete signal as heart impact The credit rating information of figure signal is the lowest class in the plurality of grade, which corresponds to the target discrete signal conduct The quality of ballistocardiography signal is worst.This second may range from the ballistocardiography signal in the target discrete signal content it is small The fluctuation range of the approximate entropy of target discrete signal when second threshold, the second threshold are less than above-mentioned first threshold.For example, First threshold is 0.8, and second threshold can be 0.2, this second may range from the target discrete signal ballistocardiography and believe Number content approximate entropy of target discrete signal when being less than 0.2 fluctuation range, if target discrete signal during less than 0.2 The fluctuation range of approximate entropy is more than 0.18, then first may range from being more than 0.18.

Wherein, the credit rating information can include multiple grades, the plurality of grade can use number, letter, symbol or The expressions such as the combination of number and letter, such as grade 1, grade 2 etc. or grade A, grade B etc. or grade A1, etc. Grade B2 etc., the embodiment of the present invention does not limit.Different brackets in the plurality of grade corresponds to the heart in the target discrete signal The different content of figure signal is impacted, and the more high corresponding content of grade is higher, the more high corresponding quality of content is better, therefore grade is got over High corresponding target discrete signal is also better as the quality of ballistocardiography signal.It is closed so that the grade in Fig. 2 is corresponding with quality Exemplified by system, credit rating information can be divided into 5 grades, be represented respectively with grade 2,3,4,5,6, and the bigger expression grade of number Higher, the more high corresponding target discrete signal of grade is also better as the quality of ballistocardiography signal, the second scope be more than 0.18, the value that terminal calculates the approximate entropy of the target discrete signal is 0.2, and the value 0.2 of the approximate entropy is in the scope more than 0.18 Interior, it is the lowest class 2 in the plurality of grade to determine the target discrete signal as the credit rating information of ballistocardiography signal, Grade 2 can represent that the target discrete signal for periodic vibration signal, has strong industrial frequency noise, i.e., the target discrete signal is as the heart Impact the poor quality of figure signal.

S305, if the value of the approximate entropy is in the 3rd scope, the initial grade of acquisition credit rating information, the 3rd model The content enclosed for ballistocardiography signal described in the target discrete signal is more than the second threshold, and less than first threshold The fluctuation range of the approximate entropy of target discrete signal during value；

S306 according to the value of the frequency spectrum proportionality coefficient and the initial grade, determines the first alternative grade；

In embodiments of the present invention, it is close can to include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate for above-mentioned fisrt feature parameter At least one of degree, signal amplitude variation coefficient.If above-mentioned fisrt feature parameter includes approximate entropy and frequency spectrum proportionality coefficient, eventually End can calculate the value of approximate entropy of the target discrete signal and the value of frequency spectrum proportionality coefficient, if the value of the approximate entropy is the 3rd Scope, the terminal can obtain the initial grade of credit rating information, and the 3rd scope is that the heart rushes in the target discrete signal The fluctuation model of the approximate entropy of target discrete signal when hitting the content of figure signal more than the second threshold, and being less than the first threshold It encloses, such as first threshold is 0.8, second threshold 0.2, the 3rd scope is the ballistocardiography signal in the target discrete signal Content be more than 0.2, and during less than 0.8 the approximate entropy of the target discrete signal fluctuation range, such as the 3rd scope for 0.15~ 0.18.The terminal can determine the first alternative grade according to the value and the initial grade of the frequency spectrum proportionality coefficient.It wherein, if should The value of frequency spectrum proportionality coefficient is more than the 3rd threshold value, which is higher than the initial grade, if the frequency spectrum proportionality coefficient Value is less than or equal to the 3rd threshold value, and the first alternative grade is identical with the initial grade.

Wherein, the credit rating information can include multiple grades, the plurality of grade can use number, letter, symbol or The expressions such as the combination of number and letter, such as grade 1, grade 2 etc. or grade A, grade B etc. or grade A1, etc. Grade B2 etc., the embodiment of the present invention does not limit.Different brackets in the plurality of grade corresponds to the heart in the target discrete signal The different content of figure signal is impacted, and the more high corresponding content of grade is higher, the more high corresponding target discrete signal conduct of grade The quality of ballistocardiography signal is also better.By taking the grade in Fig. 2 and the correspondence of quality as an example, the bigger expression grade of number is more It is high, it is assumed that grade can represent that, if score=2 can represent grade 2, score=3 can represent grade with the value of score 3 etc., the 3rd threshold value can be 0.6, initial grade score=3, if the value of frequency spectrum proportionality coefficient is more than 0.6, first is alternative etc. Grade can be that initial grade add up the first incremental change, and such as the first incremental change is 1, i.e., the first alternative grade score=score+1 =3+1=4 is higher than initial grade score=3；If the value of the frequency spectrum proportionality coefficient is less than or equal to 0.6, this is first alternative etc. Grade, i.e., first alternative grade score=score=3 identical with the initial grade.

S307 according to the value of the zero-crossing rate density and the first alternative grade, determines the second alternative grade；

In embodiments of the present invention, if above-mentioned fisrt feature parameter includes approximate entropy, frequency spectrum proportionality coefficient and zero-crossing rate Density, terminal can calculate the value of the approximate entropy of the target discrete signal, the value of frequency spectrum proportionality coefficient and zero-crossing rate density Value, the value and above-mentioned first alternative grade that terminal can be according to the zero-crossing rate density determine the second alternative grade, wherein, if should The value of zero-crossing rate density is less than or equal to the 4th threshold value, which is higher than the first alternative grade, if the zero passage The value of rate density be more than the 4th threshold value, the second alternative grade is identical with the first alternative grade, alternatively, this second alternatively etc. Grade is less than the first alternative grade.By taking the grade in Fig. 2 and the correspondence of quality as an example, number is bigger, and expression higher grade, Assuming that grade can be represented with the value of score, if score=2 can represent grade 2, score=3 can represent grade 3 etc. It is 4, the first alternative grade score=4 Deng, the 4th threshold value, if the value of zero-crossing rate density is less than or equal to 4, second is alternative etc. Grade can be that the first alternative grade add up the second incremental change, and such as the second incremental change is 1, i.e., the second alternative grade score= Score+1=4+1=5 is higher than the first alternative grade score=4；If the value of zero-crossing rate density is more than 4, the second alternative grade can With identical with the first alternative grade, i.e., the second alternative grade score=score, (score≤3), alternatively, second is alternative etc. Grade can be first alternative the first decrement of grade regressive, such as the first decrement is 1, i.e., the second alternative grade score= score-1,(score>3), because the first alternative grade score=4>3, therefore the second alternative grade score=score-1=4- 1=3.

S308 according to the value of the signal amplitude variation coefficient and the second alternative grade, determines the 3rd alternative grade；

The 3rd alternative grade is determined as the target discrete signal as quality of ballistocardiography signal etc. by S309 Grade information.

In embodiments of the present invention, if above-mentioned fisrt feature parameter include approximate entropy, frequency spectrum ratio system, zero-crossing rate density with And signal amplitude variation coefficient, terminal can calculate the value of the approximate entropy of the target discrete signal, the value of frequency spectrum proportionality coefficient, mistake The value of zero rate density and the value of signal amplitude variation coefficient, terminal can be according to the values of the signal amplitude variation coefficient and above-mentioned Second alternative grade determines the 3rd alternative grade, and the 3rd alternative grade is determined as the target discrete signal and is rushed as the heart The credit rating information of figure signal is hit, so as to the credit rating information more refined, raising judges heart activity information Accuracy.Wherein, if the value of the signal amplitude variation coefficient is less than or equal to the 5th threshold value, the 3rd alternative grade is higher than The second alternative grade, if the value of the signal amplitude variation coefficient be more than the 5th threshold value, the 3rd alternative grade with this second Alternative grade is identical, alternatively, the 3rd alternative grade is less than the second alternative grade.It is corresponding with quality with the grade in Fig. 2 Exemplified by relation, number is bigger, and expression higher grade, it is assumed that grade can be represented with the value of score, as score=2 can be with table Showing grade 2, score=3 can represent grade 3 etc., and the 5th threshold value is 0.5, and the second alternative grade is grade score=5, if The value of signal amplitude variation coefficient is less than or equal to 0.5, and the 3rd alternative grade can be that the second alternative grade add up and the 3rd pass Increment, such as the 3rd incremental change are 1, i.e., the 3rd alternative grade score=score+1=5+1=6 is higher than the second alternative grade 3rd alternative grade score=6 is determined as credit rating of the target discrete signal as ballistocardiography signal by score=5 Information；If the value of signal amplitude variation coefficient is more than 0.5, the 3rd alternative grade can be identical with the second alternative grade, i.e., Score=score, (score≤3), alternatively, the 3rd alternative grade can be second alternative the second decrement of grade regressive, such as Second decrement is 1 etc., i.e. score=score-1, (score>3), due to the second alternative grade score=5>3, the 3rd is standby Grade score=score-1=5-1=4 is selected, the 3rd alternative grade score=4 is determined as the target discrete signal as the heart Impact the credit rating information of figure signal.

The embodiment of the present invention calculates the of the target discrete signal by obtaining the target discrete signal containing wave characteristic The value of one characteristic parameter, the fisrt feature parameter include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, signal amplitude variation At least one of coefficient, if the value of the approximate entropy in the first scope, determines the target discrete signal as ballistocardiography signal Credit rating information be the highest ranking in the plurality of grade, if the value of the approximate entropy in the second scope, determine the target from Scattered signal is the lowest class in the plurality of grade as the credit rating information of ballistocardiography signal, if the value of the approximate entropy exists 3rd scope obtains the initial grade of credit rating information, according to the value of the frequency spectrum proportionality coefficient and the initial grade, determines the One alternative grade according to the value of the zero-crossing rate density and the first alternative grade, determines the second alternative grade, according to the signal width It is worth the value and the second alternative grade of variation coefficient, determines the 3rd alternative grade, the 3rd alternative grade is determined as the target Credit rating information of the discrete signal as ballistocardiography signal can obtain more specific Human Physiology information centre impact figure The credit rating information of signal, so as to more effectively judge the accuracy of heart activity information.

Fig. 4 is refer to, is a kind of schematic block diagram of signal processing apparatus provided in an embodiment of the present invention, as shown in the figure, The signal processing apparatus of the embodiment of the present invention includes：

Acquisition module 10, for obtaining target discrete signal；

Specifically, the acquisition module 10 of terminal can obtain the target discrete signal containing wave characteristic, the target discrete Signal can include at least two target discrete data, which can be used to indicate that the waveform of the target discrete signal is deposited Zero etc. is not all fluctuating amplitude of the waveform of the either target discrete signal not for straight line or the target discrete signal.It can Choosing, which can be filtered treated target discrete signal, the target after the filtering process Discrete signal can include at least two target discrete data, such as filtering process can be that cutoff frequency is 1 hertz (Hz) It is 40Hz to have limit for length's unit impulse response (Finite Impulse Response, FIR) high-pass filter and/or cutoff frequency FIR low pass filter, the FIR high-pass filters of 1Hz are mainly used for filtering out the baseline drift of below 1Hz in the target discrete signal It moves, the FIR low pass filter of 40Hz is mainly used for filtering out in the target discrete signal 40Hz with noise of super band etc..

Computing module 20, for calculating the value of the fisrt feature parameter of the target discrete signal, the fisrt feature ginseng Several values is used to represent the wave character of the target discrete signal；

Specifically, the computing module 20 of terminal can calculate the value of the fisrt feature parameter of above-mentioned target discrete signal, should The value of fisrt feature parameter is used to represent the wave character of the target discrete signal.Wherein, which can include At least one of approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, signal amplitude variation coefficient.The approximate entropy can be used for Represent the information content variation degree of the target discrete signal；The frequency spectrum proportionality coefficient can be used to indicate that the target discrete signal The spectrum energy of target frequency bands accounts for the ratio of the frequency spectrum gross energy of the target discrete signal, for example, target frequency bands can be 1Hz~ 20Hz (spectrum energy of ballistocardiography signal is concentrated mainly on 1Hz~20Hz), the frequency spectrum proportionality coefficient can be expressed as the mesh Mark the ratio of the frequency spectrum gross energy of the spectrum energy and the target discrete signal between 1Hz~20Hz of discrete signal；The zero-crossing rate Density is used to represent the fluctuation high density state of the target discrete signal；The signal amplitude variation coefficient is used to represent the target discrete The ratio of the standard deviation of the average and peak value of peak value in signal.For example, if fisrt feature parameter is approximate entropy, terminal can be counted Count stating the value of the approximate entropy of target discrete signal in；If fisrt feature parameter is approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate are close Degree, it is close that terminal can calculate the value of the approximate entropy of above-mentioned target discrete signal, the value of frequency spectrum proportionality coefficient and zero-crossing rate respectively The value of degree.

Optionally, if fisrt feature parameter is frequency spectrum proportionality coefficient, terminal can carry out above-mentioned target discrete signal fast Fast Fourier transformation, obtains the frequency spectrum of the target discrete signal, and computing module 20 calculates the target frequency bands of the target discrete signal Spectrum energy and the target discrete signal frequency spectrum gross energy ratio, which is determined as to the value of frequency spectrum proportionality coefficient. For example, the above-mentioned target discrete signal of terminal-pair carries out Fast Fourier Transform (FFT), the frequency spectrum of the target discrete signal, Frequency spectrum ratio are obtained Example coefficient accounts for the ratio of frequency spectrum gross energy for the spectrum energy between 1Hz~20Hz, can be represented with PSDr：

Optionally, if fisrt feature parameter is zero-crossing rate density, the computing module 20 of terminal can calculate above-mentioned target from The value of the zero-crossing rate density of scattered signal.Such as zero-crossing rate density can be represented with ZCRD：

Optionally, if fisrt feature parameter is signal amplitude variation coefficient, the computing module 20 of terminal can calculate above-mentioned The value of the signal amplitude variation coefficient of target discrete signal, the signal amplitude variation coefficient can be also used for reflecting above-mentioned target from The peak change degree of scattered signal.For example, signal amplitude variation coefficient can be represented with following formula：

Determining module 30 for the value according to the fisrt feature parameter, determines that the target discrete signal is rushed as the heart The credit rating information of figure signal is hit, the credit rating information is used to represent ballistocardiography described in the target discrete signal The content of signal.

Specifically, above-mentioned fisrt feature parameter can include approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, signal width It is worth at least one of variation coefficient, the determining module 30 of terminal can be according to the above-mentioned fisrt feature parameter being calculated Value, determines credit rating information of the above-mentioned target discrete signal as ballistocardiography signal, passes through the matter of the ballistocardiography signal Measure class information, it can be determined that the accuracy of heart activity information.If for example, above-mentioned fisrt feature parameter be approximate entropy, terminal According to the value of the above-mentioned approximate entropy being calculated, determine that above-mentioned target discrete signal is believed as the credit rating of ballistocardiography signal Breath；If above-mentioned fisrt feature parameter is approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, terminal is calculated according to above-mentioned The value of the value of approximate entropy, the value of frequency spectrum proportionality coefficient and zero-crossing rate density determines that above-mentioned target discrete signal is impacted as the heart The credit rating information of figure signal.

Optionally, if above-mentioned fisrt feature parameter includes approximate entropy, the computing module 20 of terminal can calculate the target from The value of the approximate entropy of scattered signal, if the value of the approximate entropy, in the first scope, determining module 30 determines the target discrete signal conduct The credit rating information of ballistocardiography signal is the highest ranking in the plurality of grade, which corresponds to target discrete letter Number as the best in quality of ballistocardiography signal.This first may range from the ballistocardiography signal in the target discrete signal The fluctuation range of the approximate entropy of target discrete signal when content is more than first threshold.For example, first threshold is 0.8, this first The approximate entropy of target discrete signal when may range from the content of the ballistocardiography signal in the target discrete signal more than 0.8 Fluctuation range, if more than the approximate entropy of target discrete signal when 0.8 fluctuation range for 0.01~0.15, then the first scope For 0.01~0.15.

Optionally, if above-mentioned fisrt feature parameter includes approximate entropy, the computing module 20 of terminal can calculate the target from The value of the approximate entropy of scattered signal, if the value of the approximate entropy, in the second scope, determining module 30 determines the target discrete signal conduct The credit rating information of ballistocardiography signal is the lowest class in the plurality of grade, which corresponds to target discrete letter Number as ballistocardiography signal quality it is worst.This second may range from the ballistocardiography signal in the target discrete signal The fluctuation range of the approximate entropy of target discrete signal when content is less than second threshold, the second threshold are less than above-mentioned first threshold Value.For example, first threshold is 0.8, second threshold can be 0.2, this second may range from the heart in the target discrete signal Impact figure signal content be less than 0.2 when the target discrete signal approximate entropy fluctuation range, if less than 0.2 when the target from The fluctuation range of the approximate entropy of scattered signal is more than 0.18, then first may range from being more than 0.18.

It is further alternative, if above-mentioned fisrt feature parameter includes approximate entropy and frequency spectrum proportionality coefficient, the calculating mould of terminal Block 20 can calculate the value of the approximate entropy of the target discrete signal and the value of frequency spectrum proportionality coefficient, if the value of the approximate entropy is Three scopes, the terminal can obtain the initial grade of credit rating information, and the 3rd scope is the heart in the target discrete signal The fluctuation of the approximate entropy of target discrete signal when impacting the content of figure signal more than the second threshold, and being less than the first threshold Scope, such as first threshold are 0.8, second threshold 0.2, and the 3rd scope is that the ballistocardiography is believed in the target discrete signal Number content be more than 0.2, and less than 0.8 when the target discrete signal approximate entropy fluctuation range, such as the 3rd scope be 0.15 ~0.18.The terminal can determine the first alternative grade according to the value and the initial grade of the frequency spectrum proportionality coefficient.Wherein, if The value of the frequency spectrum proportionality coefficient is more than the 3rd threshold value, which is higher than the initial grade, if the frequency spectrum proportionality coefficient Value be less than or equal to the 3rd threshold value, the first alternative grade is identical with the initial grade.With the grade and quality in Fig. 2 Correspondence exemplified by, number it is bigger expression higher grade, it is assumed that grade can be represented with the value of score, such as score=2 It can represent grade 2, score=3 can represent grade 3 etc., and the 3rd threshold value can be 0.6, initial grade score=3, if The value of frequency spectrum proportionality coefficient is more than 0.6, and the first alternative grade can be that initial grade add up the first incremental change, such as the first incremental change For 1 etc., i.e., the first alternative grade score=score+1=3+1=4 is higher than initial grade score=3；If the frequency spectrum ratio system Several values is less than or equal to 0.6, and the first alternative grade is identical with the initial grade, i.e., the first alternative grade score= Score=3.

It is further optional, if above-mentioned fisrt feature parameter is close including approximate entropy, frequency spectrum proportionality coefficient and zero-crossing rate Degree, the computing module 20 of terminal can calculate the value of the approximate entropy of the target discrete signal, the value and mistake of frequency spectrum proportionality coefficient The value of zero rate density, the value and above-mentioned first alternative grade that the determining module 30 of terminal can be according to the zero-crossing rate density, determines Second alternative grade, wherein, if the value of the zero-crossing rate density is less than or equal to the 4th threshold value, which is higher than should First alternative grade, if the value of the zero-crossing rate density is more than the 4th threshold value, the second alternative grade and the first alternative grade It is identical, alternatively, the second alternative grade is less than the first alternative grade.By taking the grade in Fig. 2 and the correspondence of quality as an example, Number is bigger, and expression higher grade, it is assumed that and grade can be represented with the value of score, if score=2 can represent grade 2, Score=3 can represent grade 3 etc., and the 4th threshold value is 4, the first alternative grade score=4, if the value of zero-crossing rate density is small In or equal to 4, the second alternative grade can be that the first alternative grade add up the second incremental change, and such as the second incremental change is 1, i.e., Second alternative grade score=score+1=4+1=5 is higher than the first alternative grade score=4；If the value of zero-crossing rate density is big Can be identical with the first alternative grade in 4, second alternative grades, i.e., the second alternative grade score=score, (score≤ 3), alternatively, the second alternative grade can be first alternative the first decrement of grade regressive, such as the first decrement is 1, i.e., second Alternative grade score=score-1, (score>3), because the first alternative grade score=4>3, therefore the second alternative grade Score=score-1=4-1=3.

It is further optional, if above-mentioned fisrt feature parameter include approximate entropy, frequency spectrum ratio system, zero-crossing rate density and Signal amplitude variation coefficient, the computing module 20 of terminal can calculate the value of the approximate entropy of the target discrete signal, frequency spectrum ratio The value of the value of coefficient, the value of zero-crossing rate density and signal amplitude variation coefficient, the determining module 30 of terminal can be according to the letter The value and above-mentioned second alternative grade of number amplitude variation coefficient, determine the 3rd alternative grade, and the 3rd alternative grade are determined For credit rating information of the target discrete signal as ballistocardiography signal, the credit rating so as to more be refined is believed Breath improves the accuracy for judging heart activity information.Wherein, if the value of the signal amplitude variation coefficient is less than or equal to the 5th Threshold value, the 3rd alternative grade are higher than the second alternative grade, if the value of the signal amplitude variation coefficient is more than the 5th threshold value, 3rd alternative grade is identical with the second alternative grade, alternatively, the 3rd alternative grade is less than the second alternative grade.To scheme Grade in 2 and exemplified by the correspondence of quality, number is bigger, and expression higher grade, it is assumed that grade can with the value of score come It represents, if score=2 can represent grade 2, score=3 can represent grade 3 etc., and the 5th threshold value is 0.5, and second is alternative Grade is grade score=5, if the value of signal amplitude variation coefficient is less than or equal to 0.5, the 3rd alternative grade can be the Two alternative grades add up the 3rd incremental change, such as the 3rd incremental change is 1, i.e., the 3rd alternative grade score=score+1=5+1 =6 are higher than the second alternative grade score=5, and the 3rd alternative grade score=6 is determined as the target discrete signal as the heart Impact the credit rating information of figure signal；If the value of signal amplitude variation coefficient is more than 0.5, the 3rd alternative grade can with this Two alternative grades are identical, i.e. score=score, (score≤3), alternatively, the 3rd alternative grade can be the second alternative grade The second decrement of regressive, such as the second decrement are 1, i.e. score=score-1, (score>3), due to the second alternative grade Score=5>3rd alternative grade score=4 is determined as this by 3, the 3rd alternative grade score=score-1=5-1=4 Credit rating information of the target discrete signal as ballistocardiography signal.

It is the schematic block diagram of a kind of electronic equipment provided in an embodiment of the present invention referring to Fig. 5.This implementation as depicted Electronic equipment in example can include one or more input equipments 1000, one or more output equipments 2000, one or more A processor 3000 and memory 4000.Above-mentioned processor 3000, input equipment 1000, output equipment 2000 and memory 4000 It is connected by bus 5000.For memory 4000 for storing computer program, the computer program includes program instruction, processing Device 3000 is used to perform the program instruction of the storage of memory 4000.Wherein, processor 3000 is arranged to call described program Instruction performs：

Obtain target discrete signal；

Optionally, above-mentioned fisrt feature parameter includes approximate entropy, frequency spectrum proportionality coefficient, zero-crossing rate density, signal amplitude change Change at least one of coefficient, the approximate entropy is used to represent the information content variation degree of the target discrete signal, the frequency Spectrum proportionality coefficient is used to represent that the spectrum energy of the target frequency bands of the target discrete signal to account for the frequency of the target discrete signal The ratio of gross energy is composed, the zero-crossing rate density is used to represent the fluctuation high density state of the target discrete signal, the signal Amplitude variation coefficient is used to represent the ratio of the standard deviation of the average and peak value of peak value in the target discrete signal.

Optionally, above-mentioned credit rating information includes multiple grades, described in the different brackets in the multiple grade corresponds to The different content of ballistocardiography signal described in target discrete signal, the more high corresponding content of grade are higher；

Above-mentioned fisrt feature parameter includes approximate entropy, and above-mentioned processor 3000 is specifically used for：

Wherein, the second threshold is less than the first threshold.

Optionally, above-mentioned fisrt feature parameter further includes frequency spectrum proportionality coefficient, above-mentioned processor 3000 also particularly useful for：

Optionally, above-mentioned fisrt feature parameter further includes zero-crossing rate density, and above-mentioned processor 3000 is specifically used for：

Optionally, above-mentioned fisrt feature parameter further includes signal amplitude variation coefficient, and above-mentioned processor 3000 is specifically used for：

It should be appreciated that in embodiments of the present invention, alleged processor 3000 can be central processing unit (Central Processing Unit, CPU), which can also be other general processors, digital signal processor (Digital Signal Processor, DSP), application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic Device, discrete gate or transistor logic, discrete hardware components etc..General processor can be microprocessor or this at It can also be any conventional processor etc. to manage device.

Input equipment 1000 can include piezoelectric cable, piezoresistive transducer, piezoelectric film sensor (for gathering human body The vibration information of vertebra vibration) etc., output equipment 2000 can include display (LCD etc.), loud speaker etc..

The memory 4000 can include read-only memory and random access memory, and provide instruction to processor 3000 And data.The a part of of memory 4000 can also include nonvolatile RAM.For example, memory 4000 may be used also With the information of storage device type.

In the specific implementation, input equipment 1000, output equipment 2000, processor 3000 described in the embodiment of the present invention The realization method described in signal processing method provided in an embodiment of the present invention is can perform, also can perform institute of the embodiment of the present invention The realization method of the signal processing apparatus of description, details are not described herein.

A kind of computer readable storage medium is also provided in the embodiment of the present invention, which has Computer program, the computer program include program instruction, which realizes the letter of Fig. 1 or Fig. 2 when being executed by processor Number processing method, detail refer to the description of Fig. 1 or Fig. 2 embodiments, and details are not described herein.

Above computer readable storage medium storing program for executing can be that signal processing apparatus described in foregoing any embodiment or electronics are set Standby internal storage unit, such as the hard disk or memory of electronic equipment.The computer readable storage medium can also be the electronics The plug-in type hard disk being equipped on the External memory equipment of equipment, such as the electronic equipment, intelligent memory card (Smart Media Card, SMC), secure digital (Secure Digital, SD) card, flash card (Flash Card) etc..Further, the calculating Machine readable storage medium storing program for executing can also both include the internal storage unit of the electronic equipment or including External memory equipment.The computer Readable storage medium storing program for executing is used to store the computer program and other programs and data needed for the electronic equipment.The computer can Storage medium is read to can be also used for temporarily storing the data that has exported or will export.

Those of ordinary skill in the art may realize that each exemplary moulds described with reference to the embodiments described herein Block and algorithm steps can be realized with the combination of electronic hardware, computer software or the two, in order to clearly demonstrate hardware With the interchangeability of software, each exemplary composition and step are generally described according to function in the above description.This A little functions are performed actually with hardware or software mode, specific application and design constraint depending on technical solution.Specially Industry technical staff can realize described function to each specific application using distinct methods, but this realization is not It is considered as beyond the scope of this invention.

The present invention be with reference to the embodiment of the present invention method, apparatus (equipment) and computer program product flow chart with/ Or block diagram describes.It should be understood that can by computer program instructions realize flowchart and/or the block diagram in each flow and/ Or the flow in box and flowchart and/or the block diagram and/or the combination of box.These computer program instructions can be provided To the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to generate one A machine so that the instruction generation performed by computer or the processor of other programmable data processing devices is used to implement The device for the function of being specified in one flow of flow chart or multiple flows and/or one box of block diagram or multiple boxes.

These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works so that the instruction generation being stored in the computer-readable memory includes referring to Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or The function of being specified in multiple boxes.

These computer program instructions can be also loaded into computer or other programmable data processing devices so that counted Series of operation steps is performed on calculation machine or other programmable devices to generate computer implemented processing, so as in computer or The instruction offer performed on other programmable devices is used to implement in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in a box or multiple boxes.

Although with reference to specific features and embodiment, invention has been described, it is clear that, do not departing from this hair In the case of bright spirit and scope, it can be carry out various modifications and combined.Correspondingly, the specification and drawings are only institute The exemplary illustration for the present invention that attached claim is defined, and be considered as covered in the scope of the invention arbitrary and all and repair Change, change, combining or equivalent.Obviously, those skilled in the art various changes and modifications can be made to the invention without Depart from the spirit and scope of the present invention.If in this way, these modifications and changes of the present invention belong to the claims in the present invention and its Within the scope of equivalent technologies, then the present invention is also intended to comprising including these modification and variations.

Claims

1. a kind of signal processing method, which is characterized in that including：

Obtain target discrete signal；

The value of the fisrt feature parameter of the target discrete signal is calculated, the value of the fisrt feature parameter is used to represent the mesh Mark the wave character of discrete signal；

According to the value of the fisrt feature parameter, determine that the target discrete signal is believed as the credit rating of ballistocardiography signal Breath, the credit rating information are used to represent the content of ballistocardiography signal described in the target discrete signal.

2. according to the method described in claim 1, it is characterized in that, the fisrt feature parameter includes approximate entropy, frequency spectrum ratio At least one of coefficient, zero-crossing rate density, signal amplitude variation coefficient, the approximate entropy are used to represent the target discrete letter Number information content variation degree, the frequency spectrum proportionality coefficient be used for represent the target discrete signal target frequency bands frequency spectrum energy Amount accounts for the ratio of the frequency spectrum gross energy of the target discrete signal, and the zero-crossing rate density is used to represent the target discrete signal Fluctuation high density state, the signal amplitude variation coefficient for represent in the target discrete signal average of peak value with it is described The ratio of the standard deviation of peak value.

3. according to the method described in claim 2, it is characterized in that, the credit rating information include multiple grades, it is described more Different brackets in a grade corresponds to the different content of ballistocardiography signal described in the target discrete signal, and grade is more high right The content answered is higher；

The fisrt feature parameter includes approximate entropy, and the value according to the fisrt feature parameter determines the target discrete Credit rating information of the signal as ballistocardiography signal, including：

If the value of the approximate entropy in the first scope, determines credit rating of the target discrete signal as ballistocardiography signal Information is the highest ranking in the multiple grade, and first scope is the letter of ballistocardiography described in the target discrete signal Number content be more than first threshold when the target discrete signal approximate entropy fluctuation range；

If the value of the approximate entropy in the second scope, determines credit rating of the target discrete signal as ballistocardiography signal Information is the lowest class in the multiple grade, and second scope is the letter of ballistocardiography described in the target discrete signal Number content be less than second threshold when the target discrete signal approximate entropy fluctuation range；

Wherein, the second threshold is less than the first threshold.

4. according to the method described in claim 3, it is characterized in that, the fisrt feature parameter further includes frequency spectrum proportionality coefficient, The value according to the fisrt feature parameter determines that the target discrete signal is believed as the credit rating of ballistocardiography signal Breath, further includes：

If the value of the approximate entropy in the 3rd scope, obtains the initial grade of credit rating information, the 3rd scope is described The content of ballistocardiography signal described in target discrete signal be more than the second threshold, and less than the first threshold when described in The fluctuation range of the approximate entropy of target discrete signal；

According to the value of the frequency spectrum proportionality coefficient and the initial grade, the first alternative grade is determined, wherein, if the Frequency spectrum ratio The value of example coefficient is more than the 3rd threshold value, and the first alternative grade is higher than the initial grade, if the frequency spectrum proportionality coefficient Value is less than or equal to the 3rd threshold value, and the first alternative grade is identical with the initial grade；

According to the described first alternative grade, credit rating information of the target discrete signal as ballistocardiography signal is determined.

5. method as claimed in claim 4, which is characterized in that the fisrt feature parameter further includes zero-crossing rate density, described According to the described first alternative grade, credit rating information of the target discrete signal as ballistocardiography signal is determined, including：

According to the value of the zero-crossing rate density and the first alternative grade, the second alternative grade is determined, wherein, if the zero passage The value of rate density is less than or equal to the 4th threshold value, and the second alternative grade is higher than the described first alternative grade, if the mistake The value of zero rate density is more than the 4th threshold value, and the second alternative grade is identical with the described first alternative grade, alternatively, described Second alternative grade is less than the described first alternative grade；

According to the described second alternative grade, credit rating information of the target discrete signal as ballistocardiography signal is determined.

6. method as claimed in claim 5, which is characterized in that the fisrt feature parameter further includes signal amplitude variation system Number, it is described according to the described second alternative grade, determine that the target discrete signal is believed as the credit rating of ballistocardiography signal Breath, including：

According to the value of the signal amplitude variation coefficient and the second alternative grade, the 3rd alternative grade is determined, wherein, if institute The value for stating signal amplitude variation coefficient is less than or equal to the 5th threshold value, and the 3rd alternative grade is alternative higher than described second etc. Grade, if the value of the signal amplitude variation coefficient is more than the 5th threshold value, the 3rd alternative grade is alternative with described second Grade is identical, alternatively, the 3rd alternative grade is less than the described second alternative grade；

The 3rd alternative grade is determined as credit rating information of the target discrete signal as ballistocardiography signal.

7. a kind of signal processing apparatus, which is characterized in that including：

Acquisition module, for obtaining target discrete signal；

Computing module, for calculating the value of the fisrt feature parameter of the target discrete signal, the value of the fisrt feature parameter For representing the wave character of the target discrete signal；

Determining module for the value according to the fisrt feature parameter, determines that the target discrete signal is believed as ballistocardiography Number credit rating information, the credit rating information is for representing ballistocardiography signal described in the target discrete signal Content.

8. device according to claim 7, which is characterized in that the fisrt feature parameter includes approximate entropy, frequency spectrum ratio At least one of coefficient, zero-crossing rate density, signal amplitude variation coefficient, the approximate entropy are used to represent the target discrete letter Number information content variation degree, the frequency spectrum proportionality coefficient be used for represent the target discrete signal target frequency bands frequency spectrum energy Amount accounts for the ratio of the frequency spectrum gross energy of the target discrete signal, and the zero-crossing rate density is used to represent the target discrete signal Fluctuation high density state, the signal amplitude variation coefficient for represent in the target discrete signal average of peak value with it is described The ratio of the standard deviation of peak value.

9. a kind of electronic equipment, which is characterized in that the electronic equipment includes：Processor, input equipment, output equipment and storage Device, the processor, input equipment, output equipment and memory are connected with each other, wherein, the memory is used to store computer Program, the computer program include program instruction, and the processor is arranged to call described program instruction, perform as weighed Profit requires 1-6 any one of them methods.

10. a kind of computer readable storage medium, the computer storage media is stored with computer program, the computer journey Sequence includes program instruction, and described program instruction makes the processor perform as claim 1-6 is any when being executed by a processor Method described in.