CN102238123A - Ranging code detection method and device - Google Patents

Ranging code detection method and device Download PDF

Info

Publication number
CN102238123A
CN102238123A CN2010101628838A CN201010162883A CN102238123A CN 102238123 A CN102238123 A CN 102238123A CN 2010101628838 A CN2010101628838 A CN 2010101628838A CN 201010162883 A CN201010162883 A CN 201010162883A CN 102238123 A CN102238123 A CN 102238123A
Authority
CN
China
Prior art keywords
sequence
peak
ranging
threshold value
value
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
CN2010101628838A
Other languages
Chinese (zh)
Other versions
CN102238123B (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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CN201010162883.8A priority Critical patent/CN102238123B/en
Priority to PCT/CN2010/077739 priority patent/WO2011137631A1/en
Publication of CN102238123A publication Critical patent/CN102238123A/en
Application granted granted Critical
Publication of CN102238123B publication Critical patent/CN102238123B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S11/00Systems for determining distance or velocity not using reflection or reradiation
    • G01S11/02Systems for determining distance or velocity not using reflection or reradiation using radio waves
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L7/00Arrangements for synchronising receiver with transmitter
    • H04L7/02Speed or phase control by the received code signals, the signals containing no special synchronisation information

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a ranging code detection method, which comprises the following steps of: performing fast Fourier transform (FFT) on a received time domain signal, and acquiring a ranging signal from a ranging sub-carrier of a frequency domain signal obtained by the transform, wherein the ranging signal consists of at least two symbols on a time domain; performing phase compensation and weighting combination on sequences of the other symbols except a first symbol to obtain a combined sequence; sequentially selecting local codes from a stored code table, performing a correlation operation, a differential operation and the FFT by utilizing a selected local code sequence and the combined sequence, calculating the power of the sequence generated by the FFT, and calculating the peak to average power ratio of the sequence; and when the peak to average power ratio is greater than a first threshold value, determining that the received time domain signal carries a ranging code. The invention discloses another ranging code detection method and ranging code detection devices. By the methods and the devices, the detection rate of the ranging code can be increased.

Description

Ranging code detecting method and device
Technical field
The embodiment of the invention relates to the mobile broadband wireless field, especially relates to ranging code detecting method and device in this field.
Background technology
IEEE (Institute of Electrical and Electronics Engineers, Institute of Electrical and Electronics Engineers) 802.16e be a kind of be the standard of mobile wide-band wire-less accessing system on basis with OFDMA (Orthogonal Frequency Division MultiplexingAccess, OFDM) technology.Range finding is a kind of process that is used to adjust mobile subscriber's carrier wave frequency deviation, timing offset and received power in the 802.16e agreement regulation OFDMA system.
Usually, the base station can specify time to allow the mobile subscriber find range, and this Resource Block is referred to as distance measurement zone (ranging region).802.16e agreement has been stipulated 4 kinds of distance measuring methods, one group of corresponding ranging code can be specified at different distance measuring methods in the base station, and this ranging code is a kind of CDMA (Code DivisionMultiple Access, code division multiple access inserts) sign indicating number.
When carrying out initial ranging, the mobile subscriber selects a ranging code to send at the distance measurement zone of appointment from the code table of base station appointment, is used for finishing the system parameters estimation of uplink distance measuring process and adjusting.The base station detects ranging code and the information such as calculating is inclined to one side at that time, frequency deviation that the mobile subscriber sends from the distance measuring signal that receives, feed back to the mobile subscriber who sends ranging code then.The mobile subscriber adjusts the transmission parameter of oneself again according to these information.
At present the ranging code detecting method that adopts mainly is that autocorrelation performance according to ranging code carries out, and can be divided into time domain correlation method and frequency domain correlation method.The former is meant that the base station is not FFT (FastFourier Transform, fast fourier transform) to the time domain data that receives, and directly carries out relevant, peak detection operation with local code; After the latter was meant that the base station is to the time domain data FFT conversion that receives, the value that extracts on the range finding subchannel was carried out related operation, IFFT operations such as (Inverse Fast Fourier Transform, inverse fast Fourier transforms) with local code again.
In actual practicality, problem such as amount of calculation is very big because the time domain correlation method exists, the interference that is subject to the data user and seldom being used.The thought of frequency domain correlation method is succinct, realizes that introducing simple, the IFFT computing makes amount of calculation be less than the time domain correlation method greatly again, so is used to carry out the detection of ranging code in the reality usually.
Weigh range finding detection algorithm performance two indexs are generally arranged: verification and measurement ratio and false alarm rate.False alarm rate is meant and detects the probability that the user does not have the ranging code of transmission.Along with various advanced technologies such as MIMO (Multiple-Input Multiple-Out-put, multiple-input and multiple-output), the utilization of beam forming technique beamforming etc., performance index such as the descending covering of base station, throughput have obtained significantly promoting, and the sensitivity of system's operate as normal (being signal to noise ratio) also further reduces.But the verification and measurement ratio of frequency domain correlation method under low signal-to-noise ratio is extremely low, makes systematic function be subject to up initial ranging sign indicating number and inserts detection.The ranging detection method that adopts has only been considered the application scenarios of signal to noise ratio more than 0dB mostly at present, does not take in for the situation of signal to noise ratio far below 0dB, and therefore performance is undesirable under low signal-to-noise ratio.In addition, a lot of ranging detection methods all are at AWGN (Additive white Gaussian noise, what additive white Gaussian noise) derivation was come out under the channel, and in the actual application environment, channel mostly is multidiameter fading channel, with the environment of detection method and do not match, can cause the verification and measurement ratio of ranging code is further descended.
Summary of the invention
The embodiment of the invention provides a kind of ranging code detecting method, and is extremely low in order to solve the verification and measurement ratio of ranging code detecting method under low signal-to-noise ratio that proposes in the prior art, makes systematic function be subject to up initial ranging sign indicating number and inserts the problem that detects, and comprising:
The time-domain signal that receives through single antenna is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain;
Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones;
The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol;
In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation;
The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates;
According to the power of each value in the sequence that after difference operation and conversion, generates, determine peak-to-average force ratio, described peak-to-average force ratio is the sequence that generates the respectively peak power of value and the ratio of average power after difference operation and conversion;
When peak-to-average force ratio during greater than first threshold value, carry ranging code in the time-domain signal of determining to receive, described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
The embodiment of the invention also provides a kind of ranging code detecting method, and is extremely low in order to solve the verification and measurement ratio of ranging code detecting method under low signal-to-noise ratio that proposes in the prior art, makes systematic function be subject to up initial ranging sign indicating number and inserts the problem that detects, and comprising:
During a plurality of time-domain signal of receiving via many antennas, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation: time-domain signal is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
The sequence of a plurality of sequences that generate through difference operation of obtaining respectively is worth corresponding addition, obtains the sequence after the addition, and carry out FFT;
Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
When peak-to-average force ratio during greater than first threshold value, carry ranging code in the time-domain signal of determining to receive, described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
The embodiment of the invention also provides a kind of ranging code checkout gear, and is extremely low in order to solve the verification and measurement ratio of ranging code detecting method under low signal-to-noise ratio that proposes in the prior art, makes systematic function be subject to up initial ranging sign indicating number and inserts the problem that detects, and comprising:
Compensating module is used for the time-domain signal that receives through single antenna is carried out fast fourier transform FFT, obtains distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, and described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones;
Weighting block is used for the sequence that sequence and other each symbols to first symbol obtain and is weighted merging, the sequence after obtaining to merge after phase compensation, wherein the weights of weighting are positioned at (0,1) between, each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol;
Computing module is used for selecting local code successively at the code table of having stored, utilizes the sequence of the local code of selecting and the sequence after the merging to carry out related operation; The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates; Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
First determination module is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
The embodiment of the invention provides a kind of ranging code checkout gear, and is extremely low in order to solve the verification and measurement ratio of ranging code method under low signal-to-noise ratio that proposes in the prior art, makes systematic function be subject to up initial ranging sign indicating number and inserts the problem that detects, and comprising:
Obtain module, during a plurality of time-domain signal of being used for receiving via many antennas, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation: time-domain signal is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
Summation module, the sequence that is used for a plurality of sequences that generate through difference operation that will obtain respectively is worth corresponding addition, obtains the sequence after the addition, and carries out FFT; Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
The 3rd determination module is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
The method that provides according to the embodiment of the invention, two above symbols that distance measuring signal is comprised carry out addition and merge and can improve signal to noise ratio on frequency domain, and then can solve mention in the background technology cause the ranging code verification and measurement ratio low because of low signal-to-noise ratio, thereby make system be subject to up initial ranging sign indicating number and insert the problem that detects, improve the performance of system.
Description of drawings
The flow chart of first kind of ranging code detecting method that Fig. 1 provides for the embodiment of the invention;
Fig. 2 for the embodiment of the invention provide relevant after frequency domain sequence in the flow chart of first kind of definite method of power of each frequency domain value correspondence;
The flow chart of second kind of ranging code detecting method that Fig. 3 provides for the embodiment of the invention;
Fig. 4 for the embodiment of the invention provide relevant after frequency domain sequence in the flow chart of second kind of definite method of power of each frequency domain value correspondence;
The particular flow sheet of the ranging code detecting method that Fig. 5 provides for the embodiment of the invention;
First kind of structural representation of first kind of ranging code checkout gear that Fig. 6 provides for the embodiment of the invention;
The structural representation of the computing module that Fig. 7 provides for the embodiment of the invention;
Second kind of structural representation of first kind of ranging code checkout gear that Fig. 8 provides for the embodiment of the invention;
First kind of structural representation of second kind of ranging code checkout gear that Fig. 9 provides for the embodiment of the invention;
The structural representation of the acquisition module that Figure 10 provides for the embodiment of the invention;
The structural representation of the summation module that Figure 11 provides for the embodiment of the invention;
Second kind of structural representation of second kind of ranging code checkout gear that Figure 12 provides for the embodiment of the invention.
Embodiment
Prior art mentions in order to solve, the verification and measurement ratio of frequency domain correlation method under low signal-to-noise ratio is extremely low, make systematic function be subject to up initial ranging sign indicating number and insert the problem that detects, the embodiment of the invention provides two kinds of ranging code detecting methods, a kind of application scenarios that utilizes the single antenna received signal that is applicable to, another kind is applicable to the application scenarios that utilizes many antenna receiving signals.
For solving the problems of the technologies described above, first kind of ranging code detecting method that the embodiment of the invention provides is applicable to the application scenarios that utilizes the single antenna received signal, and concrete handling process comprises as shown in Figure 1:
Step 101, obtain distance measuring signal;
Concrete, the time-domain signal that receives through single antenna is carried out fast fourier transform FFT, on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, obtain distance measuring signal, distance measuring signal is made up of at least two symbols on time domain;
Step 102, carry out phase compensation;
Concrete, the sequence of other symbols of removing first outer symbol is carried out phase compensation respectively, offset is that the equivalent phase of CP (Cyclic Prefix, Cyclic Prefix) length postpones;
Step 103, be weighted merging;
Concrete, the sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, the sequence after obtaining to merge;
Wherein, the weights of weighting are positioned between (0,1), and each weights sum equals 1, and are directly proportional with the data message amount of carrying on each symbol;
Step 104, carry out related operation, difference operation and FFT;
Concrete, in the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation;
The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates;
Step 105, determine peak-to-average force ratio;
Concrete, according to the power of each value in the sequence that after difference operation and FFT conversion, generates, determine peak-to-average force ratio, peak-to-average force ratio is the sequence that generates the respectively peak power of value and the ratio of average power after difference operation and conversion;
Step 106, when big first threshold value of peak-to-average force ratio, carry ranging code in the time-domain signal of determining to receive;
Concrete, when peak-to-average force ratio during, carry ranging code in the time-domain signal of determining to receive greater than first threshold value, first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
Step 101 is when implementing, and distance measuring signal is made of at least two symbols on time domain, and wherein, on the time domain angle, other symbols that are positioned at behind first symbol are copies of first symbol; On the frequency domain angle, the subcarrier that symbol occupies on frequency domain is identical.Therefore, distance measuring signal can be regarded as same signal and experience twice different interference through being subjected to the process of independent Gaussian white noise influence behind the channel.According to the random digital signal treatment theory, the signal to noise ratio that S different white noise interference signals of experience are carried out after addition merges is original S times, and S is a positive integer.Therefore, for distance measuring signal, two above symbols carry out addition and merge and the signal to noise ratio of distance measuring signal correspondence can be enhanced about more than once on frequency domain.For example, when distance measuring signal comprises two symbols, two symbols are carried out addition merge and its signal to noise ratio can be improved 1 times on frequency domain, be about 3dB.
The CP that removes other symbols of first outer symbol need be added on the symbol ending, for other symbols, can bring the phase place deflection on the frequency domain, therefore, can not directly plural symbol be carried out addition, for addressing this problem, step 102 is when implementing, the embodiment of the invention proposes each symbol is carried out the phase compensating method that length is CP, and concrete processing mode is as follows:
According to formula
Figure GSA00000089017200081
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure GSA00000089017200082
Be the equivalent phase delay of CP length, in other words, offset is
Figure GSA00000089017200083
In the distance measuring signal transmission course, consider that can there be the delay of certain hour in distance measuring signal when arriving receiving terminal, the information of carrying on first symbol may reduce to some extent, therefore, can weights be set respectively to plural symbol, each symbol multiply by addition again after the corresponding weight value respectively.As mentioning in the step 103, the data message amount of carrying on the weights of weighting and each symbol is directly proportional.When implementing, accuracy for the frequency domain sequence that guarantees to obtain after weighting merges can further limit the weights of weighting, for example, can be further limit according to the covering radius of sub-district, the place weights to weighting, each weights is directly proportional with the covering radius of sub-district, place.Certainly, when implementing, the weights of weighting can also limit according to other cell environments or condition, and the accuracy that can increase the sequence that obtains after weighting merges gets final product.
Be the method that clearer vivider explanation utilizes the weights weighting to merge, existing two symbols are that example describes, and according to the ratio of the data message amount of carrying on two symbols two symbols are respectively established weights, use α respectively 1, α 2Expression.The sequence of two symbol correspondences multiply by addition again after the corresponding weight value respectively, and the combination type of acquisition is as follows: Y k1Y 1, k+ α 2Y 2, k', wherein, Y kRepresent the sequence that obtains after two symbols merge.
The method that provides according to the embodiment of the invention, two above symbols that distance measuring signal is comprised carry out addition and merge and can improve signal to noise ratio on frequency domain, and then can solve mention in the background technology cause the ranging code verification and measurement ratio low because of low signal-to-noise ratio, thereby make system be subject to up initial ranging sign indicating number and insert the problem that detects, improve the performance of system.
Flow process as shown in Figure 1, step 104 is carried out difference operation and FFT when implementing, the embodiment of the invention provides a kind of preferred implementation, idiographic flow as shown in Figure 2:
Step 201, the frequency domain sequence after will being correlated with are mapped to ranging sub-carrier;
Step 202, successively the sequence of carrying on the ranging sub-carrier is carried out difference operation;
Concrete, at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
Figure GSA00000089017200091
Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side, N is the numerical value of the FFT of system window size;
Step 203, carry out zero padding and fast fourier transform;
Concrete, obtain after the difference operation, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe positioned between (N/8, N), and N tPower power with 2 goes forward one by one.
Step 203 when implementing, when local code be the user send ranging code the time, the sequence that obtains after the FFT computing is an accurate δ function waveform, peak value is positioned at l 0The place, l 0Inclined to one side during for the user.
Ranging code is being detected in definite process, channel coefficients also plays certain influence to verification and measurement ratio, if the channel coefficients influence is excessive, also can reduce verification and measurement ratio to a great extent, for addressing this problem, the embodiment of the invention utilizes the characteristic in the coherence bandwidth that a kind of solution is provided, its characteristic is: when two subcarriers that participate in difference operation should be positioned within the coherence bandwidth of channel, can eliminate the influence of channel coefficients by the difference operation between subcarrier, concrete solution is as follows: get subcarrier spacing and carry out difference operation less than two subcarriers of coherence bandwidth, at this moment, it is identical to be considered as channel coefficients approx, i.e. H k=H K+n, can get through after the difference operation: Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K, and Δ f is the subcarrier spacing of 802.16e agreement regulation.
The subcarrier of getting at interval less than coherence bandwidth calculates, and is spaced apart 1 subcarrier and it is carried out conjugate multiplication eliminates the influence of fading of channel and compare with only getting, and can obtain more performance with a spot of amount of calculation.
Flow process as shown in Figure 1, step 106 is when implementing, determine whether comprise ranging code in the time-domain signal according to the comparative result of the peak-to-average force ratio and first threshold value, mention in the background technology, weigh range finding detection algorithm performance two indexs are generally arranged, verification and measurement ratio and false alarm rate, but, in existing ranging code detecting method, verification and measurement ratio and false alarm rate are directly proportional, the verification and measurement ratio height, corresponding false alarm rate also can improve, and verification and measurement ratio is low, corresponding alert rate also can reduce, can't when guaranteeing the high detection rate, obtain lower false alarm rate, for addressing this problem, the embodiment of the invention has proposed to determine in conjunction with first threshold value and second threshold value method of ranging code, specific as follows: after confirming that peak-to-average force ratio is greater than first threshold value, carry before the ranging code in the time-domain signal of determining to receive, determine that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, second threshold value is not less than presets the pairing threshold value of ranging code police's code check lower limit.
During enforcement, the first threshold value value depends on signal to noise ratio and concrete channel circumstance, and is more excellent, and the span of first threshold value is (12,16), usually, can guarantee the higher detection rate when carrying out value in (12,16).
In addition, second threshold value depends on signal to noise ratio and concrete channel circumstance too, and is more excellent, and the span of second threshold value is (3,6), can guarantee lower false alarm rate usually when carrying out value in (3,6).
The method of ranging code is determined in employing in conjunction with first threshold value and second threshold value, overcome the verification and measurement ratio and the conflicting problem of false alarm rate that cause when simple dependence peak-to-average force ratio detects ranging code, can when improving verification and measurement ratio, guarantee to obtain lower false alarm rate.
During enforcement, remove the method for determining ranging code in conjunction with first threshold value and second threshold value, promptly utilize peak-to-average force ratio and peak-to-peak method, also may have other methods of determining ranging code, can guarantee when improving verification and measurement ratio that can obtain lower false alarm rate gets final product than definite ranging code.
For solving the problems of the technologies described above, second kind of ranging code detecting method that the embodiment of the invention provides is applicable to the application scenarios that utilizes many antenna receiving signals, and concrete handling process comprises as shown in Figure 3:
Step 301, receive via many antennas a plurality of time-domain signal the time, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation;
Concrete, time-domain signal is carried out fast fourier transform FFT, on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, obtain distance measuring signal, distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
The sequence of step 302, a plurality of sequences that generate through difference operation that will obtain respectively is worth corresponding addition, obtains the sequence after the addition, and carries out FFT;
Step 303, calculate the power of the sequence that obtains after the FFT conversion, and definite peak-to-average force ratio, peak-to-average force ratio is the peak power of each value of sequence of generating and the ratio of average power after difference operation and conversion;
Step 304, when peak-to-average force ratio during greater than first threshold value, carry ranging code in the time-domain signal of determining to receive, first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
Step 301 is when implementing, and distance measuring signal is made of at least two symbols on time domain, and wherein, on the time domain angle, other symbols that are positioned at behind first symbol are copies of first symbol; On the frequency domain angle, the subcarrier that symbol occupies on frequency domain is identical.
For the range finding symbol, the CP of second symbol is added on the symbol ending, therefore relative first symbol, can bring the phase place deflection on the frequency domain, therefore, can not directly two symbols be carried out addition, for addressing this problem, referring to step 301, the embodiment of the invention proposes each symbol is carried out the phase compensating method that length is CP, and concrete processing mode is as follows:
According to formula
Figure GSA00000089017200111
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length, Be the equivalent phase delay of CP length, in other words, offset is
Figure GSA00000089017200122
In the distance measuring signal transmission course, consider that can there be the delay of certain hour in distance measuring signal when arriving receiving terminal, the information of carrying on first symbol may reduce to some extent, therefore, can weights be set respectively to plural symbol, each symbol multiply by addition again after the corresponding weight value respectively.As mentioning in the step 103, the data message amount of carrying on the weights of weighting and each symbol is directly proportional.When implementing, accuracy for the frequency domain sequence that guarantees to obtain after weighting merges can further limit the weights of weighting, for example, can be further limit according to the covering radius of sub-district, the place weights to weighting, each weights is directly proportional with the covering radius of sub-district, place.Certainly, when implementing, the weights of weighting can also limit according to other cell environments or condition, and the accuracy that can increase the sequence that obtains after weighting merges gets final product.
Be the method that clearer vivider explanation utilizes the weights weighting to merge, existing is that example describes with two symbols, according to the ratio of the data message amount of carrying on two symbols two symbols is respectively established weights, uses α respectively 1, α 2Expression.The frequency domain sequence of two symbols multiply by addition again after the corresponding weight value respectively, and the combination type of acquisition is as follows: Y k1Y 1, k+ α 2Y 2, k', wherein, Y kRepresent the frequency domain sequence that obtains after two symbols merge.
The method that provides according to the embodiment of the invention, two above symbols that distance measuring signal is comprised carry out addition and merge and can improve signal to noise ratio on frequency domain, and then can solve mention in the background technology cause the ranging code verification and measurement ratio low because of low signal-to-noise ratio, thereby make system be subject to up initial ranging sign indicating number and insert the problem that detects, improve the performance of system.
Flow process as shown in Figure 3, step 301 and step 302 are when implementing, ranging sub-carrier after the mapping is carried out difference operation, will be through the corresponding addition of the sequential value on a plurality of antennas that difference operation generates, obtain the sequence after the addition, determine the corresponding power of each point in the sequence after the addition, calculate peak-to-average force ratio, concrete processing mode comprises as shown in Figure 4:
Step 401, the sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
Step 402, successively the sequence of carrying on the ranging sub-carrier is carried out difference operation;
Concrete, successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
Figure GSA00000089017200131
Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Inclined to one side when being the time-domain signal that receives, N is the FFT of a system window size;
Step 403, the add operation of execution phase;
Concrete, will sequence correspondence addition after the difference of a plurality of antennas that difference operation generates, obtain after the addition, length is the sequence of M;
Step 404, zero padding are also carried out fast fourier transform;
Concrete, obtain after the difference operation, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe positioned between (N/8, N), and N tPower power with 2 goes forward one by one.
Step 403 when implementing, when local code be the user send ranging code the time, the sequence that obtains after the FFT computing is an accurate δ function waveform, peak value is positioned at l 0The place, l 0Inclined to one side during for the user.
Step 402 is when implementing, the corresponding addition of the sequence of a plurality of antennas that will generate through difference operation, obtain after the addition, length is the frequency division sequence of M, avoided mentioning as in other ranging code detecting methods, the data of every antenna are all carried out a FFT computing (or IFFT computing), can reduce amount of calculation significantly, save resource.
Ranging code is being detected in definite process, channel coefficients also plays certain influence to verification and measurement ratio, if the channel coefficients influence is excessive, also can reduce verification and measurement ratio to a great extent, for addressing this problem, the embodiment of the invention utilizes the characteristic in the coherence bandwidth that a kind of solution is provided, its characteristic is: when two subcarriers that participate in difference operation should be positioned within the coherence bandwidth of channel, can eliminate the influence of channel coefficients by the difference operation between subcarrier, concrete solution is as follows: get subcarrier spacing and carry out difference operation less than two subcarriers of coherence bandwidth, at this moment, it is identical to be considered as channel coefficients approx, i.e. H k=H K+n, at this moment, Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K, and Δ f is the subcarrier spacing of 802.16e agreement regulation.
The subcarrier of getting at interval less than coherence bandwidth calculates, and is spaced apart 1 subcarrier and it is carried out conjugate multiplication eliminates the influence of fading of channel and compare with only getting, and can obtain more performance with a spot of amount of calculation.。
Flow process as shown in Figure 3, step 304 is when implementing, determine whether comprise ranging code in the time-domain signal according to the comparative result of the peak-to-average force ratio and first threshold value, mention in the background technology, weigh range finding detection algorithm performance two indexs are generally arranged, verification and measurement ratio and false alarm rate, but, in existing ranging code detecting method, verification and measurement ratio and false alarm rate are directly proportional, the verification and measurement ratio height, corresponding false alarm rate also can improve, and verification and measurement ratio is low, corresponding false alarm rate also can reduce, can't when guaranteeing the high detection rate, obtain lower false alarm rate, for addressing this problem, the embodiment of the invention has proposed to determine in conjunction with first threshold value and second threshold value method of ranging code, specific as follows: after confirming that peak-to-average force ratio is greater than first threshold value, carry before the ranging code in the time-domain signal of determining to receive, determine that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, second threshold value is not less than presets the pairing threshold value of ranging code police's code check lower limit.
During enforcement, the first threshold value value depends on signal to noise ratio and concrete channel circumstance, and is more excellent, and the span of first threshold value is (12,16), usually, can guarantee the higher detection rate when carrying out value in (12,16).
In addition, second threshold value depends on signal to noise ratio and concrete channel circumstance too, and is more excellent, and the span of second threshold value is (3,6), can guarantee lower false alarm rate usually when carrying out value in (3,6).
The method of ranging code is determined in employing in conjunction with first threshold value and second threshold value, overcome the verification and measurement ratio and the conflicting problem of false alarm rate that cause when simple dependence peak-to-average force ratio detects ranging code, can when improving verification and measurement ratio, guarantee to obtain lower false alarm rate.
During enforcement, remove the method for determining ranging code in conjunction with first threshold value and second threshold value, promptly utilize peak-to-average force ratio and peak-to-peak method, also may have other methods of determining ranging code, can guarantee when improving verification and measurement ratio that can obtain lower false alarm rate gets final product than definite ranging code.
Now describing with a specific embodiment, in this example, is example with the reception of 4 antennas, 10M system, promptly, can receive four time-domain signals synchronously during enforcement, ranging sub-carrier accounts for 144 in 1024 all subcarriers, and each distance measuring signal comprises two symbols, be respectively symbol 1, symbol 2, CP length are 128, and the big or small N of system FFT window is 1024, the local code number that has been stored in the code table is 64, and concrete implementing procedure sees also Fig. 5
Step 501, each antenna receive time-domain signal, transform to frequency domain by FFT, finish the phase compensation of symbol 2, merge with the weighting of symbol 1 and handle;
Step 502, get local code and carry out related operation; Concrete, get local code C m(0≤m<64), C m=(c M, 0, c M, 1... c M, 143) and Y I, kCarry out related operation, the sequence Y after obtaining being correlated with I, kC M, k *To Y I, kC M, k *Carry out difference operation and carry out many antennas merging; Difference sequence after being combined is carried out the FFT conversion, and sequence of calculation power is obtained its power peak and average;
Step 503, carry out the judgement of ranging code according to power peak and average, and inclined to one side when obtaining.
Concrete, step 501 specifically comprises when implementing:
The time domain data that receives is transformed to frequency domain through FFT, obtain frequency domain data;
From frequency domain data, extract the data on the ranging subcarrier, obtain the ranging frequency domain data, use Y I, j, kExpression, i wherein, j, k represent reception antenna sequence number, ranging symbol sequence number and ranging subcarrier physics sequence number respectively;
Data Y to second ranging symbol I, 2, kAccording to Y k1Y 1, k+ α 2Y 2, k' carry out 128 partial compensation for the time, use Y I, 2, kSequence behind the ' expression partial compensation for the time:
Figure GSA00000089017200151
To Y I, 2, k' multiply by symbol 2 weights α 2, add the frequency domain data Y of symbol 1 I, 1, kMultiply by symbol 1 weights α 1, obtain a new sequence Y I, k, Y I, k1Y I, 1, k+ α 2Y I, 2, k'.
Concrete, step 502 specifically comprises when implementing:
To local code C mWith Y I, kCarry out the sequence of related operation and carry out the difference operation processing, that is,
Figure GSA00000089017200161
1≤n≤M wherein, M is the maximum sub-carrier number of conjugation difference, sets according to channel condition, can obtain the sequence that a data length is M by aforesaid operations;
4 reception antenna data are handled the sequence that to obtain 4 length be M respectively according to above step;
These 4 sequences are merged, obtain sequence that a length is M and it is mended 0 to N t, N t=256, obtain sequence Z n', 0≤n<256, it is as follows to merge formula:
Figure GSA00000089017200162
To Z nCarry out 256 FFT conversion, obtain sequence after the conversion, calculate the power of this sequence, and obtain its average power P M, avg, be stored in the average memory; Find out its power peak P M, peakAnd position pos m, be stored in the peak value memory and the peak memory in;
All local code C are all carried out aforesaid operations, obtain 64 P Avg, and pos m
Concrete, step 503 specifically comprises when implementing:
For current Candidate key C m, take out its peak power P M, peakWith average power P M, avg, calculate peak-to-average force ratio PAPR m:
Figure GSA00000089017200163
Take out the peak value of every other sign indicating number in this frame, and obtain the average of these peak values
Figure GSA00000089017200164
Peak value P with Candidate key Cm M, peakDivided by
Figure GSA00000089017200165
Obtain peak-to-peak and compare PTP m, calculating formula is as follows:
P peak &OverBar; = 1 63 &Sigma; 0 &le; n < 63 n &NotEqual; m P n , peak
PTP m = P m , peak P peak &OverBar; ;
Use PAPR mCompare PTP with the first threshold T H1 mCompare with the second threshold T H2, if PAPR m>TH1 and PTP m>TH2, then this Candidate key is regarded as the ranging code that the user sends, and forwards next step processing to; Otherwise then directly handle next local code.
Based on same inventive concept, the embodiment of the invention also provides a kind of ranging code checkout gear, and concrete structure comprises as shown in Figure 6:
Compensating module 601 is used for the time-domain signal that receives through single antenna is carried out fast fourier transform FFT, obtains distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, and distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones;
Weighting block 602, be used for the sequence that sequence and other each symbols to first symbol obtain and be weighted merging after phase compensation, sequence after obtaining to merge, wherein the weights of weighting are positioned at (0,1) between, each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol;
Computing module 603 is used for selecting local code successively at the code table of having stored, utilizes the sequence of the local code of selecting and the sequence after the merging to carry out related operation; The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates; According to the power of each value in the sequence that after difference operation and conversion, generates, determine peak-to-average force ratio, peak-to-average force ratio is the sequence that generates the respectively peak power of value and the ratio of average power after difference operation and conversion;
First determination module 604 is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
In one embodiment, compensating module 601 can specifically be used for: according to formula
Figure GSA00000089017200172
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure GSA00000089017200181
Equivalent phase delay for CP length.
In one embodiment, weighting block 602 can be further used for: determine the weights of weighting according to the covering radius of sub-district, single antenna place, wherein, the weights of weighting are directly proportional with the covering radius of sub-district, single antenna place.
In one embodiment, structure as shown in Figure 7, computing module 603 specifically comprises:
First mapping submodule 701 is used for the sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
The first difference subspace module 702, be used for successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
Figure GSA00000089017200182
Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side, N is the FFT of a system window size;
First transformation submodule 703, be used for obtain after the difference operation, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
In one embodiment, the first difference subspace module 702 can be further used for: when subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure GSA00000089017200183
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
In one embodiment, structure as shown in Figure 8, the ranging code checkout gear can also comprise:
Second determination module 801, be used for determining that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that described peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, described second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
Based on same inventive concept, the embodiment of the invention also provides another ranging code checkout gear, and concrete structure comprises as shown in Figure 9:
Obtain module 901, during a plurality of time-domain signal of being used for receiving via many antennas, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation: time-domain signal is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
Summation module 902, the sequence that is used for a plurality of sequences that generate through difference operation that will obtain respectively is worth corresponding addition, obtains the sequence after the addition, and carries out FFT; Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
The 3rd determination module 903 is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
In one embodiment, obtaining module 901 can be further used for: according to formula
Figure GSA00000089017200191
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure GSA00000089017200201
Equivalent phase delay for CP length.
In one embodiment, obtaining module 901 can be further used for: determine the weights of weighting according to the covering radius of sub-district, many antennas place, wherein, the weights of weighting are directly proportional with the covering radius of sub-district, many antennas place.
In one embodiment, structure as shown in figure 10, obtaining module 901 can comprise:
Second mapping submodule 1001 is used for the sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
The second difference subspace module 1002, be used for successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
Figure GSA00000089017200202
Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side, N is the FFT of a system window size.
In one embodiment, the second difference subspace module 1002 can be further used for: when described subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure GSA00000089017200203
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
In one embodiment, structure as shown in figure 11, summation module 902 can comprise:
Addition submodule 1101, be used for after difference operation, obtaining, length is that the sequence of the sequence of M respectively is worth corresponding addition, obtains the sequence after the addition, the sequence length after the described addition is M;
Second transformation submodule 1102, be used for obtain after the addition, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
In one embodiment, structure as shown in figure 12, the ranging code checkout gear can also comprise:
The 4th determination module 1201, be used for determining that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
The method that provides according to the embodiment of the invention, two above symbols that distance measuring signal is comprised carry out addition and merge and can improve signal to noise ratio on frequency domain, and then can solve mention in the background technology cause the ranging code verification and measurement ratio low because of low signal-to-noise ratio, thereby make system be subject to up initial ranging sign indicating number and insert the problem that detects, improve the performance of system.
Further, the method of ranging code is determined in employing in conjunction with first threshold value and second threshold value, overcome the verification and measurement ratio and the conflicting problem of false alarm rate that cause when simple dependence peak-to-average force ratio detects ranging code, can when improving verification and measurement ratio, guarantee to obtain lower false alarm rate.
Further, the subcarrier of getting at interval less than coherence bandwidth calculates, be spaced apart 1 subcarrier and it is carried out conjugate multiplication eliminate the influence of fading of channel and compare with only getting, computation amount, reduced the complexity of calculating, reduced, can save resource the taking of resource, improve resource utilization, can obtain more performance with a spot of amount of calculation.
Further, the corresponding addition of the frequency domain value of a plurality of frequency division sequences that will generate through difference operation, obtain after the addition, length is the frequency division sequence of M, avoided mentioning as in other ranging code detecting methods, each sequence that obtains after difference operation is all carried out a FFT computing, can reduce amount of calculation significantly, save resource.
Obviously, those skilled in the art can carry out various changes and modification to the present invention and not break away from the spirit and scope of the present invention.Like this, if of the present invention these are revised and modification belongs within the scope of claim of the present invention and equivalent technologies thereof, then the present invention also is intended to comprise these changes and modification interior.

Claims (30)

1. a ranging code detecting method is characterized in that, comprising:
The time-domain signal that receives through single antenna is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain;
Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones;
The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol;
In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation;
The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates;
Power according to each value in the sequence that after difference operation and conversion, generates, determine peak-to-average force ratio, described peak-to-average force ratio is the ratio of peak power and average power, and described peak power is the maximum in the power of each value in the sequence, and described average power is the mean value of the power of value respectively in the sequence;
When peak-to-average force ratio during greater than first threshold value, carry ranging code in the time-domain signal of determining to receive, described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
2. the method for claim 1 is characterized in that, the sequence of other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones, and comprising:
According to formula
Figure FSA00000089017100011
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure FSA00000089017100021
Equivalent phase delay for CP length.
3. the method for claim 1 is characterized in that, the weights of described weighting are directly proportional with the covering radius of sub-district, described single antenna place.
4. as each described method of claim 1-3, it is characterized in that, the sequence that obtains behind the related operation carried out difference operation and FFT, specifically comprise:
The sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
Successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
x [ n ] = &Sigma; k &Element; K ( k + n ) &Element; K c k * Y k H k exp ( - j 2 &pi;k N l 0 ) * ( c k + n * Y k + n H k + n exp ( - j 2 &pi; ( k + n ) N l 0 ) ) *
= &Sigma; k &Element; K ( k + n ) &Element; K H k H k + n * exp ( j 2 &pi; l 0 N n ) Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side;
Obtain after the difference operation, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
5. method as claimed in claim 4 is characterized in that, when described subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure FSA00000089017100024
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
6. as each described method of claim 1-3, it is characterized in that described first threshold value is positioned between (12,16).
7. as each described method of claim 1-3, it is characterized in that, after confirming that peak-to-average force ratio is greater than first threshold value, carry before the ranging code in the time-domain signal of determining to receive, also comprise:
Determine that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that described peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, described second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
8. method as claimed in claim 7 is characterized in that, described second threshold value is positioned between (3,6).
9. a ranging code detecting method is characterized in that, comprising:
During a plurality of time-domain signal of receiving via many antennas, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation: time-domain signal is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
The a plurality of sequences that generate through difference operation that obtain respectively are worth corresponding addition, obtain the sequence after the addition, and carry out FFT;
Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
When peak-to-average force ratio during greater than first threshold value, carry ranging code in the time-domain signal of determining to receive, described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
10. method as claimed in claim 9 is characterized in that, the sequence of other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones, and comprising:
According to formula
Figure FSA00000089017100031
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure FSA00000089017100041
Equivalent phase delay for CP length.
11. method as claimed in claim 9 is characterized in that, the weights of described weighting are directly proportional with the covering radius of sub-district, described many antennas place.
12., it is characterized in that as each described method of claim 9-11, the sequence that obtains behind the related operation is carried out difference operation, specifically comprise:
The sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
Successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula
x [ n ] = &Sigma; k &Element; K ( k + n ) &Element; K c k * Y k H k exp ( - j 2 &pi;k N l 0 ) * ( c k + n * Y k + n H k + n exp ( - j 2 &pi; ( k + n ) N l 0 ) ) *
= &Sigma; k &Element; K ( k + n ) &Element; K H k H k + n * exp ( j 2 &pi; l 0 N n ) Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side.
13. method as claimed in claim 12 is characterized in that, when described subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure FSA00000089017100044
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
14. method as claimed in claim 12 is characterized in that, the sequence of a plurality of sequences that generate through difference operation of obtaining respectively is worth corresponding addition, obtains the sequence after the addition, and carries out FFT, comprising:
After difference operation, obtain, length is that the sequence of a plurality of sequences of M respectively is worth corresponding addition, obtains the sequence after the addition, the sequence length after the described addition is M;
Obtain after the addition, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
15., it is characterized in that described first threshold value is positioned between (12,16) as each described method of claim 9-11.
16., it is characterized in that as each described method of claim 9-11, after confirming that peak-to-average force ratio is greater than first threshold value, carry before the ranging code in the time-domain signal of determining to receive, also comprise:
Determine that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that described peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, described second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
17. method as claimed in claim 16 is characterized in that, described second threshold value is positioned between (3,6).
18. a ranging code checkout gear is characterized in that, comprising:
Compensating module is used for the time-domain signal that receives through single antenna is carried out fast fourier transform FFT, obtains distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, and described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones;
Weighting block is used for the sequence that sequence and other each symbols to first symbol obtain and is weighted merging, the sequence after obtaining to merge after phase compensation, wherein the weights of weighting are positioned at (0,1) between, each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol;
Computing module is used for selecting local code successively at the code table of having stored, utilizes the sequence of the local code of selecting and the sequence after the merging to carry out related operation; The sequence that obtains behind the related operation is carried out difference operation and FFT, and calculate the power of each value in the sequence that after difference operation and conversion, generates; According to the power of each value in the sequence that after difference operation and conversion, generates, determine peak-to-average force ratio, described peak-to-average force ratio is the sequence that generates the respectively peak power of value and the ratio of average power after difference operation and conversion;
First determination module is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
19. device as claimed in claim 18 is characterized in that, described compensating module specifically is used for: according to formula
Figure FSA00000089017100061
Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure FSA00000089017100062
Equivalent phase delay for CP length.
20. device as claimed in claim 18, it is characterized in that, described weighting block is further used for: determine the weights of described weighting according to the covering radius of sub-district, described single antenna place, wherein, the weights of described weighting are directly proportional with the covering radius of sub-district, described single antenna place.
21., it is characterized in that described computing module specifically comprises as each described device of claim 18-20:
First mapping submodule is used for the sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
The first difference subspace module, be used for successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula x [ n ] = &Sigma; k &Element; K ( k + n ) &Element; K c k * Y k H k exp ( - j 2 &pi;k N l 0 ) * ( c k + n * Y k + n H k + n exp ( - j 2 &pi; ( k + n ) N l 0 ) ) *
= &Sigma; k &Element; K ( k + n ) &Element; K H k H k + n * exp ( j 2 &pi; l 0 N n ) Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side;
First transformation submodule, be used for obtain after the difference operation, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
22. device as claimed in claim 21 is characterized in that, the described first difference subspace module is further used for: when described subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure FSA00000089017100071
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
23. as each described device of claim 18-20, it is characterized in that, also comprise:
Second determination module, be used for determining that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that described peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, described second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
24. a ranging code checkout gear is characterized in that, comprising:
Obtain module, during a plurality of time-domain signal of being used for receiving via many antennas, respectively each time-domain signal is all carried out the sequence of following operation to obtain to generate through difference operation: time-domain signal is carried out fast fourier transform FFT, obtain distance measuring signal on the ranging sub-carrier of the frequency-region signal that obtains after the conversion, described distance measuring signal is made up of at least two symbols on time domain; Sequence to other symbols of removing first outer symbol is carried out phase compensation respectively, and offset is that the equivalent phase of cyclic prefix CP length postpones; The sequence that sequence and other each symbols of first symbol obtained after phase compensation is weighted merging, sequence after obtaining to merge, wherein the weights of weighting are positioned between (0,1), each weights sum equals 1, and is directly proportional with the data message amount of carrying on each symbol; In the code table of having stored, select local code successively, utilize the sequence of the local code of selecting and the sequence after the merging to carry out related operation, and the sequence that obtains behind the related operation is carried out difference operation;
Summation module, the sequence that is used for a plurality of sequences that generate through difference operation that will obtain respectively is worth corresponding addition, obtains the sequence after the addition, and carries out FFT; Calculate the power of the sequence that obtains after the FFT conversion, determine peak-to-average force ratio, described peak-to-average force ratio is the peak power of each value of the sequence that generates and the ratio of average power after difference operation and conversion;
The 3rd determination module is used for carrying ranging code when peak-to-average force ratio during greater than first threshold value in the time-domain signal of determining to receive, and described first threshold value is not less than the default pairing threshold value of the ranging code verification and measurement ratio upper limit.
25. device as claimed in claim 24 is characterized in that, described acquisition module specifically is used for: according to formula Sequence to other each symbols is carried out phase compensation respectively, and wherein, i, j, k are positive integer, i 〉=2, and N is the FFT window value, Y I, kThe frequency domain value of representing k subcarrier of i symbol, Y I, kThe frequency domain value after phase compensation of k subcarrier of i symbol of ' expression, l CpBe CP length,
Figure FSA00000089017100082
Equivalent phase delay for CP length.
26. device as claimed in claim 24, it is characterized in that, described acquisition module is further used for: determine the weights of described weighting according to the covering radius of sub-district, described many antennas place, wherein, the weights of described weighting are directly proportional with the covering radius of sub-district, described many antennas place.
27., it is characterized in that described acquisition module comprises as each described device of claim 24-26:
Second mapping submodule is used for the sequence that obtains behind the related operation is shone upon back corresponding ranging sub-carrier;
The second difference subspace module, be used for successively at a distance of for the sequence of carrying on the ranging sub-carrier of n Δ f by formula x [ n ] = &Sigma; k &Element; K ( k + n ) &Element; K c k * Y k H k exp ( - j 2 &pi;k N l 0 ) * ( c k + n * Y k + n H k + n exp ( - j 2 &pi; ( k + n ) N l 0 ) ) *
= &Sigma; k &Element; K ( k + n ) &Element; K H k H k + n * exp ( j 2 &pi; l 0 N n ) Carry out difference operation, obtain the sequence that length is M, wherein, 1≤n≤M, n, M are positive integer, and M is the largest interval of subcarrier, and Δ f is for facing the interval of two sub-intercarriers mutually, and K represents the set that all ranging sub-carrier are formed, H kAnd H K+nBe channel coefficient correlation, C kAnd C K+nBe the sequence of local code C correspondence, l 0Be the time-domain signal that receives the time inclined to one side, N is the FFT window number.
28. device as claimed in claim 27 is characterized in that, the described second difference subspace module is further used for: when described subcarrier spacing during less than coherence bandwidth, and H k=H K+n, at this moment,
Figure FSA00000089017100085
Wherein, k (n) is that spacing is the right number of carrier wave of n Δ f among the ranging sub-carrier collection K.
29. device as claimed in claim 27 is characterized in that, described summation module comprises:
The addition submodule, be used for after difference operation, obtaining, length is that the sequence of the sequence of M respectively is worth corresponding addition, obtains the sequence after the addition, the sequence length after the described addition is M;
Second transformation submodule, be used for obtain after the addition, length mends N after being the sequence of M t-M 0, is N to the length that forms after the zero padding tSequence be FFT, wherein, N tBe not less than M, be positioned between (N/8, N), and N tThe value power power that equals 2.
30. as each described device of claim 24-26, it is characterized in that, also comprise:
The 4th determination module, be used for determining that the peak-to-peak ratio is greater than second threshold value, wherein, the ratio of the average of the peak power that described peak-to-peak is more corresponding with each local code than the peak power of each value of sequence that is through difference operation and generates after changing, described second threshold value are not less than presets the pairing threshold value of ranging code police's code check lower limit.
CN201010162883.8A 2010-05-04 2010-05-04 Ranging code detecting method and device Expired - Fee Related CN102238123B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201010162883.8A CN102238123B (en) 2010-05-04 2010-05-04 Ranging code detecting method and device
PCT/CN2010/077739 WO2011137631A1 (en) 2010-05-04 2010-10-14 Method and apparatus for detecting ranging code

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010162883.8A CN102238123B (en) 2010-05-04 2010-05-04 Ranging code detecting method and device

Publications (2)

Publication Number Publication Date
CN102238123A true CN102238123A (en) 2011-11-09
CN102238123B CN102238123B (en) 2016-01-20

Family

ID=44888347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010162883.8A Expired - Fee Related CN102238123B (en) 2010-05-04 2010-05-04 Ranging code detecting method and device

Country Status (2)

Country Link
CN (1) CN102238123B (en)
WO (1) WO2011137631A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103428151A (en) * 2012-05-23 2013-12-04 中兴通讯股份有限公司 Method and apparatus for removing cyclic prefix
WO2016101658A1 (en) * 2014-12-25 2016-06-30 中兴通讯股份有限公司 Method and device for estimating frequency offset of microwave communication system channel
CN106230574A (en) * 2016-09-09 2016-12-14 中国人民解放军国防科学技术大学 Signal range-measurement system that a kind of CDMA is combined with multi-subcarrier direct sequence CDMA and distance-finding method thereof
WO2017124379A1 (en) * 2016-01-21 2017-07-27 华为技术有限公司 Method and apparatus for determining time offset
CN109150780A (en) * 2018-08-13 2019-01-04 南京荣飞科技股份有限公司 A kind of WiFi ToF range positioning system based on channel state information

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060098749A1 (en) * 2004-10-12 2006-05-11 Samsung Electronics Co., Ltd. Apparatus and method for detecting ranging signal in an orthogonal frequency division multiple access mobile communication system
CN101217818A (en) * 2008-01-03 2008-07-09 北京北方烽火科技有限公司 A ranging code detecting method of base station receiver
US7586976B1 (en) * 2006-03-03 2009-09-08 Nortel Networks Limited Initial ranging detection for OFDMA systems

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060098749A1 (en) * 2004-10-12 2006-05-11 Samsung Electronics Co., Ltd. Apparatus and method for detecting ranging signal in an orthogonal frequency division multiple access mobile communication system
US7586976B1 (en) * 2006-03-03 2009-09-08 Nortel Networks Limited Initial ranging detection for OFDMA systems
CN101217818A (en) * 2008-01-03 2008-07-09 北京北方烽火科技有限公司 A ranging code detecting method of base station receiver

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103428151A (en) * 2012-05-23 2013-12-04 中兴通讯股份有限公司 Method and apparatus for removing cyclic prefix
CN103428151B (en) * 2012-05-23 2017-11-10 中兴通讯股份有限公司 Go the method and device of cyclic prefix
WO2016101658A1 (en) * 2014-12-25 2016-06-30 中兴通讯股份有限公司 Method and device for estimating frequency offset of microwave communication system channel
WO2017124379A1 (en) * 2016-01-21 2017-07-27 华为技术有限公司 Method and apparatus for determining time offset
CN108028824A (en) * 2016-01-21 2018-05-11 华为技术有限公司 A kind of method and device of definite time migration
CN108028824B (en) * 2016-01-21 2020-04-28 华为技术有限公司 Method and device for determining time offset
US10742469B2 (en) 2016-01-21 2020-08-11 Huawei Technologies Co., Ltd. Method and apparatus for determining time offset
CN106230574A (en) * 2016-09-09 2016-12-14 中国人民解放军国防科学技术大学 Signal range-measurement system that a kind of CDMA is combined with multi-subcarrier direct sequence CDMA and distance-finding method thereof
CN106230574B (en) * 2016-09-09 2019-06-04 中国人民解放军国防科学技术大学 The signal range-measurement system and its distance measuring method of a kind of CDMA in conjunction with multi-subcarrier direct sequence CDMA
CN109150780A (en) * 2018-08-13 2019-01-04 南京荣飞科技股份有限公司 A kind of WiFi ToF range positioning system based on channel state information

Also Published As

Publication number Publication date
WO2011137631A1 (en) 2011-11-10
CN102238123B (en) 2016-01-20

Similar Documents

Publication Publication Date Title
CN110290581B (en) Rapid time-frequency synchronization method and terminal in 5G system
CN101295999B (en) Detection method for accidental precursor access
CN1985451B (en) Method for accessing a wireless communication system
CN101447970B (en) Method for conducting LOFDM system timing and carrier synchronization utilizing training sequence
EP2437450B1 (en) Device and method for estimating time offset in orthogonal frequency division multiplexing (ofdm) system
CN101729465B (en) Method and device for measuring signal to noise ratio and application thereof
CN102196486B (en) Reference signal received power measurement method and device for orthogonal frequency-division multiplexing system
CN102209052B (en) Method and device for carrying out frequency offset compensation on random access signals
KR101656083B1 (en) Apparatus and method for obtaining reception synchronization in wireless communication system
CN105245479A (en) Physical frame preamble symbol receiving and processing method
JP2020519162A (en) Method and apparatus for determining uplink synchronization timing deviation
CN104918336A (en) Detection method, device and system of random access signals in interference environment
CN103237311A (en) Orthogonal frequency division multiplexing (OFDM) frequency domain interpolation pilot frequency-based cyclostationary feature spectrum sensing method
CN103220707A (en) Method for detecting quantity of antenna ports in LTE (Long Term Evolution) system
CN102238123B (en) Ranging code detecting method and device
CN101193398A (en) An uplink distance measuring method, device and base station
CN102098114B (en) Method and device for measuring signal-to-noise ratio of system
CN101582870B (en) Method and device for realizing synchronization
CN101291311B (en) Synchronization implementing method and device for multi-input multi-output orthogonal frequency division multiplexing system
CN112449403B (en) Random access channel transmission method and device in low-orbit satellite communication
CN102594737A (en) Adjacent region interference detection method and system
CN101635598A (en) Method and device for estimating noise power
CN101364846B (en) Sub-carrier noise power estimation method based on pilot
CN103415067A (en) Signal-to-noise ratio estimation method based on detection reference signal
CN101820407B (en) Serial interference cancellation based frequency domain initial ranging method and system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160120

Termination date: 20190504

CF01 Termination of patent right due to non-payment of annual fee