Randomized resonant metamaterials for single-sensor identification of elastic vibrations

Vibrations carry a wealth of useful physical information in various fields. Identifying the multi-source vibration information generally requires a large number of sensors and complex hardware. Compressive sensing has been shown to be able to bypass the traditional sensing requirements by encoding spatial physical fields, but how to encode vibration information remains unexplored. Here we propose a randomized resonant metamaterial with randomly coupled local resonators for single-sensor compressed identification of elastic vibrations. The disordered effective masses of local resonators lead to highly uncorrelated vibration transmissions, and the spatial vibration information can thus be physically encoded. We demonstrate that the spatial vibration information can be reconstructed via a compressive sensing framework, and this metamaterial can be reconfigured while maintaining desirable performance. This randomized resonant metamaterial presents a new perspective for single-sensor vibration sensing via vibration transmission encoding, and potentially offers an approach to simpler sensing devices for many other physical information. Designing efficient and flexible metamaterial with uncorrelated transmissions for spatial vibration encoding and identification remains a challenge. Here, the authors propose a randomized resonant metamaterial with randomly coupled local resonators for single-sensor identification of elastic vibrations.

V ibrations carry a wealth of physical information in the natural world, which is not only important for perceiving the surrounding environments 1 , but also useful in various fields, such as health care monitoring 2,3 , earthquake detection 4 , smart devices 5 , and Internet of Things 6 . Before being picked up by sensors, vibrations are filtered and mixed during propagation. Identifying the information carried by vibrations has great significance 7,8 . In most complex cases, directly measuring vibrations from multiple sources is impossible. The inverse methods are the main approaches to vibration identification by measuring the related quantities 9 .
One of the inverse methods is based on the filtering property of unintentionally designed vibration transfer paths to quantify which sources and paths contribute the most to the vibration issues 10 .
Other methods, such as blind source separation and array signal processing techniques, can reconstruct vibration information from mixed signals when assuming some source characteristics (e.g., statistical independence and narrowband) 11,12 . In these methods, the accuracy of vibration identification relies on the number and placements of sensors, which requires complex data acquisition systems and control circuits, and results in high power consumption.
Compressive sensing has been demonstrated to be an effective approach to circumvent the traditional sensing requirements 13 . One of the most impressive advances in compressive sensing is the design of the single-pixel camera 14 , which has inspired many studies in both the electromagnetic and acoustic fields [15][16][17][18][19][20] . In these studies, compressive sensing is combined with spatially encoded structures such as random scattering masks and metamaterials. Metamaterials are a broad family of artificially structured materials with unusual effective properties and functionalities [21][22][23][24][25] . Flexible manipulations of electromagnetic, acoustic, and elastic waves can be achieved, such as cloaking [26][27][28] , beaming 29,30 , diffusing 31 , illusion 32 , and hologram 33 . Fascinating applications, such as high-speed analog computing 34,35 , ultrasensitive detection [36][37][38] , and efficient waveguiding 39 have been demonstrated. In the field of acoustics, a properly designed metamaterial achieved encoding of acoustic waves in the spatial and frequency domains, which can be considered as a physical implementation of the measurement matrix in compressive sensing 17,20 . Unlike advances in the field of acoustics, spatial vibration encoding and identification with metamaterials remains unexplored. The core issue is how to design an eligible metamaterial with uncorrelated transmissions to encode spatial vibration information.
In this work, we propose a randomized resonant metamaterial with randomly coupled local resonators for single-sensor identification of elastic vibrations. This metamaterial is designed by developing the theory of randomly coupled resonator dynamics. The metamaterial is proved to be capable of producing highly uncorrelated transmissions for different spatial vibrations due to the disordered coupling of random effective masses. Based on encoding spatial vibration information by uncorrelated transmissions, a compressive sensing framework can be used to experimentally identify various vibration events with only a single sensor. We further demonstrate that the designed metamaterial can be used as a new type of human-machine interface. Our study not only is applicable to areas such as smart devices and Internet of Things, but also provides exciting perspectives for designing simpler vibration sensing devices and other physical sensing systems.

Results
Randomly coupled resonator system for vibration encoding. To achieve the highly uncorrelated transmissions for spatial vibration information encoding, we propose a concept of randomly coupled resonator dynamics and develop the corresponding system. The proposed system consists of multiple different coupling networks. Each coupling network is composed of N local resonators connected by springs k 0 and dampers c 0 (Fig. 1a). k n and c n are the stiffness and the damping coefficients of the nth resonator, respectively. x n and y n are displacements of matrix M and mass m. We derive the randomly coupled resonator dynamics (see Methods for details). In the dynamics, the local resonator is simplified to be an effective mass as below where k d n ¼ k n þ iωc n ; k n ¼ mω 2 n and ω n is the resonant frequency of the nth local resonator 40 . The effective mass is frequency-dependent and can be negative near the resonant frequency due to the anti-resonance effect, which corresponds to the attenuation of elastic waves ( Supplementary Fig. 1).
In the coupling network, the effective masses of all local resonators can be the same or different. For convenience, we refer to the coupling network with the same resonators as an ordinary coupling network (OCN) and that with random resonators as a randomized coupling network (RCN). By using the finite element method (FEM) and the analytical method, the transmissions of an OCN (with a local resonance frequency of 420 Hz) and a RCN can be obtained ( Fig. 1b) (see Methods for details of numerical simulation and analytical derivation). The analytical results agree well with the numerical results. For the OCN, the transmission has a wide attenuation region (gray shaded region in Fig. 1b) due to the presence of the local resonance band gap 41 . If the OCN is used for sensing, the vibration information in the band gap will be lost, which is disadvantageous for vibration encoding and identification. In contrast, the RCN can randomly encode vibration information in the frequency domain, because the randomized effective masses lead to random distributions of the amplification and attenuation bands. This randomized encoding property is useful for constructing a desirable measurement matrix for the compressive sensing. Figure 1c shows an entire randomly coupled resonator system by connecting six RCNs and a center module (without local resonators) for spatial vibration encoding. In each RCN, the spatial distribution is disordered for the effective masses of resonators. Vibrations generated from multiple sources propagate from the boundaries to the center module. The transmission of the coupling system can be obtained by measuring the vibration signals of the center mass M 0 . Figure 1d shows the transmissions of the coupling system when excitations are applied at two different locations (see Supplementary Figure 4a for transmissions corresponding to the other locations). It can be seen that substantial differences exist between the two transmissions. To evaluate the uncorrelation of the transmissions, we calculate the correlation coefficients and the correlation matrix C μ (see Supplementary Note 1). The absolute values of elements in C μ are visualized as shown in Supplementary Fig. 4b. The histogram of the absolute cross-correlation coefficients show that they narrowly distribute around zero ( Supplementary Fig. 4c). We choose the average of the absolute cross-correlation coefficients μ Ave as the metric of the uncorrelation. Here, the μ Ave is approximately equal to 0.15, indicating that the transmissions are desirably uncorrelated.
To explain the relationship between the uncorrelation and the disordered effective masses, we analyze the vibration modes of the entire coupling system. Figure 1e shows the vibration modes obtained by, respectively, applying excitations to six locations at 600 Hz. It can be seen that the displacement fields are obviously different, which reflects the complexity of the elastic vibration propagation in the coupling system. From the randomly coupled resonator dynamics, we can find that this complex vibration propagation is determined by the effective mass matrix and the coupling matrix of the coupling system. These two matrices reflect the disordered coupling relationship of random effective masses. The disordered coupling of the effective masses leads to the high uncorrelation of vibration transmissions, which is the physical basis for effective spatial vibration encoding.
We also investigate the effects of the model parameters on the vibration transmission property of the randomly coupled resonator system (see Supplementary Note 2). The relationship between the μ Ave and the parameters, including c n , c 0 , k 0 , and the number of the local resonators, can be found in Supplementary Figs. 5 and 6. The results show that the stiffness of the matrix k 0 , the damping of the resonators c n , and the number of the resonators have a great influence on the uncorrelation of the vibration transmission, which provides a guidance for the optimization model design.
Design of the randomized resonant metamaterial. The proposed randomly coupled resonator system with disordered effective masses provides the physical basis for designing an actual metamaterial system, which is feasible for uncorrelated vibration encoding. Here, we choose the spiral-based resonators to realize the metamaterial because their physical properties can be easily tuned by changing the structural parameters 30,[42][43][44] . The designed metamaterial consists of six different supercells surrounding a hexagon plate (Fig. 2a). Each supercell is carefully designed by randomly coupled resonators so that the transmissions of the metamaterial system are sensitive to the spatial locations. A single acceleration sensor is attached to the center plate to pick up the mixed modulated vibration signals generated from multiple sources. In this way, the sensing system can be expressed as y = Mx, where y is the vector form of the measured data from the single sensor (i.e., observation vector), x is the object vector containing the information of sources, and M is the measurement matrix determined by the encoding property of the metamaterial and the pre-knowledge of vibration excitations.
The supercell of the metamaterial system is composed of 19 unit cells with random geometric parameters on spiral beams (Fig. 2b). Each unit cell includes two Archimedean spiral beams and a center mass assembly consisting of a bolt and a nut. The parametric equation that determines the spiral beams is given by , and θ is a random variable from 2.6π to 4.6π. The lattice constant A 0 and the height h of the unit cell are 25 mm and 4 mm, respectively. The spiral beams and matrix of the unit cell are made of polylactic acid (PLA) with a Young's modulus E PLA = 3.8 GPa, shear modulus G PLA = 1.3 GPa, and mass density ρ PLA = 1250 kg m −3 . The center mass of the unit cell is 2.52 g. Figure 2c shows the dispersion curves calculated along the Γ-M direction by applying Bloch boundary conditions to the unit cell with θ = 2.952π. We also calculate the frequency-dependent effective mass density 45 . The out-of-plane band gap between 304.0 Hz and 493.7 Hz is consistent with the frequency band of the negative effective mass density. The out-of-plane local resonance modes at the edges of the band gap, denoted as A1 and A2, are responsible for forming the band gap. The starting frequency of the band gap can be tuned by varying the angle θ of the spiral beams, as shown in Fig. 2d. Therefore, the geometric parameter of the unit cell can be derived from the fitting function when the local resonance  Figure 2e shows the distribution of local resonance frequencies of the entire metamaterial system, where a color block corresponds to a unit cell (see Supplementary Table 3 for the parameters of unit cells). This disordered design ensures that the metamaterial system has highly uncorrelated transmissions when vibrations are excited at different locations. The highly uncorrelated transmissions can encode the information of multi-source vibrations to construct the measurement matrix M. Once M and y are determined, the object vector x can be reconstructed by solving the inverse problem using compressive sensing theory to realize vibration identification.
To investigate the encoding performance of the metamaterial system, we conduct experiments on a sample of the metamaterial as shown in Fig. 3a (see Methods). Six speakers generating vibration signals are connected to the metamaterial system at different locations. Spatial transmissions can be obtained by using the experimental modal analysis method (Fig. 3b). The correlation of the transmissions (Fig. 3c) is calculated by using the same method in Supplementary Note 1, and the average cross-correlation coefficient μ Ave is achieved as 0.15. The results demonstrate that the transmission property of the metamaterial system has desirable uncorrelation, which is beneficial for constructing an ideal measurement matrix satisfying compressive sensing theory. Single-sensor vibration identification. To verify the performance of the metamaterial system for vibration identification, we experimentally conduct the identification tasks for multi-source vibration information. Twenty different testing signals with normalized energy are collected in a signal set (see Methods for the construction of the testing signals). Because the identification needs a priori knowledge of the measurement matrix, a calibration process is experimentally performed in advance to construct the measurement matrix M by calculating the spectra of the measured signals (see Methods). After that, the testing signals are simultaneously emitted from two or three randomly selected sources out of six possible locations. This testing process is performed 40 times. The signals played each time are randomly selected from the signal set. The data measured by the single acceleration sensor form the observation vector y. A preprocessing is conducted to compress the dimensionality of the measurement matrix while maintaining the reconstruction accuracy and further improving the robustness against noises by using principal component analysis (see Methods). Then, the object vector x, which contains the information of vibration source locations and signal indices, can be reconstructed by using the two-step iterative shrinkage/thresholding (TWIST) algorithm 46 . Specifically, we take an example to demonstrate the entire identification process in detail in Supplementary Note 3 and Supplementary Fig. 8. Figure 4a shows a reconstruction result with three activated sources, where the blocks in the darkest color represent the most likely locations and signal indices. We denote this configuration of the metamaterial by Config. 1. The reconstruction results agree well with the truth (see Supplementary Fig. 9 for more results). To succinctly characterize the total reconstruction results of the 40 experiments, the reconstructed x is divided into six segments. The  Fig. 10).
Here, we define it as a correct identification when the highest reconstructed strength of the activated source is above 0.5 (whereas it is below 0.25 for an unactuated source) and matches the actual locations and signal indices. For this task, the recognition ratio is 96.7%. The results demonstrate that the metamaterial sensing system can accurately identify multi-source vibration information from the single-sensor measurement, which is of great significance for vibration sensing.
In addition to the Config. 1 in Fig. 4a, the metamaterial system can be reconfigured into other configurations according to the practical needs. A windmill-shaped metamaterial system (denoted by Config. 2) is assembled as shown in Fig. 4b. Vibrations are excited at different positions of the windmill blades. A similar recognition task is performed, and the reconstruction result is shown in Fig. 4b. The total recognition ratio of the 40 experiments is 97.9% (Supplementary Fig. 11). Figure 4c shows the configuration (Config. 3) and the reconstruction result for the case where the metamaterial is embedded in a large acrylic plate. Unlike the above configurations, excitations are not directly applied to the metamaterial. The metamaterial system still has a desirable performance, with a recognition ratio of 97.1% (Supplementary Fig. 12a). Furthermore, we investigate the noise robustness of this metamaterial sensing system. Testing signals are excited at two locations, whereas Gaussian noises (SNR = 0 dB) are excited at another location (Fig. 4d). The reconstructed result shows that the vibration information can still be correctly identified under noise interferences. The recognition ratio of the entire recognition task is 95.0% (Supplementary Fig. 12b). The above results demonstrate that the reconfigurable property of the metamaterial extends the application ranges of the metamaterial system (e.g., in quadrotor drones and airplane wings) while ensuring the accuracy of vibration identification. The proposed metamaterial system can also be used to identify impacts as shown in Fig. 4e. Different from the testing signals in the signal set, impacts are sparse in both time and spatial domains. The identification of impacts relies only on the intrinsic property of the metamaterial system. This ensures that various impact events can be identified once the impulse responses of the metamaterial system are calibrated (see Supplementary Note 4 for construction of the measurement matrix for impact identification). When impact events occur, the TWIST algorithm can be used to reconstruct the impact signals from the measured signals. Figure 4f shows the reconstruction results of an impact event with two consecutive impulses applied to Location 1. The occurrence location and time of the reconstructed impact are in good agreement with the actual ones. Figure 4g shows the successful identification of another impact event, with four impulses applied to Location 5. The results above verify the feasibility of the proposed metamaterial system for impact identification.
To demonstrate the broad application prospects of our approach, we use the proposed metamaterial system to track the trajectory of multiple vibration events. Here, eight probes (denoted by "A" to "H") are fixed on a plate where the metamaterial is embedded as shown in Fig. 5a. We tap the probes with the finger in sequence along the illustrated trajectory "A-E-F-B-H-D-C-G-A", and the vibration signals are measured as shown in Fig. 5a by the single sensor. We can clearly see the tracked trajectory according to the identification of vibration events (see the details in Supplementary Note 5). Figure 5b shows  Supplementary Fig. 13 for an example). This trajectory tracking process is dynamically visualized in Supplementary Movie 1. We also study the tracking process of the trajectories "SJTU" (i.e., the abbreviation of Shanghai Jiao Tong University). The tracked trajectories are well exhibited in Fig. 5c (see Supplementary Fig. 14 for the details). Figure 5d shows the successfully tracked trajectory "Vase" with 12 probes. The experimental details and reconstruction results are presented in Supplementary Fig. 15 and Supplementary Movie 2. The results above demonstrate the time-dependent space coding ability of the proposed metamaterial system, which can create a new type of human-machine interface for instruction, communication, and encryption without complex hardware and high power consumption. Moreover, it also has potential application prospects in fields such as robot tactile sensing and collision tracking.

Discussion
Our proposed design provides an effective strategy for the spatial encoding and single-sensor identification of elastic vibrations, which can not only reduce the complexity of traditional sensing approaches, but also flexibly meet the practical needs. The proposed randomly coupled resonator system with disordered effective masses provides the theoretical basis for designing the randomized resonant metamaterial and is easy to be extended.
The theoretical concept is a dynamics method that incorporates physical mechanism. The physical mechanism of spatial vibration encoding is that the disordered coupling of effective masses leads to uncorrelated vibration transmissions. In addition, the vibration transmission property of the metamaterial is the synergy of the resonance and anti-resonance of the local resonators. Vibration signals can thus be enhanced to some extent in certain frequency bands, which is superior to the band-stop filtering mechanism 17 and benefit the vibration identification. The designed metamaterial system achieves the identification of multiple vibration events with high recognition ratios, and has proved to be robust to noise interference. Furthermore, we have demonstrated that the designed metamaterial system can be used for trajectory tracking, which has potential application prospects in various fields.
Although the proposed metamaterial system is a prototype, it provides an attractive avenue for simpler vibration sensing, and it can be further improved in the following aspects. The metamaterial can be assembled into more configurations to identify potential spatial vibrations according to practical needs, e.g., a 3D structure configuration 47 to identify omnidirectional elastic vibrations. The operational frequency ranges can be adjusted by changing the structural parameters of unit cells. The size of the system can be further reduced with the help of advanced manufacturing techniques and microelectronics, so that it can be flexibly integrated into numerous smart devices. In addition, the field theory can be introduced to improve the analytical model for further optimization design and performance prediction of the device 28 However, there are also limitations for the proposed approach. First, although the proposed metamaterial system has achieved an effective spatial vibration encoding, the randomized design strategy is not the optimal one that can maximize the uncorrelation while minimizing the number of unit cells. This limitation can be overcome by using topological optimization 34,49 or machine learning methods 50 to achieve the on-demand design of the disordered structures in the further work. Moreover, the proposed approach needs the pre-experimental calibrations because the actual sample in experiments might not match the designed one due to the unavoidable manufacturing tolerance. The vibration identification is partly dependent on some preknowledge of vibration excitations, which is another reason for pre-experimental calibrations. This issue is hoped to be solved by utilizing artificial intelligence algorithms based on the computational and experimental data 51,52 .
Overall, we have demonstrated a randomized resonant metamaterial with disordered effective masses for single-sensor vibration identification. This work opens up avenues of vibration transmission encoding with metamaterials. We envision that the proposed metamaterial can integrate with numerous intelligent devices, platforms, and structures (e.g., wearable devices, quadrotor drones, and airplane wings) so that it can be used in broad fields such as human-machine interaction, health care monitoring, industrial field detection, information processing and communication. Moreover, we believe that the proposed metamaterial design strategy can provide the basis for designing various types of simper sensing devices for vibration and many other physical information.

Methods
Numerical simulations. Numerical simulations are conducted to analyze the randomly coupled resonator system by using the lumped mechanical system interface in the multibody dynamics module of COMSOL Multiphysics Version 5.4. Without loss of generality, the parameters of the randomly coupled resonator system are set as M = 2×10 −3 kg, m = 3×10 −3 kg, k 0 = 1.9×10 5 N m −1 , c 0 = 0.01 N s m −1 , c n = 0.25 N s m −1 , and N = 19. The local resonance frequency f n that determines k n can be found in Supplementary Table 1. The node configuration of the resonator is provided in Supplementary Fig. 2. Here, nodes "1", "2", and "3" are the internal nodes connecting the masses, springs, and dampers of the resonator. Node "a" is the input node, and nodes "b1", "b2", and "b3" are the output nodes of the subsystem. The coupling network is constructed by connecting multiple subsystems as shown in Supplementary Fig. 3a. The node connection is provided in Supplementary Table 2. Then, a unit excitation is applied to the input terminal of the coupling network. The transmission can be obtained by picking up the signal from the output terminal. The node connection of the entire coupling system is similar to that of the coupling network, and the boundary conditions are shown in Supplementary Fig. 3b. By changing the input nodes, vibration transmissions and modes of the randomly coupled resonator system can be obtained.
Derivation of the randomly coupled resonator dynamics. The dynamical equation of the nth resonator can be expressed as ÞÀk n y n À x n ð Þ¼F n ; m€ y n þ c n _ y n À _ x n ð Þþk n y n À x n ð Þ¼0; where ψ ni ¼ À1; i couples with n; ϕ n ; i ¼ n; 0; others; ϕ n is equal to the number of the resonators coupling with the nth resonator. Equation (2) can be derived as The expression of the effective mass (Eq. (1)) can thus be obtained. From Eq. (4), we can express the dynamical equation of the kth coupling network as a matrix-vector form of where M eff k ¼ diagðm eff 1 ; m eff 2 ; …; m eff N Þ is a diagonal matrix composed of effective masses given by Eq. (1), and N is the number of the resonators. The displacement vector of the resonators has the form of X k = [x 1 , x 2 , …, x N ] k T = A k e iωt , where A k is the amplitude vector. The stiffness-damping matrix K c k can be expressed as where Ψ k with elements ψ ni is a N × N coupling matrix reflecting the coupling relationship of the resonators (see Source Data file for the details). The timeharmonic excitation applied to the resonators is Thus, the displacement amplitude vector A k can be calculated by Therefore, the transmission of the coupling network can be expressed as where |x out | and |x in | are the elements in A k that correspond to the output and the input of the coupling network. The transmission of the entire coupling system can also be derived in a similar way. The dynamical equation of the coupling system is given by Here, M eff is the mass matrix of the center module of the system. X E = A E e iωt and F E = F E0 e iωt are the displacement vector and the excitation vector of the entire system. K c E = (k 0 + iωc 0 )Ψ E , where Ψ E is the coupling matrix (see Source Data file for the details). Equation (9) can thus be written as In the following analysis, let Φ E = K c E À ω 2 M eff E . Due to the fixed boundary conditions (some elements in A E are zero), the corresponding rows and columns in Φ E , A E , and F E0 must be deleted. The modified matrix and vectors are denoted bŷ Φ E ,Â E , andF E0 . Thus, the displacement amplitude vector of the entire system can be given byÂ The transmission of the entire coupling system can be obtained by calculating the ratio of the displacement amplitudes for the output and the input.
The transmissions of the coupling network and the entire coupling system obtained by the above analytical method are shown in Figs. 1b and 1d of the main text, which agrees well with the numerical results.
Experimental setup. The experimental setup to investigate the metamaterial system is shown in Fig. 3a. The main body of the metamaterial is fabricated with PLA by using fused-filament-fabrication 3D printing. The hexagon socket-head cap bolts (M4 × 8) and nuts in the metamaterial are made of 304 stainless steel. The metamaterial is fixed on the holders that are glued to the dusk caps of speakers. The speakers can be driven by a power amplifier. An impulse force hammer (Kistler, 9724 A) is used to obtain the transmissions of the metamaterial system. The vibration signals are picked up by an acceleration sensor (B&K, 4371). The measured data is acquired by a data acquisition device (NI, 9234). The transmissions of the metamaterial system are calculated as |A out (ω)|/|F in (ω)|, where A out (ω) and F in (ω) are the fast Fourier transform amplitudes of the signals from the acceleration sensor and the hammer, respectively. Then, we use a laser vibrometer (Polytec, PSV-500-3D-H) to obtain the vibration modes of the metamaterial.
For the identification tasks, we construct 20 testing signals as a signal set to verify the performance of the metamaterial system. Each signal is constructed by superposing randomized sine waves as s(t) = Σ 20 i¼1 a i sin(2πf i t + φ i ), where t ∈ [0, 1], a i ∈ [0, 1], f i ∈ [200, 800], and φ i ∈ [0, 2π]. Energy normalization is performed by dividing the root mean square of each signal to ensure that the testing signals have the same intensity. The waveforms and spectra of the 20 testing signals are provided in Supplementary Fig. 7. A laptop is used to generate the testing signals. Then, the signals are sent to an integrated acoustic card (M-Audio, M-Track Eight) and amplified by a power amplifier to drive the speakers.
Construction of the measurement matrix in frequency domain. The measurement matrix contains the spatial vibration encoding information of the metamaterial and the content of the testing signals. The construction of the measurement matrix includes two steps: experimental calibration and dimensionality compression.
Experimental calibration is currently an effective way to obtain the accurate transmissions due to the unavoidable manufacturing tolerance and the complex testing environments. Theoretically, the measurement matrix M p×q can be expressed by M q (ω p ) = H k (ω p )·S j (ω p ), where p is the index of frequency, k is the index of the location, j is the index of the testing signal, q = k × j is the column number of measurement matrix, H k (ω p ) is the transmission of the metamaterial system, and S j (ω p ) is the spectrum of the testing signal. In our experiments, the measurement matrix is calibrated by successively playing the testing signals from different locations and calculating the spectra of the measured signals. More details can be found in Supplementary Note 3.
The current calibrated measurement matrix M p×q contains redundant information as well as the noise. After the mixed vibration signal y p×1 (i.e., observation vector) is acquired, we respectively replicates the matrices y p×1 and M p×q by row for r times (r > 1 and r ∈ Z) to increase the weight of the effective information before using dimensionality compression algorithms. The principal component analysis is used to compress the dimensionality of the measurement matrix. First, M rp×q and y rp×1 are assembled into matrix A rp×(q+1) . Then, we perform zero-mean normalization on A rp×(q+1) . Next, singular value decomposition is performed on AA T to obtain the left singular vector matrix U rp×rp and the singular value matrix D. By selecting the first p′ principal singular values of D and transforming U rp×rp to U rp×p′ , the compressed measurement matrix can be obtained by M p′×q = U rp×p′ T M rp×q . In this way, the principal component analysis can compress the measurement matrix to reduce the computational complexity while maintaining the reconstruction accuracy and further improving the robustness against noises.

Code availability
Source code and processed data are available from the corresponding author upon reasonable request.