SYSTEM AND METHOD FOR DETECTING A PRESENCE IN A CLOSED ENVIRONMENT TO BE MONITORED, FOR ANTI-INTRUSION OR ANTI-THEFT PURPOSE

Abstract

A system for detecting a presence in an environment to be monitored includes an electrostatic charge variation sensor, a vibration sensor, and an environmental pressure sensor. A processing unit is configured to acquire, from the electrostatic charge variation sensor, an electrostatic charge variation signal, and detect in the electrostatic charge variation signal, first signal characteristics indicative of the presence of a subject in the environment to be monitored. The processing unit further validates the detection of presence of the subject using the vibration and pressure signals provided by the other sensors.

Claims

1. A system for detecting a presence in an environment to be monitored, comprising: a processor; an electrostatic charge variation sensor coupled to the processor and configured to detect a variation of electrostatic charge in said environment and generate an electrostatic charge variation signal; and one of a vibration sensor operatively coupled to the environment to be monitored and configured to detect an environmental vibration in the environment to be monitored and generate a vibration signal, or an environmental pressure sensor operatively coupled to the environment to be monitored and configured to detect an environmental pressure in the environment to be monitored and generate a pressure signal, wherein the processor is configured to: acquire, from the electrostatic charge variation sensor, the electrostatic charge variation signal; detect, in said electrostatic charge variation signal, first signal characteristics indicative of the presence of a subject in said environment to be monitored; acquire, from said one of the vibration sensor or the environmental pressure sensor, respectively the vibration signal or the pressure signal; detect, in said vibration signal or pressure signal acquired, respective second signal characteristics indicative of the presence of the subject in said environment to be monitored; and generate a warning signal if both the first and the second signal characteristics have been detected.

2. The system according to claim 1, further comprising the other of the vibration sensor or the environmental pressure sensor, wherein the processor is further configured to: acquire, from the other of said vibration sensor or environmental pressure sensor, respectively the vibration signal or the pressure signal; detect, in said other of the vibration signal or the pressure signal acquired, respective third signal characteristics indicative of the presence of the subject in said environment to be monitored; and generate the warning signal if all the first, the second and the third signal characteristics have been detected.

3. The system according to claim 1, wherein the detecting the first signal characteristics includes: detecting, in the electrostatic charge variation signal, the following characteristics which follow each other in temporal order: a first rising edge; a first local maximum; a first falling edge; a first local minimum; a second rising edge; and alternatively, detecting, in the electrostatic charge variation signal, the following characteristics which follow each other in temporal order: a falling edge; a first local minimum; a first rising edge; a first local maximum; a second falling edge.

4. The system according to claim 3, wherein the detecting the first signal characteristics further includes: performing a comparison of said local maximums and minimums with respective thresholds; and assessing, through comparison with respective thresholds, a value of steepness or rising rate of the first and the second rising edges and of steepness or falling rate of the falling edge.

5. The system according to claim 3, wherein the detecting, in the electrostatic charge variation signal, the characteristics that follow each other in temporal order includes: calculating the first derivative of the electrostatic charge variation signal; identifying, in the electrostatic charge variation signal and in the first derivative signal, a respective plurality of positive and negative peaks; detecting one of the following time series a) and b) with which said plurality of positive and negative peaks follow each other over time: a) a first positive peak in the first derivative signal; a second positive peak in the electrostatic charge variation signal; a first negative peak in the first derivative signal; a second negative peak in the electrostatic charge variation signal; a third positive peak in the first derivative signal, b) a third negative peak in the first derivative signal; a fourth negative peak in the electrostatic charge variation signal; a fourth positive peak in the first derivative signal; a fifth positive peak in the electrostatic charge variation signal; a fifth negative peak in the first derivative signal.

6. The system according to claim 5, wherein the detecting the first signal characteristics further includes detecting one or more of the following time relationships:
T1=T3+T4
T6=T2+T3
T7=T4+T5
T8=T6+T7 wherein: T1 is a time interval between the second positive peak and the second negative peak, T2 is a time interval between the second positive peak and the first positive peak, T3 is a time interval between the second positive peak and the first negative peak, T4 is a time interval between the second negative peak and the first negative peak, T5 is a time interval between the second negative peak and the third positive peak, T6 is a time interval between the first positive peak and the first negative peak, T7 is a time interval between the first negative peak and the third positive peak, T8 is a time interval between the first positive peak and the third positive peak.

7. The system according to claim 6, wherein said time intervals T1 through T7 are respective time distances between respective maximum or minimum points of the positive and negative peaks.

8. The system according to claim 6, wherein the detecting the first signal characteristics further includes detecting one or more of the following time relationships: T2.sub.TH_L<T2<T2.sub.TH_H, T3.sub.TH_L<T3<T3.sub.TH_H, T4.sub.TH_L<T4<T4.sub.TH_H, T5.sub.TH_L<T5<T5.sub.TH_H, where T2.sub.TH_L, T3.sub.TH_L, T4.sub.TH_L and T5.sub.TH_L may be respective lower thresholds of respective value between 30 and 70 ms, and T2.sub.TH_H, T3.sub.TH_H, T4.sub.TH_H and T5.sub.TH_H are respective higher thresholds of respective value including between 150 and 250 ms.

9. The system according to claim 1, wherein the second signal characteristics belong to the pressure signal, and wherein said detecting in said pressure signal the second signal characteristics includes: detecting a time amplitude value and maximum value of a pressure peak present in said pressure signal; calculating a first comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the pressure peak; and verifying whether said first comparison parameter is in a predetermined relationship with a first threshold.

10. The system according to claim 9, wherein detecting a time amplitude value includes calculating an integral of, or an area subtended by, the pressure peak present in said pressure signal, and wherein said first comparison parameter is calculated by dividing the result of said integral of the pressure peak, or the value of said area subtended by the pressure peak, by the maximum value of the pressure peak.

11. The system according to claim 2, wherein the third signal characteristics belong to the vibration signal, and wherein said detecting in said vibration signal the third signal characteristics includes: detecting a time amplitude value and maximum value of a vibration peak present in said vibration signal; calculating a second comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the vibration peak; and verifying whether said second comparison parameter is in a predetermined relationship with a second threshold.

12. The system according to claim 11, wherein detecting a time amplitude value includes calculating an integral of, or an area subtended by, the vibration peak present in said vibration signal, and wherein said second comparison parameter is calculated by dividing the result of said integral of the vibration peak, or the value of said area subtended by the vibration peak, by the maximum value of the vibration peak.

13. A method for detecting a presence in an environment to be monitored, comprising: detecting, by an electrostatic charge variation sensor, a variation of electrostatic charge in said environment and generating an electrostatic charge variation signal; detecting, by one of a vibration sensor or an environmental sensor operatively coupled to the environment to be monitored, respectively, an environmental vibration in the environment to be monitored and generating a vibration signal or an environmental pressure in the environment to be monitored and generating a pressure signal; acquiring, by a processor, from the electrostatic charge variation sensor, the electrostatic charge variation signal; detecting, by the processor, in said electrostatic charge variation signal, first signal characteristics indicative of the presence of a subject in said environment to be monitored; acquiring, by the processor, from said one of the vibration sensor or the environmental pressure sensor, respectively the vibration signal or the pressure signal; detecting, by the processor, in said vibration signal or pressure signal acquired, respective second signal characteristics indicative of the presence of the subject in said environment to be monitored; and generating, by the processor, a warning signal if both the first and the second signal characteristics have been detected.

14. The method according to claim 13, further comprising: detecting the environmental vibration and generating the vibration signal or detecting the environmental pressure and generating the pressure signal, by the other of the vibration sensor and the environmental pressure sensor; acquiring, from the other of said vibration sensor and environmental pressure sensor, respectively the vibration signal or the pressure signal; detecting, in said other of the vibration signal and the pressure signal acquired, respective third signal characteristics indicative of the presence of the subject in said environment to be monitored; and generating the warning signal if all the first, the second and the third signal characteristics have been detected.

15. The method according to claim 13, wherein the detecting the first signal characteristics includes: detecting, in the electrostatic charge variation signal, the following characteristics which follow each other in temporal order: a first rising edge; a first local maximum; a first falling edge; a first local minimum; a second rising edge; and alternatively, detecting, in the electrostatic charge variation signal, the following characteristics which follow each other in temporal order: a falling edge; a first local minimum; a first rising edge; a first local maximum; a second falling edge.

16. The method according to claim 15, wherein the detecting the first signal characteristics further includes: performing a comparison of said local maximums and minimums with respective thresholds; and assessing, through comparison with respective thresholds, a value of steepness or rising rate of the first and the second rising edges and of steepness or falling rate of the falling edge.

17. The method according to claim 15, wherein detecting, in the electrostatic charge variation signal, the characteristics that follow each other in temporal order includes: calculating the first derivative of the electrostatic charge variation signal; identifying, in the electrostatic charge variation signal and in the first derivative signal, a respective plurality of positive and negative peaks; detecting one of the following time series a) and b) with which said plurality of positive and negative peaks follow each other over time: a) a first positive peak in the first derivative signal; a second positive peak in the electrostatic charge variation signal; a first negative peak in the first derivative signal; a second negative peak in the electrostatic charge variation signal; a third positive peak in the first derivative signal, b) a third negative peak in the first derivative signal; a fourth negative peak in the electrostatic charge variation signal; a fourth positive peak in the first derivative signal; a fifth positive peak in the electrostatic charge variation signal; a fifth negative peak in the first derivative signal.

18. The method according to claim 17, wherein the detecting the first signal characteristics further includes detecting one or more of the following time relationships:
T1=T3+T4
T6=T2+T3
T7=T4+T5
T8=T6+T7 wherein: T1 is a time interval between the second positive peak and the second negative peak, T2 is a time interval between the second positive peak and the first positive peak, T3 is a time interval between the second positive peak and the first negative peak, T4 is a time interval between the second negative peak and the first negative peak, T5 is a time interval between the second negative peak and the third positive peak, T6 is a time interval between the first positive peak and the first negative peak, T7 is a time interval between the first negative peak and the third positive peak, T8 is a time interval between the first positive peak and the third positive peak.

19. The method according to claim 18, wherein said time intervals T1 through T7 are respective time distances between respective maximum or minimum points of the positive and negative peaks.

20. The method according to claim 18, wherein the detecting the first signal characteristics further includes detecting one or more of the following time relationships: T2.sub.TH_L<T2<T2.sub.TH_H, T3.sub.TH_L<T3<T3.sub.TH_H, T4.sub.TH_L<T4<T4.sub.TH_H, T5.sub.TH_L<T5<T5.sub.TH_H, where T2.sub.TH_L, T3.sub.TH_L, T4.sub.TH_L and T5.sub.TH_L may be respective lower thresholds of respective value between 30 and 70 ms, and T2.sub.TH_H, T3.sub.TH_H, T4.sub.TH_H and T5.sub.TH_H are respective higher thresholds of respective value including between 150 and 250 ms.

21. The method according to claim 13, wherein the second signal characteristics belong to the pressure signal, and wherein said detecting in said pressure signal the second signal characteristics includes: detecting a time amplitude value and maximum value of a pressure peak present in said pressure signal; calculating a first comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the pressure peak; and verifying whether said first comparison parameter is in a predetermined relationship with a first threshold.

22. The method according to claim 21, wherein detecting a time amplitude value includes calculating an integral of, or an area subtended by, the pressure peak present in said pressure signal, and wherein said first comparison parameter is calculated by dividing the result of said integral of the pressure peak, or the value of said area subtended by the pressure peak, by the maximum value of the pressure peak.

23. The method according to claim 14, wherein the third signal characteristics belong to the vibration signal, and wherein said detecting in said vibration signal the third signal characteristics includes: detecting a time amplitude value and maximum value of a vibration peak present in said vibration signal; calculating a second comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the vibration peak; and verifying whether said second comparison parameter is in a predetermined relationship with a second threshold.

24. The method according to claim 23, wherein detecting a time amplitude value includes calculating an integral of, or an area subtended by, the vibration peak present in said vibration signal, and wherein said second comparison parameter is calculated by dividing the result of said integral of the vibration peak, or the value of said area subtended by the vibration peak, by the maximum value of the vibration peak.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

[0019] For a better understanding of the disclosure, embodiments thereof are now described, purely by way of non-limiting example and with reference to the attached drawings, wherein:

[0020] FIG. 1 schematically illustrates a system for detecting a presence including an environmental electric charge sensor, a pressure sensor and a vibration sensor (in particular a multi-axial accelerometer) operatively coupled to a processing unit, according to an embodiment of the present disclosure;

[0021] FIG. 2 illustrates an embodiment of the environmental electrostatic charge variation sensor;

[0022] FIG. 3A illustrates an example of a pressure signal Sp generated by the pressure sensor of FIG. 1;

[0023] FIG. 3B illustrates an example of an electrostatic charge variation signal generated by the electrostatic charge variation sensor of FIG. 1;

[0024] FIG. 3C illustrates an example of a vibration signal generated by the accelerometer of FIG. 1 and partially processed by the processing unit of FIG. 1 in order to generate the modulus of the sensing axial components;

[0025] FIG. 4A illustrates the pressure signal of FIG. 3A after the removal of one of its background components, or baseline;

[0026] FIG. 4B illustrates the electrostatic charge variation signal of FIG. 3B after the removal of the relative baseline;

[0027] FIG. 4C illustrates the first derivative of the electrostatic charge variation signal of FIG. 4B;

[0028] FIG. 4D illustrates the envelope, or AC component, of the vibration signal of FIG. 3C;

[0029] FIGS. 5A and 5B illustrate, respectively, an enlarged portion of the electrostatic charge variation signal of FIG. 4B and of the first derivative of FIG. 4C;

[0030] FIG. 6 illustrates, through a flow chart, a method for detecting the human presence implemented by the system of FIG. 1, with exclusive reference to the electrostatic charge variation signal;

[0031] FIG. 7 illustrates, through a block diagram, steps of a method for analyzing the pressure signal of FIG. 3C or 4A, in order to extract or identify significant characteristics for detecting the human presence;

[0032] FIGS. 8A-8C graphically show steps for processing the pressure signal according to the method of FIG. 7;

[0033] FIG. 9 illustrates, through a block diagram, steps of a method for analyzing the vibration signal of FIG. 3C, in order to extract or identify significant characteristics for detecting the human presence;

[0034] FIGS. 10A-10C graphically show steps for processing the vibration signal according to the method of FIG. 9;

[0035] FIGS. 11 and 12 show, through a block diagram, respective methods for removing the baseline, applicable in the context of the present disclosure to generate the signals of FIGS. 4A and 4B;

[0036] FIG. 13 illustrates, through a block diagram, steps of a method for detecting peaks usable in the context of the present disclosure to recognize positive and negative peaks, applicable in the context of the methods of FIGS. 6 and 7;

[0037] FIG. 14 illustrates, through a block diagram, a method for calculating the first derivative of the electrostatic charge variation signal of FIG. 4B, to obtain the signal of FIG. 4C; and

[0038] FIG. 15 illustrates, through a block diagram, a method for extracting the envelope, or AC component, usable to generate the vibration signal of FIG. 4D from the vibration signal of FIG. 3C.

DETAILED DESCRIPTION

[0039] FIG. 1 schematically illustrates a presence detection system, or anti-intrusion system, 1. The presence detection system 1 is in particular for detecting a human presence in an environment and comprises a processing unit 2, a pressure sensor 4, coupled to the processing unit 2, an electrostatic charge variation sensor 6 coupled to the processing unit 2, and a vibrational sensor 7, in particular an accelerometer, also coupled to the processing unit 2 (hereinafter reference will explicitly be made to the accelerometer, without thereby losing generality). The pressure sensor 4, the electrostatic charge variation sensor 6 and the accelerometer 7 are arranged in the environment to be monitored; the processing unit 2 (which may be referred to herein as a processor, and which may be or include any electrical features, circuitry, or the like suitable for performing the functions described herein with respect to the processing unit, such as, for example, a computer including a microcontroller) may also be arranged in the environment to be monitored, in another adjacent environment or be of remote type, also arranged at a great distance from the environment to be monitored. The connection between the processing unit 2 and the aforementioned pressure sensor 4, electrostatic charge sensor 6 and accelerometer 7 may be implemented via cable or wireless, according to any of the available technologies.

[0040] The processing unit 2 is configured to receive, and receives during use: a signal S.sub.Q, correlated to a variation of the environmental electrical charge in the monitored environment, from the electrostatic charge variation sensor 6; a signal S.sub.A, indicative of a vibration detected in the environment monitored by the accelerometer 7; and a signal S.sub.P, indicative of a pressure (or variation of pressure) detected in the environment monitored by the accelerometer 7.

[0041] The pressure sensor 4 is arranged in the environment wherein it is desired to detect a human presence, or operatively coupled to this environment, to detect a variation of environmental pressure caused for example by the opening of a door or a window, or indicative of the entry of a foreign subject into this environment. In this case, therefore, the environment to be monitored is a closed environment, such as, for example, a room in an apartment or home. In fact, it should be remembered that the system 1 according to the present disclosure has the objective of identifying an unwanted entry into an environment to be protected, in particular for anti-theft purpose. When the system 1 is operative, the pressure detected is the environmental pressure present therein, which typically varies relatively slowly between day and night hours, due to air heating, or in conjunction with a variation in weather/climatic conditions. Any significant disturbance of this pressure may be indicative of an infringement.

[0042] Similarly, the accelerometer 7 is also arranged in the environment wherein it is desired to detect a human presence, to detect any vibrations in this environment that might be correlated to an intruder's footsteps, in particular footsteps caused by a person who has entered such an environment.

[0043] Similarly, the electrostatic charge variation sensor 6 is also arranged in the environment wherein it is desired to detect the human presence, or operatively coupled to this environment, to detect a variation of the environmental electrostatic charge caused by the entry of a foreign subject into this environment.

[0044] The analysis of the signals generated by the aforementioned sensors, and their suitable combination, allows the entry, in the environment to be monitored, of a subject or an intruder to be detected, discriminating with respect to false positives.

[0045] FIG. 2 illustrates an exemplary and non-limiting embodiment of the electrostatic charge variation sensor 6. The electrostatic charge variation sensor 6 comprises a pair of input terminals 8a, 8b, coupled to input electrodes E1, E2, respectively.

[0046] The two electrodes may be connected to a differential input (i.e., to the positive/negative “+”/“−” input pair of an amplifier stage or an ADC converter). Particular cases of this general configuration (which do not require changes to the electric diagram of FIG. 2) may provide for the use of an electrode (e.g., E1) of predominant dimensions with respect to the other electrode (e.g., E2) to the point of making this second electrode (E2), and the environmental charge variation detected thereby, completely negligible; in other cases, the second electrode (E2) may be removed.

[0047] In this embodiment, one of the electrodes E1, E2 (e.g., E2) is coupled to a reference potential, having constant value (e.g., common mode voltage, or VCM, typically half of the supply voltage of the device), while the other electrode of the electrodes E1, E2 (e.g., E1) is, for example, made of conductive material and coated with an insulating layer. The geometry of the electrode E1 determines the sensitivity which, as a first approximation, is proportional to the surface of the electrode itself. In an exemplary embodiment, the electrode E1 sensitive to the environmental charge is of square shape, with a side equal to about 2-10 cm, for example 5 cm. Other examples comprise electrodes made using conductive wires coated with insulator, having a length equal to a few centimeters or a few tens of centimeters, e.g., 10 cm-20 cm.

[0048] In particular, the input electrodes E1, E2 are arranged in the environment wherein it is desired to detect a human presence, while the rest of the electrostatic charge variation sensor 6 (e.g., the amplification stage, described hereinafter) may also be arranged outside the environment to be monitored, or inside this environment, indifferently.

[0049] The pair of input terminals 8a, 8b receives an input voltage Vd (differential signal), being supplied to an instrumentation amplifier 12. In a per se known manner, a human presence generates a variation of the environmental electrostatic charge which, in turn, after having been detected by the electrode E1, generates the input voltage Vd.

[0050] The instrumentation amplifier 12 comprises, in an exemplary embodiment, two operational amplifiers OP1 and OP2 and a biasing stage (buffer) OP3 having the function of biasing the instrumentation amplifier 12 to a common mode voltage VCM.

[0051] The inverting terminal of the amplifier OP1 is connected to the inverting terminal of the amplifier OP2 by means of a resistor R.sub.2 having a voltage equal to the input voltage

[0052] Vd at its ends; therefore, a current equal to I.sub.2=Vd/R.sub.2 will flow through this resistor R.sub.2. This current I.sub.2 does not come from the input terminals of the operational amplifiers OP1, OP2 and therefore flows through the two resistors R.sub.1 connected between the outputs of the operational amplifiers OP1, OP2, in series with the resistor R.sub.2; the current I.sub.2, therefore flowing through the series of the three resistors R.sub.1−R.sub.2-R.sub.1, produces an output voltage Vd′ given by Vd′=(2R.sub.1+R.sub.2)I.sub.2=(2R.sub.1+R.sub.2)Vd/R.sub.2. Therefore, the overall gain of the circuit of FIG. 2 is Ad=Vd′/Vd=(2R.sub.1+R.sub.2)/R.sub.2=+2R.sub.1/RR.sub.2. The differential gain depends on the value of the resistor R2 and may therefore be modified by acting on the resistor R.sub.2.

[0053] The differential output Vd′, therefore being proportional to the potential Vd between the input electrodes 8a, 8b, is provided at input to an analog-to-digital converter 14, which provides at output the charge variation signal S.sub.Q for the processing unit 2. The charge variation signal S.sub.Q is, for example, a high-resolution digital stream (16 bits or 24 bits). The analog-to-digital converter 14 is optional, since the processing unit 2 may be configured to work directly on the analog signal, or may itself comprise an analog-to-digital converter for converting the signal Vd′.

[0054] Alternatively, in presence of the analog-to-digital converter 14, the instrumentation amplifier 12 may be omitted, so that the analog-to-digital converter 14 receives the differential voltage Vd between the electrodes E1, E2 and samples this signal Vd directly.

[0055] The pressure sensor 4 is for example a pressure sensor made using MEMS technology. Examples of pressure sensors usable in the context of the present disclosure include pressure sensors with a measuring range of 200 mbar-2000 mbar and with accuracy (absolute precision) of a few mbar units; however, operating around the environmental pressure, approximately 1000 mbar, and observing relative values around it, a relevant parameter is the ability to detect variations around a working point, that is high resolution over time and amplitude and low inherent noise. Examples of sensors for this purpose include sensors with a resolution of 1 Pascal ( 1/100 of mbar), data rate equal to 200 Hz, RMS noise level equal to 0.5 Pascal (without filters applied).

[0056] In respective embodiments, other pressure sensors (other than MEMS sensors) are however usable.

[0057] The vibration sensor 7 is, as mentioned and in one embodiment, an accelerometer, for example of three-axis or six-axis type made using MEMS technology, or a sensor including the combination of accelerometer and gyroscope.

[0058] FIG. 3A illustrates an example of pressure signal Sp (raw signal) generated by the pressure sensor 4. The abscissa axis is the time, the ordinate axis is the absolute pressure values in millibars. As noted from FIG. 3A, there is a background noise and, at the time instant t=21 s, a peak 15 which differs considerably from this background noise, caused by a variation of environmental pressure, for example due to the opening of a door.

[0059] FIG. 3B illustrates an example of an electrostatic charge variation signal S.sub.Q (raw signal), generated by the electrostatic charge variation sensor 6. The abscissa axis is the time axis (in seconds, using the same time scale of FIG. 3A). The ordinate axis is the amplitude of the signal, LSB (“Least Significant Bit”), that is the minimum digital value at the output of the analog-to-digital converter, which is proportional to the voltage detected at the input electrode E1. Typically, 1 LSB corresponds to a value comprised between a few μV and a few tens of μV. The proportionality constant (or sensitivity) depends on the gain of the amplifier, on the resolution of the analog-to-digital converter and on any digital processing (e.g., oversampling, decimation, etc.). The representation in LSB is common in the art and disregards a quantification in physical units, as the aim is typically to detect relative variations, with respect to a steady or base state. The time relative to the start of the measurement is represented on the abscissa axis of the charge variation signal S.sub.Q. Being the sampling frequency, in the illustrated example, equal to 200 Hz, 200 samples correspond to each second reported on the abscissas.

[0060] As can be seen from FIG. 3B, the electrostatic charge variation sensor 6 detects the presence of the subject in the environment with a certain delay (here, about 2 seconds) with respect to the pressure sensor 4. The delay is due to the fact that, in the illustrated example, the steps in the monitored environment are not immediately successive the opening of the door; if they were, this delay would consequently be reduced up to zero in the event that the opening of the door coincides with the execution of steps in the monitored environment. The charge variation signal S.sub.Q shows a series of positive and negative peaks, which follow each other, and which identify the type of movement (here, in particular, the steps) executed by the subject in the environment. In particular, five steps may be identified, identified by a positive peak and an immediately following negative peak, delimited by respective rectangles 17 in dashed line in FIG. 3B.

[0061] FIG. 3C illustrates an example of a vibration signal SA generated by the accelerometer 7 and partially processed in order to generate the modulus of the axial sensing components. The abscissa axis is the time axis (in seconds, using the same time scale as FIGS. 3A and 3B), while the ordinate axis is the amplitude of the vibration signal S.sub.A in LSB. In this example an accelerometer having three detection axes was used, configured to detect three signals S.sub.Ax, S.sub.Ay, S.sub.Az respectively along the X-, Y-, Z-axes of a Cartesian triaxial reference system. Since the orientation of the accelerometer 7 in the environment to be monitored is not predictable in advance, according to the present disclosure the signal S.sub.A is generated by calculating the modulus of the acceleration, through the three components S.sub.Ax, S.sub.Ay, S.sub.Az, detected on each of the three axes of the accelerometer, according to the operation:

S.sub.A=√{square root over (S.sub.Ax.sup.2+S.sub.Ay.sup.2+S.sub.Az.sup.2)}

[0062] As can be seen from FIG. 3C, the accelerometer 7 detects vibrations, in the form of a plurality of close positive and negative peaks, which differ from the background noise, identifying the respective steps of the subject in the environment. In particular, five steps may be identified, substantially concurrent to the steps identified by the electrostatic charge variation sensor 6, delimited by respective rectangles 18 in dashed line in FIG. 3C. In order to be able to process the signals of FIGS. 3A-3C for the identification and extraction of the relevant characteristics to identify the presence of the subject in the environment, an aspect of the present disclosure provides for eliminating the background component, similar to the average value (DC) of the signal, also called “baseline.” Known-type algorithms may be used to remove the baseline or the background signal, for example based on the calculation of the average value of the raw signal and subtraction of this average value from the raw signal; alternatively, it is possible to use algorithms or methods specifically provided for this purpose, as for example illustrated hereinafter with reference to FIGS. 11 and 12.

[0063] FIGS. 4A, 4B and 4D respectively illustrate the signals of FIGS. 3A, 3B and 3C after respective processing aimed at removing the background component, or baseline. Hereinafter in the present description, the same references S.sub.P, S.sub.Q and S.sub.A will be used for the raw signals of FIGS. 3A-3C and for the signals subject to processing of FIGS. 4A, 4B and 4D, since the teaching of the present disclosure applies indifferently using raw or processed signals.

[0064] In particular: FIG. 4A illustrates the pressure signal Sp of FIG. 3A after the removal of the baseline; FIG. 4B illustrates the electrostatic charge variation signal S.sub.Q of FIG. 3B after the removal of the relative baseline; FIG. 4C illustrates the first derivative S.sub.Q′ of the electrostatic charge variation signal S.sub.Q of FIG. 4B; and FIG. 4D illustrates the variations of the signal of FIG. 3C with respect to the average value of this signal (i.e., the “AC component” of the vibration signal S.sub.A of FIG. 3C).

[0065] A method is now described, with reference to FIGS. 5A and 5B, for identifying significant variations in the electrostatic charge variation signal S.sub.Q and in the first derivative S.sub.Q′ thereof, that is variations of such signals that may be correlated to or associated with the human presence in the monitored environment and more particularly for verifying whether the signals generated by the sensor are similar to steps executed by a human subject. The portion of the signal S.sub.Q of FIG. 5A is an enlarged portion of a part of the signal S.sub.Q of FIG. 4B, in particular the portion 18a delimited by a dashed line in FIG. 4B. The portion of the signal S.sub.Q′ of FIG. 5B is an enlarged portion of a part of the signal S.sub.Q′ of FIG. 4B, in particular the portion 18b delimited by a dashed line in FIG. 4B.

[0066] The signal portions of FIGS. 5A and 5B have a plurality of positive and negative peaks, which follow each other with a certain periodicity. For the purpose of the present disclosure, positive and negative peaks are identified. This operation may be performed using a peak finding algorithm of known type, for example based on the comparison with predetermined or adaptive thresholds, or other algorithms specifically provided for this purpose, as described for example with reference to FIG. 13.

[0067] With reference to FIGS. 5A and 5B, the following peaks are identified, which follow each other in temporal order (both FIGS. 5A and 5B are represented with respect to the same time axis on the abscissas). The indicated values of the time instants, as well as the amplitude values of the peaks, are purely exemplarily and are not limiting of the present disclosure.

[0068] P1: is temporally the first identified peak, here a positive peak, which occurs on the first derivative signal S.sub.Q′ at about 24.3 s and has an amplitude value equal to about +30000 LSB.

[0069] P2: is temporally the second identified peak, here a positive peak, which occurs on the signal S.sub.Q at about 24.4 s and has an amplitude value equal to about +42000 LSB.

[0070] P3: is temporally the third identified peak, here a negative peak, which occurs on the first derivative signal S.sub.Q′ at about 24.55 s and has an amplitude value equal to about −38000 LSB.

[0071] P4: is temporally the fourth identified peak, here a negative peak, which occurs on the signal S.sub.Q at about 24.65 s and has an amplitude value equal to about −65000 LSB.

[0072] P5: is temporally the first identified peak, here a positive peak, which occurs on the first derivative signal S.sub.Q′ at about 24.75 s and has an amplitude value equal to about +18000 LSB.

[0073] As apparent to the skilled in the art, the positive peak P1 present on the signal S.sub.Q′ of the first derivative identifies a rising edge of the signal S.sub.Q which culminates at the positive peak P2 of the signal S.sub.Q; similarly, the negative peak P3 present on the signal S.sub.Q′ of the first derivative identifies a falling edge of the signal S.sub.Q which culminates at the negative peak P4 of the signal S.sub.Q; then, the signal S.sub.Q starts to grow again, and this new rising edge is identified by the positive peak P5 present on the signal S.sub.Q′ of the first derivative. Therefore, the aforementioned assessment of the succession of positive and negative peaks of the signals S.sub.Q and S.sub.Q′ has the function of identifying or detecting a specific trend of the signal generated by the electrostatic charge variation sensor 6, which the Applicant has identified as specific or significant of the presence (in particular, of the execution of a step) of a human subject in the monitored environment.

[0074] Summarizing, considering the time course of the electrostatic charge variation signal S.sub.Q and of the first derivative S.sub.Q′ of the electrostatic charge variation signal S.sub.Q in conjunction, the following time succession of positive and negative peaks is observed: [0075] 1. positive peak S.sub.Q′ [0076] 2. positive peak S.sub.Q [0077] 3. negative peak S.sub.Q′ [0078] 4. negative peak S.sub.Q [0079] 5. positive peak S.sub.Q′

[0080] However, the Applicant notes that, with a different arrangement of the electrodes E1, E2, the aforementioned time sequence (succession) may be inverted, that is the following time succession is observed: [0081] 1. positive peak S.sub.Q [0082] 2. positive peak S.sub.Q′ [0083] 3. negative peak S.sub.Q [0084] 4. negative peak S.sub.Q′ [0085] 5. positive peak S.sub.Q

[0086] The configuration of the electrodes may in fact have an effect on the detection of an electrostatic charge variation signal. While the geometry (firstly the surface) and the materials of the electrodes determine the sensitivity of the same, their arrangement in space and their distance affects the directionality or the cancellation ability of certain unwanted signal sources. On this last point it is noted that the two electrodes E1, E2 are coupled to the differential inputs of the differential amplifier (also referred to as instrumentation amplifier) or of the analog-to-digital converter (A/D or ADC); this stage performs the difference of the signals found at the “+” and “−” inputs of the amplifier. Therefore, by suitably dimensioning and positioning the electrodes, the common mode signals, that is, those signals that occur at both inputs with the same intensity, may be cancelled (or attenuated). Based on this, embodiments of the present disclosure include configurations with a single electrode, with two electrodes being equal but spaced apart from each other, with two electrodes having different geometries, etc. If the most stressed input is the “+” one, the signal shown in the figure is found; vice versa in the case of greater stress of the input “−.” In this context, the most stressed electrode is the one that detects potential variations (due to a variation of charge in the environment) that are more intense with respect to the other electrode; this may happen due to a different geometry and/or to a different installation point of the two electrodes.

[0087] The Applicant has verified that, when the above identified time succession (one of the two) is observed, it can be concluded that the signal portion 18a of FIG. 4B (and the first derivative 18b of FIG. 4C) is generated by a step of a human subject in the monitored environment.

[0088] In order to identify whether a variation of the signals S.sub.Q and S.sub.Q′ is one of the peaks sought, respective thresholds (positive or negative) A1.sub.TH-A5.sub.TH are provided, to be compared with the trend of the signals S.sub.Q and S.sub.Q′.

[0089] The thresholds A1.sub.TH-A5.sub.TH have a predefined/default value, identified empirically on the basis of observations of the trend of the signals S.sub.Q and SQ.sup.′, and for example are defined as follows: [0090] threshold A1.sub.TH: chosen as a fraction (e.g., between ½ and ⅙) of the maximum value reachable (or known to be such on the basis of experiments) by the first peak P1; by way of example, in this embodiment described it has a value chosen in the range 8000-12000 LSB (value expressed in modulus), in particular 10000 LSB. [0091] threshold A2.sub.TH: chosen as a fraction (e.g., between ½ and ⅙) of the maximum value reachable (or assumed to be such on the basis of experiments) by the second peak P2; by way of example, in this embodiment described it has a value chosen in the range 8000-12000 LSB (value expressed in modulus), in particular 10000 LSB. [0092] threshold A3.sub.TH: chosen as a fraction (e.g., between ½ and ⅙) of the maximum value reachable (or assumed to be such on the basis of experiments) by the third peak P3; by way of example, in this embodiment described it has a value chosen in the range 6000-8500 LSB (value expressed in modulus), here in particular 7500 LSB. [0093] threshold A4.sub.TH: chosen as a fraction (e.g., between ½ and 1/9) of the maximum value reachable (or assumed to be such on the basis of experiments) by the fourth peak P4; by way of example, in this embodiment described it has a value chosen in the range 6000-8500 LSB (value expressed in modulus), here in particular 7500 LSB. [0094] threshold A5.sub.TH: chosen as a fraction (e.g., between ½ and ⅕) of the maximum value reachable (or assumed to be such on the basis of experiments) by the fifth peak P5; by way of example, in this embodiment described it has a value chosen in the range 6000-8500 LSB (value expressed in modulus), in particular 7500 LSB.

[0095] In the example of FIGS. 5A and 5B, the thresholds have the following value: threshold A1.sub.TH: +10000 LSB; threshold A2.sub.TH: +10000 LSB; threshold A3.sub.TH: −7500 LSB; threshold A4.sub.TH: −7500 LSB; threshold A5.sub.TH: +7500 LSB.

[0096] Alternatively to what has been described, the value of the thresholds A1.sub.TH-A5.sub.TH may be chosen as a function of the background noise of the respective signals S.sub.Q and S.sub.Q′, for example equal to 8-12 times (e.g., 10 times) the average value of the noise.

[0097] The following comparisons are then performed for each threshold: [0098] the amplitude A1, in LSB, of the peak P1, exceeds, in positive, the threshold A1.sub.TH for P1 to be identified as a “positive peak”; [0099] the amplitude A2, in LSB, of the peak P2, exceeds, in positive, the threshold A2.sub.TH for P2 to be identified as a “positive peak”; [0100] the amplitude A3, in LSB, of the peak P3, exceeds, in negative, the threshold A3.sub.TH for P3 to be identified as a “negative peak”; [0101] the amplitude A4, in LSB, of the peak P4, exceeds, in negative, the threshold A4.sub.TH for P4 to be identified as a “negative peak”; and [0102] the amplitude A5, in LSB, of the peak P5, exceeds, in positive, the threshold A5.sub.TH for P5 to be identified as a “positive peak.”

[0103] To improve the robustness of the method proposed herein, by improving the discrimination between actual step and environmental noise or other perturbation, it is possible, again with reference to FIGS. 5A and 5B, to define and monitor the following additional parameters: [0104] T1: time interval between the positive peak P2 and the negative peak P4 of the electrostatic charge variation signal S.sub.Q. [0105] T2: time interval between the positive peak P2 of the electrostatic charge variation signal S.sub.Q and the positive peak P1 of the first derivative signal S.sub.Q′. [0106] T3: time interval between the positive peak P2 of the electrostatic charge variation signal S.sub.Q and the negative peak P3 of the first derivative signal S.sub.Q′. [0107] T4: time interval between the negative peak P4 of the electrostatic charge variation signal S.sub.Q and the negative peak P3 of the first derivative signal S.sub.Q′. [0108] T5: time interval between the negative peak P4 of the electrostatic charge variation signal S.sub.Q and the positive peak P5 of the first derivative signal S.sub.Q′. [0109] T6: time interval between the positive peak P1 and the negative peak P3 of the first derivative signal S.sub.Q′. [0110] T7: time interval between the negative peak P3 and the positive peak P5 of the first derivative signal S.sub.Q′. [0111] T8: time interval between the positive peak P1 and the positive peak P5 of the first derivative signal S.sub.Q′.

[0112] The existence of the following relationships is verified:

T1=T3+T4

T6=T2+T3

T7=T4+T5

T8=T6+T7

[0113] Additionally, or alternatively, the existence of the following relationships is verified, to verify that the duration of the time intervals T2-T5 complies with that expected for the signal shape that may be associated with a step of a subject: [0114] T2.sub.TH_L<T2<T2.sub.TH_H, where T2.sub.TH_L and T2.sub.TH_H represent the boundaries of a range of time values within which T2 needs to be comprised (e.g., T2.sub.TH_L=30-70 ms and T2.sub.TH_H=150-250 ms); [0115] T3.sub.TH_L<T3<T3.sub.TH_H, where T3.sub.TH_L and T3.sub.TH_H represent the boundaries of a range of time values within which T3 needs to be comprised (e.g., T3.sub.TH_L=30-70 ms and T3.sub.TH_H=150-250 ms);

[0116] T4.sub.TH_L<T4<T4.sub.TH_H, where T4.sub.TH_L and T4.sub.TH_H represent the boundaries of a range of time values within which T4 needs to be comprised (e.g., T4.sub.TH_L=30-70 ms and T4.sub.TH_H=150-250 ms); and

[0117] T5.sub.TH_L<T5<T5.sub.TH_H, where T5.sub.TH_L and T5.sub.TH_H represent the boundaries of a range of time values within which T5 needs to be comprised (e.g., T5.sub.TH_L=30-70 ms and T5.sub.TH_H=150-250 ms).

[0118] In one embodiment, the values of T1.sub.TH_L-T5.sub.TH_L are equal to each other and equal to 50 ms; and the values of T1.sub.TH_H-T5.sub.TH_H are equal to each other and equal to 200 ms.

[0119] The choice of the values of T1.sub.TH_H-T5.sub.TH_H may vary with respect to what is described herein and set on an empirical basis after experimental observations.

[0120] FIG. 6 illustrates, through a flow diagram, the method for detecting the human presence implemented by the system 1 of FIG. 1, with exclusive reference to the electrostatic charge variation signal S.sub.Q, according to what has been previously described.

[0121] At step 60 the processing unit 2 acquires the raw signal S.sub.Q from the electrostatic charge variation sensor 6.

[0122] At step 61 the raw signal S.sub.Q is processed to remove the baseline or background signal, as previously mentioned.

[0123] At step 62 the method for searching the positive/negative peaks on the electrostatic charge variation signal S.sub.Q is performed, identifying for example the time succession of the peaks P2 and P4 of FIG. 5A.

[0124] At step 63 the first derivative signal S.sub.Q′ of the electrostatic charge variation signal S.sub.Q is calculated.

[0125] Then, at step 64 the method for searching positive/negative peaks on the first derivative signal S.sub.Q′ is performed, identifying for example the time succession of the peaks P1, P3 and P5 of FIG. 5B.

[0126] The aforementioned conditions are then assessed on the amplitude A1-A5 of the detected peaks and on the time intervals T2-T5. The method proposed herein is performed in real time, that is by acquiring the samples of the raw signal S.sub.Q and assessing the conditions previously described as these samples are generated by the electrostatic charge variation sensor 6.

[0127] A counter P.sub.COUNT (initialized for example to zero) stores the number of identified peaks (five peaks P1-P5 may be utilized and, in some embodiments, may be required to confirm the identification of a step in this embodiment). At an initial instant where no peaks have yet been detected, P.sub.COUNT=0.

[0128] With reference to blocks 65-69 of FIG. 6, the value of the counter P.sub.COUNT is assessed. The increase in the value of P.sub.COUNT determines the access to the respective computational blocks 65-69, to verify the respective conditions on the peaks P1-P5 which, as previously illustrated, differ from each other in terms of amplitude thresholds and time references.

[0129] At block 65 the presence of the peak P1 in the first derivative signal S.sub.Q′ is assessed by comparing the amplitude value A1 with the respective threshold A1.sub.TH. If the comparison with the threshold determines the presence of the peak P1, then the counter P.sub.COUNT is updated (P.sub.COUNT=1) and a new datum is acquired from the raw signal S.sub.Q. Otherwise, the counter P.sub.COUNT is reset to a zero value and a new datum is acquired from the raw signal S.sub.Q.

[0130] Steps 60-64 are then performed again.

[0131] If the presence of the peak P1 has been confirmed, the assessment of the value of the counter P.sub.COUNT determines passing from step 64 to step 66, wherein the presence of the peak P2 in the electrostatic charge variation signal S.sub.Q is assessed by comparing the amplitude value A2 with the respective threshold A2.sub.TH. If the comparison with the threshold determines the presence of the peak P2, and the time conditions on the value of the interval T2 are met, such that T2.sub.TH_L<T2<T2.sub.TH_H, then the counter P.sub.COUNT is updated (P.sub.COUNT=2) and a new datum is acquired from the raw signal S.sub.Q. Otherwise, the counter P.sub.COUNT is reset to a zero value and a new datum is acquired from the raw signal S.sub.Q.

[0132] Steps 60-64 are then performed again.

[0133] If the presence of the peak P2 has been confirmed, the assessment of the value of the counter PCOUNT determines passing from step 64 to step 67, wherein the presence of the peak P3 in the first derivative signal S.sub.Q′ is assessed by comparing the amplitude value A3 with the respective threshold A3.sub.TH. If the comparison with the threshold determines the presence of the peak P3, and the time conditions on the value of the interval T3 are met, such that T3.sub.TH_L<T3<T3.sub.TH_H, then the counter P.sub.COUNT is updated (P.sub.COUNT=3) and a new datum is acquired from the raw signal S.sub.Q. Otherwise, the counter PCOUNT is reset to a zero value and a new datum is acquired from the raw signal S.sub.Q.

[0134] Steps 60-64 are then performed again.

[0135] If the presence of the peak P3 has been confirmed, the assessment of the value of the counter P.sub.COUNT determines passing from step 64 to step 68, wherein the presence of the peak P4 in the electrostatic charge variation signal S.sub.Q is assessed by comparing the amplitude value A4 with the respective threshold A4.sub.TH. If the comparison with the threshold determines the presence of the peak P4, and the time conditions on the value of the interval T4 are met, such that T4.sub.TH_L<T4<T4.sub.TH_H, then the counter P.sub.COUNT is updated (P.sub.COUNT=4) and a new datum is acquired from the raw signal S.sub.Q. Otherwise, the counter P.sub.COUNT is reset to a zero value and a new datum is acquired from the raw signal S.sub.Q.

[0136] Steps 60-64 are then performed again.

[0137] If the presence of the peak P4 has been confirmed, the assessment of the value of the counter P.sub.COUNT determines passing from step 64 to step 69, wherein the presence of the peak P5 in the first derivative signal S.sub.Q′ is assessed by comparing the amplitude value A5 with the respective threshold A5.sub.TH. If the comparison with the threshold determines the presence of the peak P5, and the time conditions on the value of the interval T5 are met, such that T5.sub.TH_L<T5<T5.sub.TH_H, then the analysis of the relative portion 18a, 18b of the signals S.sub.Q and S.sub.Q′ is concluded and a warning, or trigger signal, may be generated, which confirms the identification of a step in the signal generated by the electrostatic charge variation sensor 6.

[0138] The counter P.sub.COUNT is reset and a new datum is acquired from the raw signal S.sub.Q, to identify the presence of a successive step.

[0139] The confirmation of human presence in the environment occurs, according to an aspect of the present disclosure, after the identification of a plurality of steps, for example of five steps. However, it is apparent that, in order to speed up the detection, it is possible to report the presence of the subject even only after the identification of a single step.

[0140] As previously mentioned, to generate the actual alarm or final confirmation of the human presence, the present disclosure provides for the joint analysis of the signals S.sub.P, S.sub.A generated by the pressure sensor 4 and by the vibration sensor 7.

[0141] FIG. 7 illustrates, through a block diagram, steps of a method for analyzing the pressure signal S.sub.P.

[0142] In one embodiment, the algorithm of FIG. 7 operates in real time, similarly to the method of FIG. 6, that is the data is processed during the same acquisition step. It is assumed that the pressure signal has been converted into digital and therefore, hereinafter, the term “datum” identifies a digital value of the pressure signal S.sub.P (e.g., pressure value in mbar).

[0143] At each iteration, after the acquisition of the pressure signal S.sub.P (step 70), the i-th pressure datum Pi (amplitude value) is deducted of its baseline (step 71) and is stored in a buffer P.sub.BUFF (step 72); at the same time, or at previous or subsequent time instants, indifferently, the search for possible peaks in the pressure signal S.sub.P (step 73) is carried out, using algorithms known for this purpose, or specifically provided for this purpose. If a peak is detected (step 74, output YES), the value PK25 equal to 25% of the amplitude of the detected peak is calculated (step 75) (this percentage value is indicative and may vary for example in the range 10%-50%). Iteratively, this value PK25 is subtracted (step 76) to each pressure datum (i-th datum PK.sub.i) contained in the buffer P.sub.BUFF (operation PK.sub.i-PK25). If the value resulting from this subtraction is positive (step 77, output YES), then this value is added to a variable P.sub.AREA (representative of the area of the plane portion, between the peak and 25% of its value), to perform a calculation, in digital form, of the integral of the signal around the detected peak (step 79). The integral may be calculated as an area (variable A in step 79) subtended by the signal relating to the peak, that is in digital format by adding the amplitude values P.sub.i-PK25, only if this difference is greater than 0. At step 76 the value PK25 is in fact subtracted to each sample Pi; if the result of this operation of step 76 is positive, then this result is added to the previous area value A (where A is initialized to 0 at the beginning of the method); if the result of step 76 is negative, this result is ignored. This operation of addition is carried out for a maximum of N iterations; the count of these N iterations is carried out by increasing an index j, regardless of the value of the result of the operation of step 76 (the increase of j allows to go through the entire buffer 72).

[0144] The steps 76,77,78,79 have the function of quantifying the portion of the area subtended by the curve only in the presence of a peak, so as to be capable of performing the operation, of the successive step 80, of assessing the peak itself.

[0145] Finally, the ratio R.sub.PK between P.sub.AREA and PK.sub.i-PK25 is calculated (step 80) (obtaining a result which is greater than 1), which is indicative of the “steepness” of the peak: the smaller the value of this ratio R.sub.PK, the greater the steepness and vice versa. The value of the ratio R.sub.PK is compared with a threshold RP.sub.THRES (step 81): if R.sub.PK<R.sub.PTHRES then the peak is sufficiently steep to be similar to that generated by the opening of a door, and a signal, or trigger, indicative of this event is generated (step 82); vice versa, the method returns to step 70 by resetting the variables j and A. The choice of the threshold RP.sub.THRES includes, for example, values comprised between 10 and 30; the smaller this value, the steeper and more time-limited the detected peak.

[0146] For greater clarity of the operation of the method of FIG. 7, FIG. 8A illustrates a pressure signal Sp having the threshold PK25 and a peak PKi graphically shown thereon. Here, the peak PKi has an amplitude value equal to 0.215 mbar and therefore the threshold is PK25=0.054 mbar.

[0147] FIG. 8B illustrates the signal of FIG. 8A after step 76, wherein the operation of subtracting the value PK25 from each datum of the signal S.sub.P is performed. After this operation, the value of the peak Pki is equal to 0.16 mbar, the value of the ratio R.sub.PK is equal to 9.44, the value of the area P.sub.AREA is equal to 1.51, and the threshold RP.sub.THRES is set to 15. The assessment of step 81 gives therefore a positive result, as R.sub.PK<RP.sub.THRES.

[0148] FIG. 8C illustrates a signal wherein the door opening event is not confirmed/recognized, since R.sub.PK>RP.sub.THRES. In this example, after the operation of step 76, the amplitude of the peak PKi is equal to 0.22, the value PK25 is 0.072, the value of the area P.sub.AREA is 8.56, the ratio R.sub.PK is equal to 38.9 and the threshold RP.sub.THRES is set to 15.

[0149] FIG. 9 illustrates, through a block diagram, steps of a method for analyzing the vibration signal S.sub.A.

[0150] In one embodiment, the algorithm of FIG. 9 operates in real time, similarly to the method of FIG. 6 or FIG. 7, that is the data is processed during the same acquisition step. It is assumed that the vibration signal has been converted into digital and therefore, hereinafter, the term “datum” identifies a digital value of the vibration signal S.sub.A (e.g., signal amplitude in LSB).

[0151] At each iteration, after having acquired the vibration signal relating to the detection axes of the accelerometer (S.sub.Ax, S.sub.Ay, S.sub.Az) by the processing unit 2 (step 90), the modulus XLM (i.e., the signal S.sub.A previously discussed) of the acceleration is calculated (step 91) on the basis of the signals acquired from the three axes of the accelerometer (assuming here to use a triaxial accelerometer).

[0152] The AC component is then obtained (step 92) (i.e., amount correlated to the variation of the signal with respect to the average value of this signal, whose i-th value is indicated as XLPKi), which is stored in a buffer XLAC.sub.BUFFER; at the same time the search for possible signal peaks is carried out on such data (step 93). If a peak is detected (step 93, output YES), the value XLPK25 (step 95) is calculated equal to 25% of the peak amplitude (this percentage value may be chosen differently, for example in the range 10%-50%). Otherwise, the method returns to step 90.

[0153] Iteratively, this value XLPK25 is subtracted to all the values XLPKi contained in the buffer XLAC.sub.BUFFER (step 96). If, for each sample, the value resulting from this subtraction is positive, then (step 97, output YES) this value is added to the variable XLA (representative of the area of the plane portion, between the peak and 25% of its value), implementing an operation of calculation of the integral in digital format (step 98).

[0154] If the result of the foregoing step is negative, this result is ignored. This operation of adding and updating the variable XLA is carried out for a maximum of N iterations; the count of these N iterations is carried out by increasing an index k, regardless of the result of the assessment of step 97 (the increase of k allows to go through the entire buffer XLAC.sub.BUFFER)

[0155] The ratio R.sub.XLPK between the area A.sub.XL and XLPKi-XLPK25 (greater than 1), indicative of the steepness of the peak, is then calculated (step 100), for each i-th datum: the smaller the value of this ratio, the greater the steepness of rising of the signal, and vice versa.

[0156] R.sub.XLPK is compared (step 101) with a threshold RXLPK.sub.THRES: if R.sub.XLPK<RXLPK.sub.THRES then the identified peak is steep and similar to that generated by a subject's step (step 103) and a suitable signal, or trigger, is generated which confirms the presence of the subject in the considered environment.

[0157] In order to increase the reliability of the proposed method, so that the vibration signal is validated as generated by a subject's steps, it is optionally possible to verify (step 102) the repetition of a certain number of peaks over time (e.g., by setting a comparison threshold CountTHRES, for example equal to 2), with the condition that a time which is longer than a predefined value T.sub.THRES (for the choice of this value, considerations similar to those previously made for the pressure signal Sp are valid) does not elapse between an event of single step and the successive.

[0158] For improved clarity of the operation of the method of FIG. 9, FIG. 10A illustrates a vibration signal S.sub.A obtained by calculating the modulus of the three detection components of a triaxial accelerometer. The signal of FIG. 10A is temporally limited to the signal detected during the execution of a single step.

[0159] FIG. 10B illustrates the AC component of the signal of FIG. 10A (which actually represents the envelope of the signal of FIG. 10A). The maximum value XLPK.sub.MAX of peak amplitude, here equal to 101.9 LSB, and the calculated value of XLPK25=25.5 LSB are indicated on the signal of FIG. 10B.

[0160] FIG. 10C illustrates the signal resulting from the operation XLPKi−XLPK25 (performed for each i-th sample of the signal S.sub.A, wherein the value of the peak of FIG. 10B is now equal to 76.4, the value of the area A.sub.XL, is equal to 517.9, the value of the ratio R.sub.XLPK is equal to 6.78 and the threshold RXLPK.sub.THRES is set to the value 15. The assessment R.sub.XLPK<RXLPK.sub.THRES therefore gives a positive result, as R.sub.XLPK<RXLPK.sub.THRES.

[0161] With reference to FIGS. 11 and 12, respective methods for the removal of the baseline, applicable in the context of the present disclosure, are now described.

[0162] With reference to FIG. 11, the algorithm operates in real time, similarly to the method of FIG. 6 or FIG. 7, that is the data is processed during the same acquisition step. It is assumed that the signal received at the input (which may be any of the modulus of the vibration signal S.sub.A, the pressure signal S.sub.P and the electrostatic charge variation signal S.sub.Q) has been converted into digital and therefore, hereinafter, the term “datum” identifies a sample or digital value of the considered signal (e.g., signal amplitude in LSB or pressure value in mbar).

[0163] At each iteration, the following operations are performed.

[0164] If the input datum xi (i-th datum) is comprised between the thresholds BL.sub.THRES_H and BL.sub.THRES_L (step 110, output YES) the datum xi is accumulated (step 111) in a shift buffer having size N.sub.BLBUFF (e.g., N.sub.BLBUFF=10).

[0165] In the early iteration steps (first start of the algorithm), the thresholds BL.sub.THRES_H and BL.sub.THRES_L are ignored (i.e., the output from block 110 is “YES”), until the buffer is completely filled (a number of iterations equal to N.sub.BLBUFF are utilized and in some embodiments may be required). In other words, all the incoming samples xi will fill the buffer, as identified by the arrow 110a in dashed line.

[0166] A variable BL that stores the current baseline value is then updated with a value equal to the average value of the samples present in the buffer (step 112a) at the same time, the standard deviation value of the samples present in the buffer is calculated (step 112b). New thresholds BL.sub.THRES_H and BL.sub.THRES_L are calculated (step 113), respectively equal to the value of the variable BL increased and decreased by a multiple amount of the standard deviation of the samples present in the buffer. The parameter k adjusts the width of the band defined by the two thresholds BL.sub.THRES_H and BL.sub.THRES_L: the greater the value of k, the greater the variations of the input datum that will be absorbed in the baseline. The variable k is chosen, for example, in the range 3-6.

[0167] After calculating, for each input sample xi, the respective baseline value BL, the output datum yi=xi=BL is generated (step 114), that is the datum that will form the respective vibration S.sub.A, pressure S.sub.P or electrostatic charge signal S.sub.Q deprived of the respective baseline.

[0168] If the input datum is not comprised between the thresholds BL.sub.THRES_H and BL.sub.THRES_L (step 110, output NO), the baseline and the thresholds are not modified. The output datum yi is, in any case, equal to the input value xi deducted of the value BL calculated as the average of the samples present in the buffer. It is repeated that the operations of steps 112a, 112b are not performed until the buffer is completely filled.

[0169] FIG. 12 illustrates a further method, alternative to that of FIG. 11, for calculating the baseline and subtracting the same from the respective signal.

[0170] The algorithm operates in real time, similarly to the method of FIG. 11, that is the data is processed during the same acquisition step. It is assumed that the signal received at the input (which may be any of the modulus of the vibration signal S.sub.A, the pressure signal S.sub.P and the electrostatic charge variation signal S.sub.Q) has been converted into digital and therefore, hereinafter, the term “datum” identifies a sample or digital value of the considered signal (e.g., signal amplitude in LSB or pressure value in mbar).

[0171] At each i-th iteration, the i-th datum xi of the respective signal is acquired by the processing unit 2 (step 120). Then, the first derivative xi′ is calculated (step 121). Then, the absolute value |xi′51 of the first derivative xi′ is calculated (step 122). The absolute value |xi′| which is calculated is then entered into a buffer (step 123) having size N.sub.BLBUFF′(for example equal to 10).

[0172] If (step 124) all the values contained in the buffer are lower than a threshold BL.sub.THRES then (output YES from step 124) the input data xi is entered into a second buffer with size M.sub.BLBUFF (step 125). Being a derivative, exceeding the threshold BL.sub.THRES is indicative of the rate at which the signal increases (or decreases). This value depends on the type of quantity analyzed, on the data rate and on the “noisiness” of the environment and of the sensor itself. For example, in the case of the charge variation signal, the threshold BL.sub.THRES may be comprised between 8000 and 16000.

[0173] The baseline BL is then updated to the new value, given by the average of the elements present in this second buffer (step 126).

[0174] After calculating, for each input sample xi, the respective baseline value BL, the output datum yi=xi−BL is generated (step 127), that is the datum that will form the respective vibration S.sub.A, pressure Sp or electrostatic charge signal S.sub.Q deprived of the respective baseline.

[0175] If at least one element of the first buffer exceeds or is equal to the threshold BL.sub.THRES, the baseline variable BL is not updated (output NO from step 124).

[0176] The output value yi is, however, equal to the input value xi deducted of the value BL.

[0177] At the first start, the algorithm ignores the check of the threshold BL.sub.THRES for a number of iterations sufficient to completely fill the first buffer of size N.sub.BLBUFF′ In this starting condition, all the input samples |xi′| are used to fill this first buffer and the generation of an output datum yi is not performed.

[0178] FIG. 13 illustrates, through a block diagram, steps of a peak finding method, usable in the context of the present disclosure to recognize positive and negative peaks, for example applicable in the context of steps 62 and 73 previously described with reference to FIGS. 6 and 7, respectively.

[0179] With reference to FIG. 13, the algorithm operates in real time, similarly to the method of FIG. 6 or FIG. 7, that is the data is processed during the same acquisition step. It is assumed that the signal received at the input (which may be any of the modulus of the vibration signal S.sub.A, the pressure signal Sp and the electrostatic charge variation signal S.sub.Q) has been converted into digital and therefore, hereinafter, the term “datum” identifies a sample or digital value, or sample, of the considered signal (e.g., signal amplitude in LSB or pressure value in mbar).

[0180] With reference to the algorithm of FIG. 13, the following variables are defined and used:

[0181] xi=amplitude in LSB or pressure value in mbar of the current datum (sample) (i-th datum);

[0182] ti=time instant associated with the current datum xi;

[0183] 2N+1=width, expressed in number of samples, of a considered peak (comprehensive of the signal portion rising towards a maximum peak value, the maximum value reached, and the signal portion falling from the maximum value);

[0184] PK.sub.THRES=comparison threshold to detect the presence of a positive peak;

[0185] VL.sub.THRES=comparison threshold to detect the presence of a negative peak;

[0186] xj=local maximum reached by the positive peak;

[0187] xk=local minimum reached by the negative peak;

[0188] pka=amplitude in LSB or pressure value in mbar of the local maximum reached by the considered positive peak;

[0189] pkt=time instant associated with the local maximum pka reached by the considered positive peak;

[0190] vla=amplitude in LSB or pressure value in mbar of the local minimum reached by the considered negative peak;

[0191] vlt=time instant associated with the local minimum vla reached by the considered negative peak;

[0192] PKF=variable, or “flag,” which identifies the “positive-peak-found” event;

[0193] VLF=variable, or “flag,” which identifies the “negative-peak-found” event. At each iteration, the amplitude and the time index of the input datum are entered (steps 130a and 130b), respectively, into the two buffers X.sub.PBUFF (buffer that contains the data xi) and T.sub.PBUFF (buffer that contains the data ti). Subsequently, the maximum xj and the minimum xk of all the elements of the buffer X.sub.BUFF are calculated (steps 131a and 131b).

[0194] If the time index pkt of the local maximum xj found is not equal to the value of the index N, it means that the sample corresponding to the local maximum xj is not placed in the middle of the buffer X.sub.PBUFF; in this case no peak was found and PKF=“FALSE” (output NO from step 132a).

[0195] On the contrary, if the time index pkt of the local maximum xj found is equal to the value of the index N (output YES from step 132a), it means that the sample corresponding to the local maximum xj is placed in the middle of the buffer X.sub.PBUFF; it occurs (step 133a) if the local maximum xj found is higher than PKTHRES (e.g., PK.sub.THRES chosen of a value equal to 15000 for the electrostatic charge variation signal). If so, the presence of a peak of width 2N+1 is confirmed and the variable PKF is set to “TRUE” (step 134a).

[0196] The amplitude of the peak found and confirmed is x.sub.N, and the time index is t.sub.N.

[0197] Dual considerations may be made for searching the negative peak.

[0198] In this case, if the time index vlt of the local minimum xk found is not equal to the value of the index N, it means that the sample corresponding to the local minimum xk is not placed in the middle of the buffer TPBUFF; in this case no peak was found and VLF=“FALSE” (output NO from step 132b).

[0199] On the contrary, if the time index vlt of the local minimum xk found is equal to the value of the index N (output YES from step 132b), it means that the sample corresponding to the local minimum xk is placed in the middle of the buffer TPBUFF; it occurs (step 133b) if the local minimum xk found exceeds (for negative values) the threshold VL.sub.THRES (e.g., VL.sub.THRES chosen of a value equal to −15000 for the electrostatic charge variation signal). If so, the presence of a peak of width 2N+1 is confirmed and the variable VLF is set to “TRUE” (step 134b).

[0200] The amplitude of the negative peak found and confirmed is x.sub.N, and the time index is t.sub.N.

[0201] At the first start, the algorithm is not operative for a number of iterations equal to 2N+1, that is, until the buffers X.sub.PBUFF and T.sub.PBUFF are filled. In this step all the input samples will fill the buffers and the outputs are set to PKF=“FALSE” and VLF=“FALSE.”

[0202] FIG. 14 illustrates an algorithm or method for calculating the first derivative of the signals S.sub.P, S.sub.A and S.sub.Q, applicable in the context of the present disclosure.

[0203] The algorithm of FIG. 14 operates in real time, similarly to the methods previously described, that is the data is processed during the same acquisition step. It is assumed that the signal received at the input (which may be any of the modulus of the vibration signal S.sub.A, the pressure signal S.sub.P and the electrostatic charge variation signal S.sub.Q) has been converted into digital and therefore, hereinafter, the term “datum” identifies a sample or digital value, or sample, of the considered signal (e.g., signal amplitude in LSB or pressure value in mbar).

[0204] By definition, the output y is delayed by 2 samples with respect to the input; the first derivative of the input signal, calculated at time ti, relates to the input signal at time t(i−2). The two flows are therefore temporally aligned before calculating relative time distances.

[0205] With reference to FIG. 15, a method for extracting the envelope of the considered signal (any of the signals S.sub.P, S.sub.A and S.sub.Q), or for obtaining the aforementioned “AC component” (e.g., with reference to step 92), is illustrated through a block diagram.

[0206] With reference to FIG. 15, digital samples xi of the signal being processed are acquired and stored in the buffer 150. The buffer 150, in particular, is designed to store a plurality of samples (for example 25 samples, with a sampling rate of 50 Hz and buffers with a time window of 0.5 s). The number of samples may, in any case, vary as desired, considering that the greater this number, the smoother the signal generated at the output of the block chain of FIG. 15. For example, the number of samples in the buffer 150 is chosen in the range of 10-30.

[0207] The samples stored in the buffer 150 are sent to a first input of a subtraction block 152. The other input of the subtraction block 152 receives samples which are further processed (filtered) through the branch 154, as described hereinafter.

[0208] The branch 154 first comprises a processing block 155 which uses a Hann window 156, or Hann function, which is of a per se known type and implements the following function:

[00001]$y_i=\frac{1}{2}{x}_i\left(1-\cos \frac{2\pi i}{K}\right)$

[0209] where x.sub.i=[x.sub.0, . . . , x.sub.K−1] are the samples at the input in the processing block 155 (the subscript “i=0, . . . , K−1” identifies the i-th sample) and y.sub.i=[y.sub.0, . . . , y.sub.K−1] are the samples at the output from the processing block 155.

[0210] The branch 154 comprises an averaging block 157, which receives the samples y.sub.1=[y.sub.0, . . . , y.sub.K−1] and performs an arithmetic average operation of the value of said samples.

[0211] The branch 154 further comprises a multiplication block 158, which receives at input the average value generated at the output of block 157 and performs an operation of multiplication by 2 of said average value (since the Hann window of block 156 halves the average amplitude of the signal, the attenuation introduced is compensated with this operation), generating an output which is supplied to the second input of the subtraction block 152.

[0212] At the output of the subtraction block 152, the signal at the input in the subtraction block 152 minus its average value is obtained, therefore a signal which on average oscillates around zero, without any significant offset. The output of the subtraction block 152 is then further processed through a block 159 which implements a further Hann window, as it has been described for block 156. This further Hann window 159 has the function of smoothing the signal, smoothing the peaks and discontinuities at the ends of the analysis window.

[0213] A block 160 receives at input the array generated at the output of block 159, and performs an estimation of the variance of said array, in a per se known manner. The output from block 160 is consequently scalar.

[0214] Finally, a square root operation (block 162) of the variance value has the function of compressing the dynamic range of the output signal, as well as of bringing it back to the initial physical dimensions. In other words, the variance increases according to a power of two, and the square root restores the physical dimensions. For example, for the signal S.sub.A, if the physical dimension at the input is expressed in V, after the calculation of the variance, it is expressed in V.sup.2.

[0215] The advantages achieved by the present disclosure are apparent from the foregoing description.

[0216] In particular, the following advantages are obtained with respect to the prior art: [0217] insensitivity to environmental conditions; [0218] very low consumption when compared with other technologies (infrared, microwave, etc.); [0219] small size and ease of integration and installation; [0220] unlike common detectors, provided with “lenses” or antennas that constrain the spatial shape/arrangement thereof, the present disclosure may be physically organized on the basis of the application; [0221] reduced cost with respect to the known arts.

[0222] Further variations, with respect to what has been described, may also be provided.

[0223] For example, the present disclosure may be modified with respect to what has been described by excluding one of the pressure sensor 4 and the vibration sensor 7; in this case, the confirmation of human presence in the monitored environment is provided by the analysis steps of the electrostatic charge variation signal S.sub.Q combined with only one of the vibration signal SA and the pressure signal S.sub.P. If the sensor excluded or not present is the pressure sensor, the environment wherein the presence of the subject is detected may not be a closed environment.

[0224] Furthermore, while the present disclosure has been described with explicit reference to the processing of digital signals, what has been described applies, in a per se apparent manner, to analog signals.

[0225] A system for detecting a presence in an environment to be monitored, may be summarized as including a processing unit (2); an electrostatic charge variation sensor (6), coupled to the processing unit (2), configured to detect a variation of electrostatic charge in said environment and generate an electrostatic charge variation signal (S.sub.Q); and one of a vibration sensor (7) and an environmental pressure sensor (4), wherein the vibration sensor is operatively coupled to the environment to be monitored to detect an environmental vibration in the environment to be monitored and generate a vibration signal (S.sub.A), and wherein the environmental pressure sensor (4) is operatively coupled to the environment to be monitored to detect an environmental pressure in the environment to be monitored and generate a pressure signal (S.sub.P), wherein the processing unit (2) is configured to acquire, from the electrostatic charge variation sensor (6), the electrostatic charge variation signal (S.sub.Q); detect, in said electrostatic charge variation signal (S.sub.Q), first signal characteristics indicative of the presence of a subject in said environment to be monitored; acquire, from said one of the vibration sensor (7) and the environmental pressure sensor (4), respectively the vibration signal (S.sub.A) or the pressure signal (S.sub.P); detect, in said vibration signal (S.sub.A) or pressure signal (S.sub.P) acquired, respective second signal characteristics indicative of the presence of the subject in said environment to be monitored; and generate a warning signal if both the first and the second signal characteristics have been detected.

[0226] The system may further include the other of the vibration sensor (7) and the environmental pressure sensor (4), wherein the processing unit (2) is further configured to acquire, from the other of said vibration sensor (7) and environmental pressure sensor (4), respectively the vibration signal (S.sub.A) or the pressure signal (SP); detect, in said other of the vibration signal (S.sub.A) and the pressure signal (S.sub.P) acquired, respective third signal characteristics indicative of the presence of the subject in said environment to be monitored; and generate the warning signal if all the first, the second and the third signal characteristics have been detected.

[0227] The operation of detecting the first signal characteristics may include detecting, in the electrostatic charge variation signal (S.sub.Q), the following characteristics which follow each other in temporal order: a first rising edge; a first local maximum; a first falling edge; a first local minimum; a second rising edge; alternatively, detecting, in the electrostatic charge variation signal (S.sub.Q), the following characteristics which follow each other in temporal order: a falling edge; a first local minimum; a first rising edge; a first local maximum; a second falling edge.

[0228] The operation of detecting the first signal characteristics nay further include performing a comparison of said local maximums and minimums with respective thresholds; and assessing, through comparison with respective thresholds, a value of steepness or rising rate of the first and the second rising edges and of steepness or falling rate of the falling edge.

[0229] The operation of detecting, in the electrostatic charge variation signal (S.sub.Q), the characteristics that follow each other in temporal order may include calculating the first derivative (S.sub.Q′) of the electrostatic charge variation signal (S.sub.Q); identifying, in the electrostatic charge variation signal (S.sub.Q) and in the first derivative signal (S.sub.Q′), a respective plurality of positive and negative peaks; detecting one of the following time series a) and b) with which said plurality of positive and negative peaks follow each other over time a) a first positive peak (P1) in the first derivative signal (S.sub.Q′); a second positive peak (P2) in the electrostatic charge variation signal (S.sub.Q); a first negative peak (P3) in the first derivative signal (S.sub.Q′); a second negative peak (P4) in the electrostatic charge variation signal (S.sub.Q); a third positive peak (P5) in the first derivative signal (S.sub.Q′), b) a third negative peak in the first derivative signal (S.sub.Q′); a fourth negative peak in the electrostatic charge variation signal (S.sub.Q); a fourth positive peak in the first derivative signal (S.sub.Q′); a fifth positive peak in the electrostatic charge variation signal (S.sub.Q); a fifth negative peak in the first derivative signal (S.sub.Q′).

[0230] The operation of detecting the first signal characteristics may further include detecting one or more of the following time relationships:

T1=T3+T4

T6=T2+T3

T7=T4+T5

T8=T6+T7

wherein T1 may be a time interval between the second positive peak (P2) and the second negative peak (P4), T2 may be a time interval between the second positive peak (P2) and the first positive peak (P1), T3 may be a time interval between the second positive peak (P2) and the first negative peak (P3), T4 may be a time interval between the second negative peak (P4) and the first negative peak (P3), T5 may be a time interval between the second negative peak (P4) and the third positive peak (P5), T6 may be a time interval between the first positive peak (P1) and the first negative peak (P3), T7 may be a time interval between the first negative peak (P3) and the third positive peak (P5), T8 may be a time interval between the first positive peak (P1) and the third positive peak (P5).

[0231] Said time intervals T1-T7 may be respective time distances between respective maximum or minimum points of the positive and negative peaks.

[0232] The operation of detecting the first signal characteristics may further include detecting one or more of the following time relationships: T2.sub.TH_L<T2<T2.sub.TH_H, T3.sub.TH_L<T3<T3.sub.TH_H, T4.sub.TH_L<T4<T4.sub.TH_H, T5.sub.TH_L<T5<T5.sub.TH_H, where T2.sub.TH_L, T3.sub.TH_L, T4.sub.TH_L and T5.sub.TH_L may be respective lower thresholds of respective value including between 30 and 70 ms, and T2.sub.TH_H, T3.sub.TH_H, T4.sub.TH_H and T5.sub.TH_H are respective higher thresholds of respective value including between 150-250 ms.

[0233] The second signal characteristics may belong to the pressure signal (S.sub.P), and said operation of detecting in said pressure signal (S.sub.P) the second signal characteristics may include detecting a time amplitude value and maximum value of a pressure peak present in said pressure signal (S.sub.P); calculating a first comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the pressure peak; verifying whether said first comparison parameter is in a predetermined relationship with a first threshold.

[0234] Detecting a time amplitude value may include calculating an integral of, or an area subtended by, the pressure peak present in said pressure signal (S.sub.P), and said first comparison parameter may be calculated by dividing the result of said integral of the pressure peak, or the value of said area subtended by the pressure peak, by the maximum value of the pressure peak.

[0235] The third signal characteristics belong to the vibration signal (S.sub.A), and said operation of detecting in said vibration signal (S.sub.A) the third signal characteristics may include detecting a time amplitude value and maximum value of a vibration peak present in said vibration signal (S.sub.A); calculating a second comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the vibration peak; verifying whether said second comparison parameter is in a predetermined relationship with a second threshold.

[0236] Detecting a time amplitude value may include calculating an integral of, or an area subtended by, the vibration peak present in said vibration signal (S.sub.A), and said second comparison parameter may be calculated by dividing the result of said integral of the vibration peak, or the value of said area subtended by the vibration peak, by the maximum value of the vibration peak.

[0237] A method for detecting a presence in an environment to be monitored, may be summarized as including the steps of detecting, through an electrostatic charge variation sensor (6), a variation of electrostatic charge in said environment and generating an electrostatic charge variation signal (S.sub.Q); detecting, through a vibration sensor (7) operatively coupled to the environment to be monitored, an environmental vibration in the environment to be monitored and generating a vibration signal (S.sub.A); alternatively detecting, through an environmental pressure sensor (4) operatively coupled to the environment to be monitored, an environmental pressure in the environment to be monitored and generating a pressure signal (S.sub.P), acquiring, through a processing unit (2) from the electrostatic charge variation sensor (6), the electrostatic charge variation signal (S.sub.Q); detecting, through the processing unit (2) in said electrostatic charge variation signal (S.sub.Q), first signal characteristics indicative of the presence of a subject in said environment to be monitored; acquiring, through the processing unit (2) from said one of the vibration sensor (7) and the environmental pressure sensor (4), respectively the vibration signal (S.sub.A) or the pressure signal (S.sub.P); detecting, through the processing unit (2), in said vibration signal (S.sub.A) or pressure signal (S.sub.P) acquired, respective second signal characteristics indicative of the presence of the subject in said environment to be monitored; and generating, through the processing unit (2), a warning signal if both the first and the second signal characteristics have been detected.

[0238] The method may further include the other step of detecting the environmental vibration and generating a vibration signal (S.sub.A) or detecting the environmental pressure and generating a pressure signal (S.sub.P),

[0239] and further including the steps, performed by the processing unit (2), of acquiring, from the other of said vibration sensor (7) and environmental pressure sensor (4), respectively the vibration signal (S.sub.A) or the pressure signal (S.sub.P); detecting, in said other of the vibration signal (S.sub.A) and the pressure signal (S.sub.P) acquired, respective third signal characteristics indicative of the presence of the subject in said environment to be monitored; and generating the warning signal if all the first, the second and the third signal characteristics have been detected.

[0240] The step of detecting the first signal characteristics may include detecting, in the electrostatic charge variation signal (S.sub.Q), the following characteristics which follow each other in temporal order: a first rising edge; a first local maximum; a first falling edge; a first local minimum; a second rising edge; alternatively, detecting, in the electrostatic charge variation signal (S.sub.Q), the following characteristics which follow each other in temporal order: a falling edge; a first local minimum; a first rising edge; a first local maximum; a second falling edge.

[0241] The step of detecting the first signal characteristics may further include performing a comparison of said local maximums and minimums with respective thresholds; and assessing, through comparison with respective thresholds, a value of steepness or rising rate of the first and the second rising edges and of steepness or falling rate of the falling edge.

[0242] Detecting, in the electrostatic charge variation signal (S.sub.Q), the characteristics that follow each other in temporal order may include calculating the first derivative (S.sub.Q′) of the electrostatic charge variation signal (S.sub.Q); identifying, in the electrostatic charge variation signal (S.sub.Q) and in the first derivative signal (S.sub.Q′), a respective plurality of positive and negative peaks; detecting one of the following time series a) and b) with which said plurality of positive and negative peaks follow each other over time a) a first positive peak (P1) in the first derivative signal (S.sub.Q′); a second positive peak (P2) in the electrostatic charge variation signal (S.sub.Q); a first negative peak (P3) in the first derivative signal (S.sub.Q′); a second negative peak (P4) in the electrostatic charge variation signal (S.sub.Q); a third positive peak (P5) in the first derivative signal (S.sub.Q′), b) a third negative peak in the first derivative signal (S.sub.Q′); a fourth negative peak in the electrostatic charge variation signal (S.sub.Q); a fourth positive peak in the first derivative signal (S.sub.Q′); a fifth positive peak in the electrostatic charge variation signal (S.sub.Q); a fifth negative peak in the first derivative signal (S.sub.Q′).

[0243] The step of detecting the first signal characteristics may further include detecting one or more of the following time relationships:

T1=T3+T4

T6=T2+T3

T7=T4+T5

T8=T6+T7

wherein T1 may be a time interval between the second positive peak (P2) and the second negative peak (P4), T2 may be a time interval between the second positive peak (P2) and the first positive peak (P1), T3 may be a time interval between the second positive peak (P2) and the first negative peak (P3), T4 may be a time interval between the second negative peak (P4) and the first negative peak (P3), T5 may be a time interval between the second negative peak (P4) and the third positive peak (P5), T6 may be a time interval between the first positive peak (P1) and the first negative peak (P3), T7 may be a time interval between the first negative peak (P3) and the third positive peak (P5), T8 may be a time interval between the first positive peak (P1) and the third positive peak (P5).

[0244] Said time intervals T1-T7 may be respective time distances between respective maximum or minimum points of the positive and negative peaks.

[0245] The step of detecting the first signal characteristics may further include detecting one or more of the following time relationships: T2.sub.TH_L<T2<T2.sub.TH_H, T3.sub.TH_L<T3<T3.sub.TH_H, T4.sub.TH_L<T4<T4.sub.TH_H, T5.sub.TH_L<T5<T5.sub.TH_H, where T2.sub.TH_L, T3.sub.TH_L, T4.sub.TH_L and T5.sub.TH_L may be respective lower thresholds of respective value including between 30 and 70 ms, and T2.sub.TH_H, T3.sub.TH_H, T4.sub.TH_H and T5.sub.TH_H are respective higher thresholds of respective value including between 150-250 ms.

[0246] The second signal characteristics belong to the pressure signal (S.sub.P), and wherein said step of detecting in said pressure signal (S.sub.P) the second signal characteristics may include detecting a time amplitude value and maximum value of a pressure peak present in said pressure signal (S.sub.P); calculating a first comparison parameter which is a function of the ratio between said time amplitude value and maximum value of the pressure peak; verifying whether said first comparison parameter is in a predetermined relationship with a first threshold.

[0247] Detecting a time amplitude value may include calculating an integral of, or an area subtended by, the pressure peak present in said pressure signal (S.sub.P), and said first comparison parameter may be calculated by dividing the result of said integral of the pressure peak, or the value of said area subtended by the pressure peak, by the maximum value of the pressure peak.

[0248] The third signal characteristics belong to the vibration signal (S.sub.A), and said step of detecting in said vibration signal (S.sub.A) the third signal characteristics may include detecting a time amplitude value and maximum value of a vibration peak present in said vibration signal (S.sub.A); calculating a second comparison parameter which may be a function of the ratio between said time amplitude value and maximum value of the vibration peak; verifying whether said second comparison parameter is in a predetermined relationship with a second threshold.

[0249] Detecting a time amplitude value may include calculating an integral of, or an area subtended by, the vibration peak present in said vibration signal (S.sub.A), and said second comparison parameter may be calculated by dividing the result of said integral of the vibration peak, or the value of said area subtended by the vibration peak, by the maximum value of the vibration peak.

[0250] The various embodiments described above can be combined to provide further embodiments. All of the U.S. patents, U.S. patent application publications, U.S. patent applications, foreign patents, foreign patent applications and non-patent publications referred to in this specification and/or listed in the Application Data Sheet are incorporated herein by reference, in their entirety. Aspects of the embodiments can be modified, if necessary to employ concepts of the various patents, applications and publications to provide yet further embodiments.

[0251] These and other changes can be made to the embodiments in light of the above-detailed description. In general, in the following claims, the terms used should not be construed to limit the claims to the specific embodiments disclosed in the specification and the claims, but should be construed to include all possible embodiments along with the full scope of equivalents to which such claims are entitled. Accordingly, the claims are not limited by the disclosure.