CN102185666A - System and method for predicting multimode combined interference in complex electromagnetic environment - Google Patents

System and method for predicting multimode combined interference in complex electromagnetic environment Download PDF

Info

Publication number
CN102185666A
CN102185666A CN2011100968822A CN201110096882A CN102185666A CN 102185666 A CN102185666 A CN 102185666A CN 2011100968822 A CN2011100968822 A CN 2011100968822A CN 201110096882 A CN201110096882 A CN 201110096882A CN 102185666 A CN102185666 A CN 102185666A
Authority
CN
China
Prior art keywords
execution
interference source
interference
unit
electromagnetic environment
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
CN2011100968822A
Other languages
Chinese (zh)
Other versions
CN102185666B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN2011100968822A priority Critical patent/CN102185666B/en
Publication of CN102185666A publication Critical patent/CN102185666A/en
Application granted granted Critical
Publication of CN102185666B publication Critical patent/CN102185666B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to electromagnetic interference prediction, in particular to a system and a method for predicting multimode combined interference in a complex electromagnetic environment. The system and the method are characterized in that: multimode combined interference consisting of an interference source outside the system is analyzed in the complex electromagnetic environment by combining working radio station nonlinear receiver modeling and discrete nonlinear analysis technology by a transient-steady state equivalent method according to a multimode combined interference prediction model in the complex electromagnetic environment, and the change, caused by the multimode combined interference, of receiver output signals is evaluated by taking a bit error rate or signal-to-noise-and-distortion (SINAD) as a criterion, thereby effectively calculating the electromagnetic environmental effects of the system. The method for predicting the multimode combined interference in the complex electromagnetic environment is applied to digital communication in the complex electromagnetic environment.

Description

Multi-mode combination interference prognoses system and method under the complex electromagnetic environment
Technical field
The present invention relates to electromagnetic interference prediction, particularly multi-mode combination interference prognoses system and the method under the complex electromagnetic environment.
Background technology
In complicated and modern battlefield electromagnetic environment that state is changeable, exist a large amount of communication equipments, battlefield surveillance equipment, early warning/target-detection unit, guidance radar, weapon control system, range measurement/navigator, countermeasures set etc., they have constituted very high electromagnetic radiation source density, covered from long wave, medium wave until microwave, millimeter wave until the wave band of ultraviolet band, power is from less than 1 watt to several kilowatts even higher.So under the complex electromagnetic environments, the multi-mode combination interference that various interference source constituted causes Vehicular system to be faced with unprecedented electromagnetic interference and threatens.Multi-mode combination interference under the complex electromagnetic environment is carried out accurately quantitative prediction, is prevention and the important prerequisite that suppresses to disturb, guarantee system reliability and fail safe.
For the prediction of outside electromagnetic interference, existing Forecasting Methodology mostly is to adopt the linear discrete analytical technology, disturb at the stable state of single-mode and predict, and prediction only limits to prevention at radio-frequency port, is based on the Forecasting Methodology of energy.Obviously, numerous in interference source quantity and pattern, transient state, stable state interference source mix and exist, a large amount of modern battlefields that adopt of modes such as error correction, spread spectrum, and traditional Forecasting Methodology can not be suitable for.
Summary of the invention
It is numerous to the purpose of this invention is to provide a kind of adaptation interference source quantity and pattern, and transient state, stable state interference source mix and exist, multi-mode combination interference prognoses system and method under the complex electromagnetic environment that a large amount of modern battlefields that adopt of modes such as error correction, spread spectrum use.
The object of the present invention is achieved like this, multi-mode combination interference prognoses system under the complex electromagnetic environment, it is characterized in that: it is characterized in that: which comprises at least: work radio station unit transmitter, work unit, radio station " non-linear " receiver, multi-mode combination interference source under the complex electromagnetic environment, the Channel Elements that comprises background noise, the SINAD unit, the BER unit; The information source of work radio station unit transmitter is by output signal S1 after the source encoding, and S1 channel coding, modulating unit are handled the back and sent radiofrequency signal with certain power, enter a channel input of the Channel Elements that comprises background noise N; Simultaneously, the output in the multi-mode combination interference source under the complex electromagnetic environment disturbs D to enter the Channel Elements that comprises background noise N, background noise stack in two-way input signal and the Channel Elements, produce the output signal O2 of Channel Elements, radio frequency input as work unit, radio station " non-linear " receiver, the output signal O2 of Channel Elements enters unit, work radio station " non-linear " receiver after after the nonlinear element, channel decoding output signal be designated as O1, O1 enters output output at last through source decoding; Signal S2 after the decay of unit transmitter channel unit, described work radio station enters the input in SINAD unit, and the output signal O2 of Channel Elements enters another input of SINAD unit, by the two sides road calculated signals SINAD that enters the SINAD unit; And the output signal O1 of the channel decoding in the output signal S1 of the source encoding 5 in the work radio station unit transmitter and work unit, radio station " non-linear " receiver exports along separate routes and enters two inputs in BER unit, calculates BER by the BER unit.
Background noise in the described Channel Elements is used to describe different frequency range, the electromagnetic environment of different regions, utilize the random number sequence meet Gaussian Profile to come the simulation background noise in the unit, obtain the reasonable value of noise power by the variance of adjusting random number sequence, implementation procedure is as follows:
1) utilize a stationary random process to represent background noise
Feature according to background noise is established as the stationary random process ξ that auto-correlation function is R (τ) (t), and the power of actual noise generally all is meant its AC power, and solution procedure is as follows:
Noise average power is: R (0)=E[ξ 2(t)], E[x] expression x mathematic expectaion;
The noise direct current power is: R (∞)=E 2[ξ (t)]
So its AC power is: R (0)-R (∞)=E[ξ 2(t)]-E 2[ξ (t)]=σ 2, σ 2The variance of expression ξ (t), promptly the power of noise is represented with the variance of its process.
2) utilize the random number sequence that meets Gaussian Profile to realize this stationary random process
According in the test or the typical environment level value that provides in the standard, by adjusting the variances sigma of this random number sequence 2Obtain rational background noise, different σ 2Value characterizes different noises.
The output in described multi-mode combination interference source 3 disturbs D to be meant at the external interference source under the complex electromagnetic environment, determines its amplitude in frequency domain and time domain, the associate feature of power.
Described nonlinear element 9 refers to the non-linear partial of unit, work radio station " non-linear " receiver, and nonlinear element 9 is by the serial or parallel connection of linear filter 201,203 and memoryless non-linear components 202 or connection in series-parallel combination.
Described linear filter is meant dynamic linear with memory filter, with the systematic sampling time be the interval, the output signal of its current state is not only relevant with the input signal of current state system, and with the input signal of m state before and the output signal of n state is relevant before, the input signal of m state and the linear superposition of the output signal of n state before before being, if current state is k, be input as x (k), then the output signal y of current linear filter (k) is:
y ( k ) = Σ i = 1 n a i y ( k - i ) + Σ j = 0 m b j x ( k - j )
a i, b jParameter decision by linear filter.
Described memoryless non-linear components is meant that its output signal is irrelevant with the output signal of each state before, and only is the q order polynomial of current input state, and establish the current x of being input as (t), then current output signal y (t) can be expressed as:
y ( t ) = Σ p = 0 q a p x q ( t )
P wherein, q is by the parameter decision of memoryless nonlinear device.
It is according to probability statistics principle and signal detection means that D is disturbed in the output in described multi-mode combination interference source 3, set up the probability statistics unit of typical interference source, and according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source, with the D that is expressed as of its multi-mode combination back quantification.
The computing formula of SINAD (letter Nader) is:
SINAD = 10 lg O 2 N + D = 10 lg S 2 + N + D N + D
Described BER (error rate) computing formula is:
Figure BDA0000055939970000034
Based on the multi-mode combination interference Forecasting Methodology under the complex electromagnetic environment of the multi-mode combination interference prognoses system under the complex electromagnetic environment: it is characterized in that: which comprises at least following steps:
Execution in step 301 is carried out the multi-mode combination interference analysis under the complex electromagnetic environment;
Execution in step 302, the interference source that carries out under the complex electromagnetic environment obtains;
At first execution in step 304 is carried out the selection of interference source obtain manner, if obtain manner is a mathematical modeling, then execution in step 305 is carried out the mathematical modeling of interference source, if the interference source obtain manner detects for monitoring, then execution in step 311 utilizes measuring receiver to carry out the reception of interference source;
Step 305 is carried out the mathematical modeling of interference source;
Execution in step 306 is carried out the judgement of interference source type, if interference source is a transient state, then execution in step 307; If interference source is a stable state, then direct execution in step 308;
Interference source is a transient state, and execution in step 307 utilizes small wave converting method to carry out wink-steady-state equivalent;
Execution in step 308 is carried out the analytical calculation of equivalent stable state interference source;
Whether execution in step 309 is carried out interference source and is calculated and finish judgement, does not finish if calculate, and then repeated execution of steps 302, finish execution in step 310 if calculate;
Execution in step 310, output multi-mode combination interference source is to the radio station of working " non-linear " receiver mathematical units;
If the interference source obtain manner detects for monitoring, then execution in step 311, utilize measuring receiver to carry out the reception of interference source;
Execution in step 312, the signature analysis of the interference source that receives;
Execution in step 313 is carried out interference source type and is judged that if interference source is a transient state, then execution in step 307, if interference source is a stable state, then execution in step 308.
Execution in step 303 is carried out the receiver modeling under the complex electromagnetic environment;
Execution in step 314 is carried out the analysis of work radio station receiver nonlinear element;
Execution in step 315 then, carry out the extraction and the simplification of " non-linear " receiver mathematical units;
Execution in step 316 obtains work radio station " non-linear " receiver mathematical units, receives the multi-mode combination interference from step 310;
Execution in step 317 is obtained Discrete Nonlinear and is carried out result calculated;
Execution in step 318 calculates the output signal of " non-linear " receiver;
Execution in step 319 is carried out the analysis of signal performance;
Execution in step 320 is carried out SINAD and is calculated;
Execution in step 321 compares SINAD and the SINAD tolerance limit D0 that calculates, if SINAD more than or equal to D0 then do not disturb, prediction finishes, if SINAD is less than D0 then execution in step 322 is calculated BER;
Execution in step 322 is calculated BER;
Execution in step 323 compares BER and the BER tolerance value C0 that calculates, if BER then do not disturb smaller or equal to C0, prediction finishes, if BER greater than C0 then disturb, adjusts work radio station running parameter, repeated execution of steps 303 is predicted.
Described wink-steady-state equivalent 307, be according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source.
Described Discrete Nonlinear technology 317 is meant and utilizes discrete unit that it is simplified description at all kinds of interference sources under the complex electromagnetic environment at " non-linear " receiver front end; The linear filter of non-linear receiver in the frequency domain modeling, in the time domain modeling, is carried out time domain-frequency domain transform with FFT with its memoryless nonlinear device; Determine non-linear interference source with the multi-branched search utility, care signal is assessed through the performance behind the receiver unit.
Advantage of the present invention is: it is according to multi-mode combination interference prognoses system under the complex electromagnetic environment, utilize the wink-the steady-state equivalent analytical method, in conjunction with work radio station " non-linear " receiver modeling and Discrete Nonlinear analytical technology, the multi-mode combination interference that the interference source of system outside in the complex electromagnetic environment constitutes is analyzed, and to its receiver variation of output signals that causes, take SINAD and BER layering forecast method to carry out the electromagnetic interference prediction, thus the electromagnetic environmental effects of effective computing system.
Description of drawings
The invention will be further described below in conjunction with the embodiment accompanying drawing:
Fig. 1 is a multi-mode combination interference predicting unit under the complex electromagnetic environment;
Fig. 2 is the nonlinear element unit of " non-linear " receiver;
Fig. 3 is multi-mode combination interference analysis process figure under the complex electromagnetic environment:
Fig. 4 is the test legend of background noise unit;
Fig. 5 is the emulation legend of background noise unit.
Among the figure: 1, work radio station unit transmitter; 2, work unit, radio station " non-linear " receiver; 3, the multi-mode combination interference source under the complex electromagnetic environment; 4, information source; 5, source encoding; 6, chnnel coding; 7, modulating unit; 8, the Channel Elements of background noise; 9, nonlinear element; 10, channel decoding; 11, source decoding; 12, output; 13, SINAD unit; 14, BER unit.
Embodiment
As shown in Figure 1, multi-mode combination interference prognoses system under the complex electromagnetic environment, comprise: work radio station unit transmitter 1, work unit, radio station " non-linear " receiver 2, multi-mode combination interference source 3 under the complex electromagnetic environment, the Channel Elements 8 that comprises background noise, SINAD unit 13, BER unit 14; The information source 4 of work radio station unit transmitter 1 is by source encoding 5 back output signal S1, and S1 channel coding 6, modulating unit 7 are handled the back and sent radiofrequency signal with certain power, enter a channel input of the Channel Elements 8 that comprises background noise N; Simultaneously, the output in the multi-mode combination interference source 3 under the complex electromagnetic environment disturbs D to enter the Channel Elements 8 that comprises background noise N, background noise stack in two-way input signal and the Channel Elements, produce the output signal O2 of Channel Elements 8, radio frequency input as work unit, radio station " non-linear " receiver 2, after the output signal O2 of Channel Elements 8 enters unit, work radio station " non-linear " receiver 2, output signal is designated as O1 after nonlinear element 9, channel decoding 10, and O1 enters output 12 outputs at last through source decoding 11.
Signal S2 after 1 channel unit, 8 decay of work radio station unit transmitter enters the input in SINAD unit, and the output signal O2 of Channel Elements 8 enters another input of SINAD unit, by the two sides road calculated signals SINAD that enters the SINAD unit.
And the output signal O1 of the channel decoding 10 in the output signal S1 of the source encoding 5 in the work radio station unit transmitter 1 and work unit, radio station " non-linear " receiver 2 exports along separate routes and enters 14 two inputs in BER unit, calculates BER by BER unit 14.
As shown in Figure 4 and Figure 5, the background noise N of the Channel Elements 8 of described background noise is used to describe different frequency range, the electromagnetic environment of different regions.Utilize the random number sequence meet Gaussian Profile to come the simulation background noise in the Channel Elements 8 of background noise, obtain the reasonable value of noise power by the variance of adjusting random number sequence, implementation procedure is as follows:
1) utilize a stationary random process to represent background noise
Feature according to background noise is established as the stationary random process ξ that auto-correlation function is R (τ) (t), and the power of actual noise generally all is meant its AC power, and solution procedure is as follows:
Noise average power is: R (0)=D[ξ 2(t)], D[x] expression x mathematic expectaion.
The noise direct current power is: R (∞)=E 2[ξ (t)]
So its AC power is: R (0)-R (∞)=E[ξ 2(t)]-E 2[ξ (t)]=σ 2, σ 2The variance of expression ξ (t), promptly the power of noise is represented with the variance of its process.
2) utilize the random number sequence that meets Gaussian Profile to realize this stationary random process
According to the typical environment level value that provides in (Fig. 5) or the standard in the test, by adjusting the variances sigma of this random number sequence 2Obtain rational background noise, different σ 2Value characterizes different noises (Fig. 6).
The output in described multi-mode combination interference source 3 disturbs D to be meant at the external interference source under the complex electromagnetic environment, determines its amplitude in frequency domain and time domain, the associate feature of power.It is according to probability statistics principle and signal detection means that D is disturbed in the output in described multi-mode combination interference source 3, set up the probability statistics unit of typical interference source, and according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source, with the D that is expressed as of its multi-mode combination back quantification.
Described nonlinear element 9 refers to the non-linear partial of unit, work radio station " non-linear " receiver, and nonlinear element 9 is by the serial or parallel connection of linear filter 201,203 and memoryless non-linear components 202 or connection in series-parallel combination." non-linear " parameter can or be measured from the receiver performance index and obtain, to after the test data statistical analysis parameter default being set.
Described linear filter 201 is meant dynamic linear with memory filter, with the systematic sampling time be the interval, the output signal of its current state is not only relevant with the input signal of current state system, and with the input signal of m state before and the output signal of n state is relevant before, be before the input signal and the linear superposition of the output signal of n state before of m state.If current state is k, be input as x (k), then the output signal y of current linear filter (k) is:
y ( k ) = Σ i = 1 n a i y ( k - i ) + Σ j = 0 m b j x ( k - j )
a i, b jParameter decision by linear filter.
Described memoryless non-linear components 202 is meant that its output signal is irrelevant with the output signal of each state before, and only is the q order polynomial of current input state, and establish the current x of being input as (t), then current output signal y (t) can be expressed as:
y ( t ) = Σ p = 0 q a p x q ( t )
P wherein, q is by the parameter decision of memoryless nonlinear device.
The signal path that modulating unit 7 enters Channel Elements 8 comprises that branch road output enters the first input end of SINAD unit 13; The signal path that Channel Elements 8 enters nonlinear element 9 comprises that branch road output enters second input of SINAD unit 13; The signal path that source encoding 5 enters chnnel coding 6 comprises that branch road output enters the first input end of BER unit 14; Channel decoding 10 enters source decoding 11 signal paths and comprises that branch road output enters second input of BER unit 14; After obtaining described S2, O2, D, get according to the definition of SINAD (letter Nader):
SINAD = 10 lg O 2 N + D = 10 lg S 2 + N + D N + D
After obtaining described S1, O1, obtain BER (error rate) according to following formula:
Figure BDA0000055939970000084
After obtaining SINAD and BER, take SINAD and BER layering forecast method to carry out the electromagnetic interference prediction.
Be illustrated in figure 3 as multi-mode combination interference analysis process figure under the complex electromagnetic environment, execution in step 301 is carried out the multi-mode combination interference analysis under the complex electromagnetic environment, parallel then step 302 and 303 of carrying out;
Execution in step 302, the interference source that carries out under the complex electromagnetic environment obtains, at first execution in step 304 is carried out the selection of interference source obtain manner, if obtain manner is a mathematical modeling, then execution in step 305 is carried out the mathematical modeling of interference source, execution in step 306 is carried out the judgement of interference source type, if interference source is a transient state, then execution in step 307 utilizes small wave converting method to carry out wink-steady-state equivalent, execution in step 308 then, if interference source is a stable state, then direct execution in step 308, execution in step 308 is carried out the analytical calculation of equivalent stable state interference source; Whether execution in step 309 is carried out interference source and is calculated and finish judgement, does not finish if calculate, and then repeated execution of steps 302, finish if calculate, and then export multi-mode combination interference source to the radio station of working " non-linear " receiver mathematical units; If the interference source obtain manner detects for monitoring, then execution in step 311 utilizes measuring receiver to carry out the reception of interference source, the signature analysis of the interference source that execution in step 312 receives, execution in step 313 is carried out interference source type and is judged, if interference source is a transient state, then execution in step 307, if interference source is a stable state, then execution in step 308.
Execution in step 303, carry out the receiver modeling under the complex electromagnetic environment, at first execution in step 314, carry out the analysis of work radio station receiver nonlinear element, execution in step 315 then, carry out the extraction and the simplification of " non-linear " receiver mathematical units, execution in step 316 obtains work radio station " non-linear " receiver mathematical units, reception is from the multi-mode combination interference of step 310, execution in step 317, calculate in conjunction with the Discrete Nonlinear technology, execution in step 318, calculate the output signal of " non-linear " receiver, execution in step 319, carry out the analysis of signal performance, execution in step 320 is carried out SINAD calculating, and execution in step 321 compares SINAD and the SINAD tolerance limit D0 that calculates, if SINAD is more than or equal to D0 then do not disturb, prediction finishes, if SINAD is less than D0 then execution in step 322 is calculated BER, execution in step 323, the BER and the BER tolerance value C0 that calculate are compared, if BER then do not disturb smaller or equal to C0, prediction finishes, if BER is greater than C0 then disturb, adjustment work radio station running parameter, repeated execution of steps 303 is predicted.
Described wink-steady-state equivalent 307 is meant at the external interference source under the complex electromagnetic environment, determines its amplitude in frequency domain and time domain, the associate feature of power.According to probability statistics principle and signal detection means, set up the probability statistics unit of typical interference source, and according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source.
Described Discrete Nonlinear technology 317 is meant and utilizes discrete unit that it is simplified description at all kinds of interference sources under the complex electromagnetic environment at " non-linear " receiver front end; The linear filter of non-linear receiver in the frequency domain modeling, in the time domain modeling, is carried out time domain-frequency domain transform with FFT with its memoryless nonlinear device; Determine non-linear interference source with the multi-branched search utility, care signal is assessed through the performance behind the receiver unit.

Claims (9)

1. the multi-mode combination interference prognoses system under the complex electromagnetic environment, it is characterized in that: which comprises at least: work radio station unit transmitter, work unit, radio station " non-linear " receiver, multi-mode combination interference source under the complex electromagnetic environment, the Channel Elements that comprises background noise, the SINAD unit, the BER unit; The information source of work radio station unit transmitter is by output signal S1 after the source encoding, and S1 channel coding, modulating unit are handled the back and sent radiofrequency signal with certain power, enter a channel input of the Channel Elements that comprises background noise N; Simultaneously, the output in the multi-mode combination interference source under the complex electromagnetic environment disturbs D to enter the Channel Elements that comprises background noise N, background noise stack in two-way input signal and the Channel Elements, produce the output signal O2 of Channel Elements, radio frequency input as work unit, radio station " non-linear " receiver, the output signal O2 of Channel Elements enters unit, work radio station " non-linear " receiver after after the nonlinear element, channel decoding output signal be designated as O1, O1 enters output output at last through source decoding; Signal S2 after the decay of unit transmitter channel unit, described work radio station enters the input in SINAD unit, and the output signal O2 of Channel Elements enters another input of SINAD unit, by the two sides road calculated signals SINAD that enters the SINAD unit; And the output signal O1 of the channel decoding in the output signal S1 of the source encoding 5 in the work radio station unit transmitter and work unit, radio station " non-linear " receiver exports along separate routes and enters two inputs in BER unit, calculates BER by the BER unit.
2. the multi-mode combination interference prognoses system under the complex electromagnetic environment according to claim 1, it is characterized in that: the background noise in the described Channel Elements is used to describe different frequency range, the electromagnetic environment of different regions, utilize the random number sequence that meets Gaussian Profile to come the simulation background noise in the unit, obtain the reasonable value of noise power by the variance of adjusting random number sequence, implementation procedure is as follows:
1) utilize a stationary random process to represent background noise
Feature according to background noise is established as the stationary random process ξ that auto-correlation function is R (τ) (t), and the power of actual noise generally all is meant its AC power, and solution procedure is as follows:
Noise average power is: R (0)=E[ξ 2(t)], E[x] expression x mathematic expectaion;
The noise direct current power is: R (∞)=E 2[ξ (t)]
So its AC power is: R (0)-R (∞)=E[ξ 2(t)]-E 2[ξ (t)]=σ 2, σ 2The variance of expression ξ (t), promptly the power of noise is represented with the variance of its process.
2) utilize the random number sequence that meets Gaussian Profile to realize this stationary random process
According in the test or the typical environment level value that provides in the standard, by adjusting the variances sigma of this random number sequence 2Obtain rational background noise, different σ 2Value characterizes different noises.
3. the multi-mode combination interference prognoses system under the complex electromagnetic environment according to claim 1, it is characterized in that: the output in described multi-mode combination interference source disturbs D to be meant at the external interference source under the complex electromagnetic environment, determines its amplitude in frequency domain and time domain, the associate feature of power; It is according to probability statistics principle and signal detection means that D is disturbed in the output in described multi-mode combination interference source, set up the probability statistics unit of typical interference source, and according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source, with the D that is expressed as of its multi-mode combination back quantification.
4. the multi-mode combination interference prognoses system under the complex electromagnetic environment according to claim 1, it is characterized in that: described nonlinear element refers to the non-linear partial of unit, work radio station " non-linear " receiver, and nonlinear element is by the serial or parallel connection of linear filter 201,203 and memoryless non-linear components 202 or connection in series-parallel combination.
5. the multi-mode combination interference prognoses system under the complex electromagnetic environment according to claim 1, it is characterized in that: described linear filter is meant dynamic linear with memory filter, with the systematic sampling time be the interval, the output signal of its current state is not only relevant with the input signal of current state system, and with the input signal of m state before and the output signal of n state is relevant before, the input signal of m state and the linear superposition of the output signal of n state before before being, if current state is k, be input as x (k), then the output signal y of current linear filter (k) is:
y ( k ) = Σ i = 1 n a i y ( k - i ) + Σ j = 0 m b j x ( k - j )
a i, b jParameter decision by linear filter;
Described memoryless non-linear components is meant that its output signal is irrelevant with the output signal of each state before, and only is the q order polynomial of current input state, and establish the current x of being input as (t), then current output signal y (t) can be expressed as:
y ( t ) = Σ p = 0 q a p x q ( t )
P wherein, q is by the parameter decision of memoryless nonlinear device.
6. the multi-mode combination interference prognoses system under the complex electromagnetic environment according to claim 1 is characterized in that: the computing formula of SINAD (letter Nader) is:
SINAD = 10 lg O 2 N + D = 10 lg S 2 + N + D N + D
Described BER (error rate) computing formula is:
Figure FDA0000055939960000032
7. based on the multi-mode combination interference Forecasting Methodology under the complex electromagnetic environment of the multi-mode combination interference prognoses system under the complex electromagnetic environment: it is characterized in that: which comprises at least following steps:
Execution in step 301 is carried out the multi-mode combination interference analysis under the complex electromagnetic environment;
Execution in step 302, the interference source that carries out under the complex electromagnetic environment obtains;
At first execution in step 304 is carried out the selection of interference source obtain manner, if obtain manner is a mathematical modeling, then execution in step 305 is carried out the mathematical modeling of interference source, if the interference source obtain manner detects for monitoring, then execution in step 311 utilizes measuring receiver to carry out the reception of interference source;
Step 305 is carried out the mathematical modeling of interference source;
Execution in step 306 is carried out the judgement of interference source type, if interference source is a transient state, then execution in step 307; If interference source is a stable state, then direct execution in step 308;
Interference source is a transient state, and execution in step 307 utilizes small wave converting method to carry out wink-steady-state equivalent;
Execution in step 308 is carried out the analytical calculation of equivalent stable state interference source;
Whether execution in step 309 is carried out interference source and is calculated and finish judgement, does not finish if calculate, and then repeated execution of steps 302, finish execution in step 310 if calculate;
Execution in step 310, output multi-mode combination interference source is to the radio station of working " non-linear " receiver mathematical units;
If the interference source obtain manner detects for monitoring, then execution in step 311, utilize measuring receiver to carry out the reception of interference source;
Execution in step 312, the signature analysis of the interference source that receives;
Execution in step 313 is carried out interference source type and is judged that if interference source is a transient state, then execution in step 307, if interference source is a stable state, then execution in step 308.
Execution in step 303 is carried out the receiver modeling under the complex electromagnetic environment;
Execution in step 314 is carried out the analysis of work radio station receiver nonlinear element;
Execution in step 315 then, carry out the extraction and the simplification of " non-linear " receiver mathematical units;
Execution in step 316 obtains work radio station " non-linear " receiver mathematical units, receives the multi-mode combination interference from step 310;
Execution in step 317 is obtained Discrete Nonlinear and is carried out result calculated;
Execution in step 318 calculates the output signal of " non-linear " receiver;
Execution in step 319 is carried out the analysis of signal performance;
Execution in step 320 is carried out SINAD and is calculated;
Execution in step 321 compares SINAD and the SINAD tolerance limit D0 that calculates, if SINAD more than or equal to D0 then do not disturb, prediction finishes, if SINAD is less than D0 then execution in step 322 is calculated BER;
Execution in step 322 is calculated BER;
Execution in step 323 compares BER and the BER tolerance value C0 that calculates, if BER then do not disturb smaller or equal to C0, prediction finishes, if BER greater than C0 then disturb, adjusts work radio station running parameter, repeated execution of steps 303 is predicted.
According to claim 7 based on the multi-mode combination interference Forecasting Methodology under the complex electromagnetic environment of the multi-mode combination interference prognoses system under the complex electromagnetic environment: it is characterized in that: described wink-steady-state equivalent 307, be according to principle of wavelet analysis and minimum phase principle, interference source is carried out detecting in the time domain reproduction in--frequency domain inner analysis data characteristics--time domain, promptly realize the transient state and the analysis of stable state linked character of the transient state detection-steady-state analysis-transient state suppressing of interference source, set up the mathematical relationship of its stable state and transient state effect, form transient state-stable state conversion method, obtain the steady-state equivalent unit of interference source.
According to claim 7 based on the multi-mode combination interference Forecasting Methodology under the complex electromagnetic environment of the multi-mode combination interference prognoses system under the complex electromagnetic environment: it is characterized in that: described Discrete Nonlinear technology 317, being meant utilizes discrete unit that it is simplified description at all kinds of interference sources under the complex electromagnetic environment at " non-linear " receiver front end; The linear filter of non-linear receiver in the frequency domain modeling, in the time domain modeling, is carried out time domain-frequency domain transform with FFT with its memoryless nonlinear device; Determine non-linear interference source with the multi-branched search utility, care signal is assessed through the performance behind the receiver unit.
CN2011100968822A 2011-04-18 2011-04-18 System and method for predicting multimode combined interference in complex electromagnetic environment Expired - Fee Related CN102185666B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100968822A CN102185666B (en) 2011-04-18 2011-04-18 System and method for predicting multimode combined interference in complex electromagnetic environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100968822A CN102185666B (en) 2011-04-18 2011-04-18 System and method for predicting multimode combined interference in complex electromagnetic environment

Publications (2)

Publication Number Publication Date
CN102185666A true CN102185666A (en) 2011-09-14
CN102185666B CN102185666B (en) 2013-07-10

Family

ID=44571741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100968822A Expired - Fee Related CN102185666B (en) 2011-04-18 2011-04-18 System and method for predicting multimode combined interference in complex electromagnetic environment

Country Status (1)

Country Link
CN (1) CN102185666B (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102981064A (en) * 2012-10-09 2013-03-20 中国人民解放军63892部队 Aircraft external radio frequency electromagnetic environment prediction method and prediction system
CN103913661A (en) * 2014-04-09 2014-07-09 华北电力大学(保定) Method for estimating influences on telecommunication cable E1 wire by electromagnetic disturbance
CN103969530A (en) * 2014-05-09 2014-08-06 西安电子科技大学 Steady state test method applying time domain and frequency domain energy equivalence principle
CN104868960A (en) * 2015-04-30 2015-08-26 中国电子科技集团公司第七研究所 Method and system for predicting background interference noise level
CN106027176A (en) * 2016-06-29 2016-10-12 中国人民解放军军械工程学院 Receiver in-band dual frequency electromagnetic radiation large signal blocking interference effect prediction method
CN106209279A (en) * 2016-06-29 2016-12-07 中国人民解放军军械工程学院 Multifrequency electromagnetic radiation strong signal blocking interference effect Forecasting Methodology in receiver band
CN107465481A (en) * 2017-08-28 2017-12-12 电子科技大学 A kind of the electromagnetic spectrum umbrella cover Interference Suppression System and method of self-interference power sequence
CN109324244A (en) * 2018-09-26 2019-02-12 北京航空航天大学 A kind of Battle Field Electromagnetic generation method for unmanned plane survivability Evaluation
CN109358247A (en) * 2018-10-19 2019-02-19 中国人民解放军军事科学院国防工程研究院 Electronic equipment electromagnetic environmental effects pilot system and test method under electromagnetic pulse Composite Field environment
CN110988925A (en) * 2019-12-17 2020-04-10 北京遥测技术研究所 Pulse interference detection and parameter determination method for satellite navigation receiver
CN111624421A (en) * 2020-04-27 2020-09-04 西安电子科技大学 Electromagnetic interference monitoring and positioning method
CN112198388A (en) * 2020-10-12 2021-01-08 东南大学 Method for testing anti-interference sensitivity of artificial surface plasmon transmission line
CN114245345A (en) * 2021-11-25 2022-03-25 西安电子科技大学 Internet of vehicles power control method and system for imperfect channel state information
CN114895137A (en) * 2022-06-29 2022-08-12 广东电网有限责任公司佛山供电局 Power supply quality testing method for power supply equipment

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090238086A1 (en) * 2008-03-19 2009-09-24 Telefonaktiebolaget Lm Ericsson (Publ) Rank Dependent CQI Back-Off

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090238086A1 (en) * 2008-03-19 2009-09-24 Telefonaktiebolaget Lm Ericsson (Publ) Rank Dependent CQI Back-Off

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NANIEL J .COSTELLO: "The Road to Channel Capacity", 《PROCEEDINGS OF THE IEEE》, vol. 95, no. 6, 30 June 2007 (2007-06-30) *
周萍: "基于误码率分析的共址干扰预测研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, vol. 1, no. 10, 16 September 2010 (2010-09-16), pages 135 - 10 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102981064A (en) * 2012-10-09 2013-03-20 中国人民解放军63892部队 Aircraft external radio frequency electromagnetic environment prediction method and prediction system
CN102981064B (en) * 2012-10-09 2015-05-20 中国人民解放军63892部队 Aircraft external radio frequency electromagnetic environment prediction method and prediction system
CN103913661A (en) * 2014-04-09 2014-07-09 华北电力大学(保定) Method for estimating influences on telecommunication cable E1 wire by electromagnetic disturbance
CN103913661B (en) * 2014-04-09 2017-01-11 华北电力大学(保定) Method for estimating influences on telecommunication cable E1 wire by electromagnetic disturbance
CN103969530A (en) * 2014-05-09 2014-08-06 西安电子科技大学 Steady state test method applying time domain and frequency domain energy equivalence principle
CN103969530B (en) * 2014-05-09 2016-09-14 西安电子科技大学 A kind of steady-state method of test using time domain and the frequency domain energy principle of equal effects
CN104868960A (en) * 2015-04-30 2015-08-26 中国电子科技集团公司第七研究所 Method and system for predicting background interference noise level
CN106027176A (en) * 2016-06-29 2016-10-12 中国人民解放军军械工程学院 Receiver in-band dual frequency electromagnetic radiation large signal blocking interference effect prediction method
CN106209279A (en) * 2016-06-29 2016-12-07 中国人民解放军军械工程学院 Multifrequency electromagnetic radiation strong signal blocking interference effect Forecasting Methodology in receiver band
CN107465481A (en) * 2017-08-28 2017-12-12 电子科技大学 A kind of the electromagnetic spectrum umbrella cover Interference Suppression System and method of self-interference power sequence
CN109324244A (en) * 2018-09-26 2019-02-12 北京航空航天大学 A kind of Battle Field Electromagnetic generation method for unmanned plane survivability Evaluation
CN109358247A (en) * 2018-10-19 2019-02-19 中国人民解放军军事科学院国防工程研究院 Electronic equipment electromagnetic environmental effects pilot system and test method under electromagnetic pulse Composite Field environment
CN109358247B (en) * 2018-10-19 2021-01-12 中国人民解放军军事科学院国防工程研究院 Method for testing electromagnetic environment effect of electronic equipment in electromagnetic pulse composite field environment
CN110988925A (en) * 2019-12-17 2020-04-10 北京遥测技术研究所 Pulse interference detection and parameter determination method for satellite navigation receiver
CN111624421A (en) * 2020-04-27 2020-09-04 西安电子科技大学 Electromagnetic interference monitoring and positioning method
CN112198388A (en) * 2020-10-12 2021-01-08 东南大学 Method for testing anti-interference sensitivity of artificial surface plasmon transmission line
CN112198388B (en) * 2020-10-12 2023-09-29 东南大学 Method for testing anti-interference sensitivity of artificial surface plasmon transmission line
CN114245345A (en) * 2021-11-25 2022-03-25 西安电子科技大学 Internet of vehicles power control method and system for imperfect channel state information
CN114245345B (en) * 2021-11-25 2024-04-19 西安电子科技大学 Imperfect channel state information-oriented Internet of vehicles power control method and system
CN114895137A (en) * 2022-06-29 2022-08-12 广东电网有限责任公司佛山供电局 Power supply quality testing method for power supply equipment
CN114895137B (en) * 2022-06-29 2023-04-07 广东电网有限责任公司佛山供电局 Power supply quality testing method for power supply equipment

Also Published As

Publication number Publication date
CN102185666B (en) 2013-07-10

Similar Documents

Publication Publication Date Title
CN102185666B (en) System and method for predicting multimode combined interference in complex electromagnetic environment
CN106772268B (en) A kind of weak signal blind checking method under white Gaussian noise
CN102109553B (en) Digital spectrum analyzer and method for fast positioning of measuring signals
Altrad et al. A new mathematical analysis of the probability of detection in cognitive radio over fading channels
CN103777083B (en) Based on the capacitive equipment dielectric loss on-line monitoring method of Kalman's frequency-tracking
CN102164013B (en) Electromagnetic compatibility (EMC) comprehensive predictive analysis system and analysis method for co-site interference of radio station
Demir et al. Engine compartment UWB channel model for intravehicular wireless sensor networks
CN103941254A (en) Soil physical property classification recognition method and device based on geological radar
CN105223481A (en) Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function
CN103117822A (en) Device for receiver channel group delay measurement
CN103346847B (en) United compressing spectrum sensing method based on iteration attack detection
CN104143997A (en) Multi-access interference false locking resistance judgment method based on spread spectrum mechanism
CN115085829B (en) Method for testing sensitivity under existence of multipath fading interference and electronic equipment
CN103281141A (en) Blind spectrum sensing method and device
Li et al. Finite-state Markov wireless channel modeling for railway tunnel environments
CN102818941B (en) Electromagnetic radiation emission measurement method in outfield disturbed environment
CN104734793B (en) The energy detection method of wireless cooperation frequency spectrum perception based on p powers
CN101335567A (en) Ultra-wideband non-coherent system average bit error rate estimating method under S-V modified model fading channel of IEEE802.15.3a
CN106019250A (en) Repeating false target discriminating method based on angular glint
CN102263600B (en) Method and device for determining movement speed of terminal
CN103702357A (en) Intelligent utility network transmission packet error rate measuring and calculating method for establishing data packet collision model on basis of probability theory
CN102868403A (en) Testing system for testing main performance indexes of analog-to-digital converter
CN101373985B (en) Method and apparatus for evaluating multi-user detection interference elimination performance
KR101144527B1 (en) Apparatus for measuring short distance using overlap-based css system and distance measuring method using the same
CN104020366A (en) Method of determining receiver antenna sensitivity

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130710

Termination date: 20140418