CN102480299A - Impact-resisting method and impact-resisting system for phase-coded signal sequence - Google Patents
Impact-resisting method and impact-resisting system for phase-coded signal sequence Download PDFInfo
- Publication number
- CN102480299A CN102480299A CN201010566201XA CN201010566201A CN102480299A CN 102480299 A CN102480299 A CN 102480299A CN 201010566201X A CN201010566201X A CN 201010566201XA CN 201010566201 A CN201010566201 A CN 201010566201A CN 102480299 A CN102480299 A CN 102480299A
- Authority
- CN
- China
- Prior art keywords
- sequence
- impact
- amplitude
- 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.)
- Granted
Links
Images
Abstract
The invention discloses an impact-resisting method for a phase-coded signal sequence. According to the impact-resisting method, the impact-resisting processing of the phase-coded signal sequence is realized by processing an amplitude sequence in the phase-coded signal sequence, namely positioning and correcting impact amplitude in the amplitude sequence and recombining the amplitude sequence with a phase sequence. The invention also discloses an impact-resisting system for the phase-coded signal sequence. According to the impact-resisting method and the impact-resisting system disclosed by the invention, the problem of removal of an impact signal is solved and the loss of a useful signal is also avoided.
Description
Technical field
The present invention relates in general to the signal processing field, more specifically, relates to the shock resistance method and system of phase-coded signal sequence.
Background technology
Phase-coded signal is common in the various electronic systems, for example radar system, communication system or the like.The phase-coded signal that receives from systems such as receivers has also comprised the impact signal that various interference produce in the RX path, like the thunder and lightning signal in the atmosphere, electromagnetic interference of external environment or the like except comprising prior supposition content information.These impact signals have brought the normal process of phase-coded signal and have had a strong impact on, and have destroyed the good correlation when phase-coded signal is carried out relevant treatment.Therefore, reducing impact signal is to guarantee to comprise in the signal key that information can normal extraction to the influence of correlated process.
As far as most of electronic communication systems, the generation of impact signal generally be rendered as at random occur, the duration is extremely short, concentration of energy and far surpass the characteristic of normal signal.Through behind the narrowband systems, impact signal has certain broadening in time.The impact Signal Processing often adopted the scalping method after the position judgment in the past, promptly at first detected and impacted the moment of taking place, and then the signal in this moment was rejected.At present the impact signal processing method has two types: find and judge impact position to adopt wavelet transformation to carry out Time-Frequency Analysis Method for one type, reject then; Another kind ofly then be approximately the fractional order statistical model to impact signal, stablize distributed model, set up the signal processing criterion and handle representative the most like Alpha.
(1) small wave converting method.Build towering, Quan Taifan referring to document Lee; " higher-frequency radar suppresses to impact the research and experiment that disturbs "; Electronic letters, vol; 1999 12 phases, said method has adopted the analytical method of wavelet transformation to detect the lightning impulse signal generation moment in the radar echo signal, carries out signal compensation according to the strongest sea clutter signal period in the echo then.The detection method that this method adopts is complicated, realizes difficulty.Compensate the useful signal information of having blocked with noise signal during this external compensation, the loss of effective information is bigger.
(2) fractional order statistical model.Referring to document S.Ambike; J.Ilow; D.Hatzinakos; " Detection for Binary Transmission in a Mixture of Gaussian Noise and Impulsive Noise Modeled as an Alpha-Stable Process " .IEEE Signal Processing Letters, 1994 (the Alpha stabilization process model of mixed Gaussian noise and impact noise detects in the binary transmissions), said method is with the statistical model of the stable statistical model that distributes of fractional order statistical model Alpha as impact noise; And, carry out signal processing in this basis generation matching treatment principle.These class methods are to be based upon the matching treatment of carrying out on the fractional order statistical model, and processing procedure has been given tacit consent to the existence of impact signal, are a kind of lossy processing.
To sum up, we need a kind ofly can not produce the method for a large amount of losses of signal again when eliminating the antagonism impact signal.
Summary of the invention
To impact the problem that data align regular signal processing influence in the phase-coded signal in order solving, to the invention provides a kind of shock resistance method of phase-coded signal sequence, may further comprise the steps: step 1, with phase-coded signal sequence { s
iBe decomposed into amplitude sequence { A
iAnd phase sequence
Relation between phase-coded signal sequence and amplitude sequence and the phase sequence is:
Wherein, i is the numbering of phase-coded signal sequence; Step 2, definition minimal impact amplitude; Step 3, the sweep amplitude sequence, the position of selecting to surpass the minimal impact amplitude constitutes impact position information sequence { I
q, q is the numbering of impact position information sequence; Step 4 is according to impact position information sequence definition impact information sequence { I
p, L
p, position with impact, wherein, I
pBe impact position, L
pFor impacting width; And step 5, the impact amplitude through the location in the amplitude sequence is revised, and will have been revised the amplitude sequence and the phase sequence merging of the amplitude of impacting, to obtain phase code sequence through revising.
Wherein, comprise the steps: to calculate the first average m of amplitude sequence in the step 2
AIn the amplitude sequence, the element of selecting not surpass first average constitutes new amplitude sequence; And second average of calculating new amplitude sequence
The minimal impact amplitude is defined as
Wherein, β is given amplitude factor in advance.
Wherein, the definition impact information sequence in the step 4 comprises the steps: substep one, is 1 to begin to scan the impact position information sequence from q, makes p=1, I
p=I
q, L
p=1; Substep two makes difference DELTA=I
Q+1-I
pIf substep three is Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
Q+1And substep four, make q=q+1, and be back to substep two, up to q=Q, wherein, w is given in advance minimal impact width.
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 a natural number.
Wherein, the step of correction has been taked the linear difference processing.
In addition, also proposed a kind of shock resistance system of phase-coded signal sequence, having comprised: decomposing module is used for phase-coded signal sequence { s
iBe decomposed into amplitude sequence { A
iAnd phase sequence
Relation between phase-coded signal sequence and amplitude sequence and the phase sequence is:
Wherein, i is the numbering of phase-coded signal sequence; Definition module is used to define the minimal impact amplitude; Impact position information sequence composition module is used for through the sweep amplitude sequence, and the position of selecting to surpass the minimal impact amplitude constitutes impact position information sequence { I
q, q is the numbering of impact position information sequence; Impact locating module, be used for according to impact position information sequence definition impact information sequence { I
p, L
p, position with impact, wherein, I
pBe impact position, L
pFor impacting width; And correcting module, be used for the impact amplitude through the location of amplitude sequence is revised, and will have revised the amplitude sequence and the phase sequence merging of the amplitude of impacting, to obtain phase code sequence through revising.
Other features and advantages of the present invention will be set forth in specification subsequently, and, partly from specification, become apparent, perhaps understand by embodiment of the present invention.The object of the invention can be realized through the structure that in the specification of being write, claims and accompanying drawing, is particularly pointed out and obtained with other advantages.
Description of drawings
Accompanying drawing described herein is used to provide further understanding of the present invention, constitutes the application's a part, and illustrative examples of the present invention and explanation thereof are used to explain the present invention, do not constitute improper qualification 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 handle before with afterwards the contrast simulation figure of signal on time domain; And
Fig. 5 shows and uses the matching treatment result's of method contrast simulation figure according to an exemplary embodiment of the present invention.
Embodiment
At first, Fig. 1 shows overall thought of the present invention,, the amplitude sequence in the phase encoding pulse sequence is handled shock treatment that is, more treated amplitude sequence is combined with phase sequence afterwards, forms new coded pulse sequence.
To combine accompanying drawing to specify embodiments of the invention below.
The shock resistance processing method of the described a kind of phase encoding pulse sequence of this execution mode, as shown in Figure 2, its concrete steps are:
The first step (S201): will contain the phase-coded signal sequence { s that impacts data
i, i=0 ... N-1, the amplitude A of getting 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, the triadic relation does
Second step (S203): the average of calculating the amplitude sequence is m
A, according to given amplitude factor β in advance, the element of selecting to surpass amplitude constitutes new sequence
0≤j≤N-1, the average of calculating new sequence does
The minimal impact amplitude is defined as
The 3rd step (S205): with minimal impact amplitude sweep amplitude sequence, record surpasses the position of minimal impact amplitude, constitutes positional information sequence { I
q, q=1 ... Q.
The 4th step (S207): definition impact information sequence { I
p, L
p, p=1 ... P.According to given in advance minimal impact width w, the positional information in the scanning position information sequence.Impact width and be defined as impact position at interval.Scanning process for 1. q since 1 scanning position information sequence { I
q, q=1 ... Q, p=1, I
p=I
q, L
p=1; 2. Δ=I
Q+1-I
p3. if Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
Q+14. q=q+1; 5. repeat 2.-4. up to q=Q.
The 5th step (S209): according to impact information sequence { I
p, L
p, p=1 ... Among the P each is to information, with phase-coded signal sequence { s
i, i=0 ... Correspondence position data correction among the N-1 is following data:
Preferably, need to calculate the amplitude sequence and the phase sequence of original series in the first step.
Preferably, come from given amplitude factor in advance, thereby guaranteed the robustness of minimal impact amplitude in the minimal impact amplitude described in second step.
Preferably, in the 3rd step is standard with the minimal impact amplitude, and the scanning original data sequence produces the impact position information sequence.
Preferably, the impact information sequence comes from the judgement of impact position sequential element in the 4th step.The definition of the minimal impact width of predefined has retrained the scope that covers when taking place of impacting.
Preferably, the method for having taked amplitude and phase place to handle respectively in the data correcting method of impact position in the original data sequence in the 5th step.
Preferably, amplitude interpolation and the method for using original phase place structure new data have been taked in the 5th step in the original data sequence in the data correcting method of impact position.
In addition, also proposed a kind of shock resistance system of phase-coded signal sequence, as shown in Figure 3, comprising: decomposing module 301 is used for phase-coded signal sequence { s
iBe decomposed into amplitude sequence { A
iAnd phase sequence
Relation between phase-coded signal sequence and amplitude sequence and the phase sequence is:
Wherein, i is the numbering of phase-coded signal sequence; Definition module 303 is used to define the minimal impact amplitude; Impact position information sequence composition module 305 is used for through the sweep amplitude sequence, and the position of selecting to surpass the minimal impact amplitude constitutes impact position information sequence { I
q, q is the numbering of impact position information sequence; Impact locating module 307, be used for according to impact position information sequence definition impact information sequence { I
p, L
p, position with impact, wherein, I
pBe impact position, L
pFor impacting width; And correcting module 309, be used for the impact amplitude through the location of amplitude sequence is revised, and will have revised the amplitude sequence and the phase sequence merging of the amplitude of impacting, to obtain phase code sequence through revising.
Fig. 4 shows use and handles before and the contrast simulation figure of signal on time domain afterwards based on the method for exemplary embodiment of the present.As can be seen from Figure 4, the signal of handling amplitude on some sequences is very big, need revise.And revised waveform almost can't see any projection.
Fig. 5 shows and uses the matching treatment result's of method contrast simulation figure according to an exemplary embodiment of the present invention.Wherein, the matching result that expression contains impact signal is present in most zones of transverse axis top, this expression, and when having impact in the signal, signal need be revised.And after use method of the present invention was revised, the matching result of new method and the matching result of no impact signal were much at one.It is thus clear that the shock resistance method and system of the phase-coded signal sequence that the present invention proposes has good elimination and impacts dysgenic effect.
The described method of this execution mode; Can be generalized in the shock resistance processing procedure of all phase-coded signals; Signal Compression, Clutter suppression algorithm, the clutter that can be applied to coherent radar design, MTI and the pulse Doppler Filter Design etc. of processor that disappear mutually also can be applicable to equipment performance evaluation, denoising imaging, electronic countermeasures, electronic jamming etc.The described method of this execution mode can be applied to the Processing of Phase Coded Signals field.
The above is merely the preferred embodiments of the present invention, is not limited to the present invention, and for a person skilled in the art, the present invention can have various changes and variation.All within spirit of the present invention and principle, any modification of being done, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. the shock resistance method of a phase-coded signal sequence is characterized in that, may further comprise the steps:
Step 1 is with phase-coded signal sequence { s
iBe decomposed into amplitude sequence { A
iAnd phase sequence
Relation between said phase-coded signal sequence and said amplitude sequence and the said phase sequence is:
Wherein, i is the numbering of said phase-coded signal sequence;
Step 2, definition minimal impact amplitude;
Step 3 scans said amplitude sequence, and the position of selecting to surpass said minimal impact amplitude constitutes impact position information sequence { I
q, q is the numbering of said impact position information sequence;
Step 4 is according to said impact position information sequence definition impact information sequence { I
p, L
p, position with impact, wherein, I
pBe impact position, L
pFor impacting width; And
Step 5 is revised the impact amplitude through the location in the said amplitude sequence, and will revise the amplitude sequence and the merging of said phase sequence of the amplitude of impacting, to obtain the phase code sequence through revising.
2. method according to claim 1 is characterized in that, comprises the steps: in the said step 2
Calculate the first average m of said amplitude sequence
A
In said amplitude sequence, the element of selecting not surpass said first average constitutes new amplitude sequence; And
Calculating said new second mean amplitude sequence
The margin is defined as the minimum impact
Wherein, β is given amplitude factor in advance.
3. method according to claim 2 is characterized in that, the definition impact information sequence in the said step 4 comprises the steps:
Substep one is 1 to begin to scan said impact position information sequence from q, makes p=1, I
p=I
q, L
p=1;
Substep two makes difference DELTA=I
Q+1-I
p
If substep three is Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
Q+1And
Substep four makes q=q+1, and is back to said substep two, up to q=Q,
Wherein, w is given in advance minimal impact width.
4. method according to claim 3 is characterized in that said step 5 comprises the steps:
{ si} is modified to following sequence with the phase-coded signal sequence
5. method according to claim 1 is characterized in that the step of said correction has been taked the linear difference processing.
6. the shock resistance system of a phase-coded signal sequence is characterized in that, comprising:
Decomposing module is used for phase-coded signal sequence { s
iBe decomposed into amplitude sequence { A
iAnd phase sequence
Relation between said phase-coded signal sequence and said amplitude sequence and the said phase sequence is:
Wherein, i is the numbering of said phase-coded signal sequence;
Definition module is used to define the minimal impact amplitude;
Impact position information sequence composition module is used for through scanning said amplitude sequence, and the position of selecting to surpass said minimal impact amplitude constitutes impact position information sequence { I
q, q is the numbering of said impact position information sequence;
Impact locating module, be used for according to said impact position information sequence definition impact information sequence { I
p, L
p, position with impact, wherein, I
pBe impact position, L
pFor impacting width; And
Correcting module is used for the impact amplitude through the location of said amplitude sequence is revised, and will revise the amplitude sequence and the merging of said phase sequence of the amplitude of impacting, to obtain the phase code sequence through revising.
7. system according to claim 6 is characterized in that, said definition module comprises:
The first average module is used to calculate the first average m of said amplitude sequence
A
New amplitude sequence composition module is used for constituting new amplitude sequence at the element that said amplitude sequence is selected not surpass said first average; And
The minimal impact amplitude module, second average
that is used to calculate said new amplitude sequence is defined as said minimal impact amplitude
Wherein, β is given amplitude factor in advance.
8. system according to claim 7 is characterized in that, the definition impact information sequence in the said impact locating module comprises the steps:
Substep one is 1 to begin to scan said impact position information sequence from q, makes p=1, I
p=I
q, L
p=1;
Substep two makes difference DELTA=I
Q+1-I
p
If substep three is Δ≤w, then L
p=L
p+ Δ, otherwise p=p+1, I
p=I
Q+1And
Substep four makes q=q+1, and is back to said substep two, up to q=Q,
Wherein, w is given in advance minimal impact width.
9. system according to claim 8 is characterized in that, said correcting module is with phase-coded signal sequence { s
iBe modified to following sequence
10. system according to claim 6 is characterized in that said correcting module has been taked the linear difference 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 true CN102480299A (en) | 2012-05-30 |
| CN102480299B 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107070827A (en) * | 2016-12-12 | 2017-08-18 | 上海大学 | A kind of communication means of the linear frequency modulation keying modulation of multi-system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030022650A1 (en) * | 2001-07-26 | 2003-01-30 | Masayuki Tsuji | Multipath noise reduction method, multipath noise reducer, and FM receiver |
| 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 |
-
2010
- 2010-11-24 CN CN201010566201.XA patent/CN102480299B/en active Active
Patent Citations (3)
| 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 |
| US20030022650A1 (en) * | 2001-07-26 | 2003-01-30 | Masayuki Tsuji | Multipath noise reduction method, multipath noise reducer, and FM receiver |
| 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 |
|---|
| 蔡荣太等: "脉冲噪声检测与滤除器", 《电子器件》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107070827A (en) * | 2016-12-12 | 2017-08-18 | 上海大学 | A kind of communication means of the linear frequency modulation keying modulation of multi-system |
| CN107070827B (en) * | 2016-12-12 | 2020-06-26 | 上海大学 | Communication method for multi-system linear frequency modulation keying modulation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102480299B (en) | 2016-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105913393B (en) | A kind of adaptive wavelet threshold image de-noising method and device | |
| Qi et al. | Detection and parameter estimation of multicomponent LFM signal based on the fractional Fourier transform | |
| Bai et al. | Automatic detection and removal of high‐density impulse noises | |
| CN103135091A (en) | Adaptive impulsive noise elimination method of DOA (direction of arrival) estimation system | |
| Goring et al. | Estimation of extreme sea levels in a tide-dominated environment using short data records | |
| Chen et al. | Effective and adaptive algorithm for pepper‐and‐salt noise removal | |
| Wu et al. | Elimination of cross‐terms in the Wigner–Ville distribution of multi‐component LFM signals | |
| CN112949846B (en) | Method for constructing generated deception jamming signal suitable for direct sequence spread spectrum system | |
| CN103064063A (en) | Poly-phase code radar signal waveform automatic identification method based on continuous wave Doppler (CWD) feature | |
| Deng et al. | Hyper chaotic logistic phase coded signal and its sidelobe suppression | |
| CN107015207A (en) | Active pressing jamming classifying identification method based on FRFT domains peak value discrete feature | |
| Chunjie et al. | Synthetic algorithm for deinterleaving radar signals in a complex environment | |
| Ali et al. | An improved method of audio denoising based on wavelet transform | |
| CN107247933A (en) | FMCW laser spacings system difference frequency method for extracting signal in a kind of smoky environment | |
| CN113114408B (en) | Method and device for generating and evaluating generated deception jamming signal | |
| CN104215944A (en) | Radar sea clutter adaptive inhibition processing method based on parameter selection | |
| CN102480299B (en) | The shock resistance method and system of phase-coded signal sequence | |
| Hames et al. | Evolution of joint probability methods in coastal engineering practice in the UK | |
| CN106291470B (en) | The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar | |
| CN116184331A (en) | Anti-sorting signal design method and system based on heavy frequency sliding | |
| Shariff et al. | Modelling significant wave height data of North Sea: Rayleigh vs weibull distribution | |
| 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 | |
| Li et al. | Fractal-based weak target detection within sea clutter |
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 | ||
| 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 |