CN114384468B - Target direct positioning method and system under inconsistent impulse noise environment - Google Patents
Target direct positioning method and system under inconsistent impulse noise environment Download PDFInfo
- Publication number
- CN114384468B CN114384468B CN202210292337.9A CN202210292337A CN114384468B CN 114384468 B CN114384468 B CN 114384468B CN 202210292337 A CN202210292337 A CN 202210292337A CN 114384468 B CN114384468 B CN 114384468B
- Authority
- CN
- China
- Prior art keywords
- inconsistent
- impulse noise
- signal
- receiving stations
- cost function
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0246—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves involving frequency difference of arrival or Doppler measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a method and a system for directly positioning a target in an inconsistent impulse noise environment, and belongs to the technical field of signal processing. The method comprises the following steps: step S1, useLA receiving station, byKIntercepting the secondary signal to obtain the secondary signal from the secondary signal acquisition unitpObserving signals of a single static radiation source, and sampling the observing signals; step S2 based onLInconsistent impulse noise structure for a receiving stationLEach receiving station intercepts the cost function in the time slot; step S3 based onLWeighted coefficient sum of receiving stationsLThe global cost function is constructed by the cost function in each time slot for each of the receiving stations. Aiming at the problem of inconsistent pulse noise dispersion coefficients, the invention uses the noise dispersion coefficients to weight signals with different signal-to-noise ratios by constructing a cost function for directly positioning a target, thereby solving the problem of reduced positioning performance of a direct positioning algorithm under the condition of inconsistent pulse noise.
Description
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a method and a system for directly positioning a target in an inconsistent impulse noise environment.
Background
Passive location techniques for targets determine target location by intercepting signals transmitted or reflected by the target with a receiving station without transmitting electromagnetic signals themselves. The passive positioning has the advantages of low cost, strong anti-interference capability and the like, and has important application value in the fields of surface naval vessel positioning, sea area monitoring, illegal ground intrusion and the like.
The target direct localization algorithm is one of the important developments of passive localization. A Direct Positioning (DPD) method was first proposed in "WEISS A J. Direct Position Determination of narrow band radio transmitters [ J ]. IEEE Signal Processing Letters, IEEE, 2004, 11(5): 513-. Compared with a classical two-step passive positioning algorithm, the direct positioning algorithm directly uses an observed signal without completing the estimation of time difference/frequency difference, and generally has better positioning accuracy under the condition of low signal-to-noise ratio.
Most target passive localization algorithms, including DPD, are derived under gaussian noise conditions. However, whether in nature or due to human factors, actual noise, such as noise in power line communication systems, shallow sea acoustic channel noise, etc., often exhibits significant spike characteristics. The profile of such noise has a thicker tail than a gaussian profile. For this type of noise, it can be modeled, typically with an alpha-stationary distribution. Wherein, the dispersion coefficient gamma is one of important parameters for describing alpha-stable distribution and characterizes the dispersion degree of the distribution.
In fact, under the above impulse noise environment, the performance of many conventional parameter estimation algorithms based on gaussian noise conditions is significantly deteriorated. The article "Jinyan, Navy, Ji hong Bing. OFDM time domain parameter estimation based on correlation entropy under impulse noise [ J ] systematic engineering and electronic technology, 2015, 37(12): 2701-. The paper "Cai Rui Yan, Yang. coherent distribution source DOA estimation method based on correlation entropy under impulse noise [ J ]. electronic and informatics, 2020, 42(11): 2600-. The paper \20319j, Alpha steady distributed noise environment research [ D ]. university of big succession, 2010 "discusses the performance deterioration of the time delay estimation of the classical time delay estimation algorithm under the impulse noise condition.
Furthermore, most classical passive positioning algorithms (including DPD algorithms) typically assume that the gaussian noise at each receiving station is independent and consistent, i.e., the noise is identically distributed gaussian noise. In practice, although the noise at each receiving station is independent, the noise may not be consistent, for example, the noise at each receiving station is gaussian noise but the noise power is different; or the noise of each receiving station is pulse noise, but the dispersion parameters of the noise are different; or the noise of part of the receiving stations is gaussian noise and the noise of other receiving stations is impulse noise. It has been shown that the positioning accuracy of classical direct positioning algorithms deteriorates when the noise power of the receivers differs. In order to solve the problem, based on the maximum likelihood estimation criterion, the article "Zhonghua, Raney, Sunsmen, etc.. pulse train radiation source direct positioning algorithm [ J ] under the background of inconsistent noise, university of Western Ann traffic, 2021, 55(08): 157-. In the case that the gaussian noise power of the receiving station is inconsistent and the transmitted signal is a pulse signal, the NWO-ML-DPD has better positioning accuracy than the classical DPD algorithm. Through verification, the positioning accuracy of the algorithm is reduced to some extent under the environment that the dispersion coefficients of the impulse noise are inconsistent.
Disclosure of Invention
In order to solve the technical problem, the application provides a target direct positioning scheme under an inconsistent impulse noise environment.
The invention discloses a target direct positioning method under an inconsistent impulse noise environment in a first aspect. The method comprises the following steps:
step S1, useLA receiving station, byKIntercepting the secondary signal to obtain the secondary signal from the secondary signal acquisition unitpObserving signals of a single static radiation source, and sampling the observing signals;
wherein the observed signal is compared to the bitIn thatpThe emission signal of a single static radiation source contains the influence factors of time delay, Doppler frequency shift and noise;
step S2 based onLConstructing the non-uniform impulse noise of a receiving stationLEach receiving station intercepts the cost function in the time slot;
step S3 based onLWeighting factors for a receiving station and saidLConstructing a global cost function by the cost function of each receiving station in each time slot;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the true position of the single stationary radiation sourceLThe weighting coefficients of the receiving stations being dependent on saidLThe dispersion factor of the inconsistent impulse noise of the individual receiving stations.
According to the method of the first aspect of the present invention, in the step S1, the stepLEach of the receiving stations is in processKA secondary signal is intercepted, the time of the single interception isTFirst, oflThe receiving station is atkThe position when the signal is intercepted isAt a speed ofThe first mentionedkIn the time slot in which the secondary intercept is locatedpThe emission signal of a single stationary radiation source isWhereinIs the carrier frequency (c) of the carrier,is a bandwidth ofWA narrow band signal ofThen it is firstlAt the second receiving stationkThe secondary intercepted observation signals are:
wherein the content of the first and second substances,is the firstkIn which the transmitted signal propagates from the target to the secondlThe transmission delay of the individual receiving stations,cin order to be the speed of light,is the number of the Euclidean norm,is the firstkIn which the transmitted signal propagates from the target to the secondlThe doppler shift produced by each of the receiving stations,,is the firstkIn the secondary interception of the secondlAdditive impulse noise of a receiving station, said additive impulse noise being subject toThe distribution is stable.
According to the method of the first aspect of the invention, in said step S1, so as toSampling the observation signal for a sampling period, thenlA receiving station at the secondkThe samples of the observation signal received in the time slot in which the secondary truncation is located are:
wherein, the first and the second end of the pipe are connected with each other,
wherein the content of the first and second substances,is shown inIs a diagonal matrix of diagonal elements,to shift down an operator, saidBy passingThe acquisition of the rows of the identity matrix is cyclically shifted,indicating rounding down, useTo realizeIs performed.
According to the method of the first aspect of the present invention, in the step S2, the step S is acquiredLDispersion coefficient of inconsistent impulse noise of receiving stationThen the first mentionedkIn the time slot in which the secondary truncation is locatediThe cost function for each sample point is:
wherein the content of the first and second substances,in the form of a gaussian kernel function,in order to take out the conjugate operation,as a parameter of the length of the nucleus,representing the difference of the multiple receiver station observation samples at each sampling point in each time slot,representing a vectorTo (1) aiAnd (4) each element.
According to the method of the first aspect of the present invention, in the step S3, the first stepkThe time slot in which the secondary interception is positioned,NAdding the cost functions of the sampling points to obtain the global cost function, as follows:
according to the method of the first aspect of the present invention, in step S3, the maximum value point of the global cost function is:
the invention discloses a target direct positioning system in an inconsistent impulse noise environment in a second aspect. The system comprises:
a first processing unit configured to utilizeLA receiving station, byKIntercepting the secondary signal to obtain the secondary signal from the secondary signal acquisition unitpObserving signals of a single static radiation source, and sampling the observing signals;
wherein the observed signal is compared to the sitepThe emission signal of a single static radiation source contains the influence factors of time delay, Doppler frequency shift and noise;
a second processing unit configured to, based on theLConstructing the non-uniform impulse noise of a receiving stationLEach receiving station intercepts the cost function in the time slot;
a third processing unit configured to, based on theLWeighting factors for a receiving station and saidLConstructing a global cost function by the cost function of each receiving station in each time slot;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the true position of the single stationary radiation sourceLThe weighting coefficients of the receiving stations being dependent on saidLThe dispersion factor of the inconsistent impulse noise of the individual receiving stations.
According to the second aspect of the invention, the systemLEach of the receiving stations is in processKA secondary signal is intercepted, the time of the single interception isTOf 1 atlThe receiving station is atkThe position when the signal is intercepted secondarily isAt a speed ofThe first mentionedkIn the time slot in which the secondary intercept is locatedpThe emission signal of a single stationary radiation source isWhereinIs the carrier frequency (c) of the carrier,is a bandwidth ofWA narrow band signal ofThen it is firstlAt the second receiving stationkThe secondary intercepted observation signals are:
wherein the content of the first and second substances,is the firstkIn which the transmitted signal propagates from the target to the secondlThe transmission delay of the individual receiving stations,cin order to be the speed of light,is the number of the Euclidean norm,is the firstkIn which the transmitted signal propagates from the target to the secondlThe doppler shift produced by each of the receiving stations,,is the firstkIn the secondary interception of the secondlAdditive impulse noise of a receiving station, said additive impulse noise being subject toThe distribution is stable.
According to the system of the second aspect of the present invention, the first processing unit is specifically configured to: to be provided withSampling the observation signal for a sampling period, thenlA receiving station at the secondkThe samples of the observation signal received in the time slot where the secondary truncation is located are:
wherein, the first and the second end of the pipe are connected with each other,
wherein the content of the first and second substances,is shown inIs a diagonal matrix of diagonal elements,to shift down an operator, saidBy passingThe acquisition of the rows of the identity matrix is cyclically shifted,indicating rounding down, useTo realizeIs performed.
According to the system of the second aspect of the invention, the second processing unit is specifically configured to: obtaining theLDispersion coefficient of inconsistent impulse noise of receiving stationThen said firstkIn the time slot in which the secondary truncation is locatediThe cost function for each sample point is:
wherein, the first and the second end of the pipe are connected with each other,in the form of a gaussian kernel function,in order to take out the conjugate operation,as a parameter of the length of the nucleus,representing multiple connections at each sample point in each slotThe receiving station observes the difference in the samples,representing a vectorTo (1) aiAnd (4) each element.
According to the system of the second aspect of the present invention, the third processing unit is specifically configured to: will be the firstkThe time slot in which the secondary interception is positioned,NAdding cost functions of the sampling points to obtain the global cost function, wherein the global cost function is as follows:
according to the system of the second aspect of the present invention, the third processing unit is specifically configured to: the maximum value points of the global cost function are as follows:
a third aspect of the invention discloses an electronic device. The electronic device comprises a memory storing a computer program and a processor implementing the steps of a method for direct target localization in an environment of inconsistent impulse noise according to any one of the first aspect of the present disclosure when the computer program is executed by the processor.
A fourth aspect of the invention discloses a computer-readable storage medium. The computer readable storage medium has stored thereon a computer program which, when being executed by a processor, implements the steps of a method for direct localization of an object in an environment of inconsistent impulse noise according to any one of the first aspect of the present disclosure.
In conclusion, the technical scheme provided by the invention constructs a cost function for directly positioning the target, and weights the signals with different signal-to-noise ratios by using the noise dispersion coefficients aiming at the problem of inconsistent pulse noise dispersion coefficients, thereby solving the problem of reduced positioning performance of a direct positioning algorithm under the condition of inconsistent pulse noise.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description in the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for directly locating a target in an environment with inconsistent impulse noise according to an embodiment of the present invention;
FIG. 2 is a flow chart illustrating a direct target location algorithm according to a first embodiment of the present invention;
FIG. 3 is a diagram of the location accuracy RMSE of the ML-DPD algorithm, the NWO-ML-DPD algorithm and the NU-MCC-DPD algorithm varying with the generalized signal-to-noise ratio GSNR under the condition of inconsistent noise according to the second embodiment of the present invention;
FIG. 4 is a block diagram of a direct target location system in an environment of inconsistent impulse noise in accordance with an embodiment of the present invention;
fig. 5 is a block diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention discloses a target direct positioning method under an inconsistent impulse noise environment in a first aspect. FIG. 1 is a flow chart of a method for directly locating a target in an environment with inconsistent impulse noise according to an embodiment of the present invention; as shown in fig. 1, the method includes:
step S1, useLA receiving station, byKIntercepting the secondary signal to obtain the secondary signal from the secondary signal acquisition unitpObserving signals of a single static radiation source, and sampling the observing signals;
wherein the observed signal is compared to the sitepThe emission signal of a single static radiation source contains the influence factors of time delay, Doppler frequency shift and noise;
step S2 based onLConstructing the non-uniform impulse noise of a receiving stationLEach of the receiving stations intercepts the cost function in the time slot at each interception;
step S3 based onLWeighting factors for a receiving station and saidLConstructing a global cost function by the cost function of each receiving station in each time slot;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the true position of the single stationary radiation sourceLThe weighting coefficients of the receiving stations being dependent on saidLThe dispersion factor of the inconsistent impulse noise of the individual receiving stations.
In some embodiments, in the step S1, the stepLEach of the receiving stations is in processKA secondary signal is intercepted, the time of the single interception isTOf 1 atlThe receiving station is atkThe position when the signal is intercepted secondarily isAt a speed ofThe first mentionedkIn the time slot in which the secondary intercept is locatedpThe emission signal of a single stationary radiation source isIn whichIs the carrier frequency (c) of the carrier,is a bandwidth ofWA narrow band signal ofThen it is firstlAt the second receiving stationkThe secondary intercepted observation signals are:
wherein the content of the first and second substances,is the firstkIn which the transmitted signal propagates from the target to the secondlThe transmission delay of the individual receiving stations,cin order to be the speed of light,is the number of the Euclidean norm,is the firstkIn which the transmitted signal propagates from the target to the secondlThe doppler shift produced by each of the receiving stations,,is the firstkThe second in the secondary interceptlAdditive impulse noise of a receiving station, said additive impulse noise being subject toThe distribution is stable.
In some embodiments, in the step S1, toSampling the observation signal for a sampling period, thenlA receiving station at the secondkThe samples of the observation signal received in the time slot in which the secondary truncation is located are:
wherein the content of the first and second substances,
wherein the content of the first and second substances,is shown inIs a diagonal matrix of diagonal elements,to shift down an operator, theBy passingThe acquisition of the rows of the identity matrix is cyclically shifted,indicating rounding down, useTo realizeIs performed.
In some embodiments, in the step S2, the step S2 is obtainedLDispersion coefficient of inconsistent impulse noise of receiving stationThen said firstkIn the time slot in which the secondary truncation is locatediThe cost function for each sample point is:
wherein, the first and the second end of the pipe are connected with each other,in the form of a gaussian kernel function,in order to take out the conjugate operation,as a parameter of the length of the nucleus,representing the difference of the multiple receiver station observation samples at each sampling point in each time slot,representing a vectorTo (1) aiAnd (4) each element.
In some embodiments, in the step S3, the second stepkThe time slot where the secondary interception is located,NAdding the cost functions of the sampling points to obtain the global cost function, as follows:
in some embodiments, in the step S3, the maximum value point of the global cost function is:
first embodiment
FIG. 2 is a flow chart illustrating a direct target location algorithm according to a first embodiment of the present invention; as shown in fig. 2, the process includes:
step 1: consider the utilization ofLA receiving station pair is locatedpA single stationary radiation source target at. Assuming that each receiving station performsKA secondary signal is intercepted, the time of the single interception isT. First, thelA receiving station is atkThe position and velocity of the secondary intercepted signal are respectivelyAnd,,. Is arranged at the firstkThe transmission signal of the radiation source in a time slot isWherein, in the step (A),is the carrier frequency (c) of the carrier,is a bandwidth ofWA narrow-band signal of (2) satisfying. Taking into account the effects of time delay, Doppler shift and noiselA receiving stationkThe secondary intercepted observation signals are:
wherein the content of the first and second substances,is as followskThe sub-intercepted transmitted signal propagates from the target to the secondlThe transmission delay of the individual receiving stations,cin order to be the speed of light,is the Euclidean norm;is as followskThe sub-intercepted transmitted signal propagates from the target to the secondlA doppler shift produced by each receiving station, wherein,;is a firstkAt the time of secondary interceptionlAdditive impulse noise at a receiving station, complianceThe distribution is stable.
Step 2: to be provided withSampling the received signal for a sampling period, thenlA receiving stationkSamples of a received signal of a time slotComprises the following steps:
wherein, the first and the second end of the pipe are connected with each other,. Writing the above equation into the form of a vector
Wherein the content of the first and second substances,
wherein the content of the first and second substances,is shown inIs a diagonal matrix of diagonal elements.To shift operators downward, byThe result is obtained by circularly shifting the rows of the identity matrix,indicating rounding down, useTo realizeIs performed.
And step 3: taking into account the fact that the noise of the receiving stations is not uniform, eachThe impulse noise dispersion parameters of the receiving station are respectivelyUsing the sample value instead of expectation to obtain the secondkIn a time slotiThe cost function for each sample point is:
wherein the content of the first and second substances,is a function of a gaussian kernel, and is,in order to take out the conjugate operation,is a kernel length parameter.Representing the difference of the multiple receiver station observation samples at each sampling point in each time slot,representing a vectorTo (1) aiAnd (4) each element.
And 4, step 4: adding the cost functions of different interception moments and different sampling points to obtain a global cost function:
Second embodiment
Simulation conditions are as follows:
taking a static target position asTaking into account the number of receiving stationsIntercepting time slotsThe movement speeds are allThe initial positions of the receiving stations are respectively、. The carrier frequency of the transmitted signal is 2GHz, the signal bandwidth is 200kHz, and the interception time is 3.9ms each time.
Since impulse noise does not have finite variance, a generalized signal-to-noise ratio (GSNR) is defined:
wherein the content of the first and second substances,is the variance of the signal and is,to comply withDispersion parameters of stably distributed noise. The simulation experiment adopts Root Mean Square Error (RMSE) to measure the positioning performance of the algorithm, which is defined as follows:
wherein the content of the first and second substances,Qfor the number of Monte Carlo experiments, in the experiments described herein;Is as followsqEstimated location of the target in the submonol experiment.
Setting parameters:
given impulse noise characteristic parametersKernel length parameter of MCC criterion. The generalized signal-to-noise ratios of the signals received by the first receiving station and the third receiving station are respectively fixed to-5 dB and-2 dB, and the generalized signal-to-noise ratios of the signals received by the second receiving station, the fourth receiving station and the fifth receiving station are changed to-15, -10, -5, 0, 5, 10 and 15 dB.
FIG. 3 is a diagram of the location accuracy RMSE of the ML-DPD algorithm, the NWO-ML-DPD algorithm and the NU-MCC-DPD algorithm varying with the generalized signal-to-noise ratio GSNR under the condition of inconsistent noise according to the second embodiment of the present invention; as shown in fig. 3, the algorithm of the present invention is better than the positioning algorithm based on the maximum likelihood criterion. Compared with the NWO-ML-DPD algorithm without considering impulse noise, the NU-MCC-DPD algorithm which utilizes the noise dispersion coefficient to weight signals with different signal-to-noise ratios obviously improves the positioning precision of the target.
The invention discloses a target direct positioning system in an inconsistent impulse noise environment in a second aspect. FIG. 4 is a block diagram of a direct target location system in an environment of inconsistent impulse noise in accordance with an embodiment of the present invention; as shown in fig. 4, the system 400 includes:
a first processing unit 401 configured to utilizeLA receiving station, byKIntercepting the secondary signal to obtain the secondary signal from the secondary signalpObserving signals of a single static radiation source, and sampling the observing signals;
wherein the observed signal is compared to the sitepThe emission signal of a single static radiation source contains the influence factors of time delay, Doppler frequency shift and noise;
a second processing unit 402 configured to, based on theLConstructing the non-uniform impulse noise of a receiving stationLEach of the receiving stations intercepts the cost function in the time slot at each interception;
a third processing unit 403 configured to, based on the aboveLWeighting factors for a receiving station and saidLConstructing a global cost function by the cost function of each receiving station in each time slot;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the true position of the single stationary radiation sourceLThe weighting coefficients of the receiving stations being dependent on saidLThe dispersion factor of the inconsistent impulse noise of the individual receiving stations.
According to the second aspect of the invention, the systemLEach of the receiving stations is in processKSecondary signal interception for a single interception time ofTOf 1 atlThe receiving station is atkThe position when the signal is intercepted secondarily isAt a speed ofThe first mentionedkIn the time slot in which the secondary intercept is locatedpThe emission signal of a single stationary radiation source isWhereinIs the carrier frequency (c) of the carrier,is a bandwidth ofWA narrow band signal ofThen it is firstlAt the second receiving stationkThe secondary intercepted observation signals are:
wherein the content of the first and second substances,is the firstkIn which the transmitted signal propagates from the target to the secondlThe transmission delay of the individual receiving stations,cin order to be the speed of light,is the number of the Euclidean norm,is the firstkIn which the transmitted signal propagates from the target to the secondlThe doppler shift produced by each of the receiving stations,,is the firstkSecondary cuttingGet the center oflAdditive impulse noise of a receiving station, said additive impulse noise being subject toThe distribution is stable.
According to the system of the second aspect of the present invention, the first processing unit 401 is specifically configured to: to be provided withSampling the observation signal for a sampling period, thenlA receiving station at the secondkThe samples of the observation signal received in the time slot where the secondary truncation is located are:
wherein the content of the first and second substances,
wherein, the first and the second end of the pipe are connected with each other,is shown inIs a diagonal matrix of diagonal elements,to shift down an operator, theBy passingThe acquisition of the rows of the identity matrix is cyclically shifted,indicating rounding down, useTo realizeIs performed.
According to the system of the second aspect of the present invention, the second processing unit 402 is specifically configured to: obtaining theLDispersion coefficient of inconsistent impulse noise of receiving stationThen said firstkIn the time slot in which the secondary truncation is locatediThe cost function for each sample point is:
wherein the content of the first and second substances,in the form of a gaussian kernel function,in order to take out the conjugate operation,as a parameter of the length of the nucleus,representing the difference of the multiple receiver station observation samples at each sampling point in each time slot,representing a vectorTo (1) aiAnd (4) each element.
According to the system of the second aspect of the present invention, the third processing unit 403 is specifically configured to: will be the firstkThe time slot where the secondary interception is located,NAdding the cost functions of the sampling points to obtain the global cost function, as follows:
according to the system of the second aspect of the present invention, the third processing unit 403 is specifically configured to: the maximum value points of the global cost function are as follows:
a third aspect of the invention discloses an electronic device. The electronic device comprises a memory storing a computer program and a processor implementing the steps of a method for direct target localization in an environment of inconsistent impulse noise according to any one of the first aspect of the present disclosure when the computer program is executed by the processor.
Fig. 5 is a block diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 5, the electronic device includes a processor, a memory, a communication interface, a display screen, and an input device, which are connected by a system bus. Wherein the processor of the electronic device is configured to provide computing and control capabilities. The memory of the electronic equipment comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The communication interface of the electronic device is used for carrying out wired or wireless communication with an external terminal, and the wireless communication can be realized through WIFI, an operator network, Near Field Communication (NFC) or other technologies. The display screen of the electronic equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the electronic equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the electronic equipment, an external keyboard, a touch pad or a mouse and the like.
It will be understood by those skilled in the art that the structure shown in fig. 5 is only a partial block diagram related to the technical solution of the present disclosure, and does not constitute a limitation to the electronic device to which the solution of the present disclosure is applied, and a specific electronic device may include more or less components than those shown in the drawings, or combine some components, or have different arrangements of components.
A fourth aspect of the invention discloses a computer-readable storage medium. The computer readable storage medium has stored thereon a computer program which, when being executed by a processor, implements the steps of a method for direct localization of an object in an environment of inconsistent impulse noise according to any one of the first aspect of the present disclosure.
In summary, the technical scheme provided by the invention has the following advantages: (1) the method provided by the invention realizes the direct positioning of the moving platform to the target by utilizing the time delay and Doppler frequency shift information in the received signal; (2) the method is based on parametric modeling of the impulse noise, constructs a corresponding cost function, and can improve the positioning accuracy of the direct positioning algorithm in the impulse noise environment; (3) the method provided by the invention aims at the condition that the dispersion coefficients of the impulse noise are inconsistent, weights different signal-to-noise ratios by using the dispersion coefficients of the impulse noise, increases the positioning contribution of the signal with the high signal-to-noise ratio, and effectively improves the target positioning precision under the condition that the impulse noise is inconsistent.
It should be noted that the technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, however, as long as there is no contradiction between the combinations of the technical features, the scope of the present description should be considered. The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is specific and detailed, but not to be understood as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (9)
1. A method for directly locating an object in an environment of inconsistent impulse noise, the method comprising:
s1, acquiring observation signals from a single static radiation source at p position by utilizing L receiving stations through signal interception for K times, and sampling the observation signals;
wherein the observation signal contains time delay, Doppler shift and noise influence factors compared with the emission signal of the single stationary radiation source at p;
step S2, constructing a cost function of each of the L receiving stations in the time slot where each interception is located based on the inconsistent impulse noise of the L receiving stations;
step S3, constructing a global cost function based on the weighting coefficients of the L receiving stations and the cost function of each of the L receiving stations in each of the time slots;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the real position of the single stationary radiation source, and the weighting coefficients of the L receiving stations depend on the dispersion coefficients of the inconsistent impulse noise of the L receiving stations.
2. The method as claimed in claim 1, wherein in step S1, each of the L receiving stations is performing K signal truncations, and the single truncation time is T,the position of the receiving station at the time of the k-th signal interception is pl,kVelocity vl,kIn the time slot of the k-th interception, the emission signal of the single static radiation source at p isWherein f iscIs the carrier frequency, sk(t) is a narrow band signal with a bandwidth of W, and W < fcIf the observed signal intercepted by the ith receiving station at the kth time is:
0<t≤T,l=1,2,...,L;k=1,2,...,K
wherein the content of the first and second substances,for the transmission delay of the transmission signal from the target to the l receiving station in the k-th interception, c is the speed of light, | | | | is the Euclidean norm, fl,kA doppler shift resulting from propagation of said transmitted signal from the target to the l-th receiving station in said k-th intercept,zl,k(t) additive impulse noise for said l-th receiving station in said k-th intercept, said additive impulse noise obeying an alpha-stationary distribution.
3. The method for directly locating the target in the environment of inconsistent impulse noise of claim 2, wherein in said step S1, T is usedsSampling the observation signal for a sampling period, where samples of the observation signal received by the ith receiving station in the time slot where the kth truncation is performed are as follows:
rl,k=Al,kFl,ksk+zl,k,l=1,2,…,L;k=1,2,…,K
wherein the content of the first and second substances,
4. The method as claimed in claim 3, wherein in step S2, the dispersion coefficient γ of the inconsistent impulse noise of the L receiving stations is obtained1,γ2,...,γLIf the sampling point in the time slot where the kth truncation is located is:
wherein the content of the first and second substances,is a Gaussian kernel function, (. cndot)*To take conjugate operations, σ (σ > 0) is the kernel length parameter, el,k,i(p)=(rl,k-Al,kFl,ksk)i,l=1,2,...,L,el,k,i(p) represents the difference of the observed samples of the plurality of receiving stations at each sampling point in each time slot, (a)iRepresenting the ith element of vector a.
7. a system for direct localization of an object in an environment of inconsistent impulse noise, the system comprising:
a first processing unit, configured to acquire observation signals from a single stationary radiation source located at p by K times of signal interception using L receiving stations, and perform sampling processing on the observation signals;
wherein the observation signal contains time delay, Doppler shift and noise influence factors compared with the emission signal of the single stationary radiation source at p;
a second processing unit configured to construct a cost function of each of the L receiving stations in a time slot in which each truncation is located based on the inconsistent impulse noise of the L receiving stations;
a third processing unit configured to construct a global cost function based on the weighting coefficients of the L receiving stations and the cost function of each of the L receiving stations in the respective time slots;
wherein the maximum point of the global cost function characterizes the estimated position of the single stationary radiation source as the true position of the single stationary radiation source, and the weighting coefficients of the L receiving stations depend on the dispersion coefficients of the inconsistent impulse noise of the L receiving stations.
8. An electronic device, characterized in that the electronic device comprises a memory and a processor, the memory stores a computer program, and the processor, when executing the computer program, implements the steps in a method for direct localization of an object in an environment of inconsistent impulse noise according to any of claims 1 to 6.
9. A computer-readable storage medium, having stored thereon a computer program which, when being executed by a processor, carries out the steps of a method for direct localization of an object in an environment of inconsistent impulse noise as claimed in any one of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210292337.9A CN114384468B (en) | 2022-03-24 | 2022-03-24 | Target direct positioning method and system under inconsistent impulse noise environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210292337.9A CN114384468B (en) | 2022-03-24 | 2022-03-24 | Target direct positioning method and system under inconsistent impulse noise environment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114384468A CN114384468A (en) | 2022-04-22 |
CN114384468B true CN114384468B (en) | 2022-06-17 |
Family
ID=81205313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210292337.9A Active CN114384468B (en) | 2022-03-24 | 2022-03-24 | Target direct positioning method and system under inconsistent impulse noise environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114384468B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102210168A (en) * | 2008-11-11 | 2011-10-05 | 真实定位公司 | Use of radio access technology diversity for location |
WO2015063595A1 (en) * | 2013-11-01 | 2015-05-07 | Cgg Services Sa | Method and device for deblending seismic data using self-adapting and/or selective radon interpolation |
CN105024771A (en) * | 2015-07-20 | 2015-11-04 | 西安电子科技大学 | Frequency spectrum sensing method under Alpha stable distributed noises |
CN108872971A (en) * | 2018-07-19 | 2018-11-23 | 中国人民解放军战略支援部队信息工程大学 | A kind of object localization method and device based on the single array of movement |
CN109031196A (en) * | 2018-07-19 | 2018-12-18 | 中国人民解放军战略支援部队信息工程大学 | Based on the direct localization method of maximum likelihood of the motion view survey station to multisignal source |
CN112834981A (en) * | 2021-03-15 | 2021-05-25 | 哈尔滨工程大学 | Null array arrival direction estimation method under impulse noise background |
CN113466784A (en) * | 2021-06-28 | 2021-10-01 | 台州学院 | Self-adaptive distributed source DOA estimation method under strong pulse noise |
-
2022
- 2022-03-24 CN CN202210292337.9A patent/CN114384468B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102210168A (en) * | 2008-11-11 | 2011-10-05 | 真实定位公司 | Use of radio access technology diversity for location |
WO2015063595A1 (en) * | 2013-11-01 | 2015-05-07 | Cgg Services Sa | Method and device for deblending seismic data using self-adapting and/or selective radon interpolation |
CN105024771A (en) * | 2015-07-20 | 2015-11-04 | 西安电子科技大学 | Frequency spectrum sensing method under Alpha stable distributed noises |
CN108872971A (en) * | 2018-07-19 | 2018-11-23 | 中国人民解放军战略支援部队信息工程大学 | A kind of object localization method and device based on the single array of movement |
CN109031196A (en) * | 2018-07-19 | 2018-12-18 | 中国人民解放军战略支援部队信息工程大学 | Based on the direct localization method of maximum likelihood of the motion view survey station to multisignal source |
CN112834981A (en) * | 2021-03-15 | 2021-05-25 | 哈尔滨工程大学 | Null array arrival direction estimation method under impulse noise background |
CN113466784A (en) * | 2021-06-28 | 2021-10-01 | 台州学院 | Self-adaptive distributed source DOA estimation method under strong pulse noise |
Non-Patent Citations (2)
Title |
---|
Guosheng Yang 等.Joint Estimation of Timing and Carrier Phase Offsets for.《IEEE COMMUNICATIONS LETTERS》.2018,第22卷(第1期),全文. * |
金艳 等.α 稳定分布噪声下基于最优L-柯西加权的.《***工程与电子技术》.2016,第38卷(第7期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN114384468A (en) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guerci et al. | Joint design and operation of shared spectrum access for radar and communications | |
US10374902B2 (en) | Method for clustering wireless channel MPCs based on a KPD doctrine | |
Gante et al. | Deep learning architectures for accurate millimeter wave positioning in 5G | |
CN106970369B (en) | The broadband radar target detection method of knowledge based auxiliary | |
CN112162244B (en) | Event trigger target tracking method under related noise and random packet loss environment | |
CN108089147A (en) | Improved shortwave unit localization method | |
CN114371445A (en) | Multi-radiation source direct positioning method based on single unmanned aerial vehicle | |
CN106125059B (en) | Nonparametric Combined estimator signal and the Passive Location of position | |
Cao et al. | DNN-based peak sequence classification CFAR detection algorithm for high-resolution FMCW radar | |
CN105974362A (en) | High-precision passive positioning method for jointly estimating signal parameter and position | |
CN114384468B (en) | Target direct positioning method and system under inconsistent impulse noise environment | |
Yang et al. | Central DOA estimation of incoherently distributed noncircular sources with cross-correlation matrix | |
Liu et al. | Flexible robust adaptive beamforming method with multiple separately widened nulls | |
Al-Sadoon et al. | New and less complex approach to estimate angles of arrival | |
Xia et al. | Kalman particle filtering algorithm for symmetric alpha‐stable distribution signals with application to high frequency time difference of arrival geolocation | |
CN113938358B (en) | Time delay determination method and terminal | |
Tomikowski et al. | Acceleration of Radio Direction Finder Algorithm in FPGA Computing Platform | |
Stoilkovic et al. | A cascade‐connected neural model for improved 2D DOA estimation of an EM signal | |
CN110391820A (en) | A kind of Novel Communication method of reseptance for evading co-channel interference based on DFT | |
Xu et al. | Data‐oriented calibration method to reduce measurement bias in SFAP‐based GNSS receivers | |
CN114460534A (en) | Positioning method and system based on maximum correlation entropy under impulse noise environment | |
Kim et al. | Short‐term fading model for signals reflected by ocean surfaces in underwater acoustic communication | |
Mylonakis et al. | A Novel Three-Dimensional Direction-of-Arrival Estimation Approach Using a Deep Convolutional Neural Network | |
US10972316B1 (en) | Channel estimation using a chirp signal and the Fractional Fourier Transform | |
Berdanier et al. | Phase based 2D passive source localisation using receiver networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |