CN110247866A - A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation - Google Patents
A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation Download PDFInfo
- Publication number
- CN110247866A CN110247866A CN201910310492.7A CN201910310492A CN110247866A CN 110247866 A CN110247866 A CN 110247866A CN 201910310492 A CN201910310492 A CN 201910310492A CN 110247866 A CN110247866 A CN 110247866A
- Authority
- CN
- China
- Prior art keywords
- dmwc
- frequency spectrum
- phase difference
- spectrum perception
- signal
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
Abstract
A kind of method that the present invention proposes DMWC frequency spectrum perception phase difference calibration based on DOA estimation, the problem that phase difference causes recovery rate low is randomly generated due to spatial distribution for each sensing node of DMWC frequency spectrum perception system, this method comprises: (1) sampled signalAcquisition, obtained from the DMWC frequency spectrum perception system structure that single node is arranged as uniform linear array；(2) direction of arrival of acquisition signal is estimated by DOA, the maximum direction parameter of related coefficient is found as direction of arrival using the array signal processing TLS-ESPRIT algorithm that total least square method combines, and then obtain the phase difference that node receives signal；(3) direction of arrival is utilizedWith phase difference calibration measurement matrix.It by above method, solves in DMWC frequency spectrum perception system, influence of the phase difference to signal recovery rate between each node substantially increases the noise immunity and stability of DMWC frequency spectrum perception system.
Description
Technical field
The invention belongs to the technical fields that array signal processing is intersected with frequency spectrum perception, specifically, being estimated by DOA
Method in DMWC frequency spectrum perception system array distribution node carry out phase difference calibration.
Background technique
In recent years, with the fast development of internet and mobile communication technology, wireless communication user shows explosion type increasing
It is long, so that limited frequency spectrum resource is unable to satisfy the demand of user gradually.To improve frequency spectrum resource utilization rate, Moltila is for the first time
It is proposed understanding radio technology (CR), and in CR technology, frequency spectrum perception technology plays extensive effect.Compressed sensing (CS)
It is the Undersampling technique that a kind of sampling rate is far below Nyquist sampling rate, and can preferably restores signal, breaches
The bottleneck of Nyquist sampling thheorem.Modulation wide-band transducer (MWC) is proposed by the Elder professor of MIT, by compressed sensing
Theoretical lack sampling recovery and rebuilding technology becomes actual hardware circuit.The principle of MWC is by a sensor application to more
A channel is overlapped and is mixed with original sparse broadband signal using periodic pseudo-random signal, then uses low pass filtered
Wave device filters out low frequency signal.However, the port number of sampling can be multiplied when sparse broadband signal is continuously increased, but
It is difficult to adjust to ampling channel number in MWC, and the problem of this is also Expenses Cost in practical applications.Distribution modulation
Wide-band transducer (DMWC) solves this problem, and it is based on dynamic multinode single-pass that the thinking of DMWC, which is MWC system transition,
The networking cooperation lack sampling system in road, its each channel include pseudorandom mixing, low-pass filtering and low speed sampling, Ge Gejie
The input signal in point channel is different.Each node independently executes the local sampling task in channel all the way between each other, and
Source compression sampling data transmission to processing center.
Array signal processing is an important branch of signal processing, and one group of sensor is arranged in by it in a certain way
On the different position in space, sensor array is formed.With sensor array reception space signal, it is equivalent to the field to spatial distribution
Signal sampling obtains the spatial spreading observation data of signal source.It is the mesh of array signal processing by array received signal
It is handled, useful signal required for enhancing inhibits useless interference and noise, and extract useful signal characteristic and signal
The information for being included.Space arrival direction (DOA) estimation is a basic problem of array signal processing, it is each for determining
The deflection of signal arrival array reference array element.Each channel is distributed in the week of source signal as a sensing node in DMWC
It encloses, each sensing node collaborative perception, lack sampling data summarization, place by a data collection center all sensing nodes
It manages, make unified perception judgement, to calculate the supported collection of original signal.Each sensing node in DMWC is exactly array
Array reference array element in signal processing.
Since the spatial distribution of DMWC is random, each channel is distributed in different areas as an independent sensing node
Domain, each node difference at a distance from source signal brings phase difference problem, and phase difference problem is DMWC frequency spectrum perception system
Can system efficiently restore to receive one of the key of signal.If each node signal of sensory perceptual system is asynchronous, this is not only reduced
The stability of system, also produces certain influence to recovery rate.In order to eliminate produced by DMWC frequency spectrum perception system
Node phase difference and improve the overall stability of system, using the method for estimation direction of arrival, between sampled signal node
DOA is to estimate to obtain direction of arrival and calibration matrix to realize to calculation matrix calibration, signal is synchronous between realizing node, so as to improve
The recovery rate of DMWC and the noise immunity of system.
Summary of the invention
Present invention seek to address that each node of DMWC frequency spectrum perception system receives compression sampling signal asking there are phase difference
The method of topic and a kind of DMWC frequency spectrum perception phase alignment based on DOA estimation of proposition, its technical solution is as follows:
Step 1: the matrix after being sampled by the DMWC frequency spectrum perception system based on ULA；
Step 2: pass through sampling matrix, determine that the DOA based on array estimates the direction of arrival of obtained node signal；
Step 3: the phase difference between node is obtained according to the direction of arrival of node signal, then calculation matrix is calibrated.
The utility model has the advantages that the phase difference between different nodes can be increasing with the increase of direction of arrival, lead to the extensive of signal
Multiple rate is lower and lower.Estimated between calibration node after phase difference by DOA, the recovery with the increase of direction of arrival, after phase alignment
Rate before calibration compared to gradually increasing, and maximum difference 5% when from direction of arrival being 10 ° is until recovery rate maximum phase when direction of arrival is 80 °
Poor 65%, when direction of arrival is 0 °, calibration front and back signaling protein14-3-3 rate difference reaches maximum.Therefore by DOA estimation calibration
Better than there are the DMWC frequency spectrum perception systems of phase difference on the whole for the recovery rate of DMWC frequency spectrum perception system.
Detailed description of the invention
Flow diagram of the Fig. 1 based on the DOA DMWC frequency spectrum perception system phase calibration estimated.
Functional block diagram of the Fig. 2 based on the DOA MWC estimated.
To the evaluated error relational graph of direction of arrival under Fig. 3 difference signal-to-noise ratio.
The relational graph of Fig. 4 practical direction of arrival and evaluated error.
The relational graph of Fig. 5 difference direction of arrival and signal recovery rate.
Specific embodiment
Fig. 1 is the flow chart of the supported collection quick recovery method based on MWC, and specific embodiments of the present invention step is such as
Under, it is described below in conjunction with attached drawing.
Step 1: Fig. 2 is the functional block diagram based on the DOA DMWC estimated, it is by the uniform linear array that forms
(ULA), and the distance between each adjacent sensors are d(, c is the light velocity).All sensings
Device is consistent sampling configuration all to realize the function in the single channel MWC.It is multiplied by a cycle to original sampled signal x (t)Function p (t) (to p's (t) only requirement is that its fourier series coefficient in the Nyquist Bandwidth of signal all
Be not zero) carry out pseudo-random number sequence mixing, the spectrum information in frequency band is diffused into entire frequency band, use cutoff frequency for
Filter be filtered, and withSample rate sampled to obtain sampled signal。
Step 2: determine that the DOA based on array estimates phase difference between obtained node, the method used is estimation signal
Direction of arrival.There is identical distance between adjacent array element, this identical distance is reflected between a characteristic whistle battle array
Rotational invariance realizes the estimation to signal direction of arrival using this characteristic, the specific steps are as follows:
1) the sampled data number Q of each node, the compressed signal after sampling with DMWC are determinedIt is input with Q.According to public affairs
Formula, calculate the covariance of sampled data；
2) by carrying out generalized eigen decomposition to matrix R, it may be assumed that, wherein,；
3) it takes top n maximum eigenvalue in characteristic value to constitute signal subspace estimation S, indicates that S can according to the motion immovability of array
To be decomposed into two parts, it may be assumed that,WithFor two parts after decomposition；
4) influence for considering noise, utilizes the thought pair of total least square method (TLS)Matrix decomposes.It is overall minimum
The solution of square law meets:
It can obtain:.Then, willIt decomposes to obtain matrix according to this formula,；
5) according to formulaObtain invariable rotary matrix；
6) according to two adjacent submatrix members t moment output signal vector expression:
Wherein, matrixFor diagonal matrix, and submatrix、The rotation operator connected.Referred to as direction matrix,It indicates direction vector, is represented by.By invariable rotary matrixEigenvalues Decomposition is carried out, can be obtained special
The rotation operator of value indicative composition, wherein.；
7) rotation operator by acquiringN number of direction parameter is solved, direction parameter substitution is found out into direction matrix A, the side of finding out
To matrix column vector and the maximum column of sampled value Matrix correlation, using the maximum direction parameter of related coefficient as final DOA
The direction of arrival of estimation。
Step 3: direction of arrivalBy matrixIt indicates, wherein phase shift, formula is sampled according to DMWC frequency spectrum perception system signal:(T is former measurement
Matrix), pass through matrixPhase difference is calibrated, to improve signal reconstruction rate.In the case where direction of arrival is unknown, square is measured
Battle array is obtained according to the mixing sequence of frequency mixing stages in DMWC, results in calculation matrix there are deviation so that DMWC frequency spectrum sense
Know that the signal recovery rate of system reduces.
Fig. 3 horizontal axis is signal to noise ratio (snr), and the longitudinal axis is the evaluated error to direction of arrival, it can be seen that with mentioning for signal-to-noise ratio
Height, evaluated error will be greatly reduced.
Fig. 4 is influence of the angle of different direction of arrival to evaluated error, and when direction of arrival is 0 °, evaluated error will reduce,
It is also smaller to the influence for restoring original signal.
It is that direction of arrival compares for 10 °, 30 °, 60 °, 80 ° of lower phase alignments front and backs in tetra- width analogous diagram of Fig. 5.Upper left analogous diagram
Direction of arrival is 10 °, and upper right direction of arrival is 30 °, and lower-left direction of arrival is 60 °, and bottom right direction of arrival is 80 °, and horizontal axis is signal-to-noise ratio
(SNR), the longitudinal axis is signaling protein14-3-3 rate, and the signal (before calibration) of DOA estimation calibration is not done in star-like expression, and triangle expression is done
The signal of DOA estimation calibration (after calibration).With the increase of direction of arrival, the signal of (star-like) restores accuracy rate increasingly before calibration
Low 1 and the signal recovery rate of (triangle) can be improved greatly after calibrating.
Claims (4)
1. a kind of DMWC frequency spectrum perception system phase difference calibration method based on DOA estimation, it is characterised in that the method mistake
Journey are as follows:
Step 1: the DMWC frequency spectrum perception system by being based on uniform linear array (ULA) obtains sampled signal matrix；
Step 2: estimate that obtained direction of arrival determines that each node receives the time delay of signal according to DOA, and then between determining node
Phase difference；
Step 3: calculation matrix is calibrated according to the phase difference that node receives signal.
2. the DMWC frequency spectrum perception system phase difference calibration according to claim 1 based on DOA estimation, it is characterised in that step
In rapid one, DMWC frequency spectrum perception system uses node space random distribution, the structure of partial array homogenous linear.
3. the DMWC frequency spectrum perception system phase difference calibration method according to claim 1 or 2 based on DOA estimation, special
Sign is in step 2, using the rotational invariance between partial array, by calculating direction of arrival based on TLS-ESPRINT algorithm
To estimate phase difference and then carry out phase difference calibration.
4. the DMWC frequency spectrum perception system phase difference calibration method according to claim 1 based on DOA estimation, feature exist
In step 3, direction of arrival is obtained by the DOA estimation method in array signal, and then the phase difference between node is obtained to survey
Moment matrix is calibrated, and this method can increase substantially the recovery rate of sampled signal, improves DMWC frequency spectrum perception system
Noiseproof feature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910310492.7A CN110247866A (en) | 2019-04-17 | 2019-04-17 | A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910310492.7A CN110247866A (en) | 2019-04-17 | 2019-04-17 | A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110247866A true CN110247866A (en) | 2019-09-17 |
Family
ID=67883315
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910310492.7A Pending CN110247866A (en) | 2019-04-17 | 2019-04-17 | A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110247866A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112073130A (en) * | 2020-07-29 | 2020-12-11 | 北京邮电大学 | Frequency spectrum sensing method based on three-point shaping of phase difference distribution curve and related equipment |
| CN112073131A (en) * | 2020-07-29 | 2020-12-11 | 北京邮电大学 | Spectrum sensing method based on phase difference distribution curve analytic expression and related equipment |
| CN112213690A (en) * | 2020-09-29 | 2021-01-12 | 电子科技大学 | Time difference measuring method for phase compression sampling |
| CN112333718A (en) * | 2020-11-05 | 2021-02-05 | 哈尔滨商业大学 | Frequency and arrival angle joint estimation method based on undersampled signals |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106998306A (en) * | 2017-03-16 | 2017-08-01 | 四川大学 | Calculation matrix dynamic calibration method for distributed modulation wide-band transducer |
| CN109347774A (en) * | 2018-09-05 | 2019-02-15 | 四川大学 | The method of compensation is differed in DMWC between multichannel |
-
2019
- 2019-04-17 CN CN201910310492.7A patent/CN110247866A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106998306A (en) * | 2017-03-16 | 2017-08-01 | 四川大学 | Calculation matrix dynamic calibration method for distributed modulation wide-band transducer |
| CN109347774A (en) * | 2018-09-05 | 2019-02-15 | 四川大学 | The method of compensation is differed in DMWC between multichannel |
Non-Patent Citations (2)
| Title |
|---|
| SHAHAR STEIN IOUSHUA 等: "CaSCADE: Compressed Carrier and DOA Estimation", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| TAO CHEN 等: "A ULA-Based MWC Discrete Compressed Sampling Structure for Carrier Frequency and AOA Estimation", 《IEEE ACCESS》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112073130A (en) * | 2020-07-29 | 2020-12-11 | 北京邮电大学 | Frequency spectrum sensing method based on three-point shaping of phase difference distribution curve and related equipment |
| CN112073131A (en) * | 2020-07-29 | 2020-12-11 | 北京邮电大学 | Spectrum sensing method based on phase difference distribution curve analytic expression and related equipment |
| CN112213690A (en) * | 2020-09-29 | 2021-01-12 | 电子科技大学 | Time difference measuring method for phase compression sampling |
| CN112213690B (en) * | 2020-09-29 | 2022-05-24 | 电子科技大学 | Time difference measuring method for phase compression sampling |
| CN112333718A (en) * | 2020-11-05 | 2021-02-05 | 哈尔滨商业大学 | Frequency and arrival angle joint estimation method based on undersampled signals |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110247866A (en) | A kind of DMWC frequency spectrum perception phase alignment based on DOA estimation | |
| CN103596266B (en) | A kind of human testing and method, the apparatus and system of positioning | |
| Stoica et al. | SPICE and LIKES: Two hyperparameter-free methods for sparse-parameter estimation | |
| CN106019214B (en) | Wide-band coherent signal source DOA estimation method | |
| CN107171749A (en) | The method for determining the Doppler frequency shift for the wireless signal that moving object is directly reflected | |
| CN106802403B (en) | Acoustic vector sensors two-dimensional array MUSIC decorrelation LMS method for parameter estimation | |
| CN110412557B (en) | Method and device for measuring speed and distance based on OFDM signal | |
| CN103063300B (en) | Micro-polarization modulation array for achieving full-polarization imaging | |
| CN108955734A (en) | A kind of cavity length demodulating method of fiber F-P temperature/pressure compound sensor | |
| CN102353393A (en) | Quadrature demodulation device for interference type photo-sensor based on pi/2 phase modulation | |
| CN109633538A (en) | The maximum likelihood time difference estimation method of nonuniform sampling system | |
| CN104360305A (en) | Radiation source direction finding positioning method of uniting compressed sensing and signal cycle stationary characteristics | |
| CN109150235A (en) | Compressed sensing based multicycle direct expansion msk signal two dimension joint acquisition method | |
| CN116008976B (en) | Fusion ranging method for CHIRP signal | |
| CN107483130A (en) | One kind joint broader frequency spectrum perceives and angle-of- arrival estimation method | |
| CN112444830A (en) | Compressed sensing GNSS signal capturing method based on singular value decomposition | |
| CN113286362A (en) | Method for acquiring arrival time and arrival angle of multipath signal and related device | |
| CN102323482B (en) | Method for measuring phase frequency characteristic by using digital intermediate-frequency spectrum analyzer during network analysis and measurement | |
| CN105490682B (en) | Compression transducing signal recovery algorithms based on SVD and OMP | |
| CN102495280A (en) | Anti-noise wide-range frequency measurement method and phase locking frequency meter | |
| CN103529429A (en) | Pulse compression method for field of linear frequency modulation signals | |
| CN106533528B (en) | The two-way time compares compression and the calculation method of modem transmitted data amount | |
| CN107957579B (en) | Method and system for determining state parameters of measured object | |
| CN106154220B (en) | L-type simplifies acoustic vector-sensor array column multi-parameter Combined estimator quaternary counting method | |
| CN106998306A (en) | Calculation matrix dynamic calibration method for distributed modulation wide-band transducer |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190917 |
|
| WD01 | Invention patent application deemed withdrawn after publication |