CN110120925B - Time delay estimation method based on short wave diversity signals - Google Patents
Time delay estimation method based on short wave diversity signals Download PDFInfo
- Publication number
- CN110120925B CN110120925B CN201910358390.2A CN201910358390A CN110120925B CN 110120925 B CN110120925 B CN 110120925B CN 201910358390 A CN201910358390 A CN 201910358390A CN 110120925 B CN110120925 B CN 110120925B
- Authority
- CN
- China
- Prior art keywords
- signals
- time delay
- delay
- frame
- paths
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000005070 sampling Methods 0.000 claims description 21
- 238000009432 framing Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000005562 fading Methods 0.000 abstract description 15
- 238000004088 simulation Methods 0.000 abstract description 9
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000004891 communication Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0222—Estimation of channel variability, e.g. coherence bandwidth, coherence time, fading frequency
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
- Radio Transmission System (AREA)
Abstract
该发明公开了一种短波信道下基于多路分集信号的多径时延估计方法，属于无线通信中的信道参数估计技术领域，涉及一种基于多路分集信号、复杂度低、无需信道先验信息时延的估计方法。在非数据辅助情况下，仅靠单通道数据，利用多路携带同一信息，统计独立且衰落不相关的信号，利用分段频域匹配相关取众数的方法，可以有利对抗短波信道的不良传输特性，完成对多径时延和相对时延的估计。仿真结果表明:在较低信噪比，估计精度仅为10‑4s量级时，对多径时延和信号相对时延的估计仍能保持较高的准确率。
The invention discloses a multi-path delay estimation method based on multi-channel diversity signals under shortwave channels, which belongs to the technical field of channel parameter estimation in wireless communication, and relates to a multi-channel diversity signal-based method with low complexity and no need for channel priori Estimation method of information delay. In the case of non-data assistance, relying on single-channel data, using multiple channels to carry the same information, counting independent and fading irrelevant signals, and using the method of segmental frequency domain matching correlation to take the mode, it can be beneficial to fight against the bad transmission of shortwave channels. characteristics, complete the estimation of multipath delay and relative delay. The simulation results show that when the SNR is low and the estimation accuracy is only in the order of 10 ‑4 s, the estimation of the multipath delay and the relative signal delay can still maintain a high accuracy.
Description
Technical Field
The invention belongs to the technical field of channel parameter estimation in wireless communication, and relates to an estimation method based on multi-path diversity signals, which is low in complexity and does not need channel prior information time delay.
Background
In short-wave communication, there are many problems such as multipath delay, doppler shift and frequency selective fading, which bring great difficulty to the reception and processing of signals, so diversity technology is often adopted to combat the poor transmission characteristics of channels. The invention can utilize diversity signal, without test or pilot frequency sequence to estimate the multi-path time delay of channel parameter and relative time delay of diversity signal.
The diversity technique is a method for improving the rate of correctly determining a received signal by providing a plurality of replicas of a transmitted signal, and is used for compensating for fading channel loss and resisting adverse effects caused by fading. Diversity specifically refers to both distributed transmission and centralized reception, where distributed transmission enables a receiving end to obtain multiple statistically independent fading signals carrying the same information. Concentrated reception is where the receiver combines multiple statistically independent fading signals received to reduce the effects of fading. In accordance with the principles of information theory, copies of the original transmitted signal having other degrees of attenuation are provided to the receiver opportunity to facilitate a correct decision of the received signal.
Currently, there are three types of commonly used multipath delay estimation models. The first type: estimating the time delay between a transmission pulse and an echo based on single-channel data, wherein the auxiliary pilot frequency data is required to have good correlation performance; the second type: also based on single-channel data, the receiving end knows that the transmitted pilot frequency auxiliary data has good correlation; in the third category: based on the multi-channel data, the time difference between different position receivers is used for estimation. There are many classical algorithms under these three models: such as correlation algorithms, algorithms based on cost functions, algorithms based on signal space characteristics, etc., but these algorithms all require the receiver to know the particular test sequence or pilot signal transmitted, which is difficult to achieve in uncooperative communications. Therefore, the scheme of the invention considers a time delay estimation method under a new scene, namely the prior knowledge of unknown transmitted signals of a receiver is only based on single-channel multi-path diversity signals to complete the estimation of the multi-path time delay of the channel and the relative time delay between the signals.
Disclosure of Invention
The invention creatively utilizes the multi-path diversity signals commonly used by the short wave channel to estimate the multi-path time delay of the short wave channel and the relative time delay between the signals. And matching the multi-channel signal in a segmented frequency domain to obtain a mode, searching a reasonable secondary mode while obtaining the relative time delay of the signal, wherein the difference between the two modes is the estimated value of the multi-path time delay of the channel. The channel of the scheme is modeled by the combination of a fading process representing a direct path and a fading process representing a reflection path, a received signal comprises two path components, the method is suitable for an MFSK modulation mode of a frequency selective fading short-wave channel, the algorithm complexity is low, the required information only comprises a symbol rate and a sampling frequency, other prior information is not required to be known, and the method can also be used for a non-cooperative communication scene. Simulation shows that the scheme has accurate estimation on multipath time delay and relative time delay and can still ensure ideal performance under a lower signal-to-noise ratio.
The technical scheme of the invention is as follows:
and respectively framing the multi-path input signals in a non-overlapping way according to the length of five symbols, carrying out frequency domain matching correlation on the multi-path data frame by frame, enabling the position of the maximum correlation value of each frame to correspond to the relative time delay estimation value of the frame, and taking the mode of the relative time delay of each frame as the final relative time delay estimation value of the signals. And searching a secondary mode meeting the requirement in each frame of relative time delay, wherein the difference between the secondary mode and the final estimated value is not more than the short wave channel multipath time delay empirical value and is not less than two sampling points. And the difference between the final relative delay estimation value and the secondary mode is the multi-path delay of the channel.
The invention relates to a multipath time delay estimation method based on multipath diversity signals under a short wave channel, which comprises the following steps:
step 2, framing the two paths of signals respectively;
the number of sampling points Nsamp in a symbol is obtained from the known symbol rate Symr and the sampling frequency Fs, the signals are divided into frames without overlapping by using the frame length fl as 5, if the frame length of the last frame is not enough, the last frame is abandoned, so the total frame number is large
The m-th receiving sequence r with the length of Nm(n) becomes after framing
Wherein each element represents a frame of data after framing;
step 3, performing frequency domain matching correlation on the two paths of signals frame by frame, and recording peak positions;
respectively carrying out fast Fourier transform on each frame data of the two paths of signals to obtain the frequency spectrum Y of the two paths of signalsm(k) M is 1, 2; and (k) multiplying the k by the conjugate, performing inverse Fourier transform, and obtaining a correlation value set C after modulus value extractionk(fr), k ═ 1,2,. ann, fn; fr 1, 2.., fl, and finally the correlation peak position p (k) max (C) is recordedk(fr)),k＝1,2,...,fn；fr＝1,2,...,fl；
Step 4, taking the mode P of the peak position P (k) of the correlation valuemObtaining two paths of signals relative time delay Rdelay ═ Pm-1; if the relative delay Rdelay' is greater than
The correction delay is rfl-Rdelay', otherwise, rdlay is maintained;
step 5, searching the secondary mode P meeting the requirementsm；
Calculating the number of sampling points Ns (Fs x Mdelay ') corresponding to Mdelay ' according to the sampling frequency Fs and the empirical value Mdelay ' of the maximum multipath time delay of the short wave channel; for the correlation value peak position { P (k) ·, k ═ 1, 2.. times, fn }, the set Φ ═ P (k) |2 < | P (k) — P) satisfying the condition that the difference from the final correction delay Rdelay is greater than two sample points and less than the number Ns of sample points constitutes the set Φ ═ P (k) |2 < | P (k) > -PmNs < k 1,2, fn, P, phismThe obtained secondary mode is obtained;
step 6, the multipath time delay Mdelay is the corrected time delay Rdelay and the sub-mode PsmThe difference, i.e. Mdelay |. Rdelay-Psm|。
The invention has the beneficial effects that: under the condition of non-data assistance, the estimation method can be favorable for resisting the poor transmission characteristic of a short wave channel and finishing the estimation of multipath time delay and relative time delay by only using single-channel data, utilizing multiple channels carrying the same information, counting independent and fading uncorrelated signals and utilizing a segmented frequency domain matching correlation mode taking method. Simulation results show that the estimation precision is only 10 at a lower signal-to-noise ratio-4And when the signal is in s magnitude, the estimation of the multipath time delay and the relative time delay of the signal can still keep higher accuracy.
Drawings
Fig. 1 is a flow chart of delay estimation.
FIG. 2 is a block diagram of a short-wave ionospheric channel model.
Fig. 3 is a diagram of performance simulation of the inventive scheme versus delay estimation.
Fig. 4 is a simulation diagram of the performance of multipath delay estimation in the scheme of the present invention.
Detailed Description
The technical scheme of the invention is detailed below by combining the accompanying drawings and the embodiment. It should be understood that the scope of the present invention is not limited to the following examples, and any techniques implemented based on the present disclosure are within the scope of the present invention.
The simulation parameters are set as follows:
taking 2FSK modulation signal as an example, let symbol rate Symr be 200sps and sampling frequency Fs be 9600Hz, and the number of sampling points Nsamp in one symbol be 48. The selected simulation channel is a frequency selective fading channel with interference at the middle latitude, and the channel model is shown in fig. 2. Assume that the input signal: s (t) ═ cos (2 pi ft) or s (t) ═ ej2πft. Tau is the multipath time delay, fdIn order to maximize the doppler shift frequency,
describing the fading characteristics of a short wave channel, and outputting a signal if two variances are equal to obey a time-varying fading coefficient with a mean value of zero complex Gaussian distribution:
the steps of estimating the multipath time delay using the multipath diversity signal are as follows:
Step 2, framing the two paths of signals respectively;
the number of sampling points 48 in a symbol is obtained from the known symbol rate 200sps and sampling frequency 9600Hz, the signal is divided into frames without overlapping by sampling points with the frame length 240, if the last frame data is not enough, the last frame is discarded, so the total frame number is the same
The m-th receiving sequence r with the length of Nm(n) becomes after framing
Wherein each element represents a frame of data after framing.
Step 3, performing frequency domain matching correlation on the two paths of signals frame by frame, and recording peak positions;
respectively carrying out fast Fourier transform on each frame data of the two paths of signals to obtain the frequency spectrum Y of the two paths of signalsm(k) M is 1, 2; and (k) multiplying the k by the conjugate, performing inverse Fourier transform, and obtaining a correlation value set C after modulus value extractionk(fr), k ═ 1,2,. ann, fn; fr 1, 2.., fl, and finally the correlation peak position p (k) max (C) is recordedk(fr)),k＝1,2,...,fn；fr＝1,2,...,fl。
Step 4, taking P (k), k being 1,2mObtaining two paths of signals relative time delay Rdelay ═ Pm-1; if the relative time delay Rdelay 'is larger than 120 sampling points, the corrected time delay is Rdelay-240-Rdelay', otherwise, the Rdelay is kept;
step 5, searching the secondary mode P meeting the requirementsm；
Calculating the number 19 of sampling points corresponding to 2ms according to 9600Hz of sampling frequency and the empirical value of the maximum multipath time delay of the simulation channel for 2 ms; for { P (k), k ═ 1, 2.., fn }, the set of estimates satisfying a difference from the final relative delay greater than 2 and less than 19 samples constitutes the set Φ ═ P (k) |2 < | P (k) -PmI < 19, k ═ 1,2,.., fn }, P, ΦsmI.e. the desired secondary mode.
Step 6, the multipath time delay Mdelay is the corrected time delay Rdelay and the sub-mode PsmThe difference, i.e. Mdelay |. Rdelay-Psm|。
The performance simulation results for multipath delay and relative delay estimation are shown in fig. 3 and fig. 4, respectively, and the results show that: under the short wave frequency selective fading channel, when the signal-to-noise ratio EbN0 is 20dB and the estimation precision is 0.2ms, the method of the invention can still keep a higher precision for the estimation of the multipath time delay and the relative time delay.
Claims (1)
1. A multipath time delay estimation method based on multipath diversity signals under a short wave channel comprises the following steps:
step 1, receiving signals r of two pathsm(t), m is 1,2 is sampled with a sampling frequency Fs to obtain { r [ (r) ]m(n), m is 1, 2; n is 1,2,., N }, m represents the number of paths of the diversity signals, and N represents the data length of each path;
step 2, framing the two paths of signals respectively;
the number of sampling points Nsamp in a symbol is obtained from the known symbol rate Symr and the sampling frequency Fs, the signals are divided into frames without overlapping by using the frame length fl as 5, if the frame length of the last frame is not enough, the last frame is abandoned, so the total frame number is large
The m-th receiving sequence r with the length of Nm(n) becomes after framing
Wherein each element represents a frame of data after framing;
step 3, performing frequency domain matching correlation on the two paths of signals frame by frame, and recording peak positions;
respectively carrying out fast Fourier transform on each frame data of the two paths of signals to obtain the frequency spectrum Y of the two paths of signalsm(k) M is 1, 2; and (k) multiplying the k by the conjugate, performing inverse Fourier transform, and obtaining a correlation value set C after modulus value extractionk(fr), k ═ 1,2,. ann, fn; fr 1, 2.., fl, and finally the correlation peak position p (k) max (C) is recordedk(fr)),k＝1,2,...,fn；fr＝1,2,...,fl；
Step 4, taking the mode P of the peak position P (k) of the correlation valuemObtaining two paths of signals relative time delay Rdelay ═ Pm-1; if the relative delay Rdelay' is greater than
The correction delay is rfl-Rdelay', otherwise, rdlay is maintained;
step 5, searching the secondary mode P meeting the requirementsm；
According to the sampling frequency Fs and shortCalculating the sampling point number Ns (Fs) Mdelay ' corresponding to the Mdelay ' according to the empirical value Mdelay ' of the maximum multipath time delay of the wave channel; for the correlation value peak position { P (k) ·, k ═ 1, 2.. times, fn }, the set Φ ═ P (k) |2 < | P (k) — P) satisfying the condition that the difference from the final correction delay Rdelay is greater than two sample points and less than the number Ns of sample points constitutes the set Φ ═ P (k) |2 < | P (k) > -PmNs < k 1,2, fn, P, phismThe obtained secondary mode is obtained;
step 6, the multipath time delay Mdelay is the corrected time delay Rdelay and the sub-mode PsmThe difference, i.e. Mdelay |. Rdelay-Psm|。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910358390.2A CN110120925B (en) | 2019-04-30 | 2019-04-30 | Time delay estimation method based on short wave diversity signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910358390.2A CN110120925B (en) | 2019-04-30 | 2019-04-30 | Time delay estimation method based on short wave diversity signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110120925A CN110120925A (en) | 2019-08-13 |
| CN110120925B true CN110120925B (en) | 2021-03-30 |
Family
ID=67521671
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910358390.2A Active CN110120925B (en) | 2019-04-30 | 2019-04-30 | Time delay estimation method based on short wave diversity signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110120925B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113364541B (en) * | 2021-06-02 | 2023-04-18 | 东南大学 | Wireless signal main path transmission delay measurement method based on multi-frequency point signal intensity |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102299731A (en) * | 2010-06-28 | 2011-12-28 | 中兴通讯股份有限公司 | Transmitting device of OFDM (Orthogonal Frequency Division Multiplexing) system and method for improving cyclic delay diversity performance |
| CN104780127A (en) * | 2015-04-09 | 2015-07-15 | 浙江大学 | Multi-path channel estimation method based on time delay-Doppler R-L (Richardson-Lucy) deconvolution |
| CN107302514A (en) * | 2017-06-19 | 2017-10-27 | 武汉大学 | A kind of multi-carrier digital broadcast singal carrier smoothing high-accuracy time delay estimation method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101174868B (en) * | 2006-11-03 | 2012-01-11 | 电信科学技术研究院 | Method and device for multi-user transmission diversity |
| KR102108291B1 (en) * | 2014-04-16 | 2020-05-13 | 상하이 내셔널 엔지니어링 리서치 센터 오브 디지털 텔레비전 컴퍼니, 리미티드 | Method for generating preamble symbol, method for receiving preamble symbol, method for generating frequency domain symbol, and apparatuses |
| CN106576086B (en) * | 2014-10-23 | 2020-01-31 | 华为技术有限公司 | Channel estimation method and device |
-
2019
- 2019-04-30 CN CN201910358390.2A patent/CN110120925B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102299731A (en) * | 2010-06-28 | 2011-12-28 | 中兴通讯股份有限公司 | Transmitting device of OFDM (Orthogonal Frequency Division Multiplexing) system and method for improving cyclic delay diversity performance |
| CN104780127A (en) * | 2015-04-09 | 2015-07-15 | 浙江大学 | Multi-path channel estimation method based on time delay-Doppler R-L (Richardson-Lucy) deconvolution |
| CN107302514A (en) * | 2017-06-19 | 2017-10-27 | 武汉大学 | A kind of multi-carrier digital broadcast singal carrier smoothing high-accuracy time delay estimation method |
Non-Patent Citations (1)
| Title |
|---|
| Joint Channel Estimation and Multiuser Detection for Uplink Grant-Free NOMA;Yang Du等;《IEEE WIRELESS COMMUNICATIONS LETTERS》;20180831;第7卷(第4期);第682-685页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110120925A (en) | 2019-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101056302B (en) | Channel and Carrier Frequency Offset Estimation Method Based on UKF in OFDM System | |
| CN103716262B (en) | Channel estimation by time-domain parameter extraction | |
| CN103107969B (en) | Incremental iterative time-varying channel evaluation and inter carrier interference (ICI) elimination method of fast orthogonal frequency division multiplexing (OFDM) system | |
| JP5374508B2 (en) | Method and system for arrival time estimation of ultra-wideband multiband orthogonal frequency division multiplexed signals | |
| US8451930B2 (en) | Sparse channel estimation using order extension for frequency and time selective MIMO wireless communication networks | |
| CN101827057B (en) | Channel estimation method and channel estimator for orthogonal frequency division multiplexing (OFDM) communication system | |
| CN110048972B (en) | Underwater sound orthogonal frequency division multiplexing channel estimation method and system | |
| WO2008113216A1 (en) | A channel estimation method | |
| US7424062B2 (en) | Efficient doppler compensation method and receiver for orthogonal-frequency-division-multiplexing (OFDM) systems | |
| US7133457B2 (en) | Joint timing recovery for multiple signal channels | |
| CN110380994A (en) | Quick Bayesian matching tracks marine condition of sparse channel estimation method | |
| CN110120925B (en) | Time delay estimation method based on short wave diversity signals | |
| CN103227761B (en) | Estimation method of multi-path non-Gaussian noise channel based on empirical likelihood method | |
| US20060017613A1 (en) | High doppler channel estimation for OFD multiple antenna systems | |
| CN111628815A (en) | A Channel Estimation Method for Satellite VDES System | |
| CN102413080B (en) | Method for estimating channel in high-speed moving TDD-LTE (time division duplex-long time evolution) uplink | |
| CN103916342A (en) | Signal-to-noise-ratio estimation method | |
| CN110233687B (en) | Multi-channel data signal joint decision detection method | |
| CN110290084B (en) | Short wave channel blind symbol synchronization method based on data frequency energy peak value | |
| WO2015096716A1 (en) | Channel estimation method and device | |
| Manzoor et al. | Front-end estimation of noise power and SNR in OFDM systems | |
| Manzoor et al. | Novel SNR estimation technique in Wireless OFDM systems | |
| Kuchenbecker et al. | Estimation of the channel impulse response length and the noise variance for OFDM systems | |
| CN110336764B (en) | Short wave channel blind symbol synchronization method based on diversity signal decoding feedback | |
| US20060062333A1 (en) | Method and apparatus for channel impulse response estimation in gsm systems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |