CN108732558B - Matched filtering calculation method based on copy signal segmentation - Google Patents

Matched filtering calculation method based on copy signal segmentation Download PDF

Info

Publication number
CN108732558B
CN108732558B CN201810250668.XA CN201810250668A CN108732558B CN 108732558 B CN108732558 B CN 108732558B CN 201810250668 A CN201810250668 A CN 201810250668A CN 108732558 B CN108732558 B CN 108732558B
Authority
CN
China
Prior art keywords
signal
copy
segment
processing
copy 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
Application number
CN201810250668.XA
Other languages
Chinese (zh)
Other versions
CN108732558A (en
Inventor
韩宁
黄舒夏
方世良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201810250668.XA priority Critical patent/CN108732558B/en
Publication of CN108732558A publication Critical patent/CN108732558A/en
Application granted granted Critical
Publication of CN108732558B publication Critical patent/CN108732558B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a matched filtering calculation method based on copy signal segmentation, which divides a long copy signal into a plurality of small segments; decomposing the received signal into a plurality of processing segments; matching filtering processing and splicing of the received signal processing small section and the copy signal small section are carried out; and correcting the splicing result to obtain a matched filtering result of the undistorted receiving signal and the copy signal. The invention carries out segmentation processing and reasonable splicing on the long copy signal on the basis of carrying out segmentation processing on the received signal, thereby effectively reducing the data storage space required in the matched filtering processing process.

Description

Matched filtering calculation method based on copy signal segmentation
Technical Field
The invention relates to a matched filtering calculation method, in particular to a matched filtering calculation method based on copy signal segmentation.
Background
Signal detection is an important task for detection devices such as sonar and radar. Classical detection theory states that the best receiver to detect a known signal in a white gaussian noise background is a matched filter, which is also the most basic receiver commonly used in many detection systems such as sonar. The response function of the matched filter is the delayed conjugate mirror waveform of the matched signal u (t), which is then also referred to as the replica signal or reference signal of the matched filter. For a common input signal v (t), the output of the matched filter is y (t) ═ Rvu(t0-t), wherein RvuIs the cross-correlation value of signal v (t) to u (t). t is t0When, y (t)0)＝Rvu(0) Thus, the matched filter acts as a cross-correlator which can calculate the cross-correlation function.
In practical applications, the length of the replica signal u (t) is generally equal to the pulse width of the transmitted signal in the detection system. In order to meet the requirement of quasi-real-time matched filtering processing on long-time externally received data v (t), v (t) is generally divided into small sections, matched filtering processing is carried out one by one, and then segmented processing results are spliced. Matched filtering is one of the best filters, and its output signal-to-noise ratio is maximized when the input signal has a particular waveform. On the other hand, in order to adapt to different detection requirements of sonar, the form of the transmitted signal of the detection system is gradually increased, the pulse width of the transmitted signal is also increased, and therefore the length of the copy signal in the matched filtering process is also increased correspondingly. Under the condition of long copy signals, if the matching filtering method of only segmenting the received signals is directly adopted for processing, very high requirements are put on the storage space in the processor, and even the realization cannot be realized.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects of the prior art, the invention provides a matched filtering calculation method based on copy signal segmentation, which can effectively reduce the storage space required in the matched filtering processing process under the condition of long copy signals.
The technical scheme is as follows: the invention relates to a matched filtering calculation method based on copy signal segmentation, which comprises the following steps:
(1) the data length is LxIs divided into small segments, each segment having a length of l, and the total number of segments is NcAt the end of said copy signal, complementary lxNc-LxZero, resulting in a second copy signal and representing each said second copy signal segment;
(2) the data length is LyIs decomposed into small segments using the data length L of the first received signalyDividing by the length l of each segment of the first copy signal and rounding up to NsSupplementing l zeros at the front end of the first received signal and l x (N) at the end of the first received signals+1)-LyZero counting, obtaining a second receiving signal and representing each second receiving signal subsection;
(3) performing matched filtering processing on the same second copy signal segment, and taking the first l values to obtain a processing signal of each segment, and sequentially arranging corresponding matrixes to obtain a first processing signal corresponding to each second copy signal segment;
(4) repeating the step (3) for each second copy signal segment, and correcting each output first processing signal to obtain a second processing signal;
(5) and adding the matrixes of the second processed signals to obtain a third processed signal, and respectively removing the data of the front end and the tail end of the third processed signal to obtain the matched filtering processing result of the first copied signal and the first received signal.
Preferably, in the step (1), the total number of segments is calculated by using the data length L of the first copy signalxDivided by the length l of each segment and rounded up.
Preferably, in step (1), the second copy signal is
Figure RE-GDA0001787112180000021
The second copy signal fragment decomposed is denoted Xi＝[x(i-1)×l+1,x(i-1)×l+2,…,xi×l]，i＝1,…,Nc。
Preferably, in step (2), the second received signal is represented as
Figure BDA0001607692810000022
The total number of segments is Ns+1, each second received signal fragment is denoted Yj＝[y(j-1)×l+1,y(j-1)×l+2,…,y(j+1)×l]Wherein j is 1, …, Ns+1, each segment length being 2 × l.
Preferably, in step (3), the processed signal of each small segment is represented as zj＝[zj1,zj2,…,zjl]Wherein j is 1, …, Ns+ 1; the first processed signal for the second copy signal fragments is
Figure BDA0001607692810000023
And is represented as
Figure BDA0001607692810000024
Preferably, in step (4), the method for modifying the first processed signal is interceptionPartial matrices in the first processed signal
Figure BDA0001607692810000031
And then i × l zeros are complemented to obtain
Figure BDA0001607692810000032
Has the advantages that: the invention carries out segmentation processing and reasonable splicing on the long copy signal on the basis of carrying out segmentation processing on the received signal, thereby effectively reducing the data storage space required in the matched filtering processing process.
Drawings
FIG. 1 is a flow chart of a method according to the present invention;
FIG. 2 is a diagram illustrating a received signal being decomposed into processing fragments according to the present invention;
FIG. 3 is a graph comparing the complete results of a conventional matched filtering method with the matched filtering method of the present invention;
fig. 4 is a graph comparing local results of the conventional matched filtering method and the matched filtering method of the present invention in a time period of 1.95s to 2.2 s.
Detailed Description
Example 1
The invention discloses a matched filtering calculation method based on copy signal segmentation, which mainly solves the problem of overlarge storage space required in the matched filtering processing process under the condition of long copy signals, and firstly, the data length is L as shown in figure 1xThe first copy signal of (1) is divided into small segments, each segment has a length of l, and the total segment has a number of NcAt the end of the copy signal, complement lxnc-LxZero, resulting in a second copy signal and representing each second copy signal fragment.
The data length is LyIs divided into small segments using the data length L of the first received signalyDividing by the length l of each segment of the first copy signal and rounding up to NsAt the front of the first received signal, the first zero is complemented, and at the end, the first zero is complemented (N)s+1)-LyZero, obtaining a second received signal and representing each of the second received signalsTwo receive signal segments.
And for the same second copy signal segment, performing matched filtering processing on the second received signal segment and the second copy signal segment, and taking the first l values to obtain a processing signal of each segment and sequentially arranging corresponding matrixes of the processing signals to obtain a first processing signal corresponding to each second copy signal segment.
And (4) repeating the step (3) for each second copy signal segment, and correcting each output first processing signal to obtain a second processing signal.
And adding the matrixes of the second processed signals to obtain a third processed signal, and respectively removing the data of the front end and the tail end of the third processed signal to obtain the matched filtering processing result of the first copied signal and the first received signal. Wherein the above-mentioned "matrices" are all matrices with a row number of 1.
The method comprises the following steps:
(1) the data length is LxCopy signal of
Figure RE-GDA0001787112180000041
Dividing the obtained product into a plurality of small sections, wherein the length of each section is l; by LxDivide by l and round up to NcThen, make up LXN at the end of the copy signal xc-LxZero to obtain new copy signal
Figure RE-GDA0001787112180000042
And is represented as
Figure RE-GDA0001787112180000043
Divided NcA small segment of the copy signal is denoted Xi＝[x(i-1)×l+1,x(i-1)×l+2,…,xi×l]，i＝1,…,Nc。
(2) The data length is LyReceived signal of
Figure BDA0001607692810000044
Is broken down into several treatment fragments, where Ly>Lx(ii) a By usingLyDivide by l and round up to NsThen, the front end of the received signal y is complemented by one zero, and the tail end of the received signal y is complemented by one (N)s+1)-LyZero to obtain new received signal
Figure BDA0001607692810000045
And re-represent it as
Figure BDA0001607692810000046
Decomposed Ns+1 treatment fragments denoted Yj＝[y(j-1)×1024+1,y(j-1)×1024+2,…,y(j+1)×1024]，j＝1,…,Ns+1, each segment length being 2 × l.
(3) For the same copy signal small section Xi＝[x(i-1)×l+1,x(i-1)×l+2,…,xi×l]Carry out YjAnd XiAnd obtaining z by taking the first l valuesj＝[zj1,zj2,…,zjl]，j＝1,…,Ns+ 1; to zjCombining to obtain a corresponding XiIs/are as follows
Figure BDA0001607692810000047
And is represented as
Figure BDA0001607692810000048
(4) For each copy signal small segment Xi＝[x(i-1)×l+1,x(i-1)×l+2,…,xi×l]The operation of step (3) is performed, i is 1, …, NcFor each output
Figure BDA0001607692810000049
Making a correction, i.e. selecting ZiIn (1)
Figure BDA00016076928100000410
And then i × l zeros are complemented to obtain
Figure BDA00016076928100000411
(5) Each A isiMatrix addition to obtain
Figure BDA00016076928100000412
And represents it as
Figure BDA00016076928100000413
Respectively removing one data at the front end and the tail end of the B to obtain
Figure BDA00016076928100000414
Namely the matching filtering processing result of the original receiving signal y and the original copy signal x.
Example 2
Assuming that a transmitting signal of the active sonar is a linear frequency modulation signal, the pulse width is 1s, the lower limit frequency is 375Hz, and the upper limit frequency is 425 Hz; the length of a received signal is 5s, the signal-to-noise ratio is 0dB, and a target echo signal appears at a position with the time of 2 s; the sampling rate of the system is 5 kHz. The matched filtering processing based on the copy signal segmentation is completed according to the following steps:
(1) the copy signal x is the transmission signal with a length Lx5000, dividing the obtained product into a plurality of small segments, and selecting 1024 lengths of each small segment; in order to make the length of the new copy signal an integer multiple of L, L is usedxDivide by l and round up to NcAt the end of x, the new copy signal is obtained by complementing 120 zeros 5
Figure BDA0001607692810000051
And 5 copy signal fragments are Xi＝[x(i-1)×1024+1,x(i-1)×1024+2,…,xi×1024]，i＝1,…,5。
(2) As shown in FIG. 2, the received signal y has a length Ly25000, it is broken down into several small pieces for processing, using LyDivide by l and round up to Ns25; in order to perform the same processing operation on each small segment in the subsequent steps without losing valid data, 1024 zeros are complemented at the front end of the received signal y, and 1624 zeros are complemented at the tail end of the received signal y to obtain a new received signal
Figure BDA0001607692810000052
And re-represent it as Y ═ Y1,y2,…,y1024×27]And 26 processing fragments are denoted as Yj＝[y(j-1)×1024+1,y(j-1)×1024+2,…,y(j+1)×1024]J is 1, …,26, each segment length being 2048.
(3) For the same copy signal small section XiCarry out YjAnd XiThe first 1024 values are taken to obtain zj＝[zj1,zj2,…,zjl]J is 1, …, 26; to zjAre combined to obtain a compound corresponding to XiZ of (A)i＝[z1,z2,…,z26]And is represented by Zi＝[z1,z2,…,z26×1024]。
(4) For each copy signal small segment XiThe operation of step (3) is performed, i is 1, …,5, for each output ZiMaking a correction, i.e. selecting ZiElement [ z ] of (1)(i-1)×1024+1,z(i-1)×1024+2,…,z26×1024]And then i × 1024 zeros are complemented to obtain
Figure BDA0001607692810000053
(5) Each A isiMatrix addition, i is 1, …,5, resulting in B being a1+A2+…+A5And it is represented as B ═ B1,b2,…,b27×1024](ii) a Respectively removing 1024 data at the front end and the tail end of the B to obtain
Figure BDA0001607692810000061
Namely the matching filtering processing result of the original receiving signal y and the original copy signal x.
As can be seen from fig. 3 and 4, the matched filtering method based on the copy signal segmentation can obtain the matched filtering result of the received signal and the copy signal without distortion. In the matched filtering calculation process based on copy signal segmentation, only every small segment of copy signal needs to be stored, and the storage space can be reused and is 1024 data lengths; only every small segment of the received signal needs to be stored and the storage space can also be multiplexed, 2048 data length. When the conventional matched filtering method is used for calculation, the storage space of the copy signal needs 5000 data lengths, and the storage space of the segmented received signal is generally 10000 data lengths. Therefore, under the condition of a long copy signal, the matched filtering method based on the copy signal segmentation can effectively reduce the required storage space in the processing process.

Claims (6)

1. A method for matched filter computation based on replica signal segmentation, the method comprising the steps of:
(1) the data length is LxIs divided into small segments, each segment having a length of l, and the total number of segments is NcAt the end of said copy signal, complementary lxNc-LxZero, resulting in a second copy signal and representing each said second copy signal segment;
(2) the data length is LyIs decomposed into small segments using the data length L of the first received signalyDividing by the length l of each segment of the first copy signal and rounding up to NsSupplementing l zeros at the front end of the first received signal and l x (N) at the end of the first received signals+1)-LyZero, resulting in a second received signal and representing each of said second received signal segments;
(3) performing matched filtering processing on the same second copy signal segment, and taking the first l numerical values to obtain a processing signal of each segment, and sequentially arranging corresponding matrixes to obtain a first processing signal corresponding to each second copy signal segment;
(4) repeating the step (3) for each second copy signal segment, and correcting each output first processing signal to obtain a second processing signal;
(5) and adding the matrixes of the second processed signals to obtain a third processed signal, and respectively removing the data of the front end and the tail end of the third processed signal to obtain the matched filtering processing result of the first copied signal and the first received signal.
2. The method of claim 1, wherein in step (1), the total number of segments is calculated by using the data length L of the first copy signalxDivided by the length l of each segment and rounded up.
3. The method of claim 1, wherein in step (1), the second copy signal is
Figure FDA0003005241920000011
The second copy signal fragment decomposed is denoted Xi＝[x(i-1)×l+1,x(i-1)×l+2,…,xi×l]，i＝1,…,Nc。
4. The method of claim 1, wherein in step (2), the second received signal is represented as
Figure FDA0003005241920000012
The total number of stages is Ns+1, each second received signal fragment is denoted Yj＝[y(j-1)×l+1,y(j-1)×l+2,…,y(j+1)×l]Wherein j is 1, …, Ns+1, each segment length being 2 × l.
5. The method of claim 1, wherein in step (3), the processed signal of each segment is represented as zj＝[zj1,zj2,…,zjl]Wherein j is 1, …, Ns+ 1; the first processed signal for the second copy signal fragments is
Figure FDA0003005241920000021
And is represented as
Figure FDA0003005241920000022
6. The method of claim 1, wherein in step (4), the modifying the first processed signal is performed by truncating a partial matrix in the first processed signal
Figure FDA0003005241920000023
And then i × l zeros are complemented to obtain
Figure FDA0003005241920000024
CN201810250668.XA 2018-03-26 2018-03-26 Matched filtering calculation method based on copy signal segmentation Active CN108732558B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810250668.XA CN108732558B (en) 2018-03-26 2018-03-26 Matched filtering calculation method based on copy signal segmentation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810250668.XA CN108732558B (en) 2018-03-26 2018-03-26 Matched filtering calculation method based on copy signal segmentation

Publications (2)

Publication Number Publication Date
CN108732558A CN108732558A (en) 2018-11-02
CN108732558B true CN108732558B (en) 2021-08-24

Family

ID=63941003

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810250668.XA Active CN108732558B (en) 2018-03-26 2018-03-26 Matched filtering calculation method based on copy signal segmentation

Country Status (1)

Country Link
CN (1) CN108732558B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109567788B (en) * 2018-11-29 2021-08-20 武汉中旗生物医疗电子有限公司 Electrocardiosignal filtering method for removing ringing

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213759A (en) * 2011-04-08 2011-10-12 东南大学 Characteristic matching method of underground water target based on power spectrum
CN102353952A (en) * 2011-06-03 2012-02-15 哈尔滨工程大学 Line spectrum detection method by coherent accumulation of frequency domains
CN103733583A (en) * 2011-08-05 2014-04-16 瑞典爱立信有限公司 Reference signal generation technique
CN104113472A (en) * 2013-04-18 2014-10-22 晨星半导体股份有限公司 Wireless receiving system and signal processing method thereof
CN104462088A (en) * 2013-09-13 2015-03-25 中国银联股份有限公司 Mass text matching method
CN106383669A (en) * 2016-09-19 2017-02-08 渤海大学 Distributed storage method and system based on (n,k,m) coding
CN106772471A (en) * 2016-11-16 2017-05-31 中国人民解放军国防科学技术大学 A kind of long code segmentation overlay local correlation catching method based on GPU
CN107271980A (en) * 2017-07-17 2017-10-20 中国人民解放军国防科学技术大学 A kind of segmented matched filter processing method to interval modulated signal

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5965715B2 (en) * 2012-04-27 2016-08-10 古野電気株式会社 Ultrasonic transmission / reception apparatus, ultrasonic transmission / reception method, and ultrasonic transmission / reception program
US9678975B2 (en) * 2013-03-15 2017-06-13 International Business Machines Corporation Reducing digest storage consumption in a data deduplication system

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213759A (en) * 2011-04-08 2011-10-12 东南大学 Characteristic matching method of underground water target based on power spectrum
CN102353952A (en) * 2011-06-03 2012-02-15 哈尔滨工程大学 Line spectrum detection method by coherent accumulation of frequency domains
CN103733583A (en) * 2011-08-05 2014-04-16 瑞典爱立信有限公司 Reference signal generation technique
CN104113472A (en) * 2013-04-18 2014-10-22 晨星半导体股份有限公司 Wireless receiving system and signal processing method thereof
CN104462088A (en) * 2013-09-13 2015-03-25 中国银联股份有限公司 Mass text matching method
CN106383669A (en) * 2016-09-19 2017-02-08 渤海大学 Distributed storage method and system based on (n,k,m) coding
CN106772471A (en) * 2016-11-16 2017-05-31 中国人民解放军国防科学技术大学 A kind of long code segmentation overlay local correlation catching method based on GPU
CN107271980A (en) * 2017-07-17 2017-10-20 中国人民解放军国防科学技术大学 A kind of segmented matched filter processing method to interval modulated signal

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WSNs Microseismic Signal Subsection Compression Algorithm;Zhouzhou Liu等;《Mathematical Problems in Engineering》;20150531;第2015卷;1-9页 *
基于FPGA的高速高阶FIR滤波器的频域改进方法;郭晓伟等;《现代电子技术》;20160601;第39卷(第11期);55-62页 *

Also Published As

Publication number Publication date
CN108732558A (en) 2018-11-02

Similar Documents

Publication Publication Date Title
EP1514347B1 (en) Digital filtering method
CN111044999B (en) Method, equipment and system for reducing interference in frequency modulation continuous wave radar unit
US8891708B2 (en) Obtaining a desired non-zero phase shift using forward-backward filtering
CN101061744B (en) Audio signal enhancement
EP1640973A3 (en) Audio signal processing apparatus and method
DE102013000380A1 (en) Method and sensor system for the acoustic measurement of the characteristics of a transmission path of a measuring system between loudspeaker and microphone
JPH08505741A (en) Signal processing system
JP2951619B2 (en) Noise reduction method and circuit
CN108549078B (en) Cross-channel combination and detection method for radar pulse signals
CN102165709B (en) Signal processing method, signal processing device, and signal processing program
CN106898361A (en) Single channel blind source separation method based on feedback variation Mode Decomposition
CN108732558B (en) Matched filtering calculation method based on copy signal segmentation
EP0814920A1 (en) Adaptive optimisation process for ultrasonic measurement signals
US5900006A (en) Median filtering method and apparatus using a plurality of processing elements
EP1879292B1 (en) Partitioned fast convolution
KR20070029176A (en) Method and system for enhancing the sharpness of a video signal
US6816597B1 (en) Pseudo stereophonic device
KR20210112714A (en) System for background music
WO2013039595A1 (en) Methods and systems for determining an enhanced rank order value of a data set
Tesei et al. The asymmetric generalized gaussian function: A new hos-based model for generic noise pdfs
US6760372B1 (en) Adaptive signal processor using an eye-diagram metric
CN1684454A (en) Apparatus and method for echo indicator generation
US5528555A (en) System and method for compensating for towed array motion induced errors
CN111028857B (en) Method and system for reducing noise of multichannel audio-video conference based on deep learning
US10825443B2 (en) Method and system for implementing a modal processor

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