CN105721075B - Data are aided in frequency reuse satellite interference signal intensity measuring method and device - Google Patents
Data are aided in frequency reuse satellite interference signal intensity measuring method and device Download PDFInfo
- Publication number
- CN105721075B CN105721075B CN201610083961.2A CN201610083961A CN105721075B CN 105721075 B CN105721075 B CN 105721075B CN 201610083961 A CN201610083961 A CN 201610083961A CN 105721075 B CN105721075 B CN 105721075B
- Authority
- CN
- China
- Prior art keywords
- signal
- noise
- receiving sequence
- power
- sequence
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/345—Interference values
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Radio Relay Systems (AREA)
Abstract
The present invention discloses a kind of data auxiliary with frequency reuse satellite interference signal intensity measuring method and device.Method includes:Matched filtering is carried out to receiving sequence, makes filter output signal signal to noise ratio maximum;According to useful signal in receiving sequence and the separate characteristic of noise, the joint probability density function of whole receiving sequence is constructed;Local derviation is asked to signal power and noise power using log-likelihood function, obtains signal power estimate and noise power estimation value;Auto-correlation computation is carried out to signal, obtains the signal power estimate and noise power estimation value represented with auto-correlation amount;Signal-to-Noise information is obtained using signal and noise power estimation value.Device includes matched filtering module, joint probability density function constructing module, power estimation module, signal autocorrelation computing module and interference signal intensity computing module.The present invention is simple and easy to interference signal intensity measurement in the case where there is offset frequency situation, accuracy rate is high.
Description
Technical field
The invention belongs to satellite communication digital signal processing technique field, and particularly one kind is in the presence of having frequency deviation
It is simple to the algorithm for estimating of interference signal intensity, be easily achieved, accuracy rate is high is done based on data-aided same frequency reuse satellite
Disturb signal strength measurement method and device.
Background technology
In communication environment now, effectively measuring interference signal intensity, have in civil and military field very wide
General application value.By taking satellite communication as an example, in civil area, interference strength measure for monitor in satellite network user whether work
Make in legal running parameter, while detect disabled user's occupying and usurping to satellite repeater resource;In its military neck
Domain, accurate interference strength measurement are the preconditions realized space electromagnetic countermeasure and obtain information superiority.
Interference strength measurement is premised on estimating signal-tonoise information, and the relevant signal to noise ratio delivered both at home and abroad at present is estimated
In the document of meter, method of estimation can substantially be divided into the method for estimation in the case of non-data aided and have data-aided signal to noise ratio
Method of estimation.The signal-noise ratio estimation method of non-data aided is that signal is estimated in the case of no prior information, real
Existing complexity is high, and Feasible GLS estimation scope is smaller;There is data-aided signal-noise ratio estimation method to utilize obtained priori
Information estimates that signal to noise ratio and interference signal intensity the degree of accuracy is high, and estimation range is wide.
In the satellite communication of signal and noise with frequency multiplexing, signal is molded through ovennodulation, raised cosine, then mutual with noise
Aliasing, there is certain frequency error in receiving terminal reception signal -- frequency deviation.Under this condition, existing interference signal intensity measurement side
Method comes with some shortcomings:1st, can not accurately estimate under raised cosine condition of molding;2nd, require that reception signal has perfect synchronization,
Frequency error can not reach good estimation effect in the case of existing;3rd, estimation range is not broad enough, it is difficult to meets that small-signal is estimated
Meter requires.
The content of the invention
It is an object of the invention to provide a kind of data auxiliary with frequency reuse satellite interference signal intensity measuring method, having
It is simple to the algorithm for estimating of interference signal intensity in the presence of frequency deviation, be easily achieved, accuracy rate it is high.
Another object of the present invention is to provide a kind of data auxiliary with frequency reuse satellite interference signal intensity measurement apparatus.
The technical solution for realizing the object of the invention is:A kind of data auxiliary is surveyed with frequency reuse satellite interference signal intensity
Amount method, comprises the following steps:
(10) matched filtering:Matched filtering is carried out to receiving sequence, makes filter output signal signal to noise ratio maximum;
(20) tectonic syntaxis probability density function:According to useful signal in receiving sequence and the separate characteristic of noise,
Construct the joint probability density function of whole receiving sequence;
(30) power is estimated:According to receiving sequence joint probability density function, log-likelihood function is asked for, and utilize logarithm
Likelihood function seeks local derviation to signal power and noise power, obtains signal power estimate and noise power estimation value;
(40) signal autocorrelation computing:According to the characteristic that signal in receiving sequence is related, noise is unrelated, signal is carried out certainly
Related operation, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
(50) interference signal intensity calculates:Signal-to-Noise information is obtained using signal and noise power estimation value.
The technical solution for realizing another object of the present invention is:A kind of data auxiliary is strong with frequency reuse satellite interference signal
Measurement apparatus is spent, it includes:
Matched filtering module, is arranged on receiving terminal, for carrying out matched filtering to receiving sequence, makes filter output signal
Signal to noise ratio is maximum;
Joint probability density function constructing module, for according to useful signal in receiving sequence and the separate spy of noise
Property, construct the joint probability density function of whole receiving sequence;
Power estimation module, for according to receiving sequence joint probability density function, asking for log-likelihood function, and utilize
Log-likelihood function seeks local derviation to signal power and noise power, obtains signal power estimate and noise power estimation value;
Signal autocorrelation computing module, for according to the characteristic that signal in receiving sequence is related, noise is unrelated, entering to signal
Row auto-correlation computation, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
Interference signal intensity computing module, for obtaining Signal-to-Noise information using signal and noise power estimation value.
The present invention compared with prior art, its remarkable advantage:
1st, algorithm is simple, is easily achieved:The present invention overcomes signal frequency error using signal autocorrelation, it is allowed to receives
There is certain frequency error in signal, can still reach good estimation effect in the presence of error, algorithm is simple,
It is easily achieved;
2nd, accuracy rate is high:This algorithm estimation range can be extended to -12dB, still there is preferably estimation effect for small-signal
Fruit, robustness are good;
3rd, it is adapted to satellite communication channel:Satellite communication channel is awgn channel, signal shaping generally using raised cosine into
Shape, the present invention have taken into full account that raised cosine shapes influence to signal characteristic, the characteristics of being adapted to satellite communication channel;
The present invention is done with reference to the accompanying drawings and detailed description and further illustrated.
Brief description of the drawings
Fig. 1 is the main flow chart of the invention based on data-aided same frequency reuse satellite interference signal intensity measuring method.
Fig. 2 is that interference signal intensity measures signal-to-noise ratio (SNR) estimation Performance Simulation Results figure
Fig. 3 is that interference signal intensity measures signal-to-noise ratio (SNR) estimation NMSE simulation result figures.
Embodiment
As shown in figure 1, data auxiliary of the present invention, with frequency reuse satellite interference signal intensity measuring method, it includes following step
Suddenly:
(10) matched filtering:Receiving terminal sets matched filter, carries out matched filtering to receiving sequence, exports wave filter
Signal-to-Noise is maximum;
(20) tectonic syntaxis probability density function:According to useful signal in receiving sequence and the separate characteristic of noise,
Construct the joint probability density function of whole receiving sequence;
Tectonic syntaxis probability density function (20) step is specially:
It is r={ r by sampled sequence after matched filtering0,r1,....rk, the joint probability density function of receiving sequence
For,
Wherein,
In formula, S represents signal power, and N represents noise power,The in-phase component of receiving sequence is represented,The quadrature component of receiving sequence is represented, k=1,2 ... K represent k-th of sample value in sequence, and K grows to receive data
Degree, rkRepresent k-th of sample value in sequence.
(30) power is estimated:According to receiving sequence joint probability density function, log-likelihood function is asked for, and utilize logarithm
Likelihood function seeks local derviation to signal power and noise power, obtains signal power estimate and noise power estimation value;The work(
Rate estimation (30) step includes:
(31) log-likelihood function is asked for:Log-likelihood function is asked for as the following formula,
(32) signal power is estimated:Utilize log-likelihood function pairWithSeek local derviation so thatWhen the logarithm
The partial derivative of likelihood function is equal to 0, obtains signal power estimate,
(33) noise power estimation:Utilize log-likelihood function pairWithSeek local derviation so thatWhen this is right
The partial derivative of number likelihood function is equal to 0, obtains noise power estimation value
In formula,For reception signal general power estimate,WithThe estimation to signal and noise power is represented respectively
Value, K represent receiving sequence length, rkK-th of sample value in sequence is represented, | | the absolute value of signal is represented, * is then represented to signal
Take complex conjugate.
(40) signal autocorrelation computing:According to the characteristic that signal in receiving sequence is related, noise is unrelated, signal is carried out certainly
Related operation, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
(40) the signal autocorrelation calculation step is specially:
Receiving sequence general power estimate is,
Signal power estimate is,
Wherein, RcFor sequence auto-correlation amount after filtering,
HereinRepresent receiving sequence, wherein akIt is independent identically distributed MPSK constellations letter
Number, nkRepresent the white complex gaussian noise signal with zero mean unit variance after sampling, θk=2 π kv+ θ0, k=1,2 ... K,
V is normalized frequency difference, θ0It is phase compensation constant.
(50) interference signal intensity calculates:Signal-to-Noise information is obtained using signal and noise power estimation value.
(50) the interference signal intensity calculation procedure is specially:
Interference signal intensity is calculated using signal power estimate in interference signal and noise power estimation value,
The present invention takes into full account satellite channel feature, by improved based on most for this application scenarios of satellite communication
The method of maximum-likelihood estimation, realize the interference signal intensity estimation of relative broad range.Simultaneously there is frequency mistake in the algorithm in signal
Still good estimation performance can be realized in the case of difference.
Fig. 2, Fig. 3 give is multiplexed interference signal intensity measurement simulation result using based on data-aided with frequency.Emulation
Using QPSK signals as test signal, symbol lengths 32, character rate 200K, sampling rate 5M, signal is using liter
Cosine molding filtration, rolloff-factor 0.5, normalized frequency difference v=0.01 is set, makes phase compensation constant θ0=π/12.Test
SNR ranges are -30~10dB.Emulated 100 times with the conditions of differing without frequency difference.The simulating, verifying feasibility of this method.
Realize that the data of the inventive method are aided in frequency reuse satellite interference signal intensity measurement apparatus, it includes:
Matched filtering module, is arranged on receiving terminal, for carrying out matched filtering to receiving sequence, makes filter output signal
Signal to noise ratio is maximum;
Joint probability density function constructing module, for according to useful signal in receiving sequence and the separate spy of noise
Property, construct the joint probability density function of whole receiving sequence;
Power estimation module, for according to receiving sequence joint probability density function, asking for log-likelihood function, and utilize
Log-likelihood function seeks local derviation to signal power and noise power, obtains signal power estimate and noise power estimation value;
Signal autocorrelation computing module, for according to the characteristic that signal in receiving sequence is related, noise is unrelated, entering to signal
Row auto-correlation computation, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
Interference signal intensity computing module, for obtaining Signal-to-Noise information using signal and noise power estimation value.
Claims (5)
1. a kind of data auxiliary is the same as frequency reuse satellite interference signal intensity measuring method, it is characterised in that comprises the following steps:
(10) matched filtering:Matched filtering is carried out to receiving sequence, makes filter output signal signal to noise ratio maximum;
(20) tectonic syntaxis probability density function:According to useful signal in receiving sequence and the separate characteristic of noise, construction
The joint probability density function of whole receiving sequence;
(30) power is estimated:According to receiving sequence joint probability density function, log-likelihood function is asked for, and utilize log-likelihood
Function pair signal power seeks local derviation with noise power, obtains signal power estimate and noise power estimation value;
(40) signal autocorrelation computing:According to the characteristic that signal in receiving sequence is related, noise is unrelated, auto-correlation is carried out to signal
Computing, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
(50) interference signal intensity calculates:Signal-to-Noise information is obtained using signal and noise power estimation value
Tectonic syntaxis probability density function (20) step is specially:
It is r={ r by sampled sequence after matched filtering0,r1,....rk, the joint probability density function of receiving sequence is,
Wherein,
In formula, S represents signal power, and N represents noise power,The in-phase component of receiving sequence is represented,The quadrature component of receiving sequence is represented, k=1,2 ... K represent k-th of sample value in sequence, and K grows to receive data
Degree, rkRepresent k-th of sample value in sequence.
2. interference signal intensity measuring method according to claim 1, it is characterised in that the power estimates (30) step
Including:
(31) log-likelihood function is asked for:Log-likelihood function is asked for as the following formula,
(32) signal power is estimated:Utilize log-likelihood function pairWithSeek local derviation so thatWhen the log-likelihood
The partial derivative of function is equal to 0, obtains signal power estimate,
(33) noise power estimation:Utilize log-likelihood function pairWithSeek local derviation so thatWhen the log-likelihood
The partial derivative of function is equal to 0, obtains noise power estimation value
In formula,For reception signal general power estimate,WithThe estimate to signal and noise power, K are represented respectively
Represent receiving sequence length, rkK-th of sample value in sequence is represented, | | the absolute value of signal is represented, * then represents to take signal
Complex conjugate, akIt is independent identically distributed MPSK constellation signals.
3. interference signal intensity measuring method according to claim 2, it is characterised in that (40) the signal autocorrelation fortune
Calculating step is specially:
Receiving sequence general power estimate is,
Signal power estimate is,
Wherein, RcFor sequence auto-correlation amount after filtering,
HereinRepresent receiving sequence, wherein akIt is independent identically distributed MPSK constellation signals,
nkRepresent the white complex gaussian noise signal with zero mean unit variance after sampling, θk=2 π kv+ θ0, k=1,2 ... K, v are
Normalized frequency difference, θ0It is phase compensation constant.
4. interference signal intensity measuring method according to claim 3, it is characterised in that (50) interference signal intensity
Calculation procedure is specially:
Interference signal intensity is calculated using signal power estimate in interference signal and noise power estimation value,
5. a kind of data auxiliary is the same as frequency reuse satellite interference signal intensity measurement apparatus, it is characterised in that including:
Matched filtering module, is arranged on receiving terminal, for carrying out matched filtering to receiving sequence, makes filter output signal noise
Than maximum;
Joint probability density function constructing module, for according to useful signal in receiving sequence and the separate characteristic of noise,
Construct the joint probability density function of whole receiving sequence;
Power estimation module, for according to receiving sequence joint probability density function, asking for log-likelihood function, and utilize logarithm
Likelihood function seeks local derviation to signal power and noise power, obtains signal power estimate and noise power estimation value;
Signal autocorrelation computing module, for according to the characteristic that signal in receiving sequence is related, noise is unrelated, being carried out certainly to signal
Related operation, obtain the signal power estimate and noise power estimation value represented with auto-correlation amount;
Interference signal intensity computing module, for obtaining Signal-to-Noise information using signal and noise power estimation value;
The tectonic syntaxis probability density function is specially:
It is r={ r by sampled sequence after matched filtering0,r1,....rk, the joint probability density function of receiving sequence is,
Wherein,
In formula, S represents signal power, and N represents noise power,The in-phase component of receiving sequence is represented,The quadrature component of receiving sequence is represented, k=1,2 ... K represent k-th of sample value in sequence, and K grows to receive data
Degree, rkRepresent k-th of sample value in sequence.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610083961.2A CN105721075B (en) | 2016-02-06 | 2016-02-06 | Data are aided in frequency reuse satellite interference signal intensity measuring method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610083961.2A CN105721075B (en) | 2016-02-06 | 2016-02-06 | Data are aided in frequency reuse satellite interference signal intensity measuring method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105721075A CN105721075A (en) | 2016-06-29 |
CN105721075B true CN105721075B (en) | 2017-12-08 |
Family
ID=56155895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610083961.2A Active CN105721075B (en) | 2016-02-06 | 2016-02-06 | Data are aided in frequency reuse satellite interference signal intensity measuring method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105721075B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106330362A (en) * | 2016-08-25 | 2017-01-11 | 中国电子科技集团公司第十研究所 | Data assisted signal to noise ratio estimation method |
CN112152739B (en) * | 2020-09-24 | 2021-05-28 | 清华大学 | Method and device for calculating interference probability distribution of satellite constellation |
CN112904374B (en) * | 2021-01-29 | 2024-03-19 | 湖南国科微电子股份有限公司 | Satellite signal strength evaluation method and device, GNSS receiver and medium |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1809046A (en) * | 2006-02-14 | 2006-07-26 | 郭更生 | Method and apparatus of implementing multi-symbol OFDM synchronization |
CN100553166C (en) * | 2006-05-11 | 2009-10-21 | 上海交通大学 | Channel estimation methods based on particle filter |
CN102075476A (en) * | 2009-11-25 | 2011-05-25 | 电子科技大学 | Iteration carrier synchronization method used under condition of extremely low signal to noise ratio |
CN102983922B (en) * | 2012-11-29 | 2014-12-24 | 华为技术有限公司 | Method for measuring co-channel signal interference and receiver |
CN103869311B (en) * | 2014-03-18 | 2016-05-25 | 电子科技大学 | Real beam scanning radar super-resolution imaging method |
-
2016
- 2016-02-06 CN CN201610083961.2A patent/CN105721075B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105721075A (en) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108199994B (en) | OFDM signal frequency offset estimation method and system for wireless comprehensive measurement instrument | |
CN101087281B (en) | A measurement method and device of orthogonal range modulation N/S ratio and N/R ratio | |
CN107911329B (en) | OFDM signal demodulation method of signal analyzer | |
CN112235215B (en) | Wireless channel detection method, storage medium and terminal equipment | |
CN105721075B (en) | Data are aided in frequency reuse satellite interference signal intensity measuring method and device | |
CN107809398B (en) | MSK signal modulation parameter estimation method and communication system under impulse noise environment | |
CN108183841B (en) | Base band data processing method and system based on IEEE802.11ah in comprehensive tester | |
Li et al. | Noncircular measurement and mitigation of $ I/Q $ imbalance for OFDM-based WLAN transmitters | |
CN113037663B (en) | Improved code element rate estimation method suitable for non-constant envelope signal | |
CN107342960B (en) | Non-data-aided frequency offset estimation method suitable for amplitude phase shift keying | |
CN104363194A (en) | PSK (phase shift keying) modulation recognition method based on wave form transformation | |
CN107911323A (en) | The joint frequency phase method of estimation closed based on partly decoupled | |
WO2009007024A1 (en) | Robust channel estimation for wireless system | |
TW201406097A (en) | Method and device for estimating channel and noise | |
CN109583575B (en) | Processing method for improving instrument vector signal analysis performance based on deep learning | |
CN103188067B (en) | A kind of chip clock frequency departure estimation error of spread spectrum system and the method for correction | |
CN104901907B (en) | Data-aided steady signal-noise ratio estimation method is based under a kind of dynamic environment | |
CN111343113B (en) | Phase synchronization improvement method and device based on digital oscilloscope | |
CN108900445A (en) | A kind of method and device of signal code rate estimation | |
CN113078966A (en) | High-precision 5G large-bandwidth signal testing method | |
CN102307166B (en) | SNR (signal to noise ratio) estimation method | |
CN103516642B (en) | The method and apparatus of Combined estimator interference signal physical parameter | |
CN104270328B (en) | A kind of signal to noise ratio real-time estimation method | |
CN109257311B (en) | Method and system for determining error vector magnitude | |
CN113612711B (en) | Frequency offset estimation method for short burst modulation signal under low signal-to-noise ratio |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |