## Individual Identification Using Radar-Measured Respiratory and Heartbeat Features

Haruto Kobayashi, Yuji Tanaka, and Takuya Sakamoto

Abstract -This study proposes a method for radar-based identification of individuals using a combination of their respiratory and heartbeat features. In the proposed method, the target individual's respiratory features are extracted using the modified raised-cosine-waveform model and their heartbeat features are extracted using the mel-frequency cepstral analysis technique. To identify a suitable combination of features and a classifier, we compare the performances of nine methods based on various combinations of three feature vectors with three classifiers. The accuracy of the proposed method in performing individual identification is evaluated using a 79-GHz millimeter-wave radar system with an antenna array in two experimental scenarios and we demonstrate the importance of use of the combination of the respiratory and heartbeat features in achieving accurate identification of individuals. The proposed method achieves accuracy of 96.33% when applied to a five-day dataset of six participants and 99.39% when applied to a public one-day dataset of thirty participants.

## I. INTRODUCTION

At present, with the ongoing advancements in Internet of Things (IoT) technology, techniques for sensing of various information types in real environments have been attracting attention. IoT technologies are easing the digitalization of society and enabling the construction of sensing networks for use in daily life [1], [2]. In particular, research into identification of individuals using sensors has drawn attention in a variety of applications. Currently, facial recognition using cameras is primarily being researched and commercialized [3], [4]. Although facial recognition can achieve high-precision identification, it also raises privacy concerns and is unsuitable for sensing applications in private spaces. Additionally, external factors, such as low lighting, can lead to reduced recognition accuracy.

Rather than use cameras, contact sensor-based methods for individual identification that sense fingerprints [5], veins [6], and electrocardiograms (ECGs) [7]-[9] have been reported. Furthermore, individual identification techniques using respiratory sounds measured with microphones have also been reported [10], [11]. Han et al. [10] proposed an individual identification method that used respiratory sounds acquired from bone conduction microphones embedded in earphones and achieved identification accuracy of 95.17% for 20 participants. Although these studies focused on use of contact sensors to obtain the physiological signals required for individual

- H. Kobayashi and T. Sakamoto are with the Department of Electrical Engineering, Graduate School of Engineering, Kyoto University, Kyoto 6158510, Japan.

Y. Tanaka is with Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555, Japan.

identification, it can be difficult to attach the sensors or place the required sensors close to the body, and this factor will also hinder the identification of uncooperative individuals.

Recently, radar has been used widely for measurement of human bodies because both limb movements and physiological signals can be acquired in a noncontact manner and privacy concerns are reduced when compared with use of camerabased sensors [12]-[14]. Radar techniques have been used for respiration measurement [15]-[17], heartbeat [18]-[23] and also for individual identification applications [24], [25].

Limb motions such as walking and sitting movements are often used to identify individuals, e.g., in individual identification of walking and sitting persons using Wi-Fi signals [26]-[30] and radar techniques [31]-[42]. To enable individual identification of a person that is remaining still, Wi-Fi signals and radar have also been used to measure their physiological signals [43]-[51]. For example, Rahman et al. [46] extracted features from radar-measured respiration characteristics and identified individuals using a combination of the k -nearest neighbor ( k -NN) algorithm and majority voting. Although respiration measured via radar is useful for individual identification, respiration is also known to have both voluntary and involuntary aspects [52], which raises questions about its effectiveness for use in long-term identification. In contrast, Lin et al. [51] conducted individual identification using heartbeat features that were measured using a pair of radar systems, with accuracy of 98.61% being achieved for 78 participants.

In this paper, we propose an individual identification method that uses both respiratory and heartbeat features for the first time. A preprint of this manuscript has been posted online [53]. The contributions of this paper are listed as follows:

- · We propose a method to extract respiratory features in addition to extraction of respiratory intervals, and confirm the effectiveness of using these features to perform individual identification through experiments.
- · We propose a method to extract features of heartbeat harmonic components and demonstrate the effectiveness of use of these features in individual identification.
- · Through a series of experiments conducted over five days on six participants, we evaluate the long-term effectiveness of the proposed method. We also verify the method's accuracy when using the respiratory and heartbeat features separately and discuss the advantages of combining these features. Additionally, we validate the identification accuracy for larger groups using a public dataset with 30 participants.

## II. RESPIRATORY FEATURES FOR INDIVIDUAL IDENTIFICATION

## A. Waveform Feature Extraction Method

If the radar system receives a single dominant echo from a human body, the body's displacement can be estimated to be

$$d ( t ) = \frac { \lambda } { 4 \pi } \text{unwrap} \{ \not \leqslant s ( t ) \}, \text{ \quad \ \ } ( 1 )$$

where λ is the wavelength, s t ( ) is the signal reflected from the target, unwrap {·} is the phase unwrapping operator, and ∠ represents the phase of a complex number. When the target person is stationary, d t ( ) mainly contains the displacement components that result from the person's respiration and heartbeat.

The respiratory displacement caused by normal breathing can be expressed as a semi-periodic function of time with three phases: the expiratory phase, the expiratory plateau, and the inspiratory phase [54]. Fig. 1 shows an example of a radar-measured displacement waveform that highlights the three phases described above. To extract individual features related to the three respiratory phases, we use the modified raised-cosine-waveform (MRCW) model proposed by Hsieh et al. [55]. In the MRCW model, in addition to the respiratory period and amplitude, the user can also specify the time duration ratio of the expiratory, expiratory plateau, and inspiratory phases.

The MRCW model for one cycle T = 1 /f can be expressed using six parameters as follows:

$$\begin{array} { l } \text{using six parameters as follows: } & \text{$-0.5\bigsqcup$} \\ d _ { \text{MRCW} ( t, A, f, \beta _ { 1 }, \beta _ { 2 }, D ) = } & \text{$-1\bigsqcup$} \\ \left \{ A & \left ( \frac { T } { 2 } - T _ { \text{a1} } < t \leq \frac { T } { 2 } + T _ { \text{a2} } \right ) \\ A \cos \left \{ \frac { 2 \pi f } { \beta _ { 1 } } \left ( \left | t - \frac { T } { 2 } \right | - T _ { \text{a1} } \right ) \right \} \left ( T _ { \text{b1} } < t \leq \frac { T } { 2 } - T _ { \text{a1} } \right ) & \text{Fig. 2. Exa} \\ A \cos \left \{ \frac { 2 \pi f } { \beta _ { 2 } } \left ( \left | t - \frac { T } { 2 } \right | - T _ { \text{a2} } \right ) \right \} \left ( \frac { T } { 2 } + T _ { \text{a2} } < t \leq T - T _ { \text{a1} } \right ) \\ - A & \text{otherwise, } & \text{To det} \\ ( 2 ) & \text{sured sig} \end{array}$$

where f is the respiratory frequency, A is the amplitude, β 1 and β 2 are the roll-off rates, and D is the duty ratio. Note that β 1 , β 2 , and D correspond to the temporal ratios of the expiration, inspiration, and expiratory plateau phases, respectively. In addition, T a1 , T a2 , T b1 , and T b2 in (2) are determined as follows:

$$T _ { \text{a1} } = \frac { D ( 1 - \beta _ { 1 } ) } { 2 f }, \, T _ { \text{b1} } = \frac { ( 1 - D ) ( 1 - \beta _ { 1 } ) } { 2 f }, \quad ( 3 ) \quad \text{where} \quad \text{where} \quad \text{$\text{bits} \ \neq$}$$

$$T _ { \text{a2} } = \frac { D ( 1 - \beta _ { 2 } ) } { 2 f }, \, T _ { \text{b2} } = \frac { ( 1 - D ) \overset { \mu } { ( 1 - \beta _ { 2 } ) } } { 2 f }. \quad ( 4 ) \quad \text{shift a} \\ \quad \text{when a}$$

Example MRCW models are shown in Fig. 2 for f = 1 Hz and A = 1 . As the figure shows, the waveform for the MRCW model is dependent on the parameters ( f, β 1 , β 2 , D ) . Hereafter, the parameters p = ( f, β 1 , β 2 , D ) T are used as the respiratory features (a superscript T signifies the transpose operator).

Fig. 1. Example of respiratory displacement when measured using a radar system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-1cbe2eff-595c-4fb8-b874-d682e975934d.jpg)

Fig. 2. Examples of the respiratory displacement MRCW model with parameters (a) ( β , β 1 2 , D ) = (0 5 . , 0 5 . , 0 5) . , (b) ( β , β 1 2 , D ) = (0 1 . , 0 5 . , 0 5) . , (c) ( β , β 1 2 , D ) = (0 5 . , 0 1 . , 0 5) . , and (d) ( β , β 1 2 , D ) = (0 5 . , 0 5 . , 0 8) . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-2ef1a2ea-aeb9-49df-a4d7-17ad82c9a6fb.jpg)

To determine the MRCW model parameters from the measured signal, the following operation is performed:

$$\lambda _ { \nu } \quad \sim \cdot \quad \sim \cdot \quad \,. \\ \text{plitude}, \quad & \left ( \hat { \tau } ( t ), \hat { A } ( t ), \hat { p } ( t ) \right ) = \\ \text{at ratios} \quad & \quad \arg \min _ { ( \tau, A, p ) } \int _ { 0 } ^ { T _ { 0 } } \left | d ( t ^ { \prime } ) w ( t ^ { \prime } - t ) - d _ { \text{MRCW} } ( t ^ { \prime } - \tau, A, p ) \right | ^ { 2 } \, \text{d} t ^ { \prime }, \\ ( 2 ) \text{ are} \quad & \quad \, \hat { \tau } \quad \cdot \quad \, \hat { \nu } \quad \dots \quad \hat { \tau } \dots \quad \dots \quad \,.$$

$$\hat { d } _ { \text{MRCW} } ( t ) = d _ { \text{MRCW} } \left ( t - \hat { \tau } ( t ), \hat { A } ( t ), \hat { \mathfrak { p } } ( t ) \right ), \quad ( 6 )$$

where T 0 is the measurement time length, τ is the time shift and ˆ d MRCW ( ) t is the MRCW model that provides the best approximation of d t ( ) . In addition, w t ( ) is a rectangular window with a length of 8 s. The parameters estimated in (5) are defined here as the estimated waveform features ˆ( ) = ( ˆ ( p t f t ) , β ˆ ( ) 1 t , β ˆ ( ) 2 t , D t ˆ ( )) T at a time t .

The measured displacement waveform d t ( ) and the estimated waveform from the MRCW model ˆ d MRCW ( ) t are shown in Fig. 3. This comparison indicates that ˆ d MRCW ( ) t approximates d t ( ) . Note that the MRCW model parameters in

Fig. 3. Displacement waveform measured using radar (black solid line) and calculated using the MRCW model (red dashed line).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-a5ee1c89-a06a-45c1-978f-83447916a124.jpg)

Fig. 4. Displacement waveform measured using radar (black solid line) and calculated using the parabolic model (red solid line).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-63b29dcb-5f32-4568-9bc5-50eff359a837.jpg)

$$\text{Fig. } 3 \text{ are } \hat { f } = 0. 4 \text{ Hz and } ( \hat { \beta } _ { 1 }, \hat { \beta } _ { 2 }, \hat { D } ) = ( 0. 9, 0. 6, 0. 9 ). \quad \mathbb { R }$$

## B. Features of the Respiratory Expiratory Plateau

Although the MRCW model can approximate the respiratory displacement, the expiratory plateau waveform is modeled as a constant straight line, which is inaccurate, as illustrated in Fig. 3. In this study, the displacement waveform during the expiratory plateau time interval is used as an additional feature.

For t , parabolic fitting is performed using the parabolic curve d par ( ; t c 0 , c 1 , c 2 ) = c 2 t 2 + c 1 t + c 0 as shown:

$$( \hat { c } _ { 0 } ( t ), \hat { c } _ { 1 } ( t ), \hat { c } _ { 2 } ( t ) ) = & \quad \text{respirati} \\ \underset { ( c _ { 0 }, c _ { 1 }, c _ { 2 } ) } { \arg \min } \int _ { 0 } ^ { T _ { 0 } } \left | w _ { \text{E} } ^ { ( t ) } ( t ^ { \prime } ) \{ d ( t ^ { \prime } ) - d _ { \text{par} } ( t ^ { \prime } ; c _ { 0 }, c _ { 1 }, c _ { 2 } ) \} \right | ^ { 2 } d t ^ { \prime }, & \text{fore, o} \\ \underset { ( 7 ) } { \text{previous} } & \text{previous} \\ s ^ { \prime \prime } ( t ) \, = }$$

where w ( ) t E ( t ′ ) is a window function that is defined as follows:

$$w _ { \text{E} } ^ { ( t ) } ( t ^ { \prime } ) = \begin{cases} w ( t ^ { \prime } - t ) & \left ( \hat { d } _ { \text{MRCW} } ( t ^ { \prime } ) \geq \gamma \hat { A } ( t ) \right ), & \text{angular} \\ 0 & ( \text{otherwise} ), \end{cases}$$

where γ = 0 6 . is a threshold that was set empirically. In this paper, ˆ ( ) c 2 t is used as a feature that represents the expiratory plateau's shape. Fig. 4 shows the time interval for the expiratory plateau (gray area) and the corresponding estimated parabolic curve (red solid line).

## C. Features of the Respiratory Statistics

In Section II-A, a method to estimate the respiratory features at each instant t was presented. However, respiration includes both voluntary and involuntary components that result in temporal fluctuations in the respiratory feature. These temporal variations can have an adverse effect on the accuracy of individual identification. To mitigate this variability, we extract the statistical properties of the respiratory features.

First, the MRCW parameters ( ˆ ( f t ) , β ˆ ( ) 1 t , β ˆ ( ) 2 t , D t ˆ ( )) are transformed into a six-dimensional vector q ( ) t , as follows:

$$q ( t ) = \begin{bmatrix} q _ { 1 } \\ q _ { 2 } \\ q _ { 3 } \\ q _ { 4 } \\ q _ { 5 } \\ q _ { 6 } \end{bmatrix} = \begin{bmatrix} \hat { f } ( t ) \\ \hat { D } ( t ) \\ \hat { \beta } _ { 1 } ( t ) + \hat { \beta } _ { 2 } ( t ) \\ | \hat { \beta } _ { 1 } ( t ) - \hat { \beta } _ { 2 } ( t ) | \\ \left \{ \hat { \beta } _ { 1 } ( t ) + \hat { \beta } _ { 2 } ( t ) \right \} / \hat { f } ( t ) \\ \hat { c } _ { 2 } ( t ) \end{bmatrix}. \quad ( 9 )$$

In (9), q 1 and q 2 represent the respiratory frequency and the expiratory plateau duration, respectively. Furthermore, q 3 and q 4 represent the sum and the difference, respectively, of the inspiratory and expiratory rates. Additionally, q 5 represents the ratio of the inspiratory and expiratory rates to the sum of the inspiratory and expiratory rates relative to the respiratory frequency. Finally, q 6 is a coefficient of the parabolic curve estimated using (7).

Next, statistics are extracted for each element of the sixdimensional feature vector q ( ) t and transformed into a 24dimensional feature vector r resp as r resp = [ r 1 , r 2 , . . . , r 6 ] T ∈ R 24 × 1 , where r n = [ µ , σ n n , γ n , κ n ] ∈ R 1 × 4 is defined for n = 1 2 , , . . . , 6 ; here, µ n , σ n , γ n , and κ n represent the mean, the standard deviation, the skewness, and the kurtosis of q n , respectively.

## III. HEARTBEAT FEATURES FOR INDIVIDUAL IDENTIFICATION

## A. High-Frequency Enhancement Using Differentiation

In this section, we propose a method to extract the features of the heartbeat component required for individual identification. Because the fundamental frequency component of the heartbeat is often masked by harmonic components of the respiration, it has proven to be advantageous to use the higher harmonics of the heartbeat components [18]-[20]. Therefore, our proposed method emphasizes the high-frequency components by using second-order differentiation, as in our previous work [23]. Because the second-order differentiation s ′′ ( ) t = (d 2 / d t 2 ) s t ( ) is expressed as -ω S ω 2 ( ) in the frequency domain, the high-frequency components, including the harmonics of the heartbeat, are emphasized; here, ω is the angular frequency and S ω ( ) represents the Fourier transform of s t ( ) .

## B. Feature Extraction Using Mel-Frequency Cepstrum

Next, we use the mel-frequency cepstral coefficients (MFCCs) to generate features. Historically, MFCCs have been used to convert the short-term spectra of speech signals into feature vectors [56]. Because radar-measured heartbeat signals

are periodic within a short time, in exactly the same manner as speech signals, use of MFCCs to generate heartbeat features is appropriate when using radar [49]. In [49], MFCCs were applied to real signals that were estimated from the phases of reflected signals. However, our previous work [23] has shown that signal processing using complex reflected signals directly improves the heartbeat estimation accuracy. Therefore, in this section, we propose a method to extract heartbeat features by applying MFCCs to complex signals. The calculation of the MFCCs is performed using the following steps.

1) Short-time Fourier transform: For the differentiated signal s ′′ ( ) t , the short-time Fourier transform (STFT) is applied as follows:

$$S ( t, f ) = \int _ { - \infty } ^ { \infty } s ^ { \prime \prime } ( \tau ) w _ { S T F T } ( \tau - t ) \exp { ( - j 2 \pi f \tau ) } \, d \tau, \ ( 1 0 ) \quad 3 ) \ D _ { i }$$

where w STFT ( ) t is a rectangular window function with a width of T STFT . This window width T STFT should be greater than the heartbeat interval and we set T STFT = 2 0 . s empirically.

2) Mel filters: Next, we apply a mel filter bank to S t, f ( ) to calculate the mel-frequency spectrum. The ℓ th mel filter H f ℓ ( ) is expressed as:

$$\begin{matrix} \imath \vee \vee \infty \text{expression as}. & C _ { - k } \\ H _ { \ell } ( f ) = \begin{cases} \frac { 2 ( f - f _ { \ell } ) } { ( f _ { \ell + 1 } - f _ { \ell } ) ( f _ { \ell + 2 } - f _ { \ell } ) } & ( f _ { \ell } \leq f < f _ { \ell + 1 } ), & \text{for $k=$} \\ \frac { 2 ( f _ { \ell + 2 } - f ) } { ( f _ { \ell + 2 } - f _ { \ell + 1 } ) ( f _ { \ell + 2 } - f _ { \ell } ) } & ( f _ { \ell + 1 } \leq f < f _ { \ell + 2 } ), & C _ { + 0 } \text{ ar} \\ 0 & ( \text{otherwise)} & \text{lower} \\ & \text{(11)} & \text{extract} \end{matrix}$$

(11)

for ℓ = 0 1 , , . . . , L -1 , where

$$f _ { \ell } = \tilde { f } \left \{ \exp \left ( \frac { m _ { \ell } } { \tilde { m } } \right ) - 1 \right \}, \quad \quad ( 1 2 ) \stackrel { \iota \lambda - ( 1 ) } { \text{ and } i } \\ \text{sectity}.$$

is the center frequency of the ℓ th mel filter corresponding to the ℓ th mel scale. Here, the mel scale m ℓ is defined as

$$m _ { \ell } = \tilde { m } \frac { \ell } { L + 1 } \log \left ( 1 + \frac { f _ { s } } { 2 \tilde { f } } \right ), \quad \quad ( 1 3 ) \quad \ln _ { \text{near} \epsilon }$$

for ℓ = 0 1 , , . . . , L + 1 , where f s is the sampling frequency and ˜ m is defined as:

$$\tilde { m } = \frac { f ^ { \prime } } { \log \left ( f ^ { \prime } / \tilde { f } + 1 \right ) }. \quad \quad \ \ ( 1 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

An example of a mel filter bank and the power spectrum of s ′′ ( ) t are shown in Fig. 5 for the case where ˜ = 5 0 f . Hz, f ′ = 1 0 . kHz, and L = 64 , where these parameters were set following the convention in [49].

In this study, we apply the mel filters to a complex radar signal, where the mel filters are applied to the spectrum for both positive and negative frequencies as follows:

$$\dots \sup _ { r } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } c }sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } }sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c  sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ }sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack \\sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \\sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c c c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \subpop { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \subpop < c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup. c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } \sup _ { \substack { c } {$$

for 0 ≤ ℓ ≤ L -1 , where T 0 is the measurement time. Note that in this case, S +0 and S -0 are regarded as different variables.

Fig. 5. Examples of the mel filters and the power spectrum for s ′′ ( ) t .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-b6c1ef67-ea5f-43f3-ba8d-392691074302.jpg)

3) Discrete cosine transform: Next, we obtain the MFCC C k by applying a discrete cosine transform (DCT) to log S ℓ as shown here:

$$\frac { 1 } { 2 } \right ) ^ { 2 } C _ { k } = & \frac { 2 } { L + \delta _ { k, 0 } } \sum _ { \ell = 0 } ^ { L - 1 } \log S _ { \ell } \cos \left \{ \frac { ( 2 \ell + 1 ) k \pi } { K } \right \}, \quad \text{ }$$

$$\text{liter} \quad C _ { - k } = & \frac { 2 } { L + \delta _ { - k, 0 } } \sum _ { \ell = 0 } ^ { L - 1 } \log S _ { - \ell } \cos \left \{ \frac { - ( 2 \ell + 1 ) k \pi } { K } \right \}, \quad ( 1 8 )$$

$$& \frac { 2 } { + \delta _ { k, 0 } } \sum _ { \ell = 0 } ^ { L - 1 } \log S _ { \ell } \cos \left \{ \frac { ( 2 \ell + 1 ) k \pi } { K } \right \}, \quad ( 1 7 ) \\ & \frac { 2 } { + \delta _ { - k, 0 } } \sum _ { \ell = 0 } ^ { L - 1 } \log S _ { - \ell } \cos \left \{ \frac { - ( 2 \ell + 1 ) k \pi } { K } \right \}, \quad ( 1 8 ) \\ & + \tau - 1 \dots \dots \tilde { \kappa } \quad \kappa \quad \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \, \colon \end{array}$$

for k = 0 1 , , . . . , K -1 , where δ i,j is the Kronecker delta, K is the cepstrum dimension, and we set K = 64 . Note that C +0 and C -0 are also regarded as different variables here.

Finally, from the 2 K -dimensional MFCC C k , the lower order coefficients up to 2 K ′ dimensions are extracted, where the condition 0 &lt; K ′ &lt; K holds; the feature vector of the heartbeat r hb = [ C -( K ′ -1) , C -( K ′ -2) , . . . , C -0 , C +0 , . . . , C K ′ -1 ] T is defined and is used to perform individual identification in the following section, where we set K ′ = 24 empirically.

## IV. PERFORMANCE EVALUATION OF THE PROPOSED METHOD

In this section, we evaluate the individual identification performance of the proposed method via comparison with conventional methods. To identify an effective set of features, we compare the performance metrics obtained when the proposed feature vectors are used separately and together. We compare three cases: (1) use of the respiratory feature vector r resp only; (2) use of the heartbeat feature vector r hb only; and (3) use of both vectors as r prop = [ r T resp , r T hb ] T . We also compare three classifiers: the support vector machine (SVM), k -nearest neighbors ( k -NN), and multilayer perceptron (MLP) classifiers. For these three classifiers, we use the three feature vector types, i.e., r resp , r hb , and r prop , to compare a total of nine different approaches. Table I summarizes the definitions of nine methods for comparison; for example, method A1 indicates the identification process when using r resp and the SVM.

## A. Experiment 1: Five-Day Measurement of Six Participants

In this experiment, measurements were conducted over five days with six healthy adult participants. Each day, measurements were taken five times per participant in the morning

TABLE I DEFINITIONS OF METHODS FOR INDIVIDUAL IDENTIFICATION

| Method                        | Signal                        | Vital sign              | Feature   | Classifier    |
|-------------------------------|-------------------------------|-------------------------|-----------|---------------|
| Method A1 Method A2 Method A3 | ∠ s ( t )                     | Respiration             | r resp    | SVM k -NN MLP |
| Method B1 Method B2 Method B3 | d 2 d t 2 s ( t )             | Heartbeat               | r hb      | SVM k -NN MLP |
| Method C1 Method C2 Method C3 | ∠ s ( t ) & d 2 d t 2 s ( t ) | Respiration & heartbeat | r prop    | SVM k -NN MLP |

TABLE II PHYSICAL CHARACTERISTICS OF EACH PARTICIPANT

|   Participant |   Age (yr) |   Height (m) |   Weight (kg) |
|---------------|------------|--------------|---------------|
|             1 |         26 |         1.67 |          58.9 |
|             2 |         22 |         1.76 |          69.7 |
|             3 |         23 |         1.75 |          75.2 |
|             4 |         23 |         1.8  |          62   |
|             5 |         23 |         1.7  |          68   |
|             6 |         23 |         1.69 |          54   |

and another five times per participant in the afternoon, with each session lasting for 60 0 s . . Therefore, the total number of measurements performed was (5+5) × 5 = 50 times per each person, giving a total of 300 measurements. In other words, a total of 300 samples were taken. All participants were in good health with no history of respiratory problems or heart disease. The physical characteristics of each of the participants are listed in Table II. All measurements were conducted with each participant seated in front of the radar system at a distance of 1 5 m . . The experimental setup is pictured in Fig. 6 (a), and a schematic that indicates the physical dimensions is shown in Fig. 6 (b).

The radar specifications are listed in Table III. We used a frequency-modulated continuous-wave (FMCW) radar system with a center frequency of 79 0 . GHz and a bandwidth of 3 6 GHz . . The antenna array consists of three transmitting and four receiving elements, which form a MIMO array that can be approximated virtually as a 12-element linear array. The range and the direction of arrival of the echo from the target are identified by performing a Fourier transform in terms of the fast time [57] and the virtual element number corresponding to the beamformer method [58].

The feature vectors r resp , r hb , and r prop are visualized by compressing them into two dimensions using the t-distributed stochastic neighbor embedding (t-SNE) method, [59] as shown in Fig. 7, which shows that the feature vectors form clusters corresponding to the six participants. In particular, when compared with the other feature vectors, the clusters are separated clearly when using r prop . In the r resp case, the cluster for participant 6 is erroneously split into two separate clusters and also overlaps with other clusters.

The accuracy during individual identification was evaluated using a stratified k -fold cross-validation approach with k = 10 , thus ensuring that each subset maintained the class distribution

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-8e14d6b3-d380-4186-8c4f-d52eb79a3d38.jpg)

(a)

Fig. 6. (a) Photograph and (b) schematic of the experimental setup showing the radar system and a seated participant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-39600394-b0ce-46e3-b0c2-80d48e3b75c2.jpg)

TABLE III SPECIFICATIONS OF THE RADAR SYSTEM

| Parameter                               | Value          |
|-----------------------------------------|----------------|
| Modulation                              | FMCW           |
| Center frequency                        | 79 . 0 GHz     |
| Center wavelength                       | 3 . 8 mm       |
| Bandwidth                               | 3 . 6 GHz      |
| No. of transmitting (Tx) antennas       | 3              |
| No. of receiving (Rx) antennas          | 4              |
| Tx element spacing                      | 7 . 6 mm       |
| Rx element spacing                      | 1 . 9 mm       |
| Beamwidths of Tx elements (E-/H-planes) | ± 4 ◦ / ± 33 ◦ |
| Beamwidths of Rx elements (E-/H-planes) | ± 4 ◦ / ± 45 ◦ |
| Range resolution                        | 44 . 7 mm      |
| Sampling frequency (slow-time)          | 100 Hz         |

of the original dataset; each cross-validation iteration used 270 samples for training and 30 samples for validation. As a singlefold validation, only data acquired from the same period can be used as validation data. For example, in the single-fold method, data measured on the morning of the first day are used as the validation data and the remaining data serve as training data.

The three feature vectors r resp , r hb , and r prop are input into the three classifiers to compare their identification accuracies. Here, because both SVM and k -NN are binary classifiers, the one-versus-rest method is used to perform multi-class classification [60]. In addition, the activation function used for the MLP output layer is the softmax function. The parameters for each classifier are optimized via a grid search. The optimization processes are performed for each of the three feature

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-97ef373f-a855-4d25-ba55-ea4844f677d8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-c66947e5-dfab-4bb4-a043-de2be45e1bc0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-27cae0b3-80a1-494d-8901-84b52f7b8b11.jpg)

t-SNE 2

Fig. 7. Two-dimensional t-SNE visualizations of the feature vectors (a) r resp , (b) r hb , and (c) r prop for the 60 s data samples acquired in experiment 1.

vectors r resp , r hb , and r prop individually. The parameters to be optimized and the values selected are given in Table IV.

The confusion matrices obtained for each classifier are shown in Fig. 8. The accuracy indicates the proportion of correctly identified cases with respect to the total number of cases. When using r resp , r hb , and r prop , these three classifiers achieved accuracies of approximately 80%, 90%, and 95%, respectively. These results suggest that the heartbeat feature is superior to the respiratory feature when used in individual identification, and the combination of the respiratory and heartbeat components is superior to both the respiratory feature and the heartbeat feature when they are used separately. In Fig. 8, we see that method C1 realizes the highest accuracy of 96.33% from among the nine methods.

Fig. 9 presents the receiver operating characteristic (ROC) curves and the corresponding area under the curve (AUC) results for each classifier. Note that the ROC curves shown in Fig. 9 are the macro-averages of the one-versus-rest ROC curves for each class. Fig. 9 shows that when r prop is used, the AUC for all classifiers reaches 0.99 and is higher than that realized when using r resp and r hb separately.

Table V presents the F 1 scores and the AUCs for each participant when r prop is used as the input, along with their respective macro-averages. From Table V, we see that the SVM achieves the highest F 1 scores and AUCs among the classifiers. The AUC exceeds a macro average of 0.99 for all classifiers, thus indicating that each classifier achieves comparable accuracy to the others. This indicates that the proposed features enable high-accuracy identification, regardless of the classifier that is used.

Next, to investigate the individual identification performance for short data samples, we divide each 60 s data sample into 12 short data samples with a length of 5.0 s, thus generating 3,600 samples in total. For these shorter samples, we evaluate the performance for the heartbeat feature r hb only (corresponding to methods B1, B2, and B3) because the data sample length of 5.0 s is not sufficient to capture respiratory features, which have a typical interval of between 3 s and 5 s.

The feature vector r hb for the short data samples is visualized in two-dimensional space using the t-SNE method in Fig. 10. Although the figure shows six clusters corresponding to the six participants, these clusters overlap severely when compared with those shown in Fig. 7 (b); this is an indication of the performance degradation caused by the limited data sample length. The confusion matrices obtained using the stratified k -fold cross-validation approach are shown in Fig. 11 for k = 10 for the 5 s data samples. The accuracies obtained for methods B1, B2, and B3 with the 5 s data samples were 85.64%, 83.11%, and 84.19%, respectively, showing that the performance has been reduced by 7.69, 9.22, and 3.81 points when compared with the accuracies of 93.33%, 92.33%, and 88.00% for B1, B2, and B3, respectively, obtained with the 60 s data samples. We also confirmed that the AUCs obtained for methods B1, B2, and B3 with the 5 s data samples were 0.979, 0.968, and 0.962, respectively, which were lower than the AUCs obtained with the 60 s data samples by 0.017, 0.024, and 0.022, respectively.

## B. Experiment 2: Public Dataset of Thirty Participants

Next, we evaluated the accuracy of the proposed method using a public dataset provided by Schellenberger et al. [61], where they used a continuous-wave (CW) radar system with an operating frequency of 24 GHz to collect echoes from a total of 30 individuals, comprising 14 male and 16 female participants. The echoes were recorded with a sampling frequency of 2 0 kHz . , which was then downsampled to 100 Hz to make the conditions consistent with those of experiment 1. The measurement time varied for each participant and ranged from approximately 600 s to 700 s . In this study, we divided the data into segments with a length of T 0 and examined two cases, where T 0 = 60 0 s . and T 0 = 5 0 s . ; therefore, totals of 306 and 3,795 samples were taken, respectively. A feature vector was then generated from each sample.

Fig. 12 visualizes two-dimensional plots of the feature vector r hb as obtained from data segments with T 0 = 60 0 s . and T 0 = 5 0 s . when compressed using t-SNE. In both cases, it was observed that almost all plots were clearly separated into distinct clusters corresponding to the 30 participants. Next, the feature vector r hb obtained was input into the SVM and the identification accuracy was evaluated by performing a stratified k -fold cross-validation with k = 5 ; the resulting confusion matrices are shown in Fig. 13. First, when T 0 = 60 0 s . , 100% accuracy was achieved for all but four participants, and the overall accuracy was 98.69%, thus indicating that

TABLE IV OPTIMIZED PARAMETERS FOR EACH METHOD

| Classifier      | Parameters                                                              | Options                                                  | Feature vector       | Feature vector     | Feature vector     |
|-----------------|-------------------------------------------------------------------------|----------------------------------------------------------|----------------------|--------------------|--------------------|
|                 |                                                                         |                                                          | A) r resp            | B) r hb            | C) r prop          |
| 1) SVM 2) k -NN | Kernel function No. of neighbors Distance                               | Linear / Gaussian [1 , 300] Euclidean, Cityblock, Cosine | Gaussian 9 Cityblock | Gaussian 28 Cosine | Gaussian 28 Cosine |
| 3) MLP          | No. of layers Size of 1 st layer Size of 2 nd layer Activation function | [1 , 3] [1 , 300] [1 , 300] ReLU, † Sigmoid              | 2 39 19 ReLU †       | 2 47 49 ReLU †     | 2 15 15 ReLU †     |

- † : Rectified linear unit

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-a395cbef-4ac4-4d62-a235-4b7d1d4d4fdc.jpg)

(a)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-fe61e394-ce2b-47e6-b6a1-bc8885dd73ad.jpg)

(b)

Fig. 8. Confusion matrices for each method (a) SVM, (b) k -NN, and (c) MLP for the 60 s data samples acquired in experiment 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-b6ba5ff9-3968-44ef-a6b5-64c6dc16fd37.jpg)

individual identification was achieved successfully. Next, when T 0 = 5 0 s . , the overall accuracy realized was 99.39%, which

Fig. 9. ROC curves for (a) SVM, (b) k -NN, and (c) MLP for the 60 s data samples acquired in experiment 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-e3400099-0e67-4618-b610-7c4672add0cf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-590719dd-880f-42ea-8f64-87825afe1ed1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-c8e2da0a-6974-4d61-a171-486f09e70f3e.jpg)

TABLE V F 1 SCORE AND AUC FOR EACH PARTICIPANT USING r prop IN EXP. 1

| Classifier      | Participant   |   F 1 score |   AUC |
|-----------------|---------------|-------------|-------|
| SVM (method C1) | 1             |       0.959 | 0.994 |
|                 | 2             |       0.99  | 0.999 |
|                 | 3             |       0.904 | 0.995 |
|                 | 4             |       0.969 | 1     |
|                 | 5             |       0.961 | 0.999 |
|                 | 6             |       1     | 1     |
|                 | macro-ave.    |       0.964 | 0.998 |
| -NN (method C2) | 1             |       0.938 | 0.998 |
| -NN (method C2) | 2             |       0.96  | 0.997 |
| -NN (method C2) | 3             |       0.898 | 0.996 |
|                 | 4             |       0.961 | 0.998 |
|                 | 5             |       0.952 | 0.999 |
|                 | 6             |       0.99  | 1     |
|                 | macro-ave.    |       0.95  | 0.998 |
| MLP (method C3) | 1             |       0.951 | 0.999 |
|                 | 2             |       0.96  | 0.997 |
|                 | 3             |       0.86  | 0.979 |
|                 | 4             |       0.97  | 0.999 |
|                 | 5             |       0.923 | 0.997 |
|                 | 6             |       0.99  | 0.999 |
|                 | macro-ave.    |       0.942 | 0.995 |

also represents high-precision identification.

In experiment 1, the accuracy of method B2 was only 93.33%, whereas in experiment 2, the same method B2 achieved accuracy of 98.69%, partly because the public dataset was collected within a day and was thus unaffected by daily variations in the physiological signals. We should also note that the accuracy obtained for T 0 = 5 0 s . was almost the same as the accuracy for T 0 = 60 0 s . , which may indicate that long measurements are not necessary during individual identification if the daily variations in the physiological signals can be neglected.

## V. DISCUSSION

## A. Effect of Data Length

In experiment 1, the accuracy for the 5 s data samples was lower than that for the 60 s data samples. In contrast,

Fig. 10. Two-dimensional t-SNE visualizations of the feature vector r hb for the 5 s data samples acquired in experiment 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-e2d11704-99ec-4dfb-90be-6016cb51c7cb.jpg)

in experiment 2, the accuracy remained high even when the 5 s data samples were used. The difference between the two experiments is that the measurements in experiment 1 were taken repeatedly, twice a day over five days, and thus the respiratory and heartbeat features were not completely consistent over the datasets that were acquired on different days. In contrast, the measurements in experiment 2 were all performed over a single day, and thus the respiratory and heartbeat features remained relatively consistent over the datasets that were used for training and testing. Note that in experiment 1, when datasets acquired in the morning on day 1 were used for testing, we then selected training datasets acquired in the afternoon on day 1 in addition to those from days 2, 3, 4, and 5; we thus avoided using datasets that had been acquired in the same half-day session for both training and testing, which made the evaluation conditions used in experiment 1 difficult but practical.

In addition, the measurement setups used in experiments 1 and 2 were different; in experiment 1, the participants were seated, whereas in experiment 2, the participants were lying down on a bed. As a result, there were larger participant body movements in experiment 1 than in experiment 2, which made

Fig. 11. Confusion matrices for each method: (a) SVM, (b) k -NN, and (c) MLP for the 5 s data samples acquired in experiment 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-60698d5d-f144-44eb-a899-eaa37d54e800.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-5b35826f-23b4-47a5-8185-638167abd07f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-cdf3aa2e-c6c8-40e7-9ffe-9302279a2b7b.jpg)

Fig. 12. Two-dimensional plots of feature vectors r hb compressed by t-SNE when using the public dataset with data lengths (a) T 0 = 60 0 . s and (b) T 0 = 5 0 . s in experiment 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-1eb34320-b260-4551-9d13-f5be7bb61f91.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-87d736f8-75d8-45b4-99c9-e757280a2c7e.jpg)

problem setting in experiment 1 more challenging. For these reasons, individual identifications in experiment 1 required longer data lengths than in experiment 2. In future work, it will be important to develop a method for individual identification when using short data samples, even when the participants are seated and the datasets for training and testing are acquired on different days.

and 80.67% for the SNRs of 10.0 dB, 20.0 dB, and 30.0 dB, respectively.

## B. Effect of Background Noise

Next, we investigate the effect of background noise on the identification performance by applying the methods to data with a low signal-to-noise power ratio (SNR). The low-SNR data required were generated by adding white Gaussian noise to the measured data from in experiment 1, where this noise was approximated using a random complex-valued sequence following a normal distribution. Note here that the SNR of the data acquired in experiment 1 was originally 36.9 dB.

Table VI lists the accuracy and the AUC for each method with SNRs of 10.0, 20.0, and 30.0 dB. Using the original data with the high SNR of 36.9 dB, the average accuracy obtained for methods A1, A2, and A3 was 81.67%; when using the data with the low SNR of 10.0 dB, the average accuracy for methods A1, A2, and A3 was 80.56%, which indicates that a similar level of accuracy is maintained when using the respiratory features, even with the low-SNR data. For example, the accuracies for method A1 were 81.67%, 82.67%,

In contrast, when using the original data with the SNR of 36.9 dB, the average accuracy for methods B1, B2, and B3 was 91.22%; when using the data with the SNR of 10.0 dB, the average accuracy for methods B1, B2, and B3 was 78.33%, indicating that the averaged accuracy decreased by 12.89 points. For example, the accuracies for method B1 when using the data with the SNRs of 10.0 dB, 20.0 dB, and 30.0 dB were 79.33%, 85.00%, and 87.33%, respectively; these values were lower than the accuracies obtained for the data with the high SNR of 36.9 dB by 14.00, 8.33, and 6.00 points, respectively. These results indicate that the methods that use the heartbeat features (i.e., methods B1, B2, and B3) are affected more severely by the background noise when compared with the the methods that use respiratory features (i.e., methods A1, A2, and A3). This is attributed to the fact that the body displacement caused by the heartbeat is much smaller than that caused by respiration; thus, the heartbeat features are more susceptible to noise.

## C. Effects of Clothing and Temperature

When the measurements were conducted in experiments 1 and 2, all participants were dressed. It has been reported previously that the reflection from clothing is small enough to be considered negligible in measurements performed using

Fig. 13. Confusion matrices for experiment 2 with a data length (a) T 0 = 60 s and (b) T 0 = 5 s .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-75d56339-0a0c-44c6-8b60-283753431dd0.jpg)

TABLE VI ACCURACY AND AUC VALUES FOR LOW SNR CONDITIONS

| Method   | Accuracy (%)   | Accuracy (%)   | Accuracy (%)   | AUC   | AUC   | AUC   |
|----------|----------------|----------------|----------------|-------|-------|-------|
| Method   | 10 dB          | 20 dB          | 30 dB          | 10 dB | 20 dB | 30 dB |
| A1       | 81.67          | 82.67          | 80.67          | 0.962 | 0.965 | 0.968 |
| A2       | 78.67          | 79.00          | 78.67          | 0.948 | 0.958 | 0.944 |
| A3       | 81.33          | 78.67          | 80.33          | 0.935 | 0.930 | 0.947 |
| B1       | 79.33          | 85.00          | 87.33          | 0.951 | 0.972 | 0.981 |
| B2       | 75.33          | 81.67          | 86.33          | 0.931 | 0.968 | 0.978 |
| B3       | 80.33          | 83.33          | 83.67          | 0.938 | 0.963 | 0.965 |
| C1       | 88.67          | 92.00          | 92.33          | 0.990 | 0.992 | 0.993 |
| C2       | 86.33          | 90.33          | 93.00          | 0.975 | 0.988 | 0.993 |
| C3       | 87.67          | 91.33          | 92.33          | 0.978 | 0.985 | 0.992 |

millimeter-wave radar [62]-[64]. Therefore, it is unnecessary to consider the effects of vibration of clothes in processing the radar data. Even if there is an effect from the participant's clothes, our proposed method achieved individual identification successfully for the dressed participants, which indicates the effectiveness of the radar-based approach.

The radar-based physiological measurement principle is based on the detection of small displacements caused by the subject's heartbeat, respiration, and body motion. Therefore, the radar measurements are affected by neither the sample temperature nor the room temperature. However, as the participant's body temperature increases, their respiratory rate and heart rate also increase in general [65], which may then affect the method's accuracy when identifying multiple individuals.

## D. Performance Under Different Conditions

In this section, we evaluate the individual identification performance under eight different conditions, designated α 1, α 2, β 1, β 2, γ 1, γ 2, δ , and ε , which are defined in Table VII. We performed additional radar measurements five times for both participants 3 and 4 from experiment 1, and the additional data were used for testing only, and not for training; the classifiers were trained using the original data from experiment 1.

Fig. 14 shows the feature vectors r resp , r hb , and r prop that were obtained from the additional measurements and then compressed into two dimensions using the t-SNE method; Table VIII presents the confusion matrices for each of the eight conditions. The additional radar measurements were

performed nine months after the original radar measurements described in Section IV-A.

1) Condition α : Target Distance: First, we evaluated the performance at various distances between the radar and the target person. Conditions α 1 and α 2 correspond to distances of 1.0 m and 3.0 m, respectively, although the distance was originally set at 1.5 m in experiment 1.

Table VIII shows that although the accuracy under condition α 1 was marginally low, the accuracy under condition α 2 was high, and Fig. 14 shows that the feature vectors r resp are located far away from the corresponding cluster under condition α 1. The performance of method A1 in particular is affected negatively under condition α 1. These results may be caused by the shift in the radar reflection points on the target person's body, depending on the distance to the target, which leads to discrepancies in the respiratory feature vectors that are used for training and testing.

2) Condition β : Radar Height: Next, we evaluated the influence of the radar height from the floor; under conditions β 1 and β 2, the radar was installed at heights of 0.8 m and 1.0 m, respectively, whereas the radar was originally installed at a height of 0.9 m from the floor. As Table VIII shows, the accuracy of method A1 is low under condition β 1, whereas the accuracy of method B1 remains high under both conditions β 1 and β 2. In addition, as illustrated in Fig. 14, the respiratory feature vector r resp under condition β 1 is located near the original cluster, which indicates that the individual characteristics of each participant were captured under both conditions. Therefore, it can be deduced that the radar height does not affect the individual identification performance significantly.

3) Condition γ : Seating Direction: Next, we evaluated the effect of the target person's seating direction on the identification performance. Under conditions γ 1 and γ 2, we performed measurements for participants when seated facing at angles of 45 ◦ and 90 ◦ , respectively, relative to the radar system; the corresponding experimental setups are illustrated in Fig. 15. The results in Table VIII show that the accuracy decreases under both conditions γ 1 and γ 2, with none of the methods achieving high performances. Fig. 14 also indicates that the feature vectors determined under conditions γ 1 and γ 2 do not capture the individual features of each participant.

The results above can be explained as follows. The body displacements caused by the respiratory and heartbeat motions are particularly large at the front of the torso, whereas the displacements measured under conditions γ 1 and γ 2 are smaller because the echo is likely to come from a different side of the torso or from the shoulder/arm area. The performance evaluation results under conditions γ 1 and γ 2 indicate that the seating direction can affect the individual identification accuracy significantly.

4) Condition δ : Wall Interference: In condition δ , we evaluated the effects of clutter reflected from the mortar wall behind the participant. As shown in Fig. 6, the wall was located comparatively far behind the target person in the original measurement setup used in experiment 1. We therefore performed radar measurements where each participant's back is in contact with the wall, corresponding to a distance of 0 m from the wall. Under this condition, although the accuracy

TABLE VII DEFINITIONS OF MEASUREMENT CONDITIONS

| Condition                   | Parameter               | Value       |
|-----------------------------|-------------------------|-------------|
| Condition α 1 Condition α 2 | Target distance         | 1.0 m 3.0 m |
| Condition β 1 Condition β 2 | Radar height            | 0.8 m 1.0 m |
| Condition γ 1 Condition γ 2 | Seating direction       | 45 ◦ 90 ◦   |
| Condition δ                 | Wall interference       | -           |
| Condition ε                 | Exhaustion and sweating | -           |

decreases in the manner shown in Table VIII, Fig. 14 (c) demonstrates that the feature vector r prop remains close to the original cluster because the radar system used in experiment 1 has a wide bandwidth of 3.6 GHz, which provides a range resolution of 41.7 mm that is sufficient to prevent the wall clutter from contaminating the target person's echo.

5) Condition ε : Exhaustion and Sweating: In condition ε , we evaluated the performance when the participant is exhausted and sweating; immediately before the measurements were taken, the participants were instructed to run up 8 flights of stairs. During the measurements, the participants were out of breath and sweating as a result. Under condition ε , the accuracy was low for all methods, as shown in Table VIII. Fig. 14 (c) shows that the feature vectors r prop are all far from the original clusters, which indicates that the unique individual features were not maintained after the intense exercise. The significant increases observed in both the respiratory rate and the heart rate after exercise are likely to have a considerable impact on the performance of the proposed method. It can be concluded from these results that the proposed method is most suitable for use when the target person is at rest.

## E. Comparison with Related Works

Here we compare the results from this study with those from existing publications, as shown in Table IX. Although these studies assumed different measurement conditions (e.g., sensor types, the number of antennas used, transmitting power, distance between the antenna and the target person), we can see that the proposed method achieves an accuracy that is comparable to that reported in the existing studies.

In particular, the work of Peng et al. [49] is similar to that of this study because they also used the respiratory and heartbeat features for individual identification; they used the respiratory features for individual identification, and they also used the heartbeat features only during apnea events because the respiratory signals cannot be obtained during apnea events. In contrast, this study improves the overall identification accuracy by using a combination of the respiratory and heartbeat features simultaneously, which differentiates this study from that of Peng et al. [49].

## VI. CONCLUSION

In this study, we developed a radar-based individual identification method using both respiratory and heartbeat features.

Fig. 14. Two-dimensional t-SNE visualizations of the feature vectors (a) r resp , (b) r hb , and (c) r prop under each condition (Conditions: α 1: +; α 2: ∗ ; β 1: | ; β 2: x; γ 1: △ ; γ 2: □ ; δ : ♢ ; ε : ✩ ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-990c1d79-3486-439c-a5af-09234104456d.jpg)

TABLE VIII

## CONFUSION MATRICES UNDER DIFFERENT MEASUREMENT CONDITIONS

|            |            | Predicted   | Predicted   | Predicted   | Predicted   | Predicted   | Predicted                 | Predicted                 | Predicted                 | Predicted                 | Predicted                 | Predicted                 | Predicted              | Predicted Label        | Predicted Label        | Predicted Label        | Predicted Label        | Predicted Label        | Predicted Label        | Predicted Label        | Predicted Label       | Predicted Label       | Predicted Label       | Predicted Label       | Predicted Label       | Predicted Label       |
|------------|------------|-------------|-------------|-------------|-------------|-------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|------------------------|------------------------|------------------------|------------------------|------------------------|------------------------|------------------------|------------------------|-----------------------|-----------------------|-----------------------|-----------------------|-----------------------|-----------------------|
|            |            | Condition α | Condition α | Condition α | Condition α | Condition α | Label 1 (1.0 m) Condition | Label 1 (1.0 m) Condition | Label 1 (1.0 m) Condition | Label 1 (1.0 m) Condition | Label 1 (1.0 m) Condition | Label 1 (1.0 m) Condition | α 2 (3.0 m) Condition  | α 2 (3.0 m) Condition  | α 2 (3.0 m) Condition  | α 2 (3.0 m) Condition  | α 2 (3.0 m) Condition  | α 2 (3.0 m) Condition  | β 1 (0.8 m) Condition  | β 1 (0.8 m) Condition  | β 1 (0.8 m) Condition | β 1 (0.8 m) Condition | β 1 (0.8 m) Condition | β 1 (0.8 m) Condition | β 1 (0.8 m) Condition | β 1 (0.8 m) Condition |
| True Label |            | 1           | 2           | 3           | 4           | 5           | 6                         | 1                         | 2                         | 3                         | 4                         | 5                         | 6                      | 1                      | 2                      | 3                      | 4                      | 5                      |                        | 6                      | 1                     | 2                     | 3 4                   |                       | 5                     | 6                     |
| Method     | 3          | 1           | -           | 2           | 1           | -           | 1                         | 1                         | -                         | 4                         | -                         | -                         | -                      | 2                      | 1                      | 2                      | -                      | -                      |                        | -                      | -                     | -                     | 5                     | -                     | -                     | -                     |
| A1         | 4          | -           | 1           | -           | 3           | -           | 1                         | -                         | 1                         | -                         | 4                         | -                         | -                      | -                      | -                      | -                      | 5                      | -                      |                        | -                      | -                     | -                     | -                     | 5                     | -                     | -                     |
| Method     | 3          | -           | -           | 5           | -           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      | -                      |                        | -                      | 1                      | -                      |                        | -                     | -                     | 5                     | -                     | -                     | -                     |
| B1         | 4          | -           | -           | -           | 5           | -           | -                         | -                         | -                         | -                         | 4                         | -                         | 1                      | -                      | -                      |                        | 4 -                    | 4                      | -                      | 1                      | -                     | -                     | -                     | 4                     | -                     | 1                     |
| Method     | 3          | -           | -           | 5           | -           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      |                        | -                      | -                      | 1                      | -                      |                        | -                     | -                     | 5                     | -                     | -                     | -                     |
| C1         | 4          | -           | 1           | -           | 3           | -           | 1                         | -                         | -                         | -                         | 5                         | -                         | -                      | -                      | -                      |                        | 4 -                    | 5 -                    | -                      |                        | -                     | -                     | -                     | 5                     | -                     | -                     |
|            | Condition  | Condition   | Condition   | Condition   | Condition   | Condition   | γ 1 ( 45 ◦ ) Condition    | γ 1 ( 45 ◦ ) Condition    | γ 1 ( 45 ◦ ) Condition    | γ 1 ( 45 ◦ ) Condition    | γ 1 ( 45 ◦ ) Condition    | γ 1 ( 45 ◦ ) Condition    | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | γ 2 ( 90 ◦ ) Condition | Condition             | Condition             | Condition             | Condition             | Condition             | Condition             |
| True Label | True Label |             | 1 2         | 3           | 4           | 5           | 6                         | 1                         | 2                         | 3                         | 4                         | 5                         | 6                      | 1                      |                        | 2                      | 3                      | 4                      | 5                      | 6                      |                       | 1 2                   | 3                     | 4                     | 5                     | 6                     |
| Method     | 3          | -           | -           | 4           | -           | 1           | -                         | 2                         | -                         | -                         | 1                         | 2                         | -                      | -                      |                        | -                      | 5                      | -                      | -                      | -                      | -                     |                       | -                     | 1                     | -                     | 3                     |
| A1         | 4          | 3           | -           | -           | 2           | -           | -                         | 4                         | -                         | -                         | 1                         | -                         | -                      | -                      | -                      |                        | 2                      | 1 -                    | 2                      |                        | -                     | 1 1                   | -                     | 3                     | -                     | 1                     |
| Method     | 3          | -           | -           | 5           | -           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      | -                      |                        | 4                      | -                      | -                      |                        | -                     | -                     | 4                     | -                     | 1                     | -                     |
| B1         | 4          | -           | -           | 3           | 2           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      | 1                      |                        | -                      | 3                      | 1 1                    | -                      | -                     | 2                     | 2                     | -                     | -                     | 1                     |
| Method     | 3          | -           | -           | 5           | -           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      | -                      | 5                      | -                      | -                      |                        |                        | -                     | -                     | -                     |                       | -                     | 3                     |
| C1         | 4          | -           | 1           | 2           | 2           | -           | -                         | -                         | -                         | 5                         | -                         | -                         | -                      | -                      | 3                      | -                      |                        | 2                      | -                      | -                      | -                     | -                     | 2 3                   | - -                   | -                     | 2                     |

TABLE IX COMPARISON WITH RELATED WORKS

| Study              | Vital sign    | System   | No. of subjects   | Accuracy    |
|--------------------|---------------|----------|-------------------|-------------|
| Hwang et al. [48]  | HB            | CW       | 30                | 88.9%       |
| Peng et al. [49]   | Resp. & HB    | FMCW     | 5                 | 94.4%       |
| Korany et al. [30] | Resp.         | Wi-Fi    | 20                | 95.0%       |
| Wang et al. [44]   | Resp.         | FMCW     | 37                | 96.0%       |
| Lin et al. [51]    | HB            | CW       | 78                | 98.6%       |
| This study         | Resp. & HB HB | FMCW CW  | 6 30              | 96.3% 99.4% |

The proposed method extracts the subject's respiratory features using a modified raised-cosine-waveform model. In addition, a parabolic model is used to express the waveform during the expiratory plateau phase. The proposed method extracts the heartbeat features by applying a time-frequency analysis and a mel-frequency cepstral analysis to the second time-derivative of the complex radar signal. To identify a suitable combination of features and a classifier, we compared the performances of

Fig. 15. Schematics of the experimental setup with the radar system and seated participants under (a) condition γ 1 and (b) condition γ 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-b3cf7625-2bc6-4177-a1d9-79454a323cb2.jpg)

nine methods based on combinations of three feature vectors, comprising the respiratory feature, the heartbeat feature, and the respiratory and heartbeat feature, and three classifiers,

comprising the support vector machine, k -nearest neighbors, and multilayer perceptron methods. The accuracy of the proposed method during individual identification was evaluated in two scenarios; one scenario involved five-day radar measurements of six participants conducted by the authors, and the other scenario involved one-day radar measurement of thirty participants, which was provided as a public dataset. Through performance evaluation, we demonstrated that the combination of the respiratory and heartbeat features achieved accuracy of 96.33% for the five-day radar data from six participants, and accuracy of 99.39% was realized for a public dataset composed of thirty participants using a short data segment that lasted for only 5.0 s; these results indicate the effectiveness of the proposed approach for radar-based individual identification applications.

Although our study evaluated the effectiveness of the proposed method using controlled experiments involving six participants (experiment 1) and publicly available data (experiment 2), it is important to note that these experiments were conducted under highly specific and controlled conditions. Therefore, there are potential limitations in terms of both the validity of the proposed method and in the general applicability of the findings reported in this study. For further evaluation, it will be essential to assess the performance of the proposed method in more realistic scenarios, i.e., where the participants are unaware of the radar measurements and thus behave naturally. This will form a critical part of our future research.

## ACKNOWLEDGEMENT

This work was supported in part by SECOM Science and Technology Foundation, in part by JST under Grants JPMJMI22J2 and JPMJMS2296, in part by JSPS KAKENHI under Grants 21H03427, 23H01420, 23K19119, 23K26115 and 24K17286, and in part by the New Energy and Industrial Technology Development Organization (NEDO). This work involved human subjects in its research. Approval of all ethical and experimental procedures and protocols was granted by the Ethics Committee of the Graduate School of Engineering, Kyoto University, under Approval No. 202223.

## REFERENCES

- [1] M. Ghafurian, K. Wang, I. Dhode, M. Kapoor, P. P. Morita, and K. Dautenhahn, 'Smart home devices for supporting older adults: A systematic review,' IEEE Access , vol. 11, pp. 47137-47158, Apr. 2023, DOI: 10.1109/ACCESS.2023.3266647.
- [2] S. C. Mukhopadhyay, S. K. S. Tyagi, N. K. Suryadevara, V. Piuri, F. Scotti, and S. Zeadally, 'Artificial intelligence-based sensors for next generation IoT applications: A review,' IEEE Sensors J. , vol. 21, no. 22, pp. 24920-24932, Jan. 2021, DOI: 10.1109/JSEN.2021.3055618.
- [3] F. Schroff, D. Kalenichenko, and J. Philbin, 'FaceNet: A unified embedding for face recognition and clustering,' in Proc. IEEE Conf. Computer Vis. Pattern Recognit. (CVPR) , Boston, MA, USA, Jun. 2015, pp. 815-823, DOI: 10.1109/CVPR.2015.7298682.
- [4] R. He, J. Cao, L. Song, Z. Sun, and T. Tan, 'Adversarial cross-spectral face completion for NIR-VIS face recognition,' IEEE Trans. Pattern Anal. Mach. Intell. , vol. 42, no. 5, pp. 1025-1037, May 2020, DOI: 10.1109/TPAMI.2019.2961900.
- [5] M. M. H. Ali, V. H. Mahale, P. Yannawar, and A. T. Gaikwad, 'Fingerprint recognition for person identification and verification based on minutiae matching,' in Proc. IEEE 6th Int. Conf. Advanced Comput. (IACC) , Bhimavaram, India, Feb. 2016, pp. 332-339, DOI: 10.1109/IACC.2016.69.
- [6] A. Kumar and K. V. Prathyusha, 'Personal authentication using hand vein triangulation and knuckle shape,' IEEE Trans. Image Process. , vol. 18, no. 9, pp. 2127-2136, Sep. 2009, DOI: 10.1109/TIP.2009.2023153.
- [7] S. S. Abdeldayem and T. Bourlai, 'A novel approach for ECG-based human identification using spectral correlation and deep learning,' IEEE Trans. Biom. Behav. Identity Sci. , vol. 2, no. 1, pp. 1-14, Jan. 2020, DOI: 10.1109/TBIOM.2019.2947434.
- [8] B. Fatimah, P. Singh, A. Singhal, and R. B. Pachori, 'Biometric identification from ECG signals using Fourier decomposition and machine learning,' IEEE Trans. Instrum. Meas. , vol. 71, Art. no. 4008209, Aug. 2022, DOI: 10.1109/TIM.2022.3199260.
- [9] Y. Zhao, L. Zhao, Z. Xiao, J. Li, and C. Liu, 'Enhancing electrocardiogram identity recognition using convolutional neural networks with a multisimilarity loss model,' IEEE Trans. Instrum. Meas. , vol. 73, Art. no. 2512008, Feb. 2024, DOI: 10.1109/TIM.2024.3368481.
- [10] F. Han, P. Yang, S. Yan, H. Du, and Y. Feng, 'BreathSign: Transparent and continuous in-ear authentication using bone-conducted breathing biometrics,' in Proc. IEEE Conf. Computer Commun. (INFOCOM) , New York, NY, USA, May 2023, DOI: 10.1109/INFOCOM53939.2023.10229037.
- [11] J. Chauhan, Y. Hu, S. Seneviratne, A. Misra, A. Seneviratne, and Y. Lee, 'BreathPrint: Breathing acoustics-based user authentication,' in Proc. 15th Ann. Int. Conf. Mob. Syst. Appl. Services (MOBISYS) , New York, NY, USA, pp. 278-291, Jun. 2017, DOI: 10.1145/3081333.3081355.
- [12] A. Singh, S. U. Rehman, S. Yongchareon, and P. H. J. Chong, 'Multiresident non-contact vital sign monitoring using radar: A review,' IEEE Sensors J. , vol. 21, no. 4, pp. 4061-4084, Feb. 2021, DOI: 10.1109/JSEN.2020.3036039.
- [13] Z. Fang et al. , 'A review of emerging electromagnetic-acoustic sensing techniques for healthcare monitoring,' IEEE Trans. Biomed. Circuits Syst. , vol. 16, no. 6, pp. 1075-1094, Dec. 2022, DOI: 10.1109/TBCAS.2022.3226290.
- [14] C. Massaroni, A. Nicol` o, M. Sacchetti, and E. Schena, 'Contactless methods for measuring respiratory rate: A review,' IEEE Sensors J. , vol. 21, no. 11, pp. 12821-12839, Jun. 2021, DOI: 10.1109/JSEN.2020.3023486.
- [15] S. Iwata, T. Koda, and T. Sakamoto, 'Multiradar data fusion for respiratory measurement of multiple people,' IEEE Sensors J. , vol. 21, no. 22, pp. 25870-25879, Nov. 2021, DOI: 10.1109/JSEN.2021.3117707.
- [16] X. Liu, J. Cao, S. Tang, J. Wen, and P. Guo, 'Contactless respiration monitoring via off-the-shelf WiFi devices,' IEEE Trans. Mob. Comput. , vol. 15, no. 10, pp. 2466-2479, Oct. 2016, DOI: 10.1109/TMC.2015.2504935.
- [17] Y. S. Lee, P. N. Pathirana, C. L. Steinfort, and T. Caelli, 'Monitoring and analysis of respiratory patterns using microwave Doppler radar,' IEEE J. Transl. Eng. Health Med. , vol. 2, Art. no. 1800912, Oct. 2014, DOI: 10.1109/JTEHM.2014.2365776.
- [18] I. Iwata, T. Sakamoto, T. Matsumoto, and S. Hirata, 'Noncontact measurement of heartbeat of humans and chimpanzees using millimeterwave radar with topology method,' IEEE Sensors Lett. , vol. 7, no. 11, Art. no. 7006104, Nov. 2023, DOI: 10.1109/LSENS.2023.3322287.
- [19] P. Zheng, C. Zheng, X. Li, H. Chen, A. Wang, and Y. Luo, 'Second harmonic weighted reconstruction for non-contact monitoring heart rate,' IEEE Sensors J. , vol. 22, no. 6, pp. 5815-5823, Mar. 2022, DOI: 10.1109/JSEN.2022.3148003.
- [20] M. Le and B. Van Nguyen, 'Multivariate Correlation of Higher Harmonics for Heart Rate Remote Measurement Using UWB Impulse Radar,' IEEE Sensors J. , vol. 20, no. 4, pp. 1859-1866, Feb. 2020, DOI: 10.1109/JSEN.2019.2950635.
- [21] M. Alizadeh, G. Shaker, J. C. M. D. Almeida, P. P. Morita, and S. Safavi-Naeini, 'Remote monitoring of human vital signs using mmwave FMCW radar,' IEEE Access , vol. 7, pp. 54958-54968, Apr. 2019, DOI: 10.1109/ACCESS.2019.2912956.
- [22] T. Sakamoto et al. , 'Feature-based correlation and topological similarity for interbeat interval estimation using ultrawideband radar,' IEEE Trans. Biomed. Eng. , vol. 63, no. 4, pp. 747-757, Apr. 2016, DOI: 10.1109/TBME.2015.2470077.
- [23] I. Iwata, K. Sumi, Y. Tanaka, and T. Sakamoto, 'Accurate radarbased heartbeat measurement using higher harmonic components,' arXiv:2407.07380 , [eess.SP], Jul. 2024. [Online]. Available: http://arxiv. org/abs/2407.07380
- [24] Q. He, E. Yang, and S. Fang, 'A survey on human profile information inference via wireless signals,' IEEE Commun. Surveys Tuts. , early access, 2024, DOI: 10.1109/COMST.2024.3373397.
- [25] G. Paterniani et al. , 'Radar-based monitoring of vital signs: A tutorial overview,' Proc. IEEE , vol. 111, no. 3, pp. 277-317, Mar. 2023, DOI: 10.1109/JPROC.2023.3244362.
- [26] T. Xin, B. Guo, Z. Wang, M. Li, Z. Yu, and X. Zhou, 'FreeSense: Indoor human identification with Wi-Fi signals,' in Proc. IEEE Glob.

Commun. Conf. (GLOBECOM) , Washington, DC, USA, Dec. 2016, DOI: 10.1109/GLOCOM.2016.7841847.

- [27] Y. Zeng, P. H. Pathak, and P. Mohapatra, 'WiWho: WiFi-based person identification in smart spaces,' in Proc. 15th ACM/IEEE Int. Conf. Inf. Process. Sens. Netw. (IPSN) , Apr. 2016, DOI: 10.1109/IPSN.2016.7460727.
- [28] L. Cheng and J. Wang, 'How can I guard my AP? Non-intrusive user identification for mobile devices using WiFi signals,' in Proc. 17th ACM Int. Symp. Mob. Ad Hoc Netw. Comput. (MobiHoc) , New York, NY, USA, Jul. 2016, pp. 91-100, DOI: 10.1145/2942358.2942373.
- [29] B. Korany, H. Cai, and Y. Mostofi, 'Multiple people identification through walls using off-the-shelf WiFi,' IEEE Internet Things J. , vol. 8, no. 8, pp. 6963-6974, Apr. 2021, DOI: 10.1109/JIOT.2020.3037945.
- [30] J. Liu, Y. Chen, Y. Dong, Y. Wang, T. Zhao, and Y.-D. Yao, 'Continuous user verification via respiratory biometrics,' in Proc. IEEE Conf. Computer Commun. (INFOCOM) , Toronto, Canada, Jul. 2020, DOI: 10.1109/INFOCOM41043.2020.9155258.
- [31] P. Cao, W. Xia, M. Ye, J. Zhang, and J. Zhou, 'Radar-ID: human identification based on radar micro-Doppler signatures using deep convolutional neural networks,' IET Radar Sonar Navig. , vol. 12, no. 7, pp. 729-734, Jul. 2018, DOI: 10.1049/iet-rsn.2017.0511.
- [32] B. Vandersmissen et al. , 'Indoor person identification using a low-power FMCW radar,' IEEE Trans. Geosci. Remote Sens. , vol. 56, no. 7, pp. 3941-3952, Jul. 2018, DOI: 10.1109/TGRS.2018.2816812.
- [33] T. Sakamoto, 'Personal identification using a convolutional neural network and ultra-wideband radar measurement of walking and sitting motion,' IEICE Trans. Electron. , vol. J103-C, no. 7, pp. 321-330, Jul. 2020. [Online]. Available: https://arxiv.org/abs/2008.02182.
- [34] K. Saho, K. Inuzuka, and K. Shioiri, 'Person identification based on micro-Doppler signatures of sit-to-stand and stand-to-sit movements using a convolutional neural network,' IEEE Sensors Lett. , vol. 4, no. 3, Art. no. 3500304, Mar. 2020, DOI: 10.1109/LSENS.2020.2975219.
- [35] S. Dong, W. Xia, Y. Li, Q. Zhang, and D. Tu, 'Radar-based human identification using deep neural network for long-term stability,' IET Radar Sonar Navig. , vol. 14, no. 10, pp. 1521-1527, Oct. 2020, DOI: 10.1049/iet-rsn.2019.0618.
- [36] Y. He, H. Guo, X. Zhang, R. Li, Y. Lang, and Y. Yang, 'Person identification based on fine-grained micro-Doppler signatures and UWB radar,' IEEE Sens. J. , vol. 23, no. 18, pp. 21421-21432, Sep. 2023, DOI: 10.1109/JSEN.2023.3299558.
- [37] K. Saho, K. Shioiri, and K. Inuzuka, 'Accurate person identification based on combined sit-to-stand and stand-to-sit movements measured using Doppler radars,' IEEE Sens. J. , vol. 21, no. 4, pp. 4563-4570, Feb. 2021, DOI: 10.1109/JSEN.2020.3032960.
- [38] Z. Ni and B. Huang, 'Gait-based person identification and intruder detection using mm-wave sensing in multi-person scenario,' IEEE Sensors J. , vol. 22, no. 10, pp. 9713-9723, May 2022, DOI: 10.1109/JSEN.2022.3165207.
- [39] J. Pegoraro, F. Meneghello, and M. Rossi, 'Multiperson continuous tracking and identification from mm-wave micro-Doppler signatures,' IEEE Trans. Geosci. Remote Sens. , vol. 59, no. 4, pp. 2994-3009, Apr. 2021, DOI: 10.1109/TGRS.2020.3019915.
- [40] Y. Yang, X. Yang, T. Sakamoto, F. Fioranelli, B. Li, and Y. Lang, 'Unsupervised domain adaptation for disguised-gait-based person identification on micro-Doppler signatures,' IEEE Trans. Circuits Syst. Video Technol. , vol. 32, no. 9, pp. 6448-6460, Sep. 2022, DOI: 10.1109/TCSVT.2022.3161515.
- [41] Y. Lang, Q. Wang, Y. Yang, C. Hou, H. Liu, and Y. He, 'Joint motion classification and person identification via multitask learning for smart homes,' IEEE Internet Things J. , vol. 6, no. 6, pp. 9596-9605, Feb. 2019, DOI: 10.1109/JIOT.2019.2929833.
- [42] Y. Cheng and Y. Liu, 'Person reidentification based on automotive radar point clouds,' IEEE Trans. Geosci. Remote Sens. , vol. 60, Art. no. 5101913, 2022, DOI: 10.1109/TGRS.2021.3073664.
- [43] S. M. M. Islam et al. , 'Cross-modality continuous user authentication and device pairing with respiratory patterns,' IEEE Internet Things J. , vol. 10, no. 16, pp. 14197-14211, Aug. 2023, DOI: 10.1109/JIOT.2023.3275099.
- [44] Y. Wang, T. Gu, and H. Zhang, 'Simultaneous authentication of multiple users using a single mmwave radar,' IEEE Internet Things J. , vol. 11, no. 10, pp. 17797-17811, May 2024, DOI: 10.1109/JIOT.2024.3358548.
- [45] H. Zhao, Y. Ma, and X. Wang, 'MN-UIV: Multimodal neural network enabling user identity verification based on millimeter wave radar,' IEEE Internet Things J. , early access, 2024, DOI: 10.1109/JIOT.2024.3395287.
- [46] A. Rahman, V. M. Lubecke, O. Boric-Lubecke, J. H. Prins, and T. Sakamoto, 'Doppler radar techniques for accurate respiration char-
- acterization and subject identification,' IEEE J. Emerging Sel. Top. Circuits Syst. , vol. 8, no. 2, pp. 350-359, Jun. 2018, DOI: 10.1109/JETCAS.2018.2818181.
- [47] Y. Wang, T. Gu, T. H. Luan, and Y. Yu, 'Your breath doesn't lie: Multiuser authentication by sensing respiration using mmwave radar,' in Proc. 19th Ann. IEEE Int. Conf. Sens. Commun. Netw. (SECON) , Sep. 2022, pp. 64-72. DOI: 10.1109/SECON55815.2022.9918606.
- [48] T.-H. Hwang, J. Shi, and K. Lee, 'Enhancing privacy-preserving personal identification through federated learning with multimodal vital signs data,' IEEE Access , vol. 11, pp. 12155-121566, Jan. 2023, DOI: 10.1109/ACCESS.2023.3328641.
- [49] J. Peng, Z. Hao, Z. Zhang, R. Wang, M. Li, and X. Dang, 'Utter innocence: Contactless authentication method based on physiological signals,' IEEE Sensors J. , vol. 24, no. 6, pp. 7740-7751, Mar. 2024, DOI: 10.1109/JSEN.2024.3353256.
- [50] Y. Wang, T. Gu, T. H. Luan, M. Lyu, and Y. Li, 'HeartPrint: Exploring a heartbeat-based multiuser authentication with single mmwave radar,' IEEE Internet Things J. , vol. 9, no. 24, pp. 25324-25336, Feb. 2022, DOI: 10.1109/JIOT.2022.3196143.
- [51] F. Lin, C. Song, Y. Zhuang, W. Xu, C. Li, and K. Ren, 'Cardiac scan: A non-contact and continuous heart-based user authentication system,' in Proc. 23rd Ann. Int. Conf. Mob. Comput. Netw. (MobiCom) . New York, NY, USA, Oct. 2017, pp. 315-328. DOI: 10.1145/3117811.3117839.
- [52] D. W. Richter, 'Generation and maintenance of the respiratory rhythm,' J. Exp. Biol. , vol. 100, no. 1, pp. 93-107, Oct. 1982, DOI: 10.1242/jeb.100.1.93.
- [53] H. Kobayashi, Y. Tanaka, and T. Sakamoto, 'Individual identification using radar-measured respiratory and heartbeat features,' arXiv:2408.00972 [eess.SP], Aug. 2024. [Online]. Available: https:// arxiv.org/abs/2408.00972
- [54] I. Kerslake and F. Kelly, 'Uses of capnography in the critical care unit,' BJA Education , vol. 17, no. 5, pp. 178-183, May 2017, DOI: 10.1093/bjaed/mkw062.
- [55] C.-H. Hsieh, Y.-H. Shen, Y.-F. Chiu, T.-S. Chu, and Y.-H. Huang, 'Human respiratory feature extraction on an UWB radar signal processing platform,' in Proc. IEEE Int. Symp. Circuits Syst. (ISCAS) , Beijing, China, May, 2013, pp. 1079-1082, DOI: 10.1109/ISCAS.2013.6572037.
- [56] D. O'Shaughnessy, 'Acoustic analysis for automatic speech recognition,' Proc. IEEE , vol. 101, no. 5, pp. 1038-1053, May 2013, DOI: 10.1109/JPROC.2013.2251592.
- [57] M. Jankiraman, FMCW Radar Design , Norwood, MA, USA: Artech House, 2018, pp. 45-48.
- [58] R. A. Monzingo, R. L. Haupt, and T. W. Miller, Introduction to Adaptive Arrays , 2 nd Ed., Rijeka, Croatia: SciTech, 2011, DOI: 10.1049/SBEW046E.
- [59] L. Van der Maaten and G. Hinton, 'Visualizing data using t-SNE,' J. Mach. Learn. Res. , vol. 9, no. 11, pp. 2579-2605, 2008.
- [60] Y. Xue and M. Hauskrecht, 'Active learning of multi-class classification models from ordered class sets,' in Proc. AAAI Conf. Art. Intelligence , vol. 33, no. 1, Art. no. 1, Jul. 2019, DOI: 10.1609/aaai.v33i01.33015589.
- [61] S. Schellenberger et al. , 'A dataset of clinically recorded radar vital signs with synchronised reference sensor signals,' Sci. Data , vol. 7, Art. no. 291, Sep. 2020, DOI: 10.1038/s41597-020-00629-5.
- [62] Z. Xiao, J. Xu, and T. Hu, 'Research on the transmissivity of some clothing materials at millimeter-wave band,' in Proc. Int. Conf. Microw. Mil. Wave Technol. , Nanjing, China, Apr. 2008, pp. 1750-1753, DOI: 10.1109/ICMMT.2008.4540812.
- [63] A. Brizzi and Y. Hao, 'Experimental investigation of the effect of clothes on the propagation of millimeter waves over the human trunk,' in Proc. IEEE Int. Wireless Symp. (IWS) , Beijing, China, Apr. 2013, DOI: 10.1109/IEEE-IWS.2013.6616739.
- [64] N. Alexander, C. Callejero, and R. Gonzalo, 'Millimetre-wave material properties,' in Proc. Joint Int. Conf. Infrared Mil. Waves and Int. Conf. Terahertz Electron. , Cardiff, UK, Sep. 2007, pp. 264-266, DOI: 10.1109/ICIMW.2007.4516490.
- [65] G. W. Kirschen, D. D. Singer, H. C. Thode Jr, and A. J. Singer, 'Relationship between body temperature and heart rate in adults and children: A local and national study,' Am. J. Emerg. Med. , vol. 38, no. 5, pp. 929-933, May 2020, DOI: 10.1016/j.ajem.2019.158355.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-e79913aa-1ceb-485a-ae7e-180acf42412d.jpg)

Haruto Kobayashi (Graduate Student Member, IEEE) received the B.E. degree in electrical and electronic engineering from Kyoto University, Kyoto, Japan, in 2023, where he is currently pursuing the M.E. degree in electrical engineering at the Graduate School of Engineering.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-3d819562-fe76-435e-9afb-c99e4622cffa.jpg)

Yuji Tanaka (Member, IEEE) received the B.E., M.E., and Ph.D. degrees in engineering from Kanazawa University, Ishikawa, Japan, in 2015, 2017, and 2023, respectively. From 2017 to 2020, he worked at the Information Technology R&amp;D center, Mitsubishi Electric Corporation. He is currently an assistant professor at the Nagoya Institute of Technology. His research interests include radar signal processing, radio science, and radar measurement of physiological signals. He was awarded the Second Prize in the URSI-JRSM 2022 Student Paper

Competition and the Young Researcher's Award from the Institute of Electrical and Electronics Engineers (IEICE) Technical Committee on Electronics Simulation Technology. He is a member of the IEICE and the Society of Geomagnetism and Earth, Planetary and Space Sciences (SGEPSS).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00972v2-figure-4f4a444d-40fc-4482-a228-844907074c17.jpg)

Takuya Sakamoto (Senior Member, IEEE) received the B.E. degree in electrical and electronic engineering from Kyoto University, Kyoto, Japan, in 2000 and M.I. and Ph.D. degrees in communications and computer engineering from the Graduate School of Informatics, Kyoto University, in 2002 and 2005, respectively. From 2006 through 2015, he was an Assistant Professor at the Graduate School of Informatics, Kyoto University. From 2011 through 2013, he was also a Visiting Researcher at Delft University of Technology, Delft, the Netherlands.

From 2015 until 2019, he was an Associate Professor at the Graduate School of Engineering, University of Hyogo, Himeji, Japan. In 2017, he was also a Visiting Scholar at the University of Hawaii at Manoa, Honolulu, HI, USA. From 2019 until 2022, he was an Associate Professor at the Graduate School of Engineering, Kyoto University. From 2018 through 2022, he was also a PRESTO researcher of the Japan Science and Technology Agency, Japan. Since 2022, he has been a Professor at the Graduate School of Engineering, Kyoto University. His current research interests are wireless human sensing, radar signal processing, and radar measurement of physiological signals.

Prof. Sakamoto was a recipient of the Best Paper Award from the International Symposium on Antennas and Propagation (ISAP) in 2004, the Young Researcher's Award from the Institute of Electronics, Information and Communication Engineers of Japan (IEICE) in 2007, the Best Presentation Award from the Institute of Electrical Engineers of Japan in 2007, the Best Paper Award from the ISAP in 2012, the Achievement Award from the IEICE Communications Society in 2015, 2018, and 2023, the Achievement Award from the IEICE Electronics Society in 2019, the Masao Horiba Award in 2016, the Best Presentation Award from the IEICE Technical Committee on Electronics Simulation Technology in 2022, the Telecom System Technology Award from the Telecommunications Advancement Foundation in 2022, and the Best Paper Award from the IEICE Communication Society in 2007 and 2023.
Haruto Kobayashi Yuji Tanaka Takuya Sakamoto
Individual Identification Using Radar-Measured Respiratory and Heartbeat Features
This study proposes a method for radar-based identification of individuals
using a combination of their respiratory and heartbeat features. In the
proposed method, the target individual's respiratory features are extracted
using the modified raised-cosine-waveform model and their heartbeat features
are extracted using the mel-frequency cepstral analysis technique. To identify
a suitable combination of features and a classifier, we compare the
performances of nine methods based on various combinations of three feature
vectors with three classifiers. The accuracy of the proposed method in
performing individual identification is evaluated using a 79-GHz
millimeter-wave radar system with an antenna array in two experimental
scenarios and we demonstrate the importance of use of the combination of the
respiratory and heartbeat features in achieving accurate identification of
individuals. The proposed method achieves accuracy of 96.33% when applied to a
five-day dataset of six participants and 99.39% when applied to a public
one-day dataset of thirty participants.