CN103743403A - Calculating method of single-frequency signal reverberation center - Google Patents
Calculating method of single-frequency signal reverberation center Download PDFInfo
- Publication number
- CN103743403A CN103743403A CN201410029611.9A CN201410029611A CN103743403A CN 103743403 A CN103743403 A CN 103743403A CN 201410029611 A CN201410029611 A CN 201410029611A CN 103743403 A CN103743403 A CN 103743403A
- Authority
- CN
- China
- Prior art keywords
- reverberation
- center
- frame
- frequency
- calculation
- 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.)
- Pending
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Abstract
The invention discloses a calculating method of a single-frequency signal reverberation center. The method comprises the processing processes: (a), calculating the frequency of a reverberation center by adopting a gravity center method according to a frequency response result; (b), smoothening the frequency of the reverberation center of the current frame by adopting a time smoothening technology according to the calculation result of the frequency of the reverberation center of the former frame, wherein the smoothened frequency is used as a final calculation result. The calculating method has the beneficial effects that the calculation of the single-frequency signal reverberation center of a homing system is completed, the single-frequency signal reverberation center can be estimated accurately by adopting the method, and the navigation speed deviation, the movement posture change and the ocean current influence can be adapted, so that the false alarm probability is reduced and the detection performance of single-frequency signals is improved; by the adoption of the gravity center method and the time smoothening technology, according to the method, the convenience in calculation is achieved, and the requirement on the detection on the single-frequency signals of an XX homing system can be met.
Description
Technical field
The computing method that the present invention relates to a kind of simple signal reverberation center, belong to self-conductance field of engineering technology, be mainly used in calculating simple signal reverberation center and improve detection performance under real boat condition.
Background technology
In homing system, simple signal is one of topmost signal form.In single-frequency signal detection process, reverberation is the key factor that impact detects performance.At present, the interval that homing system mainly arranges certain Doppler spread is based on experience value reverberation district, under real boat state, due to the deviation of the speed of a ship or plane and setting speed and the impact of ocean current and state of motion, there is fluctuation in reverberation center, thereby has increased the false-alarm probability of system and reduced detection performance.There is no at present the reverberation center calculation method based on gravity model appoach.
Summary of the invention
The object of the invention is: propose first a kind of reverberation center calculation method based on gravity model appoach, this scheme have calculate easy, result is accurate, adaptable advantage, and by reality navigate checking can accurately estimate reverberation center, there is very high practical value and application prospect.
Technical scheme of the present invention is: the computing method at a kind of simple signal reverberation center, and main methods is: a. adopts gravity model appoach to calculate reverberation centre frequency according to frequency response results; B. utilize the reverberation centre frequency result of calculation of previous frame, adopt time smoothing technology to carry out smoothly the reverberation centre frequency of present frame, and as final calculation result.
The method specifically comprises the steps:
Step 1: gravity model appoach is determined reverberation centre frequency
According to the frequency response P (k, i) of i frame data, utilize gravity model appoach to calculate the centre frequency of reverberation, k=0 wherein ..., N-1.I frame reverberation centre frequency Φ (i) computing formula is:
Step 2: reverberation centre frequency is carried out to time smoothing
Since L frame, reverberation center is carried out smoothly, can improving like this stability of reverberation center calculation, specific formula for calculation is:
In formula, the weight that w is i frame, should meet 0<w<1.By regulating weight can reduce the impact of target echo on reverberation center calculation.F (i) is the final calculation result of i frame reverberation centre frequency.
The computing method at the simple signal reverberation center that the present invention proposes have following significant advantage and the innovation that is different from traditional reverberation district method:
Feature 1: proposed first the simple signal reverberation center calculation method based on gravity model appoach
According to the frequency response of reverberation, adopt gravity model appoach to calculate reverberation center, the reverberation center of causing for thunder speed deviation, state of motion and ocean current fluctuation has certain adaptability;
Feature 2: adopt time smoothing technology to improve stability and the accuracy of reverberation center calculation
In real boat, in the single work period of homing system, do not carry out motor-driven, therefore in the single work period of homing system, in the short time, reverberation center is basicly stable, adopt time smoothing technology can improve stability and the accuracy of reverberation center calculation, and can reduce echoed signal to calculating the impact at reverberation center.
The present invention relates to the computing method at a kind of simple signal reverberation center, the method is on simple signal frequency domain detection basis, adopt gravity model appoach and time smoothing technique computes reverberation center, and there is certain adaptability for thunder speed deviation, state of motion variation and ocean current in real boat, not only can accurately estimate reverberation center, and there is higher stability.The method, can practical requirement through repeatedly real boat checking.
The invention has the beneficial effects as follows: the present invention proposes a kind of computing method of the simple signal reverberation center based on gravity model appoach and time smoothing technology:
1) the present invention has completed the calculating at homing system simple signal reverberation center, adopt the method can accurately estimate single-frequency signal reverberation center, and can adapt to that speed of a ship or plane deviation, state of motion change and ocean current impact, thereby reduce false-alarm probability and improve the detection performance of simple signal;
2) the present invention adopts gravity model appoach and time smoothing technology, and the method is calculated easy, can meet the needs of homing system single-frequency signal detection.
Accompanying drawing explanation
Fig. 1: the processing flow chart of the computing method at a kind of simple signal reverberation center
Embodiment
Now in conjunction with the embodiments, the invention will be further described for accompanying drawing:
The computing method at the simple signal reverberation center in the present embodiment comprise the steps:
Step 1: gravity model appoach is determined reverberation centre frequency
According to the frequency response P (k, i) of i frame data, utilize gravity model appoach to calculate the centre frequency of reverberation, k=0 wherein ..., N-1.I frame reverberation centre frequency Φ (i) computing formula is:
In the present embodiment, N=20, the value of each frame rate response P is as shown in table 1, according to above formula, can be calculated the value of reverberation centre frequency Φ (i) as shown in table 2.
The value of table 1 frequency response P
The 1st frame | The 2nd frame | The 3rd frame | The 4th frame | The 5th frame | The 6th frame | The 7th frame | The 8th frame | The 9th frame | The 10th frame |
4 | 4 | 5 | 23 | 15 | 7 | 8 | 2 | 17 | 2 |
9 | 1 | 13 | 6 | 17 | 12 | 14 | 20 | 12 | 0 |
17 | 10 | 11 | 2 | 20 | 3 | 4 | 0 | 23 | 3 |
6 | 1 | 10 | 9 | 26 | 13 | 5 | 3 | 15 | 35 |
21 | 3 | 19 | 0 | 9 | 6 | 15 | 10 | 1 | 2 |
5 | 5 | 13 | 12 | 2 | 4 | 17 | 16 | 4 | 16 |
16 | 8 | 1 | 3 | 10 | 6 | 2 | 1 | 7 | 19 |
3 | 0 | 4 | 12 | 12 | 1 | 4 | 5 | 8 | 6 |
5 | 19 | 4 | 10 | 5 | 11 | 8 | 10 | 7 | 7 |
3 | 4 | 11 | 7 | 15 | 1 | 4 | 8 | 5 | 26 |
4 | 6 | 31 | 13 | 17 | 11 | 7 | 3 | 8 | 5 |
4 | 10 | 6 | 3 | 5 | 9 | 7 | 12 | 4 | 2 |
4 | 6 | 3 | 9 | 1 | 2 | 15 | 0 | 8 | 8 |
6 | 8 | 2 | 3 | 20 | 2 | 5 | 5 | 5 | 11 |
6 | 5 | 4 | 5 | 8 | 9 | 1 | 1 | 5 | 18 |
8 | 9 | 15 | 5 | 5 | 1 | 8 | 3 | 6 | 5 |
19 | 2 | 8 | 6 | 19 | 12 | 7 | 2 | 7 | 11 |
3 | 2 | 11 | 12 | 4 | 3 | 12 | 11 | 7 | 6 |
2 | 6 | 5 | 6 | 4 | 3 | 9 | 4 | 11 | 4 |
7 | 1 | 13 | 0 | 12 | 23 | 13 | 7 | 7 | 17 |
Table 2 reverberation centre frequency Φ (i) result of calculation
The 1st frame | The 2nd frame | The 3rd frame | The 4th frame | The 5th frame | The 6th frame | The 7th frame | The 8th frame | The 9th frame | The 10th frame |
9.5 | 10.4636 | 10.2169 | 9.2534 | 9.2124 | 10.8921 | 10.5333 | 9.5528 | 8.982 | 10.6355 |
Step 2: reverberation centre frequency is carried out to time smoothing
Since L frame, reverberation center is carried out smoothly, can improving like this stability of reverberation center calculation, specific formula for calculation is:
In formula, the weight that w is i frame, should meet 0<w<1.By regulating weight can reduce the impact of target echo on reverberation center calculation.F (i) is the final calculation result of i frame reverberation centre frequency.
In the present embodiment, L=5, w=0.9.Finally calculate to such an extent that F (i) is worth as shown in table 3.
The final calculation result of table 3 reverberation centre frequency
The 1st frame | The 2nd frame | The 3rd frame | The 4th frame | The 5th frame | The 6th frame | The 7th frame | The 8th frame | The 9th frame | The 10th frame |
9.5000 | 10.4636 | 10.2169 | 9.2534 | 9.7398 | 10.0973 | 9.9566 | 9.7234 | 9.7145 | 10.0130 |
Claims (1)
1. the computing method at simple signal reverberation center, is characterized in that, comprise the steps:
Step 1: gravity model appoach is determined reverberation centre frequency
According to the frequency response P (k, i) of i frame data, utilize gravity model appoach to calculate the centre frequency of reverberation, k=0 wherein ..., N-1.I frame reverberation centre frequency Φ (i) computing formula is:
Step 2: reverberation centre frequency is carried out to time smoothing
Since L frame, reverberation center is carried out smoothly, specific formula for calculation is:
In formula, the weight that w is i frame, should meet 0<w<1.F (i) is the final calculation result of i frame reverberation centre frequency.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410029611.9A CN103743403A (en) | 2014-01-22 | 2014-01-22 | Calculating method of single-frequency signal reverberation center |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410029611.9A CN103743403A (en) | 2014-01-22 | 2014-01-22 | Calculating method of single-frequency signal reverberation center |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103743403A true CN103743403A (en) | 2014-04-23 |
Family
ID=50500446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410029611.9A Pending CN103743403A (en) | 2014-01-22 | 2014-01-22 | Calculating method of single-frequency signal reverberation center |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103743403A (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1343120A (en) * | 1970-12-23 | 1974-01-10 | North American Rockwell | Focusing control of synthetic aperture processing for side- looking radar |
US4053886A (en) * | 1976-11-19 | 1977-10-11 | The United States Of America As Represented By The Secretary Of The Navy | Stepped dual-frequency, ocean-wave spectrometer |
US5418536A (en) * | 1981-12-21 | 1995-05-23 | Westinghouse Electric Corporation | Bandwidth and amplitude insensitive frequency discriminator |
JP3455188B2 (en) * | 2001-02-06 | 2003-10-14 | 株式会社エム・シー・シー | Target detection method and target detection device |
CN101738601A (en) * | 2009-11-27 | 2010-06-16 | 西安电子科技大学 | System and method for measuring speed of locomotive based on radar near field echo power spectrum characteristics |
CN101833035A (en) * | 2010-04-19 | 2010-09-15 | 天津大学 | Linear frequency-modulated parameter estimating method and implementing device thereof |
CN102624418A (en) * | 2012-03-14 | 2012-08-01 | 东南大学 | Hydroacoustic biphase modulation direct sequence spread spectrum signal carrier frequency estimation method |
CN103248591A (en) * | 2013-05-27 | 2013-08-14 | 清华大学 | Coarse frequency offset estimation method based on frequency spectrum barycenter |
-
2014
- 2014-01-22 CN CN201410029611.9A patent/CN103743403A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1343120A (en) * | 1970-12-23 | 1974-01-10 | North American Rockwell | Focusing control of synthetic aperture processing for side- looking radar |
US4053886A (en) * | 1976-11-19 | 1977-10-11 | The United States Of America As Represented By The Secretary Of The Navy | Stepped dual-frequency, ocean-wave spectrometer |
US5418536A (en) * | 1981-12-21 | 1995-05-23 | Westinghouse Electric Corporation | Bandwidth and amplitude insensitive frequency discriminator |
JP3455188B2 (en) * | 2001-02-06 | 2003-10-14 | 株式会社エム・シー・シー | Target detection method and target detection device |
CN101738601A (en) * | 2009-11-27 | 2010-06-16 | 西安电子科技大学 | System and method for measuring speed of locomotive based on radar near field echo power spectrum characteristics |
CN101833035A (en) * | 2010-04-19 | 2010-09-15 | 天津大学 | Linear frequency-modulated parameter estimating method and implementing device thereof |
CN102624418A (en) * | 2012-03-14 | 2012-08-01 | 东南大学 | Hydroacoustic biphase modulation direct sequence spread spectrum signal carrier frequency estimation method |
CN103248591A (en) * | 2013-05-27 | 2013-08-14 | 清华大学 | Coarse frequency offset estimation method based on frequency spectrum barycenter |
Non-Patent Citations (1)
Title |
---|
裴益轩等: "滑动平均法的基本原理及应用", 《火炮发射与控制学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
PH12018550059A1 (en) | Systems and methods for recommending an estimated time of arrival | |
CN106844935B (en) | Large-damping engineering structure modal parameter identification method | |
WO2020131187A3 (en) | Ocean weather forecasting system | |
GB201305116D0 (en) | Earth model estimation through an acoustic full waveform inversion of seismic data | |
AU2017419266A1 (en) | Methods and systems for estimating time of arrival | |
CN104076355A (en) | Method for conducting before-detection tracking on weak and small target in strong-clutter environment based on dynamic planning | |
CN105738880A (en) | Moment estimation method for reverse inverse gauss texture composite compound gauss sea clutter amplitude distributed parameters | |
CN103324083A (en) | Non-linear ship motion control method based on robust observer | |
CN104182609A (en) | Decorrelation based three-dimensional target tracking method for unbiased converted measurement | |
CN103913736A (en) | Laser micro Doppler parameter estimation method based on spectrogram rearrangement | |
CN107607092A (en) | A kind of wave based on unmanned aerial vehicle remote sensing and floating motion forecast system | |
CN103323815A (en) | Underwater acoustic locating method based on equivalent sound velocity | |
CN105607042A (en) | Method for locating sound source through microphone array time delay estimation | |
CN106896363A (en) | A kind of submarine target active tracing track initiation method | |
WO2016198093A9 (en) | Determining of model parameters for positioning purposes | |
CN102455423A (en) | Method for eliminating sound reflection interference in ultrasonic location | |
WO2019130246A3 (en) | Location determination using acoustic models | |
CN104215939A (en) | Knowledge assisted space-time adaptive processing method integrating generalized symmetrical structure information | |
CN104614554A (en) | Self-amending method of standard error of ship-based wind speed and direction transducer | |
CN101876585B (en) | ICA (Independent Component Analysis) shrinkage de-noising method evaluating noise variance based on wavelet packet | |
CN103837859A (en) | Method for calculating position of sound emission source through wavelet neural network of mixed leapfrog | |
CN104865602A (en) | Method of determining epicenter location and earthquake start time through multiple data sources | |
CN110677140A (en) | Random system filter containing unknown input and non-Gaussian measurement noise | |
CN110398775A (en) | Tunnel is dashed forward discharge disaster microseismic event signal fluctuation first break pickup method and system | |
Shi et al. | Active heave compensation prediction research for deep sea homework crane based on KPSO-SVR |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140423 |