This application claims priority to European Application 5 No. 20161699.2 filed on Mar. 9, 2020 and European Application No. 20178419.6 filed on Jun. 5, 2020, incorporated herein by reference.
The invention relates to a method for automatic multi-object localization and/or vital sign monitoring using radar, a system for automatic multi-object localization and/or vital sign monitoring using radar, and a computer program product comprising computer program code means adapted for automatic multi-object localization and/or vital sign monitoring using radar. The object can, for example, be people.
Generally, in times of an increasing number of older people, there is a growing need of contactless health monitoring in order to assist clinicians and caregivers in their daily work.
EP 3 425 419 A1 relates to a method for localization and monitoring of living being targets in an environment. Said method comprises the steps of transmitting a sequence of radio frequency waveforms, the waveforms being a continuous-wave waveform modulated in frequency and/or phase, and detecting a sequence of reflected waveforms being reflected by a target and Doppler-shifted due to a movement of the target, and forming a sequence of waveform transforms. The waveform transform comprises discretized information in a plurality of range bins. The information in a single range bin corresponds to reflections occurring at a specific sector in the environment. The information for a single specific sector in a sub-sequence of the sequence of waveform transforms is analyzed. Finally, the movement of a target in the specific sector based on the waveform transform information for that specific sector during a time period corresponding to the sub-sequence is detected. However, said method does not allow for automatic multi-object localization and/or vital sign monitoring in a particularly accurate and efficient manner.
Accordingly, there is an object to provide a method for automatic multi-object localization and/or vital sign monitoring, a system for automatic multi-object localization and/or vital sign monitoring using radar, and a computer program product comprising computer program code means adapted for automatic multi-object localization and/or vital sign monitoring.
This problem is solved by the features of claim 1 for a method for automatic multi-object localization and/or vital sign monitoring, the features of claim 14 for a system for automatic multi-object localization and/or vital sign monitoring using radar, and the features of claim 15 for a computer program product comprising computer program code means adapted for automatic multi-object localization and/or vital sign monitoring. The dependent claims contain further developments.
According to a first aspect of the invention, a method for automatic multi-object localization and/or vital sign monitoring is provided. The method comprises the steps of transmitting a radar signal, receiving a corresponding radar signal in order to form a corresponding observation matrix, reducing noise by applying singular value decomposition to the observation matrix, processing the result of the singular value decomposition by an independent component analysis in order to estimate the corresponding sources, and estimating propagation channels of the estimated sources by minimizing the corresponding residual error based on the observation matrix and the estimated sources. Advantageously, multiple objects, such as people, can automatically be localized and respective vital signs can be monitored in a highly accurate and efficient manner.
According to a first preferred implementation form of the first aspect of the invention, before the independent component analysis, the method further comprises the step of estimating the number of respective targets, preferably the number of respective persons, especially by calculating the signal-to-noise ratio of the respective uncorrelated sources in the result of the singular value decomposition. Advantageously, for example, inaccuracies can further be reduced.
According to a second preferred implementation form of the first aspect of the invention, the method further comprises the step of locating the respective objects and/or removing the respective order ambiguity on the basis of the estimated propagation channels. Advantageously, for instance, accuracy can further be increased.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of performing a phase demodulation with respect to the estimated sources. Advantageously, for example, inaccuracies can further be reduced. Further advantageously, the phase demodulation may preferably comprise or be a linear demodulation.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of extracting the respective vital signs information or signal from the estimated sources after the phase demodulation. Advantageously, for instance, both efficiency and accuracy can further be increased. Further advantageously, the respective vital sings information or signal may preferably comprise or be a time domain signal.
According to a further preferred implementation form of the first aspect of the invention, the step of extracting is performed on the basis of the following equation:
ŷ(t)=ŝ(t)λ0/4π,
wherein λ0 denotes the corresponding wavelength at the first frequency of a respective chirp signal, and
wherein π denotes the constant Pi.
Advantageously, for example, complexity can be reduced, which leads to an increased efficiency.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of obtaining the corresponding respiration and heartbeat signals by performing a filtering operation with respect to the respective vital signs information or signal. Advantageously, for instance, accuracy can further be increased.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of obtaining the corresponding respiration and heartbeat signals by performing a wavelet decomposition with respect to the respective vital signs information or signal. Advantageously, for example, inaccuracies can further be reduced.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of obtaining the corresponding respiration and heartbeat signals by performing a Hilbert transform followed by a convex optimization with respect to the respective vital signs information or signal. Advantageously, for instance, accuracy can further be increased.
According to a further preferred implementation form of the first aspect of the invention, the wavelet decomposition comprises the usage of a discrete Meyer mother-wavelet. Advantageously, for example, complexity can further be reduced, thereby increasing efficiency.
According to a further preferred implementation form of the first aspect of the invention, the method further comprises the step of estimating the corresponding heart rates by performing a frequency transform with respect to the respective respiration and heartbeat signals. Advantageously, for instance, inefficiencies can further be reduced.
According to a further preferred implementation form of the first aspect of the invention, the frequency transform comprises or is a fast Fourier transform. Advantageously, for example, by further reducing complexity, efficiency can be increased.
According to a further preferred implementation form of the first aspect of the invention, the radar signal comprises or is a chirp signal. Advantageously, for instance, accuracy can further be increased.
According to a further preferred implementation form of the first aspect of the invention, the radar signal originates from a frequency-modulated continuous wave radar, preferably a single input single output frequency-modulated continuous wave radar. Advantageously, for example, both efficiency and accuracy can further be increased. Further advantageously, instead of single input single output, it may be used multiple input multiple output or multiple input single output or single input multiple output. As a further advantage, any other type of modulated radar may be used.
According to a second aspect of the invention, a system for automatic multi-people localization and vital sign monitoring using radar is provided. The system comprises a transmitter unit, a receiving unit and a processing unit. In this context, the transmitting unit is configured to transmit a radar signal, the receiving unit is configured to receive a radar signal in order to form a corresponding observation matrix. In addition to this, the processing unit is configured to reduce noise by applying singular value decomposition to the observation matrix, to process the result of the singular value decomposition by an independent component analysis in order to estimate the corresponding sources, and to estimate propagation channels of the estimated sources by minimizing the corresponding residual error based on the observation matrix and the estimated sources. Advantageously, multiple objects, such as people, can automatically be localized and respective vital signs can be monitored in a highly accurate and efficient manner. Further advantageously, the processing unit may additionally be configured to estimate the number of respective targets, preferably the number of respective persons, especially by calculating the signal-to-noise ratio of the respective uncorrelated sources in the result of the singular value decomposition before the independent component analysis.
According to a third aspect of the invention, a computer program product is provided. The computer program product comprising computer program code means adapted for automatic multi-object localization and/or vital sign monitoring according to the steps of the method and of any of the preferred implementation forms thereof when said program is run on a computer or any electronic system. Advantageously, for example, multiple people can automatically be localized and respective vital signs can be monitored in a highly accurate and efficient manner.
Exemplary embodiments of the invention are now further explained with respect to the drawings by way of example only, and not for limitation. In the drawings:
Firstly,
It is noted that it might be particularly advantageous if the method further comprises the step of locating the respective targets and/or removing the respective order ambiguity on the basis of the estimated propagation channels Ĥ.
In addition to this or as an alternative, the method may further comprise the step of performing a phase demodulation with respect to the estimated sources Ŝ.
In this context, it might be particularly advantageous if the method further comprises the step of extracting the respective vital signs information or signal ŷ(t) from the estimated time domain sources ŝ(t) after the phase demodulation.
Further advantageously, the step of extracting may especially be performed on the basis of the following equation:
In this context, λ0 denotes the corresponding wavelength at the first frequency of a respective chirp signal, and π denotes the constant Pi. Moreover, the method may further comprise the step of obtaining the corresponding respiration and heartbeat signals by performing a filtering operation with respect to the respective vital signs information or signal ŷ(t).
It is further noted that the method may further comprise the step of obtaining the corresponding respiration and heartbeat signals by performing a wavelet decomposition with respect to the respective vital signs information or signal ŷ(t). It might be particularly advantageous if the method further comprises the step of obtaining the corresponding respiration and heartbeat signals by performing a Hilbert transform followed by a convex optimization with respect to the respective vital signs information or signal ŷ(t).
With respect to the above-mentioned wavelet decomposition, it is noted that the wavelet decomposition may preferably comprise the usage of a discrete Meyer mother-wavelet.
Furthermore, the method may additionally or alternatively comprise the step of estimating the corresponding heart rates by performing a frequency transform with respect to the respective respiration and heartbeat signals. With respect to said frequency transform, it is noted that the frequency transform may preferably comprise or be a fast Fourier transform.
Moreover, with respect to the above mentioned radar signal, it is noted that the radar signal may preferably comprise or be a chirp signal. In addition to this or as an alternative, the radar signal may especially originate from a frequency-modulated continuous wave radar, preferably a single input single output frequency-modulated continuous wave radar.
Furthermore,
In this context, as a first step 99 (before the above-mentioned step 100), a radar signal is transmitted. In addition to this, especially in the sense of a target existence probability, between the above-mentioned steps 101 and 102, an additional step 104 has been inserted. As it can be seen from
With respect to
Advantageously, especially in the sense of a target existence probability, the processing unit 12 may additionally be configured to estimate the number of respective targets, preferably the number of respective persons, P especially by calculating the signal-to-noise ratio of the respective uncorrelated sources in the result of the singular value decomposition before the independent component analysis.
It might be particularly advantageous if the processing unit 12 is further configured to locate the respective targets and/or removing the respective order ambiguity on the basis of the estimated propagation channels Ĥ. In addition to this or as an alternative, the processing unit 12 may preferably be configured to perform a phase demodulation with respect to the estimated sources Ŝ. In this context, the processing unit 12 of the system 10 may especially be configured to extract the respective vital signs information or signal ŷ(t) from the estimated time domain sources ŝ(t) after the phase demodulation.
It might be particularly advantageous if said extraction is performed on the basis of the following equation:
wherein λ0 denotes the corresponding wavelength at the first frequency of a respective chirp signal, and
wherein π denotes the constant Pi.
In addition to this or as an alternative, the processing unit 12 may preferably be configured to obtain the corresponding respiration and heartbeat signals by performing a filtering operation with respect to the respective vital signs information or signal ŷ(t). Further additionally or further alternatively, the processing unit 12 may especially be configured to obtain the corresponding respiration and heartbeat signals by performing a wavelet decomposition with respect to the respective vital signs information or signal ŷ(t).
As a further alternative or in further addition to this, the processing unit 12 may especially be configured to obtain the corresponding respiration and heartbeat signals by performing a Hilbert transform followed by a convex optimization with respect to the respective vital signs information or signal ŷ(t).
With respect to the above-mentioned wavelet decomposition, it is noted that it might be particularly advantageous if the wavelet decomposition comprises the usage of a discrete Meyer mother-wavelet.
Again, with respect to the processing unit 12 of the system 10, it is noted that the processing unit 12 may additionally or alternatively be configured to estimate the corresponding heart rates by performing a frequency transform with respect to the respective respiration and heartbeat signals. In this context, with respect to said frequency transform, it is noted that it might be particularly advantageous if the frequency transform comprises or is a fast Fourier transform.
Furthermore, with respect to the radar signal, it is noted that the radar signal may preferably comprise or be a chirp signal.
In addition to this or as an alternative, the radar signal may especially originate from a frequency-modulated continuous wave radar, preferably a single input single output frequency-modulated continuous wave radar.
With respect to
In this context, whereas said system 20 comprises the already above-mentioned transmitter unit 9 and the receiving unit 11, the processing unit 12 is located outside the system 20.
Accordingly,
In addition to the explanations above, the invention will be explained in more detail especially considering a frequency-modulated continuous wave (FMCW) radar in the following. However, the invention can be used with any radar with range-Doppler capability, for instance, stepped-frequency continuous wave (SFCW), phase-modulated continuous wave (PMCW), ultra-wideband impulse radio (UWB-IR) etc. A linear FMCW radar transmits a series of signals, called chirps, whose instantaneous frequency increases linearly over time as illustrated by
It is noted that said time interval can be avoided so that the respective chirps are continuously transmitted without interruption.
The baseband signals (called beat signals), resulting by mixing the copies of the transmitted chirps with the received chirps, are digitized and arranged in a K×M data matrix, where K is the number of samples acquired per beat signal and depends both on T and on the sampling time Tf of the analog-to-digital converter (ADC) while M is the number of beat signals that are used in the data processing, collected each Ts. In radar terminology, Tf determines the fast time while Ts determines the slow time. This matrix is preprocessed by applying the fast Fourier transform (FFT) per column. The result is a K×M complex data matrix whose elements are denoted by X(k,m), with k=0, . . . , K−1 especially being the range bin index, K=N, and m=0, . . . , M−1 is the index in slow time.
The columns contain the range profiles, determined each Ts, while the rows contain the Doppler signals (i.e., phase shifts or vital signs information) from the range bins. This corresponds to ideally divide the whole radar detectable space (in radar terminology known as radar unambiguous range) into K evenly spaced concentric circles, namely range bins, whose size (i.e., range resolution) depends on the radar waveform. For example, it is to be considered an ideal situation with a single still subject at the k-th range bin in an empty room with no multipath. If it is scanned each row, and if we neglect the spreading effect of the FFT, there will be phase variations only at the k-th range bin. Exploiting especially this feature, we can determine the location of the target multiplying k by the range resolution and extract the vital signs from the complex signal in the k-th raw.
In indoor environments, the multipath interference cannot be ignored while designing a radar-based system. It occurs when a signal takes two or more paths from the transmitting antenna to the receiving antenna and depends on the room structure and the object's locations. This involves that the multipath signals generated by a subject may have identical delays (i.e., time of flight) of the direct path signals of one or more persons, causing false localizations and incorrect vital signs extraction. This is especially equivalent to saying that the multipaths and the direct paths signals comes from the same range bins.
The situation becomes even more problematic in presence of static reflectors (i.e., clutter, objects). Especially based on the limited room size, we consider only the first L range bins of the K total available ones. Within a range bin, differences in delays may be small and translated into phase shifts. Correspondingly, we assume that the number of possible path delays may preferably also be L. As a result, the (slowly time-varying) multipath propagation channel impulse response h(t,m) may be modeled as:
h(t,m)=Σl=0L-1βl·δ(t−τl/(m)), (1)
wherein δ(⋅) denotes the Dirac delta function, l denotes the path index, τ denotes the propagation of the l-th path, and β denotes the complex path gain which indicates the overall attenuation and phase shift.
In this context, it is assumed that there are P people and Q static clutter in the room, which are regarded as point scatters and, therefore, the summation of all reflections from each infinitesimal point of a target's surface which falls into the same range bin can be regarded as one overall reflection through one propagation path. The only difference between the two types of targets is that the subjects have physiological activities (i.e., variable phase in slow time or Doppler signal), leading to a time varying τl(m), while the static clutter involves a constant τl (i.e., constant phase in slow time or Doppler signal). A sketch of this data model is shown in
The complex elements S(k,m) form a two-dimensional (2D) observation matrix X that can be decomposed as:
X=HS+C+N, (2)
wherein H denotes an L×P complex mixing matrix, S denotes a P×M complex matrix containing the Doppler shifts (i.e., vital signs information) caused by the P subjects at each Ts, C denotes an L×M matrix with identical columns containing the DC (direct current) information in slow-time resulting from static reflections, and N denotes an L×M matrix containing zero mean additive noise.
The next step is to perform AC (alternate current) coupling to foregoing equation in order to remove the DC components (i.e., the mean values) while still preserving the physiological (i.e., variable) motions. After that, the observation matrix X becomes:
wherein the noise remains the same because of the zero-mean nature. The statistic properties of this data model especially are: (1)
It is noted that the model of the foregoing equation is especially used as starting point for the vital-sign monitoring and automatic localization in accordance with the invention.
Another exemplary block diagram of the inventive method is shown in
In an ideal situation, with a range resolution of 20 cm, the Doppler information would be expected only in range bins 6 and 11 respectively for subject 1 and subject 2. However, as it can be seen in
Similar considerations can be done for the signals of bin 11 (which can be considered as three sines in 20 seconds), which are replicated in bins 12-14, 16-19. In addition, the two subjects interfere each other in bins 14, 17-19 where the I signal corresponds to one subject while the Q signal corresponds to the other subject.
It is noted that some solutions have been proposed for automatic target localization based on vital signs according to the state of the art. Those approaches rely on the variation of the range spectrum to differentiate humans from stationary targets. The idea is that the physiological movements, involve a larger standard deviation (std) than static objects. Therefore, the range bins with large variations indicate the location of human targets. However, this approach fails in presence of multipaths, which result both in copies of the same target at wrong distances and in new (therefore false) targets when direct paths and multipaths combine in certain range bins.
This is demonstrated by the standard deviation profile depicted in
Especially looking at bin 13 in
With the aid of the invention, first the vital signs signals of the subjects are determined, and then they are located. In accordance with
At this point, especially an order ambiguity issue has to be faced: it is not yet possible to indicate which source (i.e., vital signs signal) corresponds to which target. From the equation
Another approach can be generating the analytic signals of ŷ(t) by using the Hilbert transform:
ŷ*(t)=(ŷ(t))=ar sin(2πfrt)+ah sin(2πfht), (4)
wherein denotes the Hilbert transform operator.
The analytic signal ŷA(t) is:
ŷ
A(t)=ŷ(t)+jŷ*(t)=arej2πf
The second last equation can be rewritten as:
ŷ
A=ψHf, (6)
wherein
f=[ej2πf
Ψ is an unknown vector which indicates the complex weights of the sinusoidal candidates in f, and Δf is the resolution of the discretized dictionary.
The estimation of Ψ from ŷA(t) is a convex optimization problem of the form:
wherein ζ denotes a positive real value and the l−1 norm term is added to preserve the sparsity of the weight vector.
Since the respiration rate is normally within the 0.1-0.8 Hz range and the respiration signal is stronger than the heartbeat signal, it can be assumed:
y(t)≈yr(t)=arej2πf
Therefore, one can temporally focus on frequencies between 0.1 and 0.8 Hz to reconstruct ŷA(t) and to estimate the respiration rate. The heartbeat signal ŷh(t) can be estimated by subtracting the reconstructed respiration signal ŷr(t) from ŷA(t). The corresponding heartbeat can be obtained with the same methodology of the second last equation.
Again, with respect to the above-mentioned singular value decomposition, the target existence probability, the independent component analysis, the linear demodulation, and the automatic localization, the corresponding steps thereof should be discussed in greater detail in the following.
In this context, an array signal model is introduced: The received signal sR(t,m) over a multipath channel can be modelled as the convolution between the channel impulse response h(t,m) and the transmitted signal sT(t).
Considering a single target, it can be expressed as:
wherein aT is the complex amplitude indicating the amplitude of the chirp and its initial phase, ρ=B/T is the sweeping rate, m is the slow time index, i is the index corresponding to the i-th target/object, d is the path distance, y(m) is the chest surface displacement, c0 is the speed of light, and * is the convolution operator.
The baseband signal sB(t,m) can be modelled as:
where the superscript * indicates the complex conjugate operation. The quadratic term is neglected because τi,j(m) is close to zero. The digitized baseband signal can be expressed as:
s
B(n,m)≈ΣiΣlβi,l*aT2ej2πf
where n is the fast time index and Tf is the fast time sampling time (i.e., ADC sampling rate). After performing FFT in fast time, the frequency domain signal X(k,m) becomes:
wherein F is the Fourier transform operator, w(n) is a rectangular window function in fast time, w(n) and W(k) are a Fourier pair. Since the rectangular window in frequency domain is a sinc function with gradients close to zero around 2πρτi,l the frequency domain window function Wi,j(k) can be considered as a fixed one in slow time.
Assuming P subjects and Q static clutter in a room, equation (14) can be written as:
wherein ϕi(m) is the Doppler shift caused by the vital signs on the FMCW signal. The observation signal X(k,m) is a dual-variable function of k and m and forms a two-dimensional observation matrix X with factorization
Here, 1 is a length M all-ones column vector. In presence of additive noise, the data model becomes as in equation (2) mentioned above. The next step is to remove C in X using a projection (or centering matrix) P such that CP=0. This requires 1T P=0, where 1 is a vector with all entries equal to 1, so that
P=I−1(1T1)−11T, (24)
wherein I is an identity matrix. The result is:
XP=HSP+CP+NP=HSP+N, (25)
wherein the noise statistics remain the same due to their zero mean nature. Equation (25) is especially equivalent to the mentioned above equation (3) with:
Now, with respect to the singular value decomposition, it is noted that the noise is reduced by applying singular value decomposition (SVD) to equation (3) as:
wherein U denotes an L×L unitary matrix containing left singular vectors, V is an M×M unitary matrix containing right singular vectors, and τ is a diagonal matrix containing all the singular values.
The first P columns in U and V are denoted respectively as Us and Vs while the rest columns as Un and Vn respectively. As noise is assumed to be uncorrelated to the sources in
Therefore, Vs can be expressed as the following linear transform:
V
s
=A
wherein A is a P×P square mixing matrix. Therefore, by taking the signal space, one cannot only remove the noise falling into the null space, but also pre-whiten the data based on the second order statistic.
With respect to the target existence probability, it is noted that in order to estimate P, and therefore to determine Vs, we calculate the signal-to-noise ratio (SNR) of the uncorrelated sources in V. An exemplary calculation of the SNR is explained in the following.
A spectrum of a canonical radar-based vital signs signal consists essentially of the respiration fundamental, which is the dominant component of the signal, of one or two decreasing in magnitude respiration harmonics, and of the very small heartbeat fundamental.
We consider as signal power the power within the fundamental and its first harmonic while the rest of the spectrum is considered as signal noise. The first P sources of V (i.e., first P columns) produce a high SNR and indicate Vs while the remaining sources have very low SNR and indicate Vn.
Furthermore, also other checks may be performed on the spectrum's local maxima:
We consider as noise any source producing a ratio less than 2. In fact, in the canonical spectrum, the two highest peaks indicate the respiration fundamental and its first harmonic and their ration is always greater than a factor 2.
In those two situations, we fix the SNR to 0 dB. The last operation is to scan the obtained SNR profile starting from the first estimation and stopping when the first descending order uncorrelated source produced an SNR below a threshold. The latter source indicates the starting of Vn while the previous ones are the Vs sources corresponding to the P subjects.
It is further noted that with the aid of the target existence probability, we especially extract the Doppler signals in order to first determine the number of targets and then to estimate their vital signs and location. The first step is to reduce the noise by applying the singular value decomposition to
Again, it is noted that the foregoing explanations show an exemplary way to calculate the SNR. For instance, the signal power can be calculated from the fundamental to the second harmonic of the respiration or other way. Further exemplarily, the checks on the peaks could not be necessary may be a preferred implementation form.
With respect to the independent component analysis, it is noted that in equation (29), the rows of
Since the first step of the independent component analysis, i.e., the pre-whitening, has been done by the singular value decomposition, the estimation of
Ŝ=W
H
V
s. (30)
Based on the central limit theorem, the sum of independent random variables tends to a Gaussian distribution. This gives the inspiration that the estimated independent sources should tend to be as non-Gaussian distributed as possible. There are various ways to measure non-Gaussianity.
With respect to the phase demodulation or the linear demodulation, respectively, it is noted that especially based on a small angle approximation, preferably valid for sub-10 GHz radar, a source can be approximated as:
wherein ϕ(t) is the phase (Doppler) shift caused by the vital signs. Especially after DC removal, we observe that the imaginary part is more powerful than the real part:
(t)≈jØ(t). (32)
The vital signs information or signal ŷ(t) can be extracted applying the linear demodulation to the estimated source ŝ(t) and it can be expressed as:
wherein ŝLIN(t) is the estimated source after the linear demodulation and λ0 is the wavelength corresponding with the start frequency of the chirp. If the demodulation is performed correctly, equation (33) is equivalent to the motion of the chest surface x(t) caused by the vital signs. This movement can be modeled as the sum of two sinusoidal functions as:
x(t)=xr(t)+xh(t)=ar cos(2πfrt)+ah cos(2πfht), (34)
wherein xr(t) and xh(t) indicate respectively the mechanical chest surface displacements caused by respiration (expansion of thorax and lungs) and heart contractions; ar and ah are the maximum mechanical displacements caused by the lungs and heart on the chest surface (with typical amplitudes around 4-12 mm and 0.1-0.5 mm, respectively), and fr and fh are the vital signs frequencies which represent the information to be extracted.
Finally, with respect to the automatic localization, it is noted that in order to locate the targets and to remove the order ambiguity, knowing the observation matrix
However, due to the path loss, the observation vectors (rows) in
Assuming
In order to preserve the sparsity of the multipath propagation channel, the final cost function should be
wherein ζ is the penalty coefficient which represents a trade-off between the residual error and the sparsity. As the scaling matrix Ā is known, it is obtained
Ĥ=
Again, with respect to the target existence probability, two further exemplary embodiments being based on the embodiment according to
As it can been from said
Moreover,
In this context, as it can be seen, fixing a threshold thr (dashed line in
Finally,
In this context, it is noted that a visual check of the time domain signals clearly reveals that the first two sources resemble the typical radar-based vital signs signals.
In addition to this, these sources result less noisy than the other three. This is also confirmed by the spectra plots. In fact, the first two spectra consist essentially of the respiration fundamental, which is the dominant component of the signal, and of its harmonic, while the other three spectra contain several components.
We can notice that, although noisy, the fourth source could be misleading as it might resemble the second one. However, a check of the fundamental-harmonic ratio resolves the ambiguity. In fact, in a normal radar-based vital signs signals this ratio is always greater than a factor 2 (see the spectrum of the second source).
While various embodiments of the present invention have been described above, it should be understood that they have been presented by way of example only, and not limitation. Numerous changes to the disclosed embodiments can be made in accordance with the disclosure herein without departing from the spirit or scope of the invention. Thus, the breadth and scope of the present invention should not be limited by any of the above described embodiments. Rather, the scope of the invention should be defined in accordance with the following claims and their equivalents.
Although the invention has been illustrated and described with respect to one or more implementations, equivalent alterations and modifications will occur to others skilled in the art upon the reading and understanding of this specification and the annexed drawings. In addition, while a particular feature of the invention may have been disclosed with respect to only one of several implementations, such feature may be combined with one or more other features of the other implementations as may be desired and advantageous for any given or particular application.
Number | Date | Country | Kind |
---|---|---|---|
20161699.2 | Mar 2020 | EP | regional |
20178419.6 | Jun 2020 | EP | regional |