CN113466838B - Iterative compensation target radiation noise data simulation method and system - Google Patents
Iterative compensation target radiation noise data simulation method and system Download PDFInfo
- Publication number
- CN113466838B CN113466838B CN202110583185.3A CN202110583185A CN113466838B CN 113466838 B CN113466838 B CN 113466838B CN 202110583185 A CN202110583185 A CN 202110583185A CN 113466838 B CN113466838 B CN 113466838B
- Authority
- CN
- China
- Prior art keywords
- data
- radiation noise
- target radiation
- simulation data
- spectrum
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 168
- 238000004088 simulation Methods 0.000 title claims abstract description 135
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000001228 spectrum Methods 0.000 claims abstract description 98
- 238000012545 processing Methods 0.000 claims abstract description 51
- 238000012937 correction Methods 0.000 claims abstract description 21
- 238000010606 normalization Methods 0.000 claims abstract description 21
- 230000009466 transformation Effects 0.000 claims description 20
- 238000005070 sampling Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 241000251729 Elasmobranchii Species 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
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/52001—Auxiliary means for detecting or identifying sonar signals or the like, e.g. sonar jamming signals
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- 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/52004—Means for monitoring or calibrating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention relates to the field of sonar signal processing, in particular to an iterative compensation target radiation noise data simulation method and system, wherein the method comprises the following steps: in the processing frequency band, performing spectrum structure correction on the noise data subjected to the power spectrum normalization processing according to a target radiation noise spectrum structure to be simulated to obtain target radiation noise simulation data; and transmitting and receiving the target radiation noise simulation data for multiple times in sequence, compensating the currently transmitted target radiation noise simulation data by adopting the data received each time, and obtaining the target radiation noise data consistent with the target radiation noise spectrum structure to be simulated through multiple iterative compensation. The method of the invention well solves the problem that the spectrum structure of the analog data output by the transducer has distortion.
Description
Technical Field
The invention relates to the field of sonar signal processing, in particular to an iterative compensation target radiation noise data simulation method and system.
Background
In the field of underwater acoustic signal processing, radiation noise signals of targets in seawater are always the hot points of research of people, the radiation noise signals are the information sources of passive sonar equipment, and the information can be used for detecting, orienting, positioning, tracking, identifying and the like of the targets. In the research of fields such as passive detection and identification of targets, analog simulation research of noise signals has important practical significance and operational significance, such as scheme demonstration of a command system, simulation experiments, analog training, testing of a sonar system in a laboratory stage and the like.
Especially in underwater sound countermeasure technology, the target sound characteristic realistic simulation technology is the key technology of the sound bait. The simulation of the target radiation noise mainly aims to research and analyze the acoustic characteristics of the target, and can simulate the characteristics with high fidelity for the evaluation of various acoustic systems and the purpose of luring torpedoes and sonars in underwater sound confrontation.
In the process of simulating the underwater acoustic target radiation noise data, the spectrum structure of the simulated data output by the transducer is affected by different frequency point responses of the transmitting system, and if the simulated target radiation noise data is compensated by only adopting the response data measured by the transducer 1/3 octave, the compensated data spectrum structure can not meet the requirement when the frequency band response difference of the transmitted data of the transducer is large.
Disclosure of Invention
The invention aims to solve the problem of distortion of a spectrum structure of analog data output by a transducer, and provides an iterative compensation target radiation noise data simulation method and system.
In order to achieve the above object, the present invention proposes an iterative compensation target radiation noise data simulation method, including:
in the processing frequency band, performing spectrum structure correction on the noise data subjected to the power spectrum normalization processing according to a target radiation noise spectrum structure to be simulated to obtain target radiation noise simulation data;
and transmitting and receiving the target radiation noise simulation data for multiple times in sequence, compensating the currently transmitted target radiation noise simulation data by adopting the data received each time, and obtaining the target radiation noise data consistent with the target radiation noise spectrum structure to be simulated through multiple iterative compensation.
As an improvement of the above method, the method specifically comprises:
step 1) calculating a spectrum filter H (f) of target radiation noise simulation data according to the shape of a target radiation noise spectrum SL (f) to be simulated in a processing frequency band, wherein f is 1, …, f issH (f) is in dB, and f represents the frequency domain;
H(f)＝[0,…,0,SL(fl),SL(fl+1),…,SL(fh),0,…,0,SL(fs-fh+1),SL(fs-fh+2),…,SL(fs-fl+1),0,…,0]
where SL (f) has the unit of dB,fl,fhrespectively SL (f) continuum start-stop frequency, fsThe system sampling frequency;
step 2), generating random noise data noise (t) by a white noise generator, and carrying out power spectrum normalization processing on the noise (t) to obtain target radiation noise simulation data reference data consig (t) as follows:
wherein std (-) is a standard deviation solving function;
step 3), performing time-frequency transformation on the target radiation noise simulation data reference data consig (t) to obtain corresponding frequency domain data fconsig (f), where f is 1, …, fs；
fconsig(f)＝FFT[consig(t)]
Wherein, FFT [. cndot ] is a fast Fourier transform function;
step 4), performing spectrum structure correction on the frequency domain data fconsig (f) of the target radiation noise simulation data through a spectrum filter H (f), and obtaining the frequency domain data mfconsig (f) of the target radiation noise simulation data after spectrum correction, wherein the frequency domain data mfconsig (f) is as follows:
step 5), performing time-frequency inverse transformation on mfconsig (f), and then taking a real part to obtain target radiation noise simulation data mfconsig (t) as follows:
mfconsig(t)＝IFFT(mfconsig(f))
wherein IFFT [ ] is an inverse fast Fourier transform function;
step 6) standard deviation normalization processing is carried out on the target radiation noise simulation data mfconcesig (t) to obtain target radiation noise simulation data targetsig with power of 00(t) is:
step 7) simulating a total sound source SL according to the target radiation noiseoutputFor targetsig0(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig (t) output by the signal source, wherein the target radiation noise simulation data targetsig (t) is as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
step 8) setting a compensation filter H0(f),f＝1,…,fsComprises the following steps:
H0(f)＝[0,…,0,0,0,…,0,0,…,0,0,0,…,0,0,…,0]；
step 9), distributing a standard hydrophone at a fixed distance from the transmitting transducer, and collecting data received by the standard hydrophone as sample data sig (t);
step 10), performing time-frequency transformation on sample data sig (t) to obtain frequency domain data fsig (f) of the sample data, where f is 1, …, and f iss：
fsig(f)＝FFT[sig(t)]；
Step 11) using frequency domain data fsig (f) of sample data, f is 1, …, fsAnd target radiation noise simulation data frequency domain data mfconsig (f), calculating a compensation filter H0(f),f＝1,…,fs，H0(f) Has a unit of dB, and satisfies the following formula:
wherein lg (·) represents a base logarithm of 10;
step 12) applying a compensation filter H0(f),f＝1,…,fsPerforming spectrum structure compensation on the target radiation noise simulation data frequency domain data mfconsig (f), and obtaining target radiation noise simulation data frequency domain data mfconsig1(f) after the spectrum compensation at this time:
step 13) performing time-frequency inverse transformation on mfconsig1(f), and taking a real part to obtain time domain data mfconsig1(t) of the target radiation noise simulation data, wherein the time domain data mfconsig1(t) is as follows:
mfconsig1(t)＝IFFT(mfconsig1(f))；
step 14) carrying out standard deviation normalization processing on time domain data mfconsig1(t) of the target radiation noise simulation data to obtain equal target radiation noise simulation data targetsig1 with power of 00(t) is:
step 15) simulating data total sound source SL according to target radiation noiseoutputFor targetsig10(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig1(t) output by the signal source as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
step 16) repeatedly executing the step 9) to the step 15) to obtain target radiation noise simulation data targetsig1(t) compensated by the iteration;
step 17) comparing whether the spectrum structure of the targetsig1(t) is consistent with the structure of the target radiation noise spectrum SL (f) to be simulated; if not, continue to step 16) until the two structures are consistent.
An iteratively compensated target radiated noise data modeling system, the system comprising: the system comprises a spectrum structure correction module, an iteration output module and a compensation module; wherein,
the spectrum structure correction module is used for correcting the spectrum structure of the noise data subjected to the power spectrum normalization processing according to the target radiation noise spectrum structure to be simulated in the processing frequency band to obtain target radiation noise simulation data;
the iterative output module is used for calling the compensation module for multiple times to carry out iterative compensation until target radiation noise data consistent with a target radiation noise spectrum structure to be simulated are obtained;
and the compensation module is used for transmitting and receiving the target radiation noise simulation data and compensating the currently transmitted target radiation noise simulation data by adopting the received data.
As an improvement of the above system, the specific processing procedure of the spectrum structure modification module includes:
within the processing frequency band, a spectrum filter H (f) for calculating target radiation noise simulation data according to the shape of a target radiation noise spectrum SL (f) to be simulated, wherein f is 1, … and fsH (f) is in dB, and f represents the frequency domain;
H(f)＝[0,…,0,SL(fl),SL(fl+1),…,SL(fh),0,…,0,SL(fs-fh+1),SL(fs-fh+2),…,SL(fs-fl+1),0,…,0]
wherein the unit of SL (f) is dB, fl,fhRespectively SL (f) continuum start-stop frequency, fsThe system sampling frequency;
generating random noise data noise (t) by a white noise generator, and carrying out power spectrum normalization processing on the noise (t) to obtain target radiation noise simulation data reference data consig (t) as follows:
wherein std (-) is a standard deviation solving function;
performing time-frequency transformation on the target radiation noise simulation data reference data consig (t) to obtain corresponding frequency domain data fconsig (f), wherein f is 1, …, fs；
fconsig(f)＝FFT[consig(t)]
Wherein, FFT [. cndot ] is a fast Fourier transform function;
performing spectrum structure correction on the frequency domain data fconsig (f) of the target radiation noise simulation data through a spectrum filter H (f), and obtaining the frequency domain data mfconsig (f) of the target radiation noise simulation data after spectrum correction, wherein the frequency domain data mfconsig (f) is as follows:
performing time-frequency inverse transformation on mfconsig (f), and then taking a real part to obtain target radiation noise simulation data mfconsig (t) as follows:
mfconsig(t)＝IFFT(mfconsig(f))
where IFFT [ ] is an inverse fast Fourier transform function.
As an improvement of the above system, the specific processing procedure of the compensation module includes:
carrying out standard deviation normalization processing on the target radiation noise simulation data mfconsig (t) to obtain target radiation noise simulation data targetsig with power of 00(t) is:
simulating data total sound source SL according to target radiation noiseoutputFor targetsig0(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig (t) output by the signal source, wherein the target radiation noise simulation data targetsig (t) is as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
setting a compensation filter H0(f),f＝1,…,fsComprises the following steps:
H0(f)＝[0,…,0,0,0,…,0,0,…,0,0,0,…,0,0,…,0]；
distributing a standard hydrophone at a fixed distance from the transmitting transducer, and acquiring data received by the standard hydrophone as sample data sig (t);
performing time-frequency transformation on sample data sig (t) to obtain frequency domain data fsig (f) of the sample data, where f is 1, …, and f iss：
fsig(f)＝FFT[sig(t)]；
Frequency domain data fsig (f) using sample data, f 1, …, fsAnd target radiation noise simulation data frequency domain data mfconsig (f), calculating a compensation filter H0(f),f＝1,…,fs，H0(f) Has a unit of dB, and satisfies the following formula:
wherein lg (·) represents a base logarithm of 10;
using compensating filters H0(f),f＝1,…,fsPerforming spectrum structure compensation on the target radiation noise simulation data frequency domain data mfconsig (f), and obtaining target radiation noise simulation data frequency domain data mfconsig1(f) after the spectrum compensation at this time:
and (3) performing time-frequency inverse transformation on the mfconsig1(f), and taking a real part to obtain time domain data mfconsig1(t) of the target radiation noise simulation data, wherein the time domain data mfconsig1(t) is as follows:
mfconsig1(t)＝IFFT(mfconsig1(f))；
the time domain data mfconsig1(t) of the target radiation noise simulation data is subjected to standard deviation normalization processing to obtain equal target radiation noise simulation data targetsig1 with the power of 00(t) is:
simulating data total sound source SL according to target radiation noiseoutputFor targetsig10(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig1(t) output by the signal source：
Wherein SLnormThe unimodal 1V signal corresponds to the sound source level of the system output signal.
Obtaining target radiation noise simulation data targetsig1(t) after the iteration compensation;
as an improvement of the above system, the specific processing procedure of the iterative output module includes:
comparing whether the spectrum structure of the target radiation noise simulation data targetsig1(t) output by the current compensation module is consistent with the structure of a target radiation noise spectrum SL (f) to be simulated or not; if not, the compensation module is called to continue the iterative compensation process until targetsig1(t) is consistent with the SL (f) structure
Compared with the prior art, the invention has the advantages that:
the method of the invention utilizes the correlation between the transducer transmitting data sample and the known data sample, realizes the analog data compensation by transmitting and receiving the required target radiation noise analog data for multiple times and adopting the compensation mode of the data received each time, improves the consistency of the output analog data spectrum structure of the transducer after compensation and the required spectrum structure, and well solves the problem of distortion of the output analog data spectrum structure of the transducer.
Drawings
FIG. 1 is a schematic structural diagram of an example embodiment of the present invention;
FIG. 2 is a block diagram of a method implementation of the present invention;
FIG. 3 is target radiated noise simulation data for a desired spectral structure;
FIG. 4 is data after the transducer transmits analog data that is received by a standard hydrophone;
FIG. 5 is the data after the transducer transmits 1 compensated analog data received by a standard hydrophone;
FIG. 6 is data after 2 times of compensated analog data transmission by the transducer is received by a standard hydrophone;
fig. 7 is data after 3 times of compensated analog data transmission by the transducer is received by a standard hydrophone.
Detailed Description
The problem of distortion of the spectrum structure of the analog data output by the transducer due to different influences of different frequency point responses of a transmitting system is solved. According to the invention, the analog data compensation is realized by transmitting and receiving the required target radiation noise analog data twice and adopting the received data to compensate according to the correlation between the transducer transmission data sample and the known data sample, so that the consistency of the spectrum structure of the analog data output by the transducer after compensation and the required spectrum structure is improved. Finally, the effectiveness and feasibility of the method are further verified through example analysis.
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and examples.
Example 1
Before describing the process of the present invention in detail, a description will first be given of an apparatus to which the process of the present invention is applied. Fig. 1 is a schematic structural diagram of a target radiation noise data simulation device, which comprises 6 parts, namely a program control computer 1, a digital-to-analog conversion device 2, a power amplification device 3, a transducer 4, a standard hydrophone 5 and a collection device 6.
Target radiation noise data simulation technology
The simulation process of the broadband continuous spectrum is mainly realized by a white noise generator, a correction filter, a compensation filter and the like, and the realization block diagram is shown in fig. 2.
White gaussian noise simulation
The simulation method of the target radiation noise is to utilize the characteristics of the noise, firstly generate white noise and then realize the white noise through a spectrum filter. In order to simulate the target radiated noise data, white gaussian noise is first simulated.
According to the central limit theorem, if R0,R1,......,RN-1N random variables which are mutually independent and distributed in the same way are adopted, so that when N is larger, the random variables are larger
Wherein μ ═ E [ R ]k]，σ2＝Var[Rk]K is 0, 1. At this time, N independent [0,1 ] s are added]And adding the uniformly distributed random variables to obtain the random variables which are approximately in Gaussian distribution.
Let RN (k) be [0,1]In a uniformly distributed sequence of independent random variables, then μ ═ E [ rn (k)]＝0.5，σ2＝Var[RN(k)]1/12, making a random variable sequence according to the above formula
Wherein k is 0, 1. The sequence is a Gaussian white noise sequence with the length M, the mean value of 0 and the variance of 1. Here, the Gaussian distribution is guaranteed by N (N ≧ 30), and the independence is guaranteed by the combined sum in the intercept RN (k) that does not overlap each other.
Spectral filter design
The design of the spectrum filter can adopt two methods of FIR filter and linear prediction, which is the key of the target radiation noise.
FIR frequency sampling method
And designing a corresponding FIR filter by adopting a frequency sampling method according to the specific noise frequency characteristic. Setting a desired frequency spectrum with certain frequency characteristics as H (f), then according to the frequency sequence f1,f2,...,fi(i ═ 1,2,. N) and the amplitude sequence H (f)i) The coefficients of the spectral filter can be obtained.
Linear prediction method
The least square error forward linear prediction method is one of the methods for designing the spectrum filter, and the input signal sequence is set as x (n), the order of the linear prediction filter is p, and the coefficient is ai1,2, p, the desired output of linear prediction is
The filter coefficients can be given based on the output error e (n) with the smallest variance.
Target radiation noise data simulation method
Within the processing frequency band, according to the simulated target radiation noise spectrum shape (continuous spectrum start and stop frequency f)l,fhAnd) calculating a target radiation noise simulation data spectrum filter H (f), wherein f is 1, … and fs(in dB).
H(f)＝[0,…,0,SL(fl),SL(fl+1),…,SL(fh),0,…,0,SL(fs-fh+1),SL(fs-fh+2),…,SL(fs-fl+1),0,…,0] (3)
Where SL (f) (in dB), fsThe system sampling frequency.
Generating random noise data noise (t) by adopting a formula (2), and carrying out power spectrum normalization processing on the noise (t) to obtain target radiation noise simulation data reference data consig (t).
In the formula, std (. cndot.) is a function of the standard deviation.
Performing time-frequency transformation on the target radiation noise simulation data reference data consig (t) to obtain corresponding frequency domain data fconsig (f), wherein f is 1, …, fs。
fconsig(f)＝FFT[consig(t)] (5)
Wherein FFT [. cndot. ] is a fast Fourier transform function.
And (3) carrying out spectrum structure correction on the target radiation noise simulation data frequency domain data fconsig (f) to obtain the target radiation noise simulation data frequency domain data mfconsig (f) after spectrum correction.
And (f) performing time-frequency inverse transformation on the mfconsig (f), and taking a real part of the mfconsig (f) to obtain target radiation noise simulation data mfconsig (t).
mfconsig(t)＝IFFT(mfconsig(f)) (7)
In the formula, IFFT [. cndot. ] is an inverse fast Fourier transform function.
Performing standard deviation normalization processing on the target radiation noise simulation data mfconsig (t) to obtain target radiation noise simulation data targetsig with power of 0 and the like0(t)。
Simulating data total sound source SL according to target radiation noiseoutputAnd (5) carrying out amplitude compensation processing on the targetsig (t) to obtain target radiation noise simulation data output by the signal source.
In the formula, SLnormThe unimodal 1V signal corresponds to the sound source level of the system output signal.
Setting a compensation filter H0(f),f＝1,…,fs(in dB) is:
H0(f)＝[0,…,0,0,0,…,0,0,…,0,0,0,…,0,0,…,0] (10)
and a standard hydrophone is arranged at the position 1m of the transmitting transducer, and the acquisition hydrophone receives the data transmitted by the data transducer as sample data sig (t).
Performing time-frequency transformation on the sample data sig (t) to obtain corresponding frequency domain data fsig (f), where f is 1, …, and fs。
fsig(f)＝FFT[sig(t)] (11)
Using sample data frequency domain data fsig (f), f is 1, …, fsAnd target radiation noise simulation data frequency domain data mfconsig (f), calculating a target radiation noise simulation data spectrum compensation filter H0(f),f＝1,…,fs(in dB).
To obtain H0(f) The unit is dB:
lg (-) is base logarithm of 10.
Using compensating filters H0(f),f＝1,…,fs(unit is dB), spectrum structure compensation is carried out on the target radiation noise simulation data frequency domain data mfconsig (f), and target radiation noise simulation data frequency domain data mfconsig1(f) after spectrum compensation is obtained.
And (3) performing time-frequency inverse transformation on mfconsig1(f), and taking a real part of the mfconsig1(f) to obtain target radiation noise simulation data mfconsig (1 t).
mfconsig1(t)＝IFFT(mfconsig1(f)) (15)
And (3) performing standard deviation normalization processing on the target radiation noise simulation data mfconsig1(t) to obtain target radiation noise simulation data targetsig1(t) with the power of 0 and the like.
Simulating data total sound source SL according to target radiation noiseoutputAnd (5) carrying out amplitude compensation processing on the targetsig (t) to obtain target radiation noise simulation data targetsig1(t) output by the signal source.
In the formula, SLnormThe unimodal 1V signal corresponds to the sound source level of the system output signal.
And updating the sample data sig (t), and repeatedly executing the equations (11) to (17), wherein the target radiation noise simulation data targetsig1(t) can be obtained through iterative compensation.
Example 2
Embodiment 2 of the present invention proposes an iterative compensation target radiation noise data simulation system, which includes: the system comprises a spectrum structure correction module, an iteration output module and a compensation module; wherein,
the spectrum structure correction module is used for correcting the spectrum structure of the noise data subjected to the power spectrum normalization processing according to the target radiation noise spectrum structure to be simulated in the processing frequency band to obtain target radiation noise simulation data;
the iterative output module is used for calling the compensation module for multiple times to carry out iterative compensation until target radiation noise data consistent with a target radiation noise spectrum structure to be simulated are obtained;
and the compensation module is used for transmitting and receiving the target radiation noise simulation data and compensating the currently transmitted target radiation noise simulation data by adopting the received data.
Example analysis
The effects of the method of the present invention are compared below with examples.
In order to further verify that the method can effectively obtain the target radiation noise simulation data of the required spectrum structure, the following numerical simulation analysis is carried out.
FIG. 3 is target radiated noise simulation data for a desired spectral structure;
FIG. 4 is data after the transducer transmits analog data that is received by a standard hydrophone;
FIG. 5 is the data after the transducer transmits 1 compensated analog data received by a standard hydrophone;
FIG. 6 is data after 2 times of compensated analog data transmission by the transducer is received by a standard hydrophone;
fig. 7 is data after 3 times of compensated analog data transmission by the transducer is received by a standard hydrophone.
According to simulation results, after 3 times of compensation, the simulation technology emitted by the transducer is target radiation noise simulation data of a required spectrum structure, and therefore the target radiation noise data simulation can be achieved.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention and are not limited. Although the present invention has been described in detail with reference to the embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (2)
1. A method of iteratively compensated target radiated noise data simulation, the method comprising:
in the processing frequency band, performing spectrum structure correction on the noise data subjected to the power spectrum normalization processing according to a target radiation noise spectrum structure to be simulated to obtain target radiation noise simulation data;
the target radiation noise simulation data are transmitted for multiple times and received respectively, the currently transmitted target radiation noise simulation data are compensated by adopting the data received each time, and the target radiation noise data consistent with the target radiation noise spectrum structure to be simulated are obtained through multiple iterative compensation;
the method specifically comprises the following steps:
step 1) calculating a spectrum filter H (f) of target radiation noise simulation data according to the shape of a target radiation noise spectrum SL (f) to be simulated in a processing frequency band, wherein f is 1, …, f issH (f) is in dB, and f represents the frequency domain;
H(f)＝[0,…,0,SL(fl),SL(fl+1),…,SL(fh),0,…,0,SL(fs-fh+1),SL(fs-fh+2),…,SL(fs-fl+1),0,…,0]
wherein the unit of SL (f) is dB, fl,fhRespectively SL (f) continuum start-stop frequency, fsThe system sampling frequency;
step 2), generating random noise data noise (t) by a white noise generator, and carrying out power spectrum normalization processing on the noise (t) to obtain target radiation noise simulation data reference data consig (t) as follows:
wherein std (-) is a standard deviation solving function;
step 3), performing time-frequency transformation on the target radiation noise simulation data reference data consig (t) to obtain corresponding frequency domain data fconsig (f), where f is 1, …, fs；
fconsig(f)＝FFT[consig(t)]
Wherein, FFT [. cndot ] is a fast Fourier transform function;
step 4), performing spectrum structure correction on the frequency domain data fconsig (f) of the target radiation noise simulation data through a spectrum filter H (f), and obtaining the frequency domain data mfconsig (f) of the target radiation noise simulation data after spectrum correction, wherein the frequency domain data mfconsig (f) is as follows:
step 5), performing time-frequency inverse transformation on mfconsig (f), and then taking a real part to obtain target radiation noise simulation data mfconsig (t) as follows:
mfconsig(t)＝IFFT(mfconsig(f))
wherein IFFT [ ] is an inverse fast Fourier transform function;
step 6) standard deviation normalization processing is carried out on the target radiation noise simulation data mfconcesig (t) to obtain target radiation noise simulation data targetsig with power of 00(t) is:
step 7) simulating a total sound source SL according to the target radiation noiseoutputFor targetsig0(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig (t) output by the signal source, wherein the target radiation noise simulation data targetsig (t) is as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
step 8) setting a compensation filter H0(f),f＝1,…,fsComprises the following steps:
H0(f)＝[0,…,0,0,0,…,0,0,…,0,0,0,…,0,0,…,0]；
step 9), distributing a standard hydrophone at a fixed distance from the transmitting transducer, and collecting data received by the standard hydrophone as sample data sig (t);
step 10), performing time-frequency transformation on sample data sig (t) to obtain frequency domain data fsig (f) of the sample data, where f is 1, …, and f iss：
fsig(f)＝FFT[sig(t)]；
Step 11) using frequency domain data fsig (f) of sample data, f is 1, …, fsAnd target radiation noise simulation data frequency domain data mfconsig (f), calculating a compensation filter H0(f),f＝1,…,fs，H0(f) Has a unit of dB, and satisfies the following formula:
wherein lg (·) represents a base logarithm of 10;
step 12) applying a compensation filter H0(f),f＝1,…,fsPerforming spectrum structure compensation on the target radiation noise simulation data frequency domain data mfconsig (f), and obtaining target radiation noise simulation data frequency domain data mfconsig1(f) after the spectrum compensation at this time:
step 13) performing time-frequency inverse transformation on mfconsig1(f), and taking a real part to obtain time domain data mfconsig1(t) of the target radiation noise simulation data, wherein the time domain data mfconsig1(t) is as follows:
mfconsig1(t)＝IFFT(mfconsig1(f))；
step 14) carrying out standard deviation normalization processing on time domain data mfconsig1(t) of the target radiation noise simulation data to obtain equal target radiation noise simulation data targetsig1 with power of 00(t) is:
step 15) simulating data total sound source SL according to target radiation noiseoutputFor targetsig10(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig1(t) output by the signal source as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
step 16) repeatedly executing the step 9) to the step 15) to obtain target radiation noise simulation data targetsig1(t) compensated by the iteration;
step 17) comparing whether the spectrum structure of the targetsig1(t) is consistent with the structure of the target radiation noise spectrum SL (f) to be simulated; if not, continue to step 16) until the two structures are consistent.
2. An iteratively compensated target radiated noise data modeling system, the system comprising: the system comprises a spectrum structure correction module, an iteration output module and a compensation module; wherein,
the spectrum structure correction module is used for correcting the spectrum structure of the noise data subjected to the power spectrum normalization processing according to the target radiation noise spectrum structure to be simulated in the processing frequency band to obtain target radiation noise simulation data;
the iterative output module is used for calling the compensation module for multiple times to carry out iterative compensation until target radiation noise data consistent with a target radiation noise spectrum structure to be simulated are obtained;
the compensation module is used for transmitting and receiving the target radiation noise simulation data and compensating the currently transmitted target radiation noise simulation data by adopting the received data;
the specific processing procedure of the spectrum structure correction module comprises the following steps:
within the processing frequency band, a spectrum filter H (f) for calculating target radiation noise simulation data according to the shape of a target radiation noise spectrum SL (f) to be simulated, wherein f is 1, … and fsH (f) is in dB, and f represents the frequency domain;
H(f)＝[0,…,0,SL(fl),SL(fl+1),…,SL(fh),0,…,0,SL(fs-fh+1),SL(fs-fh+2),…,SL(fs-fl+1),0,…,0]
wherein the unit of SL (f) is dB, fl,fhRespectively SL (f) continuum start-stop frequency, fsThe system sampling frequency;
generating random noise data noise (t) by a white noise generator, and carrying out power spectrum normalization processing on the noise (t) to obtain target radiation noise simulation data reference data consig (t) as follows:
wherein std (-) is a standard deviation solving function;
performing time-frequency transformation on the target radiation noise simulation data reference data consig (t) to obtain corresponding frequency domain data fconsig (f), wherein f is 1, …, fs；
fconsig(f)＝FFT[consig(t)]
Wherein, FFT [. cndot ] is a fast Fourier transform function;
performing spectrum structure correction on the frequency domain data fconsig (f) of the target radiation noise simulation data through a spectrum filter H (f), and obtaining the frequency domain data mfconsig (f) of the target radiation noise simulation data after spectrum correction, wherein the frequency domain data mfconsig (f) is as follows:
performing time-frequency inverse transformation on mfconsig (f), and then taking a real part to obtain target radiation noise simulation data mfconsig (t) as follows:
mfconsig(t)＝IFFT(mfconsig(f))
wherein IFFT [ ] is an inverse fast Fourier transform function;
the specific processing procedure of the compensation module comprises the following steps:
carrying out standard deviation normalization processing on the target radiation noise simulation data mfconsig (t) to obtain target radiation noise simulation data targetsig with power of 00(t) is:
simulating data total sound source SL according to target radiation noiseoutputFor targetsig0(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig (t) output by the signal source, wherein the target radiation noise simulation data targetsig (t) is as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
setting a compensation filter H0(f),f＝1,…,fsComprises the following steps:
H0(f)＝[0,…,0,0,0,…,0,0,…,0,0,0,…,0,0,…,0]；
distributing a standard hydrophone at a fixed distance from the transmitting transducer, and acquiring data received by the standard hydrophone as sample data sig (t);
performing time-frequency transformation on sample data sig (t) to obtain frequency domain data fsig (f) of the sample data, where f is 1, …, and f iss：
fsig(f)＝FFT[sig(t)]；
Using sample dataF 1, …, fsAnd target radiation noise simulation data frequency domain data mfconsig (f), calculating a compensation filter H0(f),f＝1,…,fs，H0(f) Has a unit of dB, and satisfies the following formula:
wherein lg (·) represents a base logarithm of 10;
using compensating filters H0(f),f＝1,…,fsPerforming spectrum structure compensation on the target radiation noise simulation data frequency domain data mfconsig (f), and obtaining target radiation noise simulation data frequency domain data mfconsig1(f) after the spectrum compensation at this time:
and (3) performing time-frequency inverse transformation on the mfconsig1(f), and taking a real part to obtain time domain data mfconsig1(t) of the target radiation noise simulation data, wherein the time domain data mfconsig1(t) is as follows:
mfconsig1(t)＝IFFT(mfconsig1(f))；
the time domain data mfconsig1(t) of the target radiation noise simulation data is subjected to standard deviation normalization processing to obtain equal target radiation noise simulation data targetsig1 with the power of 00(t) is:
simulating data total sound source SL according to target radiation noiseoutputFor targetsig10(t) carrying out amplitude compensation processing to obtain target radiation noise simulation data targetsig1(t) output by the signal source as follows:
wherein SLnormThe sound source level of the output signal of the system corresponding to the unimodal 1V signal;
the specific processing procedure of the iteration output module comprises the following steps:
comparing whether the spectrum structure of the target radiation noise simulation data targetsig1(t) output by the current compensation module is consistent with the structure of a target radiation noise spectrum SL (f) to be simulated or not; if not, the call compensation module continues the iterative compensation process until targetsig1(t) is consistent with the SL (f) structure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110583185.3A CN113466838B (en) | 2021-05-27 | 2021-05-27 | Iterative compensation target radiation noise data simulation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110583185.3A CN113466838B (en) | 2021-05-27 | 2021-05-27 | Iterative compensation target radiation noise data simulation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113466838A CN113466838A (en) | 2021-10-01 |
| CN113466838B true CN113466838B (en) | 2022-03-08 |
Family
ID=77871577
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110583185.3A Active CN113466838B (en) | 2021-05-27 | 2021-05-27 | Iterative compensation target radiation noise data simulation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113466838B (en) |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2651082B1 (en) * | 1989-08-16 | 1991-12-06 | Safare Crouzet | SUBMARINE BROADBAND TRANSMITTER |
| CN103900682B (en) * | 2012-12-26 | 2016-02-03 | 中国科学院声学研究所 | A kind of time-frequency combination formula signal fluctuation alignment detection method and system |
| CN104101871A (en) * | 2013-04-15 | 2014-10-15 | 中国科学院声学研究所 | Narrowband interference suppression method and narrowband interference suppression system used for passive synthetic aperture |
| CN104931955B (en) * | 2015-06-08 | 2017-06-20 | 哈尔滨工程大学 | A kind of method of compensating wide band Sonar system transmission signal amplitude |
| CN106358108B (en) * | 2016-08-31 | 2019-11-12 | 菁音电子科技（上海）有限公司 | Compensating filter is fitted system, sound equipment compensation system and method |
| CN107179535A (en) * | 2017-06-01 | 2017-09-19 | 东南大学 | A kind of fidelity based on distortion towed array strengthens the method for Wave beam forming |
| CN108917907A (en) * | 2018-05-17 | 2018-11-30 | 中国船舶科学研究中心（中国船舶重工集团公司第七0二研究所） | The conversion method of characteristics of ship underwater radiation noise test result in different water environments |
| CN110673100A (en) * | 2019-09-02 | 2020-01-10 | 江泽林 | Pulse compression method based on real-time spectrum estimation |
| CN111024209B (en) * | 2019-12-25 | 2022-03-29 | 深圳金信诺高新技术股份有限公司 | Line spectrum detection method suitable for vector hydrophone |
-
2021
- 2021-05-27 CN CN202110583185.3A patent/CN113466838B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113466838A (en) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ramírez et al. | Modeling nonlinear audio effects with end-to-end deep neural networks | |
| Feder et al. | Maximum likelihood noise cancellation using the EM algorithm | |
| Boya et al. | Identification of multiple partial discharge sources using acoustic emission technique and blind source separation | |
| CN109270525B (en) | Through-wall radar imaging method and system based on deep learning | |
| CN110109058A (en) | A kind of planar array deconvolution identification of sound source method | |
| CN112146751B (en) | Real-time sound field separation method based on time domain equivalent source method | |
| CN106954136A (en) | A kind of ultrasonic directional transmissions parametric array of integrated microphone receiving array | |
| CN107462317B (en) | A kind of water surface target radiation noise simulation method under complex scene | |
| CN109764956B (en) | Near-field acoustic holography method based on combinatorial optimization regularization method | |
| JP4293986B2 (en) | Method and system for representing a sound field | |
| CN113109794B (en) | Deep sea sound source depth setting method based on deep neural network in strong noise environment | |
| Sarabia et al. | Spatial LibriSpeech: An Augmented Dataset for Spatial Audio Learning | |
| CN113466838B (en) | Iterative compensation target radiation noise data simulation method and system | |
| CN106017837B (en) | A kind of analogy method of equivalent sound simulation source | |
| Sheng et al. | Intelligent adaptive filtering algorithm for electromagnetic-radiation field testing | |
| CN111965599B (en) | Sound source identification method for two-dimensional dynamic grid compression beam forming | |
| CN114358046B (en) | Multi-complexity-level complex electromagnetic interference environment simulation generation method and system | |
| US20200166608A1 (en) | System and method for modeling environmental data | |
| Delabie et al. | An acoustic simulation framework to support indoor positioning and data driven signal processing assessments | |
| CN113109795B (en) | Deep sea direct sound zone target depth estimation method based on deep neural network | |
| CN114252148A (en) | Sound field reconstruction method based on prolate ellipsoid wave superposition | |
| Okubo et al. | Three-Dimensional Sound Wave Propagation Reproduction by CE-FDTD Simulation Applying Actual Radiation Characteristics | |
| CN118295032B (en) | High-precision joint estimation method and system for instantaneous frequency and azimuth angle of underwater target | |
| Zhu et al. | A validation study for a bat-inspired sonar sensing simulator | |
| CN114216464B (en) | Intelligent underwater target positioning method |
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 |