CN113810326A - Method and device for estimating time offset, electronic equipment and storage medium - Google Patents

Method and device for estimating time offset, electronic equipment and storage medium Download PDF

Info

Publication number
CN113810326A
CN113810326A CN202111138548.9A CN202111138548A CN113810326A CN 113810326 A CN113810326 A CN 113810326A CN 202111138548 A CN202111138548 A CN 202111138548A CN 113810326 A CN113810326 A CN 113810326A
Authority
CN
China
Prior art keywords
value
channel
time offset
time
deviation
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.)
Granted
Application number
CN202111138548.9A
Other languages
Chinese (zh)
Other versions
CN113810326B (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.)
New H3C Technologies Co Ltd
Original Assignee
New H3C Technologies Co Ltd
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 New H3C Technologies Co Ltd filed Critical New H3C Technologies Co Ltd
Priority to CN202111138548.9A priority Critical patent/CN113810326B/en
Publication of CN113810326A publication Critical patent/CN113810326A/en
Application granted granted Critical
Publication of CN113810326B publication Critical patent/CN113810326B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2668—Details of algorithms
    • H04L27/2669—Details of algorithms characterised by the domain of operation
    • H04L27/2671—Time domain
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2656—Frame synchronisation, e.g. packet synchronisation, time division duplex [TDD] switching point detection or subframe synchronisation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W56/00—Synchronisation arrangements
    • H04W56/004—Synchronisation arrangements compensating for timing error of reception due to propagation delay
    • H04W56/005—Synchronisation arrangements compensating for timing error of reception due to propagation delay compensating for timing error by adjustment in the receiver
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The application discloses a time offset estimation method and device. A method for estimating time offset includes: obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot; obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel; when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS; determining a target time offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals. By applying the technical scheme, the time offset value is calibrated through the deviation value, the accuracy of the time offset estimation value is improved, and N is not increasedifftNamely, the accuracy of the time offset estimation is improved without increasing the amount of calculation.

Description

Method and device for estimating time offset, electronic equipment and storage medium
Technical Field
The present invention relates to the field of communications technologies, and in particular, to a method and an apparatus for estimating a time offset, an electronic device, and a storage medium.
Background
In wireless communication, a certain time offset is generated by transmission of signals in space and transmission of signals in an interface. In order to obtain accurate sampling results, the requirement for time-offset synchronization is high. Although a part of the time offset is eliminated in the random access process, a certain time offset error still remains. Therefore, at the receiving end, further time offset synchronization is required.
Sounding Reference Signals (SRS) are used to assess the channel quality of the uplink and measure related parameters, which are then reported to higher layersAnd flexible scheduling is realized. In practical application, since the SRS is a frequency domain signal, Inverse Fast Fourier Transform (IFFT) is required before time domain time offset estimation. Due to the number of transform points NifftThere is a slight deviation between the actual peak position of the SRS and the measured peak position. Which in turn causes a bias in the strategy for time-offset synchronization. Theoretically NifftThe smaller the deviation of the time offset synchronization is; in order to ensure the accuracy of time-offset synchronization, N is increasedifft. But N isifftThis increase in the number of operations also increases.
Disclosure of Invention
In order to solve the above technical problems, embodiments of the present application provide a time offset estimation method, device, electronic device, and storage medium, so as to improve the estimation accuracy of a time offset value without increasing the required computational power. The technical scheme is as follows:
a method of time offset estimation, comprising:
obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot;
obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel;
when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS;
determining a target time offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
A time offset estimation apparatus, comprising:
a time offset estimation unit for obtaining a time offset value obtained when performing time offset estimation based on a channel Sounding Reference Signal (SRS) in a first time slot;
a receiving unit, configured to obtain a channel parameter of an uplink channel in a second time slot, where the channel parameter is used to characterize channel quality of the uplink channel;
a deviation determining unit configured to determine a deviation value corresponding to the time deviation value according to a time domain channel function when performing time deviation estimation based on the SRS, when the channel quality of the uplink channel satisfies a preset condition;
the synchronization unit determines a target offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
An electronic device, comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor implements a method of time offset estimation by executing the executable instructions;
the time offset estimation comprises the following steps:
obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot;
obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel;
when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS;
determining a target time offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
A computer readable storage medium having stored thereon computer instructions which, when executed by a processor, implement a method of time bias estimation;
the time offset estimation comprises the following steps:
obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot;
obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel;
when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS;
determining a target time offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
According to the technical scheme provided by the embodiment of the application, a time offset value obtained by estimating the time offset of the SRS is obtained in a first time slot; obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel; when the channel quality of an uplink channel meets a preset condition, determining a deviation value corresponding to the time deviation value according to a time domain channel function obtained by SRS time deviation estimation; summing the time offset value and the deviation value to obtain a target time offset value; for time-offset synchronization of the signals. The time offset value is calibrated by using the offset value, and N is not increased while the accuracy of the time offset estimation value is improvedifftNamely, the accuracy of the time offset estimation is improved without increasing the amount of calculation.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and other drawings can be obtained by those skilled in the art according to the drawings.
FIG. 1 is a flow chart illustrating a method of time offset estimation according to an embodiment of the present application;
fig. 2 is a schematic structural diagram of a time offset estimation apparatus according to an embodiment of the present application;
fig. 3 is a hardware structure diagram of a computer device in which a time offset estimation apparatus according to an embodiment of the present application is located.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments that can be derived from the embodiments given herein by a person of ordinary skill in the art are intended to be within the scope of the present disclosure.
Some relevant concepts designed by this application are first introduced:
air interface: the base station and the mobile device such as a mobile phone are referred to as an air interface because they propagate in the air by electromagnetic waves.
SRS (Sounding Reference Signal): in wireless communication, the method is used for estimating uplink channel frequency domain information and performing frequency selective scheduling; used for estimating the downlink channel and carrying out downlink beam forming.
MCS (Modulation and Coding Scheme, Modulation and Coding strategy): the configuration of the radio frequency rate of the wireless transmission protocol 802.11n is realized by the MCS index value. The MCS modulation coding table is a representation of 802.11n proposed for characterizing the communication rate of a wireless communication technology. The MCS forms a rate table with the MCS index as a row and the columns of the table as the factors that affect the communication rate concerned. Therefore, each MCS index corresponds to a physical transmission rate under a set of parameters,
with the rapid development of the mobile internet, new services and new services are continuously emerging, the flow of the mobile data service is increased explosively, the 4G mobile communication system is difficult to meet the requirement of the sudden increase of the mobile data flow, and the 5G era is produced. 5G is a novel mobile communication system, which not only solves the communication between people, but also solves the communication problem between people and objects, and between objects and objects. The concept of a new air interface (NR) has emerged in 5G mobile communication systems. The 5G era enhances the independence of base stations, which are divided into two classes, 4G and 5G. The link with the 4G base station is the "old air interface" and the connection with the 5G base station is the so-called "new air interface". The 5G new air interface is designed based on the 4G air interface technology, and also uses the modulation mode of Orthogonal Frequency Division Multiplexing (OFDM), corrects some unreasonable parts of 4G on the frame structure, and increases the support for large connection and low time delay, thereby being more flexible and having higher frequency spectrum efficiency.
In the NR communication system, both the propagation of signals in space and the transmission of signals over interfaces cause time offsets of transmitted and received signals. In order to ensure the sampling accuracy, the signals need to be time-offset synchronized. In order to solve the problem, a random access procedure is proposed, which eliminates a part of the time offset, but still retains a certain residual time offset. Therefore, the receiving end, i.e. the base station, needs to be able to eliminate the residual time offset, so as to better achieve the time offset synchronization.
The NR system uses a Sounding Reference Signal (SRS) to evaluate the channel quality of an uplink, and measures related parameters, and then reports to a higher layer, thereby implementing flexible scheduling. The application is time domain based time offset estimation. Since the SRS is a frequency domain signal, it needs to be converted to the time domain for time offset estimation. Limited by the number of transform points NifftThe actual peak position and the measured peak position have a slight deviation, which may cause the result of the time offset estimation to be inaccurate. To improve the accuracy of the time offset estimation, the general scheme is to increase NifftBut with NifftThe amount of computation of the scheme itself is also increasing.
In summary, the above increase of N is usedifftThe technical scheme of (1) has the contradiction between the time offset estimation accuracy and the calculation complexity, and the calculation complexity is increased while the time offset estimation accuracy is pursued; while the calculation is easy, the time offset estimation generates larger errors.
The application provides a time offset estimation method which can be applied to a base station and can be used for estimating the time offset without increasing NifftThe accuracy of the time offset estimation is improved in the case of (1), see fig. 1. The method comprises the following steps:
s101, obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot;
s102, obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel;
s103, when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS;
s104, determining a target offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
The application provides a time offset estimation method, which is used for calibrating a time offset value through a deviation value to realize that N is not increasedifftThe time offset estimation accuracy is improved.
Firstly, a time offset value is obtained in a first time slot, an SRS is received, and the SRS with time division, frequency division, distributed frequency division and code division is separated by analyzing information such as symbol position, frequency domain position, dressing interval and the like of the SRS. As an example, a channel initial estimation may be performed, and the received SRS and a corresponding code base sequence may be conjugate-multiplied based on the SRS and a local channel, and the channel initial estimation is shown in formula (1).
Figure BDA0003283147760000061
Wherein the content of the first and second substances,
Figure BDA0003283147760000062
representing the channel function, YSRSDenotes a received SRS, X* SRSRepresenting the conjugate of the corresponding code base sequence.
FFT point number (N) of OFDM with 100M bandwidth and 30KHz subcarrierOFDM)4096 example channel, time domain channel estimation algorithm selects IFFT points NifftSelecting a sampling and extracting interval Ncomb，NcombThe value is 1, 2 and 4. The time domain time bias estimation algorithm is undersampling, order
Figure BDA0003283147760000063
nw> 1, wherein nwIs a multiple of undersampling, NifftSelection of IFFT points N representing time domain channel estimation algorithmcombRepresenting the sampling interval, and obtaining H after initial estimation of channelest(k) As shown in equation (2).
Figure BDA0003283147760000064
Where k is the sampled value of the frequency, Hest(k) Representing the frequency domain function of the channel, H (k) representing the channel power, NTARepresenting the actual value of the time-domain peak, nHRepresents the maximum series of fourier transforms, and w (k) represents the frequency domain noise function.
NifftThe point IFFT transformed into the time domain results in the following equation (3).
Figure BDA0003283147760000065
Where N is the sample value of time, h (N) represents the time domain function of the channel, NifftThe time domain channel estimation algorithm selects IFFT points, k is a sampling value of frequency, H (k) represents channel power, nHRepresenting the maximum number of Fourier transforms, NTARepresenting the actual value of the time-domain peak, nwThe number of times of undersampling is expressed,
Figure BDA0003283147760000066
Ncombrepresents the sample decimation interval, and w (n) represents the time domain noise function.
In the first time slot, searching the peak position of the time domain of the formula (3), calculating the peak position, and setting the peak position as a first sampling point; analyzing an ideal peak position of the time domain of the formula (3), wherein the ideal peak position is the peak position of the time domain without time offset and is set as a second sampling point; and (4) calculating the difference value of the two sampling points by using a formula (4) to obtain the time offset value.
Figure BDA0003283147760000071
Wherein TO represents a time offset value, NOFDMNumber of FFT points, N, representing OFDMifftThe time domain channel estimation algorithm is represented by selecting IFFT points, delta N represents the difference value of the first sampling point and the second sampling point, NcombRepresenting the sample decimation interval and Ts the sampling period.
According to the technical scheme, the sampling point positions of the two peak values are obtained through calculation, and the time offset value is calculated by using the difference value of the sampling point positions. Of course, it can be understood that there are many schemes for calculating the time offset value, for example, using impulse signal test acquisition, etc., which do not affect the implementation of the technical solution of the present embodiment, and the present embodiment is not limited herein.
Then acquiring an offset value in a second time slot; the second time slot and the first time slot are the same time slot, or the second time slot and the first time slot are different time slots separated by a specified number of time slots. First, a channel parameter of an uplink channel is obtained in a second time slot to characterize channel quality of the uplink channel, and as an example, the received channel parameter may be an MCS order. In the field of communications technologies, there are many parameters for determining channel quality, such as reference signal received power (RSPR), signal to interference plus noise ratio (SINR), and reference signal received quality (RSPQ). In this embodiment, the parameter that the channel quality of the uplink channel satisfies the preset condition is a signal-to-noise ratio of the uplink channel. The selection of these parameters does not affect the implementation of the solution of this embodiment, and the solution of this embodiment does not limit this.
And after the channel parameters of the uplink channel representing the channel quality are obtained, judging whether the channel quality of the uplink channel meets the preset conditions. When the channel quality of the uplink channel meets a preset condition, for example, in the embodiment of the present application, when the MCS order is greater than a preset threshold. Determining that the channel quality of the uplink channel meets a preset condition, for example, the uplink channel has a high signal-to-noise ratio; the high signal-to-noise ratio means that the signal-to-noise ratio of the uplink channel is greater than a preset signal-to-noise ratio threshold.
And after the judgment is established, determining a deviation value corresponding to the time deviation value according to a time domain channel function for performing time deviation estimation based on the SRS. As an example, a peak deviation is calculated based on a peak value in the time domain channel function, a first value and a second value respectively adjacent to the peak value, and a deviation function is further obtained; the independent variable in the deviation function is n, the n is used for representing the deviation variable, the value range of the n is between a third value and a fourth value, the fourth value is an undersampling factor, and the undersampling factor is determined according to IFFT points adopted during time deviation estimation; and finally, selecting a numerical value meeting a preset condition from the value range of the n according to the peak deviation and the deviation function, and determining the selected numerical value as the deviation value. One embodiment of obtaining the deviation value is detailed with the following formula as an example.
Order to
Figure BDA0003283147760000081
Then N isTA＝nwnp+nδWherein n isδ∈[0，nw). Where n is the sample value of time, h (n) represents the time domain function of the channel, npSample values representing time-domain peaks, NTARepresenting the actual value of the time-domain peak, nwThe number of times of undersampling is expressed,
Figure BDA0003283147760000082
Nifftselection of IFFT points N representing time domain channel estimation algorithmcombRepresenting the sampling interval, nδRepresenting the deviation value of the time domain peak.
The expression of the peak sample value and the first and second values respectively adjacent to the peak sample value in the time domain channel function may be as shown in equation (5).
Figure BDA0003283147760000083
Wherein n ispSample values representing time-domain peaks, h (N) representing the time-domain function of the channel, NifftThe time domain channel estimation algorithm selects IFFT points, k is a sampling value of frequency, H (k) represents channel power, nHRepresenting the maximum number of Fourier transforms, nδDeviation value, n, representing a time domain peakwThe number of times of undersampling is expressed,
Figure BDA0003283147760000084
Ncombrepresents the sample decimation interval, and w (n) represents the time domain noise function.
The channel of the present application may be considered "flat" if it meets specified conditions, e.g., if the channel has a high signal-to-noise ratio, i.e., the channel power H (k) HAWherein H isAIs a constant value for a period of time, i.e. the channel power can be regarded as a constant HA. Order to
Figure BDA0003283147760000085
nHRepresenting the maximum number of fourier transforms. The first formula reduction substituted into formula (5) results in formula (6).
Figure BDA0003283147760000086
Wherein n ispSample values representing time-domain peaks, H (n) representing the time-domain function of the channel, HARepresenting the channel power constant, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000091
nHrepresenting the maximum number of Fourier transforms, nδRepresents the deviation value of the temporal peak, and w (n) represents the temporal noise function.
Order to
Figure BDA0003283147760000092
Where x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000093
nHrepresenting the maximum number of fourier transforms. Substituting equation (5) yields equation (7).
Figure BDA0003283147760000094
Wherein n ispRepresenting the time domainSampled values of the peak, H (n) representing the time-domain function of the channel, HADenotes the channel power constant, nδThe deviation value representing the time-domain peak,
Figure BDA0003283147760000095
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000096
nHrepresenting the maximum series of fourier transforms, w (n) represents the time-domain noise function.
Let fh1＝h(np+1)-h(np-1)，
fh2＝h(np+1)+h(np-1)-2h(np)，
Wherein n ispAnd (3) representing the sampling value of the time domain peak, and h (n) representing the channel time domain function, so as to obtain a formula (8) and a formula (9).
Figure BDA0003283147760000097
Wherein n ispSample values representing time-domain peaks, h (n) representing the time-domain function of the channel, nδDeviation value, H, representing a time domain peakARepresents the channel power constant, w (n) represents the time domain noise function,
Figure BDA0003283147760000098
Figure BDA0003283147760000099
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000101
nHrepresenting the maximum number of fourier transforms.
Figure BDA0003283147760000102
Wherein n ispSample values representing time-domain peaks, h (n) representing the time-domain function of the channel, nδRepresents the deviation value of the temporal peak, w (n) represents the temporal noise function,
Figure BDA0003283147760000103
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000104
nHrepresenting the maximum number of fourier transforms.
The channel of the present application satisfies a specified condition, for example, in the case that the channel has a high signal-to-noise ratio, noise in the calculation process can be ignored, and the formula (8) and the formula (9) can be further solved as the formula (10).
Figure BDA0003283147760000105
Figure BDA0003283147760000106
Wherein HADenotes the channel power constant, nδDeviation value, f, representing the time domain peakh1＝h(np+1)-h(np-1)，fh2＝h(np+1)+h(np-1)-2h(np)，npSample values representing time domain peaks, h (n) represents a channel time domain function,
Figure BDA0003283147760000107
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000108
nHrepresenting the maximum number of fourier transforms.
Dividing the formula (10) into upper and lower formulas to obtain the peak deviation
Figure BDA0003283147760000109
Wherein f ish1＝h(np+1)-h(np-1)，fh2＝h(np+1)+h(np-1)-2h(np)，npSample values representing time domain peaks, and h (n) represents a channel time domain function.
For f (n)δ) Is obtained by Taylor expansion
Figure BDA00032831477600001010
Wherein
Figure BDA00032831477600001011
nδThe deviation value representing the time-domain peak,
Figure BDA00032831477600001012
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000111
nHrepresenting the maximum series of Fourier transforms, m being the Taylor series, cmIs a variable which varies with m, SmIs a variable that varies with m, and N is a constant. The expansion term is shown in equation (11).
f(nδ+1)-f(nδ-1)
＝2{(-2c2+4c4-6c6+8c8+…)nδ+j(c1-c3+c5+…)}+o(Δ1)f(nδ+1)+f(nδ-1)-2f(nδ)
＝2{(-c2+c4-c6+c8+…)+j(-3c1+5c3-7c5+…)nδ}+o(Δ2) (11)
Wherein n isδA deviation value representing a time domain peak, c is a variable obtained by taylor expansion,
Figure BDA0003283147760000112
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000113
nHrepresenting the maximum number of Fourier transforms, o (Δ)1)、o(Δ2) Is a taylor high order expansion, which is an infinitesimal value. And f (x) is continuously derivable when x is 1, resulting in equation (12).
f(1)＝0
Figure BDA0003283147760000114
Where N is a constant in the taylor expansion.
Thus c0＝N，c1Substituting pi (N-1) into equation (11) yields equation (13).
f(nδ+1)-f(nδ-1)＝2{πcot(π/N)nδ}+o(Δ1)
f(nδ+1)+f(nδ-1)-2f(nδ)＝2{-N-jπNnδ}+o(Δ2) (13)
Wherein n isδThe deviation value representing the time-domain peak,
Figure BDA0003283147760000115
x is an input argument, f (x) is an output function, NifftThe time domain channel estimation algorithm is expressed by selecting IFFT points, k is a sampling value of frequency,
Figure BDA0003283147760000116
nHrepresenting the maximum series of Fourier transforms, N is a constant in the Taylor expansion, o (Δ)1)、o(Δ2) Is a taylor high order expansion, which is an infinitesimal value.
The upper and lower equations of equation (13) are divided to obtain equation (14).
Figure BDA0003283147760000117
Wherein n isδDeviation value, N, representing a time domain peakifftThe time domain channel estimation algorithm is represented by selecting the number of IFFT points,
Figure BDA0003283147760000118
x is the input argument, f (x) is the output function, k is the sampled value of frequency,
Figure BDA0003283147760000121
nHrepresenting the maximum number of fourier transforms.
This is substituted into equation (10) to obtain equation (15).
Figure BDA0003283147760000122
Wherein n isδDeviation value, N, representing a time domain peakifftThe time domain channel estimation algorithm is represented by selecting the number of IFFT points,
Figure BDA0003283147760000123
denotes the peak deviation, fh1＝h(np+1)-h(np-1)，fh2＝h(np+1)+h(np-1)-2h(np)，npSample values representing time domain peaks, and h (n) represents a channel time domain function.
Let deviation function
Figure BDA0003283147760000124
Wherein n isδ∈[0，nw)，nδDeviation value, N, representing a time domain peakifftSelection of IFFT points n representing time domain channel estimation algorithmwIs undersampled by timesThe number of the first and second groups is,
Figure BDA0003283147760000125
Ncombrepresenting the sample decimation interval. g (n) represents a deviation function with an argument of n, n being used to represent a deviation variable nδThe value range of the n is between a third value and a fourth value, and the fourth value is an undersampling factor nwSaid undersampling factor nwBy estimating the number of IFFT points N used in accordance with the time offsetifftIt is determined that,
Figure BDA0003283147760000126
Ncombrepresenting the sample decimation interval.
From said n according to said peak deviation and said deviation functionδSelecting a value meeting a preset condition from the value range of (a), and determining the selected value as the deviation value of the time domain peak value. Calculating the absolute value of the difference between the peak deviation and the deviation function to obtain a reference function, i.e.
Figure BDA0003283147760000127
The argument of the reference function is nδ. From said nδIs selected to minimize the value of the reference function, and the selected value is determined to be the one satisfying the preset condition, i.e. the value satisfying the preset condition
Figure BDA0003283147760000128
The application adopts an exhaustive method to respectively calculate
Figure BDA0003283147760000129
Obtaining a correspondence of the minimum value of the reference function
Figure BDA00032831477600001210
Wherein the content of the first and second substances,
Figure BDA00032831477600001211
denotes the peak deviation, g (n)δ) Shows a deviationDifference function, nδDeviation value, n, representing a time domain peakwIs a multiple of the undersampling of the sample,
Figure BDA00032831477600001212
Nifftselection of IFFT points N representing time domain channel estimation algorithmcombRepresenting the sampling interval, fh1＝h(np+1)-h(np-1)，fh2＝h(np+1)+h(np-1)-2h(np)，npSample values representing time domain peaks, and h (n) represents a channel time domain function. The above method is only one method for obtaining the minimum value of the reference function, and there are many methods for obtaining the minimum value of the reference function, and for example, the minimum value of the structural function, the dichotomy, the plotting, and the like may be used, and the present application does not limit the method.
Of course, it should be understood that the above is only for the purpose of finding the deviation value n of the time domain peakδThe present application is not limited to such a method. For example, when the time domain channel estimation algorithm selects the number of IFFT points NiiftWhen larger, the deviation function can be used
Figure BDA0003283147760000131
Wherein n isδ∈[0，nw)，nδDeviation value, N, representing a time domain peakifftSelection of IFFT points n representing time domain channel estimation algorithmwIs a multiple of the undersampling of the sample,
Figure BDA0003283147760000132
Ncombrepresents the sampling interval, and is further decomposed into formula (16):
Figure BDA0003283147760000133
wherein g (n)δ) Representing a deviation function, nδDeviation value, N, representing a time domain peakifftThe time domain channel estimation algorithm is represented by selecting the number of IFFT points.
The time domain peak deviation value n can then be solved using equation (17)δ
Figure BDA0003283147760000134
Wherein the content of the first and second substances,
Figure BDA0003283147760000135
denotes the peak deviation, fh1＝h(np+1)-h(np-1)，fh2＝h(np+1)+h(np-1)-2h(np)，npSample values representing time-domain peaks, h (n) representing the time-domain function of the channel, nδDeviation value, N, representing a time domain peakifftThe time domain channel estimation algorithm is represented by selecting the number of IFFT points.
Thus, the deviation value n for the time domain peakδThe calculation method in (2) does not affect the implementation of the technical solution in this embodiment, and this embodiment does not limit this.
After the deviation value is obtained, calibrating the time deviation value to obtain a target time deviation value, namely the target time deviation value position is the sum of the time deviation value position and the deviation value position, namely NTA＝nwnp+nδWherein N isTARepresenting the actual value of the time-domain peak, nwThe number of times of undersampling is expressed,
Figure BDA0003283147760000141
Nifftselection of IFFT points N representing time domain channel estimation algorithmcombRepresenting the sampling interval, npSample values representing time-domain peaks, nδRepresenting the deviation value of the time domain peak.
When the channel quality of the uplink channel does not satisfy a preset condition, for example, in the embodiment of the present application, when the MCS order is smaller than a preset threshold. The channel quality of the uplink channel does not meet a preset condition, for example, the uplink channel does not have a high signal-to-noise ratio condition; the high signal-to-noise ratio means that the signal-to-noise ratio of the uplink channel is greater than a preset signal-to-noise ratio threshold. And determining the time bias value as the target time bias value, namely determining the target time bias value as the time bias value.
That is to say NTAThe calculation is shown in equation (18).
Figure BDA0003283147760000142
Wherein N isTARepresenting the actual value of the time-domain peak, nwThe number of times of undersampling is expressed,
Figure BDA0003283147760000143
Nifftselection of IFFT points N representing time domain channel estimation algorithmcombRepresenting the sampling interval, npThe sample values representing the time domain peaks,
Figure BDA0003283147760000144
deviation value, MCS, representing time domain peakindexIndicates MCS order, MCSrefIndicating a system specific value for giving the MCS order reference.
Based ultimately on the resulting NTAThe target offset value is obtained by equation (19).
Figure BDA0003283147760000145
Where TO _ final represents the target offset, NOFDMNumber of FFT points, N, representing OFDMifftSelection of IFFT points N representing time domain channel estimation algorithmTARepresenting the actual value of the time-domain peak, NcombRepresenting the sample decimation interval and Ts the sampling period.
And denoising the obtained target time bias value, and performing time bias synchronization according to the denoised target time bias value.
Corresponding to the above method embodiment, the present application further provides a time offset estimation apparatus, which can be applied to a base station, and as shown in fig. 2, the time offset estimation apparatus includes:
a time offset estimation unit 110 that obtains a time offset value obtained when performing time offset estimation based on a channel sounding reference signal SRS in a first time slot;
a receiving unit 120, configured to obtain a channel parameter of an uplink channel in a second time slot, where the channel parameter is used to characterize channel quality of the uplink channel;
a deviation determining unit 130, configured to determine, when the channel quality of the uplink channel meets a preset condition, a deviation value corresponding to the time deviation value according to a time domain channel function when performing time deviation estimation based on the SRS;
a synchronization unit 140 for determining a target offset value according to the offset value and the time offset value; and the target time offset value is used for performing time offset synchronization of the signals.
In addition, the present application also provides an electronic device, as shown in fig. 3, including:
a processor;
a memory for storing processor-executable instructions;
wherein the processor implements the steps in the various embodiments of the time offset estimation method described above by executing the executable instructions.
Furthermore, the present application also provides a computer-readable storage medium, on which computer instructions are stored, which when executed by a processor, implement the steps in the various embodiments of the time offset estimation method described above.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be essentially or partially implemented in the form of a software product, which may be stored in a storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, etc., and includes several instructions for enabling a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the embodiments or some parts of the embodiments of the present application.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the apparatus embodiment, since it is substantially similar to the method embodiment, it is relatively simple to describe, and reference may be made to some descriptions of the method embodiment for relevant points. The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The foregoing is directed to embodiments of the present application and it is noted that numerous modifications and adaptations may be made by those skilled in the art without departing from the principles of the present application and are intended to be within the scope of the present application.

Claims (10)

1. A method of time offset estimation, comprising:
obtaining a time offset value obtained when time offset estimation is carried out on the basis of a channel Sounding Reference Signal (SRS) in a first time slot;
obtaining channel parameters of an uplink channel in a second time slot, wherein the channel parameters are used for representing the channel quality of the uplink channel;
when the channel quality of the uplink channel meets a preset condition, determining an offset value corresponding to the time offset value according to a time domain channel function when time offset estimation is carried out based on the SRS;
determining a target time offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
2. The method of claim 1, wherein the channel parameter is an MCS order;
the channel quality of the uplink channel satisfying a preset condition includes:
and when the MCS order is greater than a preset threshold value, the channel quality of the uplink channel meets a preset condition.
3. The method of claim 1, wherein the channel quality of the uplink channel satisfying a predetermined condition comprises the uplink channel having a high signal-to-noise ratio;
the high signal-to-noise ratio means that the signal-to-noise ratio of the uplink channel is greater than a preset signal-to-noise ratio threshold.
4. The method of claim 1, further comprising:
and when the channel quality of the uplink channel does not meet a preset condition, determining the time offset value as the target time offset value.
5. The method of claim 1, wherein the determining the offset value corresponding to the time offset value according to a time domain channel function when estimating the time offset based on the SRS comprises:
calculating a peak deviation based on a peak value in the time domain channel function, a first value and a second value which are respectively adjacent to the peak value;
obtaining a deviation function; the independent variable in the deviation function is n, the n is used for representing the deviation variable, the value range of the n is between a third value and a fourth value, the fourth value is an undersampling factor, and the undersampling factor is determined according to IFFT points adopted during time deviation estimation;
and selecting a numerical value meeting a preset condition from the value range of the n according to the peak deviation and the deviation function, and determining the selected numerical value as the deviation value.
6. The method of claim 5, wherein selecting a value from the range of values of n that satisfies a predetermined condition based on the peak deviation and the deviation function comprises:
calculating the absolute value of the difference between the peak value deviation and the deviation function to obtain a reference function; the argument of the reference function is the n;
selecting a value from the value range of n, the selected value minimizing the value of the reference function; the selected value is determined as the one satisfying the preset condition.
7. The method of claim 1, further comprising:
and denoising the target time bias value, and performing time bias synchronization according to the denoised target time bias value.
8. A time offset estimation device based on the method of any one of claims 1 to 7, comprising:
a time offset estimation unit for obtaining a time offset value obtained when performing time offset estimation based on a channel Sounding Reference Signal (SRS) in a first time slot;
a receiving unit, configured to obtain a channel parameter of an uplink channel in a second time slot, where the channel parameter is used to characterize channel quality of the uplink channel;
a deviation determining unit configured to determine a deviation value corresponding to the time deviation value according to a time domain channel function when performing time deviation estimation based on the SRS, when the channel quality of the uplink channel satisfies a preset condition;
the synchronization unit determines a target offset value according to the time offset value and the deviation value; and the target time offset value is used for performing time offset synchronization of the signals.
9. An electronic device, comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor implements the method of any one of claims 1-7 by executing the executable instructions.
10. A computer readable storage medium having stored thereon computer instructions which, when executed by a processor, carry out the steps of the method according to any one of claims 1 to 7.
CN202111138548.9A 2021-09-27 2021-09-27 Time offset estimation method, device, electronic equipment and storage medium Active CN113810326B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111138548.9A CN113810326B (en) 2021-09-27 2021-09-27 Time offset estimation method, device, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111138548.9A CN113810326B (en) 2021-09-27 2021-09-27 Time offset estimation method, device, electronic equipment and storage medium

Publications (2)

Publication Number Publication Date
CN113810326A true CN113810326A (en) 2021-12-17
CN113810326B CN113810326B (en) 2023-07-25

Family

ID=78896878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111138548.9A Active CN113810326B (en) 2021-09-27 2021-09-27 Time offset estimation method, device, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN113810326B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114745245A (en) * 2022-06-10 2022-07-12 广州世炬网络科技有限公司 Method, device and equipment for adjusting time frequency offset and storage medium
CN115988629A (en) * 2023-03-13 2023-04-18 新华三技术有限公司 Timing estimation method, device, equipment and readable storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6177906B1 (en) * 1999-04-01 2001-01-23 Arraycomm, Inc. Multimode iterative adaptive smart antenna processing method and apparatus
CN102710562A (en) * 2012-06-01 2012-10-03 武汉邮电科学研究院 Combined time bias estimation method based on phase position and channel impact
CN106878205A (en) * 2015-12-10 2017-06-20 电信科学技术研究院 A kind of timing offset method of estimation, device and terminal
CN112217756A (en) * 2019-07-12 2021-01-12 华为技术有限公司 Doppler frequency offset estimation method and communication device
CN112422250A (en) * 2020-11-06 2021-02-26 上海擎昆信息科技有限公司 Method and system for multi-reference signal joint measurement time offset

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6177906B1 (en) * 1999-04-01 2001-01-23 Arraycomm, Inc. Multimode iterative adaptive smart antenna processing method and apparatus
CN102710562A (en) * 2012-06-01 2012-10-03 武汉邮电科学研究院 Combined time bias estimation method based on phase position and channel impact
CN106878205A (en) * 2015-12-10 2017-06-20 电信科学技术研究院 A kind of timing offset method of estimation, device and terminal
CN112217756A (en) * 2019-07-12 2021-01-12 华为技术有限公司 Doppler frequency offset estimation method and communication device
CN112422250A (en) * 2020-11-06 2021-02-26 上海擎昆信息科技有限公司 Method and system for multi-reference signal joint measurement time offset

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114745245A (en) * 2022-06-10 2022-07-12 广州世炬网络科技有限公司 Method, device and equipment for adjusting time frequency offset and storage medium
CN115988629A (en) * 2023-03-13 2023-04-18 新华三技术有限公司 Timing estimation method, device, equipment and readable storage medium
CN115988629B (en) * 2023-03-13 2023-06-27 新华三技术有限公司 Timing estimation method, apparatus, device and readable storage medium

Also Published As

Publication number Publication date
CN113810326B (en) 2023-07-25

Similar Documents

Publication Publication Date Title
CN113810326A (en) Method and device for estimating time offset, electronic equipment and storage medium
CN110191071B (en) Channel estimation-based measurement method and device in narrow-band Internet of things system
KR20090010238A (en) System and method of calculating noise variance
CN114615122B (en) Method and device for determining frequency offset of communication signal
WO2020238349A1 (en) Group delay ripple calibration method, storage medium and electronic apparatus
JP2019009493A (en) Testing device for mimo mode system, and testing method
CN104348763A (en) Channel measuring method for large-scale antennas, and user terminal
CN112449403B (en) Random access channel transmission method and device in low-orbit satellite communication
CN102075465A (en) Method and equipment for estimating noise power
CN110430156A (en) The frame synchornization method and system of the OFDM data that happens suddenly transmission
CN113132284B (en) Carrier phase tracking method and device
CN101312597B (en) Channel estimation unit and channel estimation method
CN113938358B (en) Time delay determination method and terminal
CN111726308B (en) Orthogonal matching pursuit channel estimation method based on frequency response pre-interpolation
CN111917676B (en) Linear frequency modulation interference cancellation method
JP5099124B2 (en) Wireless communication quality estimation method and apparatus
CN113709658B (en) Method, device and equipment for estimating moving speed of terminal
CN112883787A (en) Short sample low-frequency sinusoidal signal parameter estimation method based on spectrum matching
JP5129832B2 (en) Interference wave power measuring apparatus, interference wave power measuring method, SIR measuring apparatus and SIR measuring method
CN103856960A (en) Communication terminal and cell measurement method and device in long term evolution system
CN115987742B (en) OFDM communication system-based frequency domain Doppler frequency offset expansion resisting method
CN111181889B (en) Frequency offset estimation sample receiving control method, system, equipment and storage medium
CN114500192A (en) Channel estimation method and terminal equipment
CN102130874A (en) Channel estimation method and device
KR100902842B1 (en) Fading channel modeling equipment of communication system and channel modeling method thereof

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