CN114338297B - Combined timing synchronization and frequency offset estimation method under incoherent LoRa system - Google Patents

Combined timing synchronization and frequency offset estimation method under incoherent LoRa system Download PDF

Info

Publication number
CN114338297B
CN114338297B CN202111419634.7A CN202111419634A CN114338297B CN 114338297 B CN114338297 B CN 114338297B CN 202111419634 A CN202111419634 A CN 202111419634A CN 114338297 B CN114338297 B CN 114338297B
Authority
CN
China
Prior art keywords
signal
frame structure
data
chirp
data frame
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111419634.7A
Other languages
Chinese (zh)
Other versions
CN114338297A (en
Inventor
余忠洋
高继勋
李博
何栎
王禹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Henan Institute of Engineering
Original Assignee
Henan Institute of Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Henan Institute of Engineering filed Critical Henan Institute of Engineering
Priority to CN202111419634.7A priority Critical patent/CN114338297B/en
Publication of CN114338297A publication Critical patent/CN114338297A/en
Application granted granted Critical
Publication of CN114338297B publication Critical patent/CN114338297B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Radio Relay Systems (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The invention provides a combined timing synchronization and frequency offset estimation method under an incoherent LoRa system, which comprises the steps of establishing a data frame structure, sequentially carrying out de-modulation operation, de-chirp operation, discrete Fourier transform operation, modulus taking operation, maximum taking operation, amplitude taking operation and the like on a received pilot signal to obtain a combined offset estimation quantity related to time delay and Doppler frequency shift, and compensating the combined offset estimation quantity into the received data signal, thereby realizing reliable estimation of large transmission time delay and Doppler frequency shift under extremely low signal-to-noise ratio. The invention supports reliable satellite Internet of things communication under extremely low signal-to-noise ratio, large transmission delay and Doppler frequency shift, and has lower processing complexity and implementation complexity compared with the existing estimation algorithm.

Description

Combined timing synchronization and frequency offset estimation method under incoherent LoRa system
Technical Field
The invention relates to the technical field of communication of the Internet of things, in particular to a combined timing synchronization and frequency offset estimation method under an incoherent LoRa system.
Background
For satellite internet of things transmission, its transmission link generally has the following three disadvantages:
(1) The received signal power is low: in free space propagation, the amplitude fading of the received signal is proportional to the square of the transmission distance. The satellite terminals are relatively far from the user or ground control station, which can result in the signal power received at the receiving end becoming very low. This requires that the designed receiver synchronization scheme work well at very low signal to noise ratios;
(2) The Doppler shift is large: taking carrier frequency of 4GHz (S band), short distance of satellite to user or ground control station of 200km and relative speed of Mach 10 (Mach 1 is calculated according to 340 m/S) as examples, according to calculation formula f of Doppler frequency shift d =f c vcos θ/c (where f c The maximum obtainable doppler shift (i.e., θ=0) is about 45kHz for the carrier frequency, v is the relative velocity, θ is the angle with the ground horizontal. This will limit the design of the receiver synchronization scheme both in terms of estimation accuracy and estimation range;
(3) The timing and remaining carrier deviation is large: the low-rail satellite ground transmission delay is 30ms, which can lead to the closed-loop control period being far longer than that of ground mobile communication, so that the performance of timing synchronization and frequency offset estimation adopting closed-loop control is reduced. And therefore a larger timing error and residual frequency offset remain during closed loop control.
In fact, internet of things transmissions, including satellite internet of things, belong to short packet transmissions, so the available pilot resources become very limited. And the larger doppler shift and propagation delay can lead to serious challenges for a communication system based on the LoRa technology. This requires a study of the receiver synchronization scheme with respect to the LoRa technique. For studies of receiver synchronization schemes under LoRa technology, tapprel J et al, in the article "An open-source LoRa physical layer prototype on GNU radio" (2020 IEEE SPAWC,2020:1-5), describes a software defined radio implementation of LoRa transceivers, while also designing An estimation module for carrier frequency offset and sampling time offset. Colavolpe G et al, in "Reception of LoRa signals from LEO satellites" (IEEE Transactions on Aerospace and Electronic Systems,2019,55 (6): 3587-3602), herein, describe a joint parameter estimation scheme based on a fast Fourier transform (Fast Fourier transformation, FFT) for carrier parameters such as Doppler shift, doppler rate and transmission delay. Bernier C et al, in the text, "Low complexity LoRa frame synchronization for ultra-low power software-defined radios" (IEEE Transactions on Communications,2020,68 (5): 3140-3152), propose a low complexity LoRa frame synchronization scheme intended for use in the recently proposed ultra-low power consumption software defined radio systems. Xhonneux M et al, in the text, "A low-complexity synchronization scheme for LoRa end nodes" (https:// arxiv. Org/abs/1912.11344v1, 2019), states that a low complexity synchronization scheme for LoRa end nodes is proposed that can estimate and correct carrier offset versus sampling time. However, these synchronization schemes require a relatively complex process and do not take into account the characteristics of the incoherent LoRa system.
Disclosure of Invention
Aiming at the defects of low received signal power, large Doppler frequency shift, large deviation between timing and residual carrier and the like in satellite Internet of things communication and the technical problem that the processing process of a synchronous scheme based on the LoRa technology in the prior art is complex, the invention provides a combined timing synchronization and frequency offset estimation method under an incoherent LoRa system, which supports reliable satellite Internet of things communication with extremely low signal-to-noise ratio, large transmission delay and Doppler frequency shift and has lower processing complexity and implementation complexity.
In order to solve the technical problems, the invention adopts the following technical scheme: a joint timing synchronization and frequency offset estimation method under an incoherent LoRa system comprises the following steps:
step S1: establishing a data frame structure, wherein the data frame structure comprises a received pilot signal and a received data signal;
step S2: performing de-modulation operation on the received pilot signal in the data frame structure to obtain a de-modulated signal;
step S3: carrying out de-chirp operation on the de-modulated signal to obtain a de-chirp signal;
step S4: performing discrete Fourier transform on the de-chirped signal to obtain a frequency domain signal;
step S5: performing a modulus operation on the frequency domain signal to obtain a corresponding amplitude;
step S6: performing maximum value and amplitude angle operation on the obtained amplitude value to obtain a corresponding maximum value index, and taking the maximum value index as an estimated value of the joint offset of the transmission delay and the Doppler frequency shift;
step S7: and compensating the received data signals in the data frame structure by using the estimated value of the joint offset to obtain corrected data signals, and realizing reliable estimation of large transmission delay and Doppler frequency shift under extremely low signal-to-noise ratio.
The method for establishing the data frame structure in the step S1 is as follows:
step S1.1: given a length L p Pilot block and length L d Is a data block of (a);
step S1.2: will have a length of L p Insertion of pilot block of length L d Number of (2)And obtaining a data frame structure F according to the head of the block.
The method for obtaining the demodulation signal in the step S2 is as follows:
step S2.1: firstly, a received pilot signal in a data frame structure F is obtained, and the data frame structure F is traversed through sampling time k to obtain a sampling time set kappa corresponding to a pilot block p ={k:0≤k≤L p -1 and a set of sample moments k corresponding to a data block d ={k:L p ≤k≤L p +L d -1}; then the sampling time set kappa corresponding to the pilot block is collected p Extracting the first chirp received pilot signal one by one through sampling time k:
Figure SMS_1
wherein: b is the transmission bandwidth, m=2 SF Is orthogonal chirp number, SF is spreading factor, τ, f d And θ are respectively transmission delay, doppler shift and phase offset, n k (l) Is the mean value is 0 and the variance is sigma 2 Complex gaussian random variable s k (l-tau) is the LoRa modulated signal with the transmission delay tau attached,
Figure SMS_2
is an imaginary unit;
step S2.2: finally receives the pilot signal r k (l) p Performing a de-modulation operation to obtain a de-modulated signal:
Figure SMS_3
wherein: d, d k For transmitting pilot symbols.
In the step S3, the de-modulated signal r' k (l) p And performing a de-chirp operation to obtain a de-chirp signal:
Figure SMS_4
in the step S4 of the above-mentioned process,for de-chirp signal z k (l) Performing Discrete Fourier Transform (DFT) operation to obtain a frequency domain signal:
Figure SMS_5
wherein: q represents the frequency index of the DFT.
In the step S5, the frequency domain signal Z (q) is subjected to a modulus operation to obtain a corresponding amplitude value:
Figure SMS_6
in the step S6, the amplitude is calculated
Figure SMS_7
And performing maximum value and amplitude angle taking operation to obtain corresponding maximum value indexes, namely estimated values of the joint offset:
Figure SMS_8
the method for obtaining the correction data signal in the step S7 is as follows:
step S7.1: first, the received data signal in the data frame structure F is obtained, and the sampling time set kappa corresponding to the data block is collected d And extracting the first chirp received data signal one by one through the sampling time k:
Figure SMS_9
step S7.2: then use the estimate of the joint offset
Figure SMS_10
Receiving a data signal r k (l) d Performing compensation operation to obtain a corrected data signal:
Figure SMS_11
the beneficial effects of the invention are as follows:
1. compared with the existing estimation algorithm, the method has lower processing complexity. Because the estimation method of the present invention does not require separate estimation of the transmission delay and the doppler shift, instead, a joint offset is estimated for the transmission delay and the doppler shift.
2. In actual operation, the Discrete Fourier Transform (DFT) operation in the method based on the joint timing synchronization and frequency offset estimation can be replaced by high-efficiency Fast Fourier Transform (FFT) operation, so that the method has lower implementation complexity.
3. The invention finally obtains a joint offset estimation quantity related to delay and Doppler frequency shift through a series of operation processing on the received pilot signal, and then compensates the joint offset estimation quantity into the received data signal, thereby realizing reliable estimation of large transmission delay and Doppler frequency shift under extremely low signal-to-noise ratio.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a block diagram of a data frame according to the present invention;
FIG. 2 is a schematic flow chart of the present invention;
fig. 3 is an uncoded incoherent LoRa system performance based on a joint timing synchronization and frequency offset estimation method at spreading factor sf=12 according to the present invention;
fig. 4 shows the performance of an uncoded incoherent LoRa system of the present invention based on a joint timing synchronization and frequency offset estimation method at a spreading factor sf=14;
fig. 5 shows performance of a Turbo coding incoherent LoRa system based on a joint timing synchronization and frequency offset estimation method under a spreading factor sf=12;
fig. 6 shows performance of a Turbo coding incoherent LoRa system based on a joint timing synchronization and frequency offset estimation method under the spreading factor sf=14.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without any inventive effort, are intended to be within the scope of the invention.
In order to support reliable satellite Internet of things communication under extremely low signal-to-noise ratio, large transmission delay and Doppler frequency shift, the invention provides a joint timing synchronization and frequency shift estimation method under an incoherent LoRa system. The specific implementation steps comprise:
step S1: first a data frame structure is established.
The method for establishing the data frame structure comprises the following steps:
step S1.1: given a length L p Pilot block and length L d Is a data block of (a);
step S1.2: will have a length of L p Insertion of pilot block of length L d The header of the data block of (a) is the data frame structure F shown in fig. 1.
The received pilot signal in the data frame structure F is then processed accordingly to obtain an estimate of the joint offset for the propagation delay and doppler shift. The specific operation process is shown in fig. 2, and comprises the following steps:
step S2: the method comprises the following specific operations of performing a de-modulation operation on a received pilot signal in a data frame structure F to obtain a de-modulated signal:
step S2.1: firstly, a received pilot signal in a data frame structure F is obtained, and the data frame structure F is traversed through sampling time k to obtain a sampling time set kappa corresponding to a pilot block p ={k:0≤k≤L p -1 and a set of sample moments k corresponding to a data block d ={k:L p ≤k≤L p +L d -1}; then the sampling time set kappa corresponding to the pilot block is collected p The receiving pilot signal r of the first chirp is extracted one by one through the sampling time k k (l) p
Figure SMS_12
Wherein: b is the transmission bandwidth, m=2 SF Is orthogonal chirp number, SF is spreading factor, τ, f d And θ are respectively transmission delay, doppler shift and phase offset, d k To transmit pilot symbols, n k (l) Is the mean value is 0 and the variance is sigma 2 Complex gaussian random variable s k (l-tau) is the LoRa modulated signal with the transmission delay tau attached,
Figure SMS_13
is imaginary unit, ++>
Figure SMS_14
Referred to as a joint offset for transmission delay and doppler shift.
Step S2.2: finally receives the pilot signal r k (l) p Performing de-modulation operation to obtain a de-modulated signal r' k (l) p
Figure SMS_15
Wherein:
Figure SMS_16
is gaussian noise.
Step S3: for the obtained de-modulated signal r' k (l) p Performing the de-chirp operation to obtain a de-chirp signal z k (l):
Figure SMS_17
Wherein:
Figure SMS_18
is gaussian noise.
Step S4: for the obtained de-chirp signal z k (l) Performing Discrete Fourier Transform (DFT) operation to obtain a frequency domain signal Z (q):
Figure SMS_19
wherein: delta (·) is the dirac function,
Figure SMS_20
is the noise term n k (l) DFT results of->
Figure SMS_21
Is N 1 (q) random inversion result, +.>
Figure SMS_22
For the joint offset S (τ, f d ) The round (·) function represents the rounding operation of the rounding operation, and q represents the frequency index of the DFT.
Step S5: performing a modulus operation on the obtained frequency domain signal Z (q) to obtain a corresponding amplitude value
Figure SMS_23
Figure SMS_24
Wherein:
Figure SMS_25
is N 2 One rotation result of (q).
Step S6: for the obtained amplitude
Figure SMS_26
Performing maximum value taking and amplitude angle taking operations to obtain corresponding maximum value index q * I.e. estimate of joint offset +.>
Figure SMS_27
Figure SMS_28
Step S7: using estimates of joint offsets
Figure SMS_29
For the received data signal r in the data frame structure F k (l) d Performing compensation operation to obtain corrected data signal +.>
Figure SMS_30
The specific operation is as follows:
step S7.1: first, the received data signal in the data frame structure F is obtained, and the sampling time set kappa corresponding to the data block is collected d The first chirp received data signal r is obtained by extracting one by one through the sampling time k k (l) d
Figure SMS_31
Step S7.2: then use the estimate of the joint offset
Figure SMS_32
Receiving a data signal r k (l) d Performing compensation operation to obtain correctionData signal->
Figure SMS_33
Figure SMS_34
Wherein:
Figure SMS_35
also the mean value is 0 and the variance is sigma 2 Is a complex gaussian random variable of (c).
In this embodiment, the received signal is separated into a data signal and a pilot signal by a demultiplexer, and the received pilot signal is sequentially subjected to a demodulation operation, a chirp removal operation, a DFT operation, a modulus taking operation, a maximum taking operation, an amplitude angle taking operation, and the like to obtain a joint offset estimation amount related to delay and doppler shift, and then the joint offset estimation amount is compensated into the received data signal, so as to realize reliable estimation of large doppler shift and transmission delay under extremely low signal-to-noise ratio.
In order to further illustrate the beneficial effects of the present invention, in this embodiment, a comparison is performed through a simulation experiment, which is specifically as follows:
simulation 1:
1.1 simulation conditions
The transmission bandwidth b=20 MHz, the spreading factor SF is 12 and 14, and the corresponding adjacent chirp intervals Δf are 488Hz and 122Hz, respectively. Assuming again a transmission delay τ=256 chips, maximum doppler shift f d =45 kHz > Δf and random phase bias θ∈ [ -pi, pi). For LoRa modulated signals, the incoherent demodulation method proposed in "Capacity approaching codes for non-coherent orthogonal modulation" (IEEE Transactions on Wireless Communications,2007,6 (11): 4004-4013) by Fabregas AG I and Grant AJ is used.
1.2 simulation results and analysis
With pilot symbol length L p =4 and data symbol length L d And=60 is an example. The pilot overhead at this time is about 0.067. In practice, for spreading factors sf=12 and sf=14, forThe corresponding transmission data sequences are 720 bits and 840 bits in length, respectively.
Fig. 3 and fig. 4 show Bit Error Rate (BER) performance of an uncoded incoherent LoRa system based on a joint timing synchronization and frequency offset estimation method under these two spreading factors, respectively. Wherein the diamond marked curves in fig. 3 are shown in the uncoded LoRa system, sf=12, f d Under the conditions of=45 kHz and τ=256 chirp, the uncoded incoherent LoRa system performance based on the joint timing synchronization and frequency offset estimation method is adopted. The curve marked with square in fig. 3 shows the results in an uncoded LoRa system, sf=12, f d And under the conditions of=0 and τ=0, adopting an uncoded incoherent LoRa system performance based on a joint timing synchronization and frequency offset estimation method. The diamond marked curves in fig. 4 are shown in uncoded LoRa system, sf=14, f d Under the conditions of=45 kHz and τ=256 chirp, the uncoded incoherent LoRa system performance based on the joint timing synchronization and frequency offset estimation method is adopted. The curve marked with square in fig. 4 shows the results in an uncoded LoRa system, sf=14, f d And under the conditions of=0 and τ=0, adopting an uncoded incoherent LoRa system performance based on a joint timing synchronization and frequency offset estimation method.
As can be seen from the simulation results of fig. 3 and 4, the incoherent LoRa modulation system based on the proposed synchronization method achieves near-ideal conditions (i.e., f, no matter whether the spreading factor SF is 12 or 14 d =0 and τ=0). Specifically, when ber=10 -4 When for sf=12 and sf=14, the incoherent LoRa system performance based on the proposed synchronization scheme is 0.4dB and 0.2dB worse than the ideal case, respectively. These results demonstrate the effectiveness of the proposed joint timing synchronization and frequency offset estimation method. In addition, the spreading gain from a large spreading factor can be found to be very substantial.
Simulation 2:
2.1 simulation conditions
Take a Turbo code with a code rate of 0.5 and an information bit length of 384 as an example. The adopted inner and outer interleavers are both quadratic permutation polynomial interleavers, and the adopted decoding algorithm is a modified Max-Log-MAP algorithm; pilot symbol length L used p =4, for sf=12 and sf=14, the transmitted data symbol length L d 768/12=64 and respectively
Figure SMS_36
(/>
Figure SMS_37
The representation is rounded up, which may be achieved by zero padding the codeword), the corresponding pilot overhead is approximately 0.063 and 0.073, respectively.
2.2 simulation results and analysis
Fig. 5 and fig. 6 show the performance of the Turbo coding incoherent LoRa system based on the joint timing synchronization and frequency offset estimation method under the two spreading factors, respectively. Wherein, the curve marked with diamond in fig. 5 represents the non-coherent LoRa system, sf=12, f in Turbo coding d =45kHz、τ=256chirps、L p And under the condition of=4, the performance of the non-coherent LoRa system of Turbo coding based on the joint timing synchronization and frequency offset estimation method is adopted. The curve marked with square in fig. 5 shows the result of coding incoherent LoRa system in Turbo, sf=12, f d And under the conditions of=0 and τ=0, the performance of the Turbo coding incoherent LoRa system based on the joint timing synchronization and frequency offset estimation method is adopted. The curve marked with diamond in fig. 6 shows the result of coding incoherent LoRa system in Turbo, sf=14, f d =45kHz、τ=256chirps、L P And under the condition of=4, the performance of the non-coherent LoRa system of Turbo coding based on the joint timing synchronization and frequency offset estimation method is adopted. The curve marked with square in fig. 6 shows the result of coding incoherent LoRa system in Turbo, sf=14, f d And under the conditions of=0 and τ=0, the performance of the Turbo coding incoherent LoRa system based on the joint timing synchronization and frequency offset estimation method is adopted.
From the simulation results of fig. 5 and fig. 6, it can be found that the Turbo coding incoherent LoRa system based on the proposed synchronization method achieves BER performance very close to the ideal case, regardless of whether the spreading factor SF is 12 or 14. Specifically, when ber=10 -4 ~10 -5 When the range is in the range, the performance difference of the two is not more than 0.1dB.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, alternatives, and improvements that fall within the spirit and scope of the invention.

Claims (4)

1. A combined timing synchronization and frequency offset estimation method under an incoherent LoRa system is characterized by comprising the following steps:
step S1: establishing a data frame structure, wherein the data frame structure comprises a received pilot signal and a received data signal;
step S2: performing de-modulation operation on the received pilot signal in the data frame structure to obtain a de-modulated signal;
the method for obtaining the demodulation signal in the step S2 is as follows:
step S2.1: firstly, a received pilot signal in a data frame structure F is obtained, and the data frame structure F is traversed through sampling time k to obtain a sampling time set kappa corresponding to a pilot block p ={k:0≤k≤L p -1 and a set of sample moments k corresponding to a data block d ={k:L p ≤k≤L p +L d -1},L p For the length of the pilot block, L d Is the length of the data block; then the sampling time set kappa corresponding to the pilot block is collected p Extracting the first chirp received pilot signal one by one through sampling time k:
Figure FDA0004251075160000011
wherein: b is the transmission bandwidth, m=2 SF Is orthogonal chirp number, SF is spreading factor, τ, f d And θ are respectively transmission delay, doppler shift and phase offset, n k (l) Is the mean value is 0 and the variance is sigma 2 Complex gaussian random variable s k (l-tau) is the LoRa modulated signal with the transmission delay tau attached,
Figure FDA0004251075160000012
is an imaginary unit;
step S2.2: finally butt jointReceiving pilot signal r k (l) p Performing a de-modulation operation to obtain a de-modulated signal:
Figure FDA0004251075160000013
wherein: d, d k For transmitting pilot symbols;
step S3: carrying out de-chirp operation on the de-modulated signal to obtain a de-chirp signal;
in the step S3, the de-modulated signal r k ′(l) p And performing a de-chirp operation to obtain a de-chirp signal:
Figure FDA0004251075160000014
step S4: performing discrete Fourier transform on the de-chirped signal to obtain a frequency domain signal;
step S5: performing a modulus operation on the frequency domain signal to obtain a corresponding amplitude;
step S6: performing maximum value and amplitude angle operation on the obtained amplitude value to obtain a corresponding maximum value index, and taking the maximum value index as an estimated value of the joint offset of the transmission delay and the Doppler frequency shift;
in the step S6, the amplitude is calculated
Figure FDA0004251075160000015
And performing maximum value and amplitude angle taking operation to obtain corresponding maximum value indexes, namely estimated values of the joint offset:
Figure FDA0004251075160000021
step S7: performing compensation operation on the received data signals in the data frame structure by using the estimated value of the joint offset to obtain corrected data signals, so as to realize reliable estimation of large transmission delay and Doppler frequency shift under extremely low signal-to-noise ratio;
the method for obtaining the correction data signal in the step S7 is as follows:
step S7.1: first, the received data signal in the data frame structure F is obtained, and the sampling time set kappa corresponding to the data block is collected d And extracting the first chirp received data signal one by one through the sampling time k:
Figure FDA0004251075160000022
step S7.2: then use the estimate of the joint offset
Figure FDA0004251075160000023
Receiving a data signal r k (l) d Performing compensation operation to obtain a corrected data signal:
Figure FDA0004251075160000024
2. the method for jointly timing synchronization and frequency offset estimation in an incoherent LoRa system according to claim 1, wherein the method for establishing a data frame structure in step S1 is as follows:
step S1.1: given a length L p Pilot block and length L d Is a data block of (a);
step S1.2: will have a length of L p Insertion of pilot block of length L d The header of the data block of (2) to obtain a data frame structure F.
3. The method for joint timing synchronization and frequency offset estimation in incoherent LoRa system according to claim 2, wherein in step S4, the de-chirp signal z is de-chirped k (l) Performing discrete Fourier transform operation to obtain a frequency domain signal:
Figure FDA0004251075160000025
wherein: q represents the frequency index of the DFT.
4. The method for joint timing synchronization and frequency offset estimation in incoherent LoRa system according to claim 3, wherein in step S5, the frequency domain signal Z (q) is subjected to a modulo operation to obtain a corresponding amplitude value:
Figure FDA0004251075160000026
CN202111419634.7A 2021-11-26 2021-11-26 Combined timing synchronization and frequency offset estimation method under incoherent LoRa system Active CN114338297B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111419634.7A CN114338297B (en) 2021-11-26 2021-11-26 Combined timing synchronization and frequency offset estimation method under incoherent LoRa system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111419634.7A CN114338297B (en) 2021-11-26 2021-11-26 Combined timing synchronization and frequency offset estimation method under incoherent LoRa system

Publications (2)

Publication Number Publication Date
CN114338297A CN114338297A (en) 2022-04-12
CN114338297B true CN114338297B (en) 2023-06-30

Family

ID=81046477

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111419634.7A Active CN114338297B (en) 2021-11-26 2021-11-26 Combined timing synchronization and frequency offset estimation method under incoherent LoRa system

Country Status (1)

Country Link
CN (1) CN114338297B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115022143A (en) * 2022-05-16 2022-09-06 西安电子科技大学 Carrier synchronization method suitable for LFM modulation system
CN117439849B (en) * 2023-12-19 2024-03-15 深圳市国电科技通信有限公司 Signal demodulation method, device, electronic equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3438693A1 (en) * 2017-08-03 2019-02-06 Vodafone IP Licensing limited A method of signal encoding, a receiver device and associated transmission system
CN109818660A (en) * 2019-01-23 2019-05-28 南京信息工程大学 MIMO-OFDM underwater acoustic system design method based on FRFT
CN111190183A (en) * 2018-11-13 2020-05-22 通用汽车环球科技运作有限责任公司 Sliding window integration scheme for target detection in radar systems

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7221701B2 (en) * 2002-08-28 2007-05-22 Altratek, Inc. System and method for CDMA communications
CN101645865B (en) * 2009-09-15 2011-12-28 哈尔滨工业大学 Multipath time delay and Doppler frequency shift estimation method in signal channel based on fractional Fourier transform and system for implementing method
CN104168232B (en) * 2014-08-11 2018-07-10 厦门大学 The assay method of multidiameter delay and Doppler frequency shift in a kind of underwater acoustic channel
US10107895B2 (en) * 2014-09-19 2018-10-23 The Boeing Company Amplitude calibration of a stepped-chirp signal for a synthetic aperture radar
EP3264622B1 (en) * 2016-07-01 2019-03-27 Semtech Corporation Low complexity, low power and long range radio receiver
CN106597407B (en) * 2016-12-06 2019-05-07 西安电子科技大学 More stellar associations are the same as lower weak echo signal time delay and Doppler frequency shift combined estimation method
KR102120021B1 (en) * 2018-01-22 2020-06-05 라디우스 가부시키가이샤 Reception method, reception device, transmission method, transmission device, transmission and reception system
CN112787730A (en) * 2020-12-25 2021-05-11 青岛三海科技有限公司 Doppler estimation method based on joint waveform design
CN113156365B (en) * 2021-03-12 2022-04-12 复旦大学 Speed, angle and distance joint estimation method based on conjugate ZC sequence pair
CN113472705B (en) * 2021-06-16 2022-06-21 复旦大学 Estimation and prediction method of RIS auxiliary channel based on ZC sequence
CN113567946B (en) * 2021-07-20 2023-08-29 同济大学 Method for detecting real target and virtual image of millimeter wave radar

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3438693A1 (en) * 2017-08-03 2019-02-06 Vodafone IP Licensing limited A method of signal encoding, a receiver device and associated transmission system
CN111190183A (en) * 2018-11-13 2020-05-22 通用汽车环球科技运作有限责任公司 Sliding window integration scheme for target detection in radar systems
CN109818660A (en) * 2019-01-23 2019-05-28 南京信息工程大学 MIMO-OFDM underwater acoustic system design method based on FRFT

Also Published As

Publication number Publication date
CN114338297A (en) 2022-04-12

Similar Documents

Publication Publication Date Title
US7466768B2 (en) IQ imbalance compensation
CN114338297B (en) Combined timing synchronization and frequency offset estimation method under incoherent LoRa system
JP5450645B2 (en) Blind frequency offset estimation for temporally and / or spatially correlated signals
CN101710885A (en) Method for estimating frequency deviation in carrier synchronization of wireless communication system
CN102724150B (en) Method, device and system for realizing IQ (in-phase component and quadrature component) unbalance estimation
US20170373893A1 (en) Frequency shift keying (fsk) demodulator and method therefor
CN102413089B (en) Shannon limit coding GMSK demodulation method used for satellite communication system
CN101018083A (en) Dopla frequency shift compensation method in the MPSK mobile communication system
WO2012088288A1 (en) Automatic frequency offset compensation in zero-intermediate frequency receivers using minimum-shift keying (msk) signaling
US20040171398A1 (en) Channel estimation method for a mobile communication system
US7643571B2 (en) Wireless communication system and method
US20020159531A1 (en) OFDM data demodulators synchronization
US7529306B2 (en) Estimation of asymmetries between inphase and quadrature branches in multicarrier transmission systems
CN111800364B (en) Method for estimating and correcting frequency offset of coded CPM (continuous phase modulation) signal based on waveform matching
CN111901269A (en) Gaussian frequency shift keying modulation method, device and system with variable modulation index
CN113225284B (en) 8FSK incoherent soft decision demodulation method for high-performance channel coding and decoding
US7606331B2 (en) Frequency offset compensation in radio receiver
CN115022143A (en) Carrier synchronization method suitable for LFM modulation system
CN109818894B (en) GMSK signal detection method and detection device under multipath channel
CN112671684B (en) Self-adaptive demodulation method of short-time burst BPSK signal
US20100027727A1 (en) Recovering symbols in a communication receiver
CN114465865A (en) Frequency offset estimation method for pilot frequency segmentation recombination
KR100489409B1 (en) Method for improving a performance of transmitter and receiver in wireless communication systems
US20170180182A1 (en) Joint noncoherent demodulation and carrier frequency offset correction based on non-linear filtering
US20030112061A1 (en) Method for demodulating signal

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant