Research Article Sparse GLONASS Signal Acquisition Based on Compressive Sensing and Multiple Measurement Vectors

. A sparse global navigation satellite system (GLONASS) signal acquisition method based on compressive sensing and multiple measurement vectors is proposed. The nonsparse GLONASS signal can be represented sparsely on our proposed dictionary which is designed based on the signal feature. Then, the GLONASS signal is sensed by a normalized orthogonal random matrix and acquired by the improved multiple measurement vectors acquisition algorithm. There are 10 cycles of pseudorandom codes in a navigation message, and these 10 pseudorandom codes have the same row sparse structure. So, the acquisition probability can be raised by row sparse features theoretically. A large number of simulated GLONASS signal experiments show that the acquisition probability increases with the increase in the measurement vector column dimension. Finally, the practical availability of the new method is veriﬁed by acquisition experiments with the real record GLONASS signal. The new method can reduce the storage space and energy loss of data transmission. We hope that the new method can be applied to ﬁeld receivers that need to record and transmit navigation data for a long time.


Introduction
Global navigation satellite system (GLONASS), which was developed by the Soviet Union first, is maintained and developed by Russia now. e purpose of GLONASS is to provide unlimited number of air, marine, and any other type of users with all-weather three-dimensional positioning, velocity measuring, and timing anywhere in the world or near-earth space [1]. Until the writing date of this article (25/ 4/2020), the number of GLONASS satellites in orbit has reached 27. 24 satellites are working normally, such as launching navigation data and time, 1 satellite is in the flight test phase, and 2 satellites are spares. 24 satellites are evenly distributed on three near-circular orbital planes, the three orbital planes are separated by 120 degrees, and the satellites in the same plane are separated by 45 degrees. e orbit altitude is about 19100 kilometers, and the running cycle is 11 hours 15 minutes 44 seconds [2]. Due to the geographical position of Russia, the inclination angle of the GLONASS orbital plane is 64.8 degrees, which is larger than that of the American global position system (GPS), which is 55 degrees, so the coverage effect is better than GPS at high latitudes [3,4].
For the GLONASS receiver, it must acquire the satellite signal in the line of sight firstly. e GLONASS signal is a kind of a spread spectrum signal, so the spread spectrum code phase is unknown. Due to the relative motion of the satellite and the receiver, the carrier frequency is changed, called Doppler frequency. e aim of acquisition is to estimate the coarse spread spectrum code phase and Doppler frequency of the satellite in the field of view [5]. e acquisition is a two-dimensional searching process, which needs to search all the likely spread spectrum code phase and Doppler frequency at the same time. e serial acquisition method is a common method, which needs to reconstruct all possible code phase and Doppler frequency and correlate with the satellite signal. e parallel acquisition method which uses Fourier transform can improve the acquisition effect, but it increases the complexity [6]. e traditional acquisition methods are based on the Nyquist sampling theory, which is a redundant sampling method [7]. e recently proposed compressive sensing (CS) theory can greatly reduce the sampling frequency of sparse signals and reconstruct the sparse signal from the single measurement vector (SMV) through some reconstruction algorithms [8,9]. Now, CS has been applied in many fields including medical image processing [10] and narrowband spectrum sensing [11]. In recent years, several CS-based methods have been improved for global navigation satellite system (GNSS) signal acquisition. In [12], the author proposed a two-stage deterministic GNSS signal acquisition technique using the Walsh-Hadamard matrix. e proposed method makes fast acquisition possible, but the method needs two stages which may be a little complicated. In [13], the authors proposed an efficient acquisition technique for the GPS receiver. e new method sampled and detected the GPS signal using a bank of randomized correlators with much fewer chip-matched filters than the traditional acquisition method. e GPS signals can be acquired by the proposed technique, but it increased the hardware complexity of the receiver. In [14], the authors proposed a GPS signal acquisition method based on CS. e matching pursuit algorithm is used to recover the GPS sparse signal. Numerical simulations illustrate that the method can recover the GPS sparse spike signal effectively. e application of compressive sensing in GNSS is limited to GPS with SMV, and until now, the application in GLONASS with MMV is not reported.
In radar signal processing, it needs to estimate the direction of arrival (DOA), direction finding, direction of departure (DOD), and so on [15,16]. Multi-input multioutput radar (MIMO) uses the correlation between multiple signals and the sparsity of the radar signal to improve the performance of estimation and obtain better performance than the traditional radar [17,18]. A natural extension of CS is the multiple measurement vectors (MMVs) where the sparse signal has the same sparse structure and is measured by the same measurement matrix [19]. When the underlying signals have the common support set, it will help to reduce the successful recovery condition. In recent years, the MMVs have been used for off-the-grid line spectrum denoising and estimation [20], heart rate monitoring [21], classification [22], etc. At present, our research group is working on modeling and simulation of the GLONASS satellite signal. In the project, we found that it is a heavy burden for storage and transmission of the satellite signal at high sampling rate, and the every millisecond signal has the same sparse feature. So, our team studies to apply the theory of compressive sensing and MMV for acquisition. In this paper, a new GLONASS signal acquisition method based on CS and MMV is proposed. We design the sparse representation dictionary according to the same pseudorandom ranging code and the different satellite carrier frequencies. An improved MMV acquisition algorithm is used to acquire the GLONASS signal. Experiment results show that the new method can acquire the GLONASS signal efficiently. e rest of this paper is organized as follows. In Section 2, the CS and MMV are introduced briefly. Section 3 introduces the GLONASS signal and the serial acquisition method. Section 4 describes the new GLONASS signal sparse dictionary and gives the improved acquisition method based on CS and MMV. Section 5 gives the simulation results of the new method. Finally, the paper is concluded in Section 6.

Compressive Sensing and Multiple Measurement Vectors
According to the Nyquist sampling theory, the sampling rate of the signal must be more than twice the maximum frequency of the signal which will get a large number of redundant data. e CS theory is proposed by Donoho et al. [8,9], which points out that the sparse signal can be sensed by a lower sampling rate, which is much lower than that of the Nyquist theory. e sparse signal can be recovered accurately from the measurement data of CS. e CS consists of three parts: sparse representation of the signal, measurement matrix, and signal reconstruction [23].
In CS, the sparse signal x ∈ R N×1 can be measured by the measurement matrix Φ ∈ R M×N (M < N), which can be expressed as where y ∈ R M×1 is the CS measurement data, Ψ is the sparse dictionary, and n ∈ R M×1 is the interference noise. e sparse signal can be reconstructed from the measurement data through some reconstruction algorithms. e sparse signal reconstruction problem can be described as follows: where ‖ · ‖ 0 means zero-norm which counts the number of nonzero elements. We can get the best sparse solution by solving the optimization argument. e zero-norm provides a simple notion of sparsity, but is not a right notion of norm in mathematics because it does not satisfy the norm definition by simple verification. e sparse signal reconstruction is an NP-hard problem and has strict condition. Naturally, some methods are proposed to promote the recovery effect, for example, MMV where a single matrix Φ ∈ R M×N is used to obtain the measurements of multiple sparse signals X � [x 1 , x 2 , . . . , x L ] ∈ R N×L with the same support set. e measurements of MMV can be expressed as y 1 � Φx 1 , y 2 � Φx 2 , . . . , y L � Φx L . Given Y ∈ R M×L and Φ ∈ R M×N , the sparse signal reconstruction based on MMV is an optimization problem [24]: where | · | row−0 counts the nonzero row of matrix X and ‖ · ‖ F is the Frobenius norm. To solve the optimization problem (3), some SMV recovery algorithms have extended to MMV. Orthogonal matching pursuit for MMV (OMPMMV) [25], simultaneous orthogonal matching pursuit (SOMP) [26], MMV focal underdetermined system solver (MFOCUSS) [27], etc., are proposed. When the sparse signals have a common sparse feature, the MMV can lead to a reconstruction advantage. e same support set provides more information to the optimization, so the reconstruction result based on MMV is better than SMV. In [20], authors proposed a novel fast sparse direction-of-arrival (DOA) estimation algorithm based on the MMV and a smooth 0-norm algorithm for the multiple-input multiple-output radar. For the MMV-based joint smooth function, the DOA estimation can be achieved by using the corresponding gradient-based sparse signal reconstruction. e proposed method can solve the MMV problem and achieve better DOA estimation performance for both white and colored Gaussian noises. In [28], authors proposed a new sparse synthetic aperture radar (SAR) imaging method using the MMV model to reduce the computation cost and enhance the imaging result. e sparse row structure information of MMV is used to improve the imaging quality under sub-Nyquist sampling rate.
e experimental results suggest that the proposed method can obtain SAR imaging more effectively than SMV.

GLONASS Signal and Acquisition
GLONASS satellite signal consists of navigation message, pseudorandom (PR) ranging code, and carrier [29]. Navigation message contains system time, ephemeris, satellite health, ionosphere delay parameters, etc. e navigation message is transmitted by the satellite with a rate of 50 bps (20 milliseconds per data). e rate becomes 100 bps after Meander coding. In order to improve the anti-interference ability of the system, Modulo-2 addition with the PR ranging code [30] is performed on the navigation message. PR ranging code is a sequence of the maximum length shift register (m sequence) which is sampled at the output of the 7 th stage of the 9 th stage shift register. All the GLONASS satellites utilize the same PR ranging code and are distinguished by the different carrier frequencies. e sequence length is 511, and the bit rate is 511 kilobits per second, so it is repeated every 1 millisecond. e generating polynomial of the 9th stage shift register is After Modulo-2 addition with the PR ranging code, the navigation message is modulated to different frequencies by binary phase shift keying. e GLONASS is a frequency division multiple access (FDMA) system, and satellites are distinguished according to different carrier frequencies. e navigation signals transmitted by the GLONASS satellite have two carrier frequencies, L1 and L2: where k is the frequency number of the satellite. e range of k is −7 ∼ 6, which can provide 14 different carrier frequencies. A receiver on Earth cannot receive the signals coming from satellites in opposite points of an orbit plane (antipodal), so the two satellites can transmit navigation radio signals on equal frequencies [31]. at is to say, the 14 frequencies can support the usage of 28 satellites. Multipath propagation phenomenon is inevitable in the wireless communication system because of the radio signal reflection or refraction [32,33]. e GLONASS terminal receives the satellite signal with two or more paths including the line-of-sight component. Multipath disturbance is an important error source in the navigation system, which can result in position error because of the biased ranging measurements. In this paper, we give the GLONASS signal model which takes the multipath into account; however, how to mitigate the multipath is out of scope of this paper.
In the receiver, the kth GLONASS intermediate frequency signal on the L1 band can be expressed as (L1 band is discussed in this paper only) where R is the multipath number, a i is the multipath signal amplitude, C(t) is the GLONASS satellite PR ranging code, and D k (t) is the navigation message. After down converting, Serial acquisition algorithm is an often-used spread spectrum signal acquisition method for GLONASS or GPS [34,35]. e method is based on multiplication of the incoming signal and the locally generated reference signal, which is the PR ranging code with phase delay and Doppler frequency shift carrier signals. e PR ranging code has a certain code phase, from 0 to 510 chips. e incoming signal is multiplied by this locally generated PR ranging code. After the multiplication with the PR ranging code, the incoming signal is multiplied by the locally generated carrier signal. Because the phase of the incoming signal is unknown, the incoming signal is necessary to be multiplied with the in-phase (I) and quadrature (Q) signal. e two products are integrated over the period of the PR ranging code, finally squared, and added. e output exceeds the predefined threshold if the PR ranging code delay and Doppler frequency are corresponding to the signal. e receiver concludes that there is the GLONASS satellite in the line of sight. e serial search acquisition method has two different sweeps: a PR ranging code phase sweep over all 511 different code phases and a frequency sweep over all possible carrier frequencies of ±10 kHz with the searching steps of 500 Hz for the high dynamic receiver. In this case, the total number of inner product operation for one satellite is 511 × 2 10000 500 which is a very large number of combinations for the whole navigation system. is exhausting search routine is the main weakness of the acquisition method.

GLONASS Signal Acquisition Based on CS and MMV
A new GLONASS signal acquisition method based on CS and MMV has been proposed. A sparse dictionary corresponding to the satellite is designed according to the satellite PR ranging code and intermediate frequency. e satellite signal has been measured by the same random matrix first, which is compressed largely. en, a MMV reconstruction algorithm is used for the acquisition of the satellite. e new acquisition method is mainly composed of two parts: sparse dictionary and reconstruction algorithm.

Sparse Dictionary of GLONASS Satellite
Signal. e GLONASS signal can be simplified as in navigation data ignoring multipath and noise. e premise of compressive sensing application is sparsity, so the signal needs to be represented sparsely [36]. e signal is nonsparse in common transform basis, for example, discrete cosine transform (DCT), fast Fourier transform (FFT), and discrete wavelet transform (DWT). A sparse dictionary based on the signal feature must be designed. C(t) is the PR ranging code which is an orthogonal sequence, and the cosine signal is sparse based on Fourier transform. For the PR ranging code, which is a sparse signal. All the GLONASS satellites use the same PR ranging code, so the sparse dictionary of the kth satellite signal can be expressed as where C t−τ is the PR ranging code phase delay, e − j(2π(f I +k×0.5625+Δf d )t+θ k ) is the frequency with the step of Δf d , f I + k × 0.5625 stands for the kth satellite signal IF, and ⊗ is the Kronecker product. e position of the sparse value is corresponding to the PR ranging code phase and Doppler frequency. e PR ranging code length is 511, so the satellite signal possibly has 511 different code phases. For a high dynamic receiver, the Doppler frequency range is ±10 kHz with the searching step of 500 Hz, that is to say, there are 41 frequencies. e Kronecker product of code phase delay and frequency sweep is 511 × 41 � 20951. If the length of the GLONASS signal is N, the dictionary dimension is 20951 × N. According to the kth satellite sparse dictionary, the kth GLONASS signal x k can be sparsely represented as where s k is a sparse signal if the kth GLONASS satellite is in the field of view. e histogram of the k � 2 GLONASS satellite under different transform bases is shown in Figure 1. e receiver intermediate frequency is 5.5 MHz, so the satellite signal frequency is 5.5 + 2 × 0.5625 � 6.6250 MHz. e signal-tonoise ratio (SNR) is −20 dB, and the noise is additive white Gaussian noise which is added using the MATLAB function "awgn." e simulation signal has two paths, and multipath signal amplitude is half the line-of-sight signal. PR ranging code phase and Doppler frequency are set randomly. e GLONASS PR ranging code period is 1 ms, so the selected experimental signal size is 1 ms with 40 MHz sampling frequency. In the experiment, the "haar" wavelet is chosen as the analysis function. e sparse representation results are the same based on other wavelet functions. From the histogram, we can find that the signal is nonsparse and cannot be represented sparsely under DCT, FFT, and DWT dictionary. Most of the values are below normalized 0.1 based on the designed sparse dictionary of this paper. e GLONASS signal can be represented sparsely because all the values below 0.1 are considered as noise which can be set to 0.

GLONASS Signal Acquisition Based on CS and MMV.
Sparsity is the essential prerequisite for CS. e GLONASS signal is a sparse signal based on the sparse dictionary, so the CS can be used for compressing and reconstructing the signal. At the satellite, the signal x k transmitted by the satellite is a clean signal which is not polluted by the noise. But, at the receiver, the incoming signal is x k + n which is contaminated by the noise and multipath interference. According to CS theory, the polluted GLONASS signal can be compressed by a random matrix Φ which has carried out orthogonalization between the arbitrary column vectors first, and then a compressed signal y k is got: Duration of one navigation message is 10 ms after Meander coding, but the period of PR ranging code is 1 ms.
at means there are 10 repetitions of the PR ranging code in one navigation message. If the satellite signals are sensed by every millisecond, then MMV can be expressed as a matrix: where X � [x 1 , x 2 , . . . , x L ] is the signal matrix and N � [n 1 , n 2 , . . . , n L ] is the noise matrix. e sparse signal can be obtained from formula (3) according to the MMV reconstruction algorithm. But, the GLONASS signal acquisition is different from the sparse signal reconstruction problem. e purpose of CS signal reconstruction is to reconstruct the sparse signal value and position accurately. e GLONASS signal acquisition is interested in the position of the largest sparse value only. e corresponding PR ranging code phase and Doppler frequency can be obtained from the position. Secondly, the satellite sparse dictionary dimension is large, so some classical reconstruction algorithms which have inverse operation are unfit for GLONASS signal acquisition. In this paper, the greedy algorithm is used for GLONASSS signal acquisition. e greedy algorithm is a successful sparse signal reconstruction method which carries out exhaustive iteration for reducing the residual error by updating the optimal dictionary atom [37]. At iteration, one atom (one column of the sensing matrix) is chosen to add to the support set once for reducing the residual error maximally. e algorithm terminates when the iteration numbers exceed the given value or the error falls below the given threshold. Many variants are available in the family of the greedy algorithm, which hope to offer improvements either in complexity or/and in the accuracy. e orthogonal matching pursuit (OMP) can get good reconstruction effect, which is more famous than other greedy algorithms [38]. e OMPMMV [25] is a successful algorithm for MMV sparse signal reconstruction. e algorithm needs some improvements for MMV GLONASS signal acquisition [14]. e PR ranging code is an orthogonal sequence which can simplify the updating of residual error, and the error is not suitable as a termination condition because the GLO-NASS signal has strong noise. e more important reason is that the GLONASS signal acquisition is more interested in the position of the sparse value than the sparse value. Taking these into consideration, an improved MMV greedy acquisition algorithm for GLONASS signal acquisition is proposed.
e new algorithm has some similarities with the greedy algorithm. e inner product of the residual error with the sensing matrix is identical to the greedy algorithm in iteration step 1. But, there are some differences. Firstly, the PR ranging code is an orthogonal sequence, so the maximum inner product value is selected for updating the residual error which requires a least square calculation in the OMP algorithm. Secondly, the residual error is not suitable as the stopping condition of iteration because the GLONASS signal SNR can deteriorate into −20 dB or worse. In the greedy acquisition algorithm, the judgment condition is the total iteration I. A specific sparse dictionary is designed for one satellite once, so the sparsity is 1 for a dictionary in theory. To enhance the result, the total number of iterations I is set to 5 − 10. Finally, the purpose of acquisition is to estimate the coarse PR ranging code phase and Doppler frequency of the satellite signal but not reconstruct the sparse signal. e final output of the acquisition algorithm is acquiring the satellite or not (Algorithm 1).

Numerical Simulation Experiments
e implementation of the new acquisition algorithm is shown in Figure 2.
e incoming GLONASS signal is polluted by noise. e GLONASS signal is compressed by the matrix Φ ∈ R M×N which is a normalized orthogonal random matrix. e MMV is obtained by rearranging the compressed signal. In the receiver, the satellite dictionary is generated for acquisition. e PR ranging code module generates the spreading spectrum code for all satellites. e local oscillation module generates the carrier.
ere is a 90°phase-shifted version because the phase of the GLONASS signal is unknown [39]. e different GLONASS satellites have different intermediate frequencies. e dictionary module needs to adjust the intermediate frequency of the local oscillator according to the satellite number. Simulation and analysis of the new method with MATLAB is given in the following sections.

MMV Acquisition at Different Noise and Compressive
Ratio. Noise and compressive ratio are the main reasons which affect the acquisition effect of MMV. In this experiment, the influences of noise and compressive ratio at different MMV dimensions are analyzed. For the convenience of discussion, we define the compressive ratio (CR) as where N and M are the length of the GLONASS signal and the compressed data, respectively. e compression ratio intuitively describes the ratio of the compressed signal to the original signal. e third GLONASS satellite (which is selected at random) on the L1 band is selected for numerical simulations. e IF is 5.5 MHz after down-conversion, so the IF of the third satellite is 5.5 + 3 × 0.5625 � 7.1875 MHz. e sampling frequency is set to 40 MHz. e time to transmit a complete PR ranging code needs 1 ms, so the length of the 1 ms GLONASS signal is 40000. e incoming signal has two paths, and the amplitude of the multipath signal is set to half the line-of-sight component. To reduce the simulation burden without loss of generality, we assume that the statistical model for the delay and Doppler Mathematical Problems in Engineering 5 frequency is uniformly distributed. e PR ranging code delay range is 0-255 codes; and the Doppler frequency range is f I ± 5 kHz. In order to improve the acquisition efficiency, the acquisition dictionary of the PR ranging code phase and frequency are set to half chip and 500 Hz, respectively. e Kronecker product of code phase and frequency is 256 × 2 × 21 � 10752, which is the abscissa point. e noise is additive white Gaussian noise (AWGN) and added by the MATLAB function "awgn." e SNR is set to −15 dB, −20 dB, and −25 dB, respectively, and the CR is set to 0.6, 0.4, and 0.2. e combinations of SNR and CR have 9 cases at the same MMV dimension. Iteration time I is set to 5. e aim of this experiment is to analyze the influence of noise and compression ratio at different MMV signal dimensions L � 1, 5, and 10. e MMV acquisition degenerates into an SMV problem when L � 1. To save the page, the MMV acquisition results at −20 dB and −25 dB are shown in

Acquisition Probability Analysis
. At the MMV model, GLONASS signal acquisition based on CS is mainly affected by noise, compression ratio, and MMV signal dimension.
is experiment analyzes the effect of these three factors on the acquisition performance. e k � −1 GLONASS satellite (which is selected at random) on the L1 band is selected for numerical simulations. e center frequency is 5.5 MHz after down-conversion, so the IF of the k � −1 satellite is 5.5 − 1 × 0.5625 � 4.9375 MHz. e sampling frequency is set to 30 MHz. e amplitude of the multipath signal is set to half the direct signal, and the PR ranging phase and Doppler frequency are uniformly distributed. e range of noise is from −40 dB to −10 dB with different MMV signal Input: sparse dictionary Ψ, measurement data Y � ΦX + ΦN � ΦΨS + ΦN � ΘS + ΦN, total iteration times I, threshold coefficient c; Initial: residual error R 0 � Y, sparse solution S 0 � 0, support set C 0 � ∅, and k � 0; Iteration: (1) Perform the inner product operation of residual error and the sensing matrix Θ according to is the j th column of P k ; (2) Find the maximum inner product value p k j m � arg max‖p k (j)‖ 2 , add the index of the maximum into the support set C k � C k− 1 ∪ j m , and add the maximum inner product value into the corresponding position of sparse signal e experimental vector dimension is set to L � 1, 4, 7, and 10. From the figures, we can see that the acquisition probabilities are low when the MMV dimension is 1, and the probability increases greatly when the MMV dimension increases to 4. e acquisition probabilities are close to the classical serial acquisition method when the dimension becomes 10. e experiments conclude that the MMV provides more relevant information between the rows and improves the performance of acquisition.

Acquisition of Real GLONASS Signal.
e real environment is very complicated, the terminal receives the signals of multiple satellites at the same time, which may interfere with each other, and the noise is not only additive white Gaussian noise. In this experiment, we select a real GLONASS signal for testing the MMV algorithm, and the signal is collected on July 6, 2017, by the UTREK310 receiver in the Nanjing University of Aeronautics and Astronautics. e UTREK310 receiver is a GLONASS and GPS dualchannel signal collecting system with its matching software receiver SDRG310. e receiver radio frequency (RF) front end filters and amplifies the signal. en, the GLONASS signal is saved at the computer after the down-converter and the analog-to-digital converter.
e IF of the collected GLONASS signal is 4.5 MHz, and the sampling frequency is 20 MHz. All the satellites in the line of sight are shown in Figure 5, which is arranged from large to small according to the intensity of the serial acquisition correlation value. at is to say, the k � −4 satellite signal intensity is the weakest. e satellites' constellation and acquisition results are shown in Figures 5(a) and 5(b).
In this experiment, the 10 ms GLONASS signal is used for simulation, and GPS signals are removed to reduce acquisition interference. e ±5 kHz searching range of Doppler frequency is enough because the receiver is static. e PR ranging code phase is 0-510 codes. e searching step is 500 Hz and half code, respectively, so the combination of code phase and frequency sweep is 511 × 2 × 21 � 21462. According to the acquisition results of experiments 1 and 2, it is appropriate to set the CR to 0.4 for high compressive ratio. All the six satellites in the line of sight can be acquired successfully. Only the acquisition results of the k � −4 satellite with L � 1, 5, and 10 are shown in Figures 6(a), 6(b), and 6(c) because the satellite signal is the weakest, and the acquisition results of other five satellites are better. In the MMV acquisition algorithm, the judgment threshold c is set to 2 for the six satellites, but this value is not suitable for the weakest k � −4 satellites when the MMV signal dimension L � 1 because the second maximum normalized value is greater than 0.5.
For the receiver, acquisition is a blind searching process. e receiver does not know which satellites are in the line of sight. at is to say, the acquisition algorithm needs to search all the possible satellites. e new algorithm needs to establish all satellites' sparse dictionary in turn. ere are no −7, −6, etc., satellites in the recorded signal. Taking the k � 1 satellite, for example, acquisition results of the recorded signal based on the sparse dictionary of the k � 1 satellite with L � 1, 5, and 10 are shown in Figures 7(a), 7(b), and 7(c), respectively. We can find that these values are all very close to each other. We can conclude that there is no k � 1

Conclusions
In this paper, we proposed a new sparse GLONASS signal acquisition method based on CS and MMV. A new sparse dictionary is designed for representing the GLONASS signal, and an improved sparse MMV reconstruction algorithm is used for acquisition. A large number of simulation results show that the new method can effectively acquire simulated or real GLONASS signals at high compression ratio and low SNR circumstances with high MMV column dimension. e new method can greatly reduce the sampling rate and energy loss in signal processing and has important application prospects, such as low altitude aircraft with less noise and interference and position and monitoring of field animals which need to record the GLONASS data for a long time.
e acquisition performances of the new algorithm are not ideal at lower SNR circumstances, which may be improved by applying reference [40] methods which are   successful in low SNR. Our research team will continue to delve into how to improve the acquisition effect of this method.

Data Availability
All experimental data can be downloaded from the network disk: "https://pan.baidu.com/s/ 1v3uNVXDFa_s31vu0XMfsPQ," with keyword: "1234." (1) e [ Figure 1.mat] data used to support the acquisition result of this study have been shown in Figure 1. (2) e [ Figure 3.rar] data used to support the acquisition results of this study have been shown in Figures 3(a)

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.