CN102480299B - The shock resistance method and system of phase-coded signal sequence - Google Patents
The shock resistance method and system of phase-coded signal sequence Download PDFInfo
- Publication number
- CN102480299B CN102480299B CN201010566201.XA CN201010566201A CN102480299B CN 102480299 B CN102480299 B CN 102480299B CN 201010566201 A CN201010566201 A CN 201010566201A CN 102480299 B CN102480299 B CN 102480299B
- Authority
- CN
- China
- Prior art keywords
- sequence
- amplitude
- impact
- phase
- coded signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses the shock resistance method of phase-coded signal sequence, by processing the amplitude sequence in phase-coded signal sequence, namely, impact amplitude in amplitude sequence is positioned and revised, and amplitude sequence and phase sequence are reconsolidated, thus realize the shock resistance process to phase-coded signal sequence.In addition, the invention also discloses a kind of shock resistance system of phase-coded signal sequence.Method and system of the present invention not only solves the removal problem of impact signal, avoids the loss of useful signal simultaneously.
Description
Technical field
The present invention relates in general to signal transacting field, more specifically, relates to the shock resistance method and system of phase-coded signal sequence.
Background technology
Phase-coded signal is common in various electronic system, such as radar system, communication system etc.From the system acceptance such as receiver to phase-coded signal except comprising in advance except supposition content information, further comprises the impact signal that in RX path, various interference produces, as the electromagnetic interference etc. of the thunder and lightning signal in air, external environment.The normal process of these impact signals to phase-coded signal brings and has a strong impact on, and destroys good correlation when phase-coded signal carries out relevant treatment.Therefore, reduce the impact of impact signal on correlated process be ensure to comprise in signal information can the key of normal extraction.
For most of electronic communication system, the generation of impact signal is generally rendered as random appearance, duration extremely short, concentration of energy and far exceedes the characteristic of normal signal.After narrowband systems, impact signal has certain broadening in time.Scalping method after position judgment was often adopted to the process of impact signal in the past, and namely first detected the moment of impacting and occurring, then the signal in this moment is rejected.Current impact signal processing method has two classes: a class is carried out Time-Frequency Analysis Method discovery to adopt wavelet transformation and judges impact position, is then rejected; Another kind of, impact signal is approximately fractional order statistical model, as Alpha Stable distritation model, sets up signal transacting criterion and carry out processing representing the most.
(1) small wave converting method.Towering, Quan Taifan is built see document Lee, " higher-frequency radar suppresses the research and experiment of impulsive disturbance ", electronic letters, vol, 12 phases in 1999, the analytical method that described method have employed wavelet transformation detects the lightning impulse signal generation moment in radar echo signal, then carries out signal compensation according to the sea clutter signal period the strongest in echo.The detection method that the method adopts is more complicated, realizes difficulty.Compensate with noise signal during this external compensation and blocked useful signal information, the loss of effective information is larger.
(2) fractional order statistical model.See document S.Ambike, J.Ilow, D.Hatzinakos, " DetectionforBinaryTransmissioninaMixtureofGaussianNoisea ndImpulsiveNoiseModeledasanAlpha-StableProcess " .IEEESignalProcessingLetters, the 1994 Alpha stabilization process model inspection of mixed Gaussian noise and impact noise (in the binary transmissions), described method is using the statistical model of the statistical model of fractional order statistical model Alpha Stable distritation as impact noise, and produce matching treatment principle on this basis, carry out signal transacting.These class methods are based upon the matching treatment that fractional order statistical model carries out, and processing procedure has given tacit consent to the existence of impact signal, is a kind of lossy process.
To sum up, we need a kind of method that can not produce again a large amount of loss of signal while eliminating antagonism impact signal.
Summary of the invention
Impacting to solve in phase-coded signal the problem that data align regular signal process impact, the invention provides a kind of shock resistance method of phase-coded signal sequence, comprising the following steps: step one, by phase-coded signal sequence { s
ibe decomposed into amplitude sequence { A
iand phase sequence
phase-coded signal sequence and the pass between amplitude sequence and phase sequence are:
wherein, i is the numbering of phase-coded signal sequence; Step 2, definition minimal impact amplitude; Step 3, sweep amplitude sequence, selects the position exceeding minimal impact amplitude to form impact position information sequence { I
q, q is the numbering of impact position information sequence; Step 4, according to impact position information sequence definition impact information sequence { I
p, L
p, to position impact, wherein, I
pfor impact position, L
pfor impacting width; And step 5, the impact amplitude after positioning in amplitude sequence is revised, and the amplitude sequence and phase sequence that have modified impact amplitude are merged, to obtain the PSK sequence through revising.
Wherein, the first average m calculating amplitude sequence is comprised the steps: in step 2
a; In amplitude sequence, select to form new amplitude sequence more than the element of the first average; And calculate the second average of new amplitude sequence
by minimal impact amplitude definition be
wherein, β is amplitude factor given in advance.
Wherein, the definition impact information sequence in step 4 comprises the steps: sub-step one, is scan impact position information sequence 1, makes p=1, I from q
p=I
q, L
p=1; Sub-step two, makes difference DELTA=I
q+1-I
p; Sub-step three, if Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
q+1; And sub-step four, make q=q+1, and be back to sub-step two, until q=Q, wherein, w is minimal impact width given in advance.
Wherein, step 5 comprises the steps: phase-coded signal sequence { s
ibe modified to following sequence,
i
p≤ n≤I
p+ L
p-1, wherein, n is natural number.
Wherein, the step of correction takes linear difference process.
In addition, also proposed a kind of shock resistance system of phase-coded signal sequence, comprising: decomposing module, for by phase-coded signal sequence { s
ibe decomposed into amplitude sequence { A
iand phase sequence
phase-coded signal sequence and the pass between amplitude sequence and phase sequence are:
wherein, i is the numbering of phase-coded signal sequence; Definition module, for defining minimal impact amplitude; Impact position information sequence composition module, for by sweep amplitude sequence, selects the position exceeding minimal impact amplitude to form impact position information sequence { I
q, q is the numbering of impact position information sequence; Im-pact location module, for defining impact information sequence { I according to impact position information sequence
p, L
p, to position impact, wherein, I
pfor impact position, L
pfor impacting width; And correcting module, for revising the impact amplitude after positioning in amplitude sequence, and the amplitude sequence and phase sequence that have modified impact amplitude are merged, to obtain the PSK sequence through revising.
Other features and advantages of the present invention will be set forth in the following description, and, partly become apparent from specification, or understand by implementing the present invention.Object of the present invention and other advantages realize by structure specifically noted in write specification, claims and accompanying drawing and obtain.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, and form a application's part, schematic description and description of the present invention, for explaining the present invention, does not form inappropriate limitation of the present invention.In the accompanying drawings:
Fig. 1 shows the block diagram of overall thought of the present invention;
Fig. 2 shows the flow chart of the shock resistance method of phase-coded signal sequence according to an exemplary embodiment of the present invention;
Fig. 3 shows the block diagram of the shock resistance system of phase-coded signal sequence according to an exemplary embodiment of the present invention;
Fig. 4 show use method according to an exemplary embodiment of the present invention process before and after the contrast simulation figure of signal in time domain; And
Fig. 5 shows the contrast simulation figure of the matching treatment result using method according to an exemplary embodiment of the present invention.
Embodiment
First, Fig. 1 shows overall thought of the present invention, that is, carry out process to the amplitude sequence in phase encoding pulse sequence and go shock treatment, be combined by treated amplitude sequence more afterwards with phase sequence, form new coded pulse sequence.
Embodiments of the invention are described in detail below in conjunction with accompanying drawing.
The shock resistance processing method of a kind of phase encoding pulse sequence described in present embodiment, as shown in Figure 2, its concrete steps are:
The first step (S201): by the phase-coded signal sequence { s containing impact data
i, i=0 ... N-1, gets the amplitude A of each element of sequence
iand phase place
be expressed as amplitude sequence { A
i, i=0 ... N-1 and phase sequence
i=0 ... N-1, triadic relation is
Second step (S203): the average calculating amplitude sequence is m
a, according to amplitude factor β given in advance, select the element not exceeding amplitude to form new sequence
0≤j≤N-1, the average calculating new sequence is
minimal impact amplitude definition is
3rd step (S205): with minimal impact amplitude scan amplitude sequence, record exceedes the position of minimal impact amplitude, forms positional information sequence { I
q, q=1 ... Q.
4th step (S207): definition impact information sequence { I
p, L
p, p=1 ... P.According to minimal impact width w given in advance, the positional information in scanning position information sequence.Impact width and be defined as impact position interval.Scanning process is 1. q scanning position information sequence { I from 1
q, q=1 ... Q, p=1, I
p=I
q, L
p=1; 2. Δ=I
q+1-I
p; If 3. Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
q+1; 4. q=q+1; 5. repeat 2.-4. until q=Q.
5th step (S209): according to impact information sequence { I
p, L
p, p=1 ... every a pair information in P, by phase-coded signal sequence { s
i, i=0 ... correspondence position data correction in N-1 is following data:
I
p≤n≤I
p+L
p-1。
Preferably, the amplitude sequence and the phase sequence that calculate original series is needed in the first step.
Preferably, the minimal impact amplitude described in second step comes from amplitude factor given in advance, thus ensure that the robustness of minimal impact amplitude.
Preferably, in the 3rd step with minimal impact amplitude for standard, scanning original data sequence produces impact position information sequence.
Preferably, in the 4th step, impact information sequence comes from the judgement of impact position sequential element.The definition of the minimal impact width of predefined constrains the scope of impacting and covering when occurring.
Preferably, in the 5th step in original data sequence impact position data correcting method in take the method that amplitude and phase place process respectively.
Preferably, in the 5th step in original data sequence impact position data correcting method in take amplitude interpolation and use the method for original phase formation new data.
In addition, also proposed a kind of shock resistance system of phase-coded signal sequence, as shown in Figure 3, comprising: decomposing module 301, for by phase-coded signal sequence { s
ibe decomposed into amplitude sequence { A
iand phase sequence
phase-coded signal sequence and the pass between amplitude sequence and phase sequence are:
wherein, i is the numbering of phase-coded signal sequence; Definition module 303, for defining minimal impact amplitude; Impact position information sequence composition module 305, for by sweep amplitude sequence, selects the position exceeding minimal impact amplitude to form impact position information sequence { I
q, q is the numbering of impact position information sequence; Im-pact location module 307, for defining impact information sequence { I according to impact position information sequence
p, L
p, to position impact, wherein, I
pfor impact position, L
pfor impacting width; And correcting module 309, for revising the impact amplitude after positioning in amplitude sequence, and the amplitude sequence and phase sequence that have modified impact amplitude are merged, to obtain the PSK sequence through revising.
Fig. 4 show use method according to an exemplary embodiment of the present invention process before and after the contrast simulation figure of signal in time domain.As can be seen from Figure 4, the signal amplitude in some sequences not carrying out processing is very large, needs to revise.And revised waveform almost can't see any projection.
Fig. 5 shows the contrast simulation figure of the matching treatment result using method according to an exemplary embodiment of the present invention.Wherein, represent that the matching result containing impact signal is present in the most areas above transverse axis, this represents, when there is impact in signal, signal needs to revise.And after use method of the present invention is revised, the matching result of new method is almost identical with the matching result without impact signal.Visible, the shock resistance method and system of the phase-coded signal sequence that the present invention proposes has good elimination and impacts dysgenic effect.
Method described in present embodiment, can be generalized in the shock resistance processing procedure of all phase-coded signals, the Signal Compression of coherent radar, Clutter suppression algorithm, clutter can be applied to disappear mutually the design etc. of the design of processor, MTI and pulse Doppler filter, also can be applicable to equipment performance evaluation, denoising imaging, electronic countermeasures, electronic jamming etc.Method described in present embodiment can be applied to Processing of Phase Coded Signals field.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (6)
1. a shock resistance method for phase-coded signal sequence, is characterized in that, comprise the following steps:
Step one, by phase-coded signal sequence { s
ibe decomposed into amplitude sequence { A
iand phase sequence
described phase-coded signal sequence and the pass between described amplitude sequence and described phase sequence are:
wherein, i is the numbering of described phase-coded signal sequence;
Step 2, definition minimal impact amplitude, wherein, definition minimal impact amplitude comprises: the first average calculating amplitude sequence is m
a, according to amplitude factor β given in advance, select the element not exceeding amplitude to form new sequence
the second average calculating new sequence is
minimal impact amplitude definition is
Step 3, scans described amplitude sequence, selects the position exceeding described minimal impact amplitude to form impact position information sequence { I
q, q is the numbering of described impact position information sequence, q=1 ... Q;
Step 4, according to described impact position information sequence definition impact information sequence { I
p, L
p, to position impact, wherein, I
pfor impact position, L
pfor impacting width; And
Step 5, impact amplitude after positioning in described amplitude sequence is revised, and the amplitude sequence and described phase sequence that have modified impact amplitude are merged, to obtain the PSK sequence through revising, wherein, carry out correction to the impact amplitude after positioning in described amplitude sequence to comprise: according to impact information sequence { I
p, L
p, every a pair information in p=1...P, by phase-coded signal sequence { s
i, the correspondence position data correction in i=0...N-1 is following data:
i
p≤ n≤I
p+ L
p-1, wherein, n is natural number.
2. method according to claim 1, is characterized in that, the definition impact information sequence in described step 4 comprises the steps:
Sub-step one, is scan described impact position information sequence 1 from q, makes p=1, I
p=I
q, L
p=1;
Sub-step two, makes difference DELTA=I
q+1-I
p;
Sub-step three, if Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
q+1; And
Sub-step four, makes q=q+1, and is back to described sub-step two, until q=Q,
Wherein, w is minimal impact width given in advance.
3. method according to claim 1, is characterized in that, the step of described correction takes linear interpolation processing.
4. a shock resistance system for phase-coded signal sequence, is characterized in that, comprising:
Decomposing module, for by phase-coded signal sequence { s
ibe decomposed into amplitude sequence { A
iand phase sequence
described phase-coded signal sequence and the pass between described amplitude sequence and described phase sequence are:
wherein, i is the numbering of described phase-coded signal sequence;
Definition module, for defining minimal impact amplitude, wherein, definition minimal impact amplitude comprises: the first average calculating amplitude sequence is m
a, according to amplitude factor β given in advance, select the element not exceeding amplitude to form new sequence
the second average calculating new sequence is
minimal impact amplitude definition is
Impact position information sequence composition module, for by the described amplitude sequence of scanning, selects the position exceeding described minimal impact amplitude to form impact position information sequence { I
q, q is the numbering of described impact position information sequence, q=1 ... Q;
Im-pact location module, for defining impact information sequence { I according to described impact position information sequence
p, L
p, to position impact, wherein, I
pfor impact position, L
pfor impacting width; And
Correcting module, for revising the impact amplitude after positioning in described amplitude sequence, and the amplitude sequence and described phase sequence that have modified impact amplitude are merged, to obtain the PSK sequence through revising, wherein, carry out correction to the impact amplitude after positioning in described amplitude sequence to comprise: according to impact information sequence { I
p, L
p, every a pair information in p=1...P, by phase-coded signal sequence { s
i, the correspondence position data correction in i=0...N-1 is following data:
i
p≤ n≤I
p+ L
p-1, wherein, n is natural number.
5. system according to claim 4, is characterized in that, the definition impact information sequence in described im-pact location module comprises the steps:
Sub-step one, is scan described impact position information sequence 1 from q, makes p=1, I
p=I
q, L
p=1;
Sub-step two, makes difference DELTA=I
q+1-I
p;
Sub-step three, if Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
q+1; And
Sub-step four, makes q=q+1, and is back to described sub-step two, until q=Q,
Wherein, w is minimal impact width given in advance.
6. system according to claim 4, is characterized in that, described correcting module takes linear interpolation processing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010566201.XA CN102480299B (en) | 2010-11-24 | 2010-11-24 | The shock resistance method and system of phase-coded signal sequence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010566201.XA CN102480299B (en) | 2010-11-24 | 2010-11-24 | The shock resistance method and system of phase-coded signal sequence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102480299A CN102480299A (en) | 2012-05-30 |
| CN102480299B true CN102480299B (en) | 2016-04-06 |
Family
ID=46092802
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201010566201.XA Active CN102480299B (en) | 2010-11-24 | 2010-11-24 | The shock resistance method and system of phase-coded signal sequence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102480299B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107070827B (en) * | 2016-12-12 | 2020-06-26 | 上海大学 | Communication method for multi-system linear frequency modulation keying modulation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6795559B1 (en) * | 1999-12-22 | 2004-09-21 | Mitsubishi Denki Kabushiki Kaisha | Impulse noise reducer detecting impulse noise from an audio signal |
| CN101408616A (en) * | 2008-11-24 | 2009-04-15 | 江南大学 | Inverse synthetic aperture radar imaging distance aligning method applicable to low signal-noise ratio data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4151243B2 (en) * | 2001-07-26 | 2008-09-17 | 三菱電機株式会社 | Multipath noise removal method and removal apparatus, FM receiver |
-
2010
- 2010-11-24 CN CN201010566201.XA patent/CN102480299B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6795559B1 (en) * | 1999-12-22 | 2004-09-21 | Mitsubishi Denki Kabushiki Kaisha | Impulse noise reducer detecting impulse noise from an audio signal |
| CN101408616A (en) * | 2008-11-24 | 2009-04-15 | 江南大学 | Inverse synthetic aperture radar imaging distance aligning method applicable to low signal-noise ratio data |
Non-Patent Citations (1)
| Title |
|---|
| 脉冲噪声检测与滤除器;蔡荣太等;《电子器件》;20080414;第30卷(第6期);第2053-2056页,2060页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102480299A (en) | 2012-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Qi et al. | Detection and parameter estimation of multicomponent LFM signal based on the fractional Fourier transform | |
| CN105913393B (en) | A kind of adaptive wavelet threshold image de-noising method and device | |
| Presti et al. | GNSS signal acquisition in the presence of sign transition | |
| Bai et al. | Automatic detection and removal of high‐density impulse noises | |
| Wu et al. | Elimination of cross‐terms in the Wigner–Ville distribution of multi‐component LFM signals | |
| CN108267784A (en) | A kind of seismic signal random noise compression process method | |
| CN109871733A (en) | A kind of adaptive sea clutter signal antinoise method | |
| Chen et al. | Effective and adaptive algorithm for pepper‐and‐salt noise removal | |
| CN103995950A (en) | Wavelet coefficient partial discharge signal noise elimination method based on related space domain correction threshold values | |
| CN112949846A (en) | Method for constructing generated deception jamming signal suitable for direct sequence spread spectrum system | |
| Ali et al. | An improved method of audio denoising based on wavelet transform | |
| CN106019257B (en) | The interpolation method of feature when based on high-frequency ground wave radar Current Observations result sky | |
| CN102480299B (en) | The shock resistance method and system of phase-coded signal sequence | |
| CN113114408B (en) | Method and device for generating and evaluating generated deception jamming signal | |
| Bayat et al. | Loran phase code revisited for continuous wave interference cancellation | |
| Qiu et al. | Sparse and low‐rank matrix decomposition based micro‐motion target indication in through‐the‐wall radar | |
| Luo et al. | Approach for transient interference detection based on straight line extraction for high‐frequency sky‐wave radar | |
| CN114675253A (en) | Water surface micro-amplitude wave frequency estimation method and device, electronic equipment and storage medium | |
| Wang et al. | Radar emitter intrapulse signal blind sorting under modified wavelet denoising | |
| Chunrui et al. | SVD-based Method for Radio Frequency Interference Suppression Applied to SAR. | |
| Bing et al. | Evaluating the complexity degree of electromagnetic environment utilizing morphological covering based fractal dimension | |
| Zhang et al. | Image denoising using Bandelets and hidden Markov tree models | |
| Wang et al. | Denoising of chaotic signal using independent component analysis and empirical mode decomposition with circulate translating | |
| Wang et al. | Efficient TR‐TBD algorithm for slow‐moving weak multi‐targets in heavy clutter environment | |
| Yang et al. | Adaptive median threshold algorithm used in FDIS of DSSS receivers |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20170516 Address after: 602, gate 2, 150001 Building 1, Nangang District, Heilongjiang, Harbin Co-patentee after: Liu Yongtan Patentee after: Yang Qiang Co-patentee after: Zhang Ning Co-patentee after: Zhang Xinchao Co-patentee after: Wang Haihong Co-patentee after: Xi Kun Co-patentee after: Harbin Institute of Technology Asset Investment Management Co., Ltd. Address before: 338 box 150001, Harbin Institute of Technology, 92 West straight street, Heilongjiang, Harbin Patentee before: Harbin Institute of Technology |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20170802 Address after: 150001, room 1017, science and technology building, Harbin National University, 434 Post Street, Nangang District, Harbin, Heilongjiang Patentee after: Harbin science and Technology Co., Ltd. ray letter Address before: No. 2, building 602, building 1, Nangang District, Harbin Co-patentee before: Liu Yongtan Patentee before: Yang Qiang Co-patentee before: Zhang Ning Co-patentee before: Zhang Xinchao Co-patentee before: Wang Haihong Co-patentee before: Xi Kun Co-patentee before: Liu Yongtan Zhang Ning Zhang Xinchao Wang Haihong Harbin Institute of Technology Xi Kun asset investment management limited liability company |
|
| TR01 | Transfer of patent right |