673 results about "Wavenumber" patented technology

The invention relates to the field of ocean conductivity data processing, and discloses an Argo conductivity data correction method based on the influence of thermohaline internal waves, which comprises the following steps: installing a thermohaline depth sensor and an acoustic Doppler flow velocity profiler on an Argo buoy, sampling seawater in the vertical direction, and calculating the depth of the sampled seawater; acquiring temperature, salinity and flow velocity profile data of seawater; determining the upper and lower boundary depths of the thermohalocline and the thickness of the thermohalocline according to the acquired data, calculating the internal wave frequency and the internal wave number, and calculating the temperature internal wave amplitude, the salinity internal wave amplitude, the temperature internal wave phase and the salinity internal wave phase; establishing a conductivity data correction model according to the calculated parameters; and performing dynamic correction on the conductivity data by using the model. According to the method disclosed by the invention, the accuracy and the reliability of conductivity data can be remarkably improved, and support can be provided for research in many fields such as research on physical characteristics of the ocean thermocline and halocline, ocean ecological environment and ocean material circulation.

The invention relates to a two-parameter calibration-free pulverized coal temperature measurement method based on acoustic tomography, and the method comprises the following steps: 1, arranging a plurality of loudspeaker and microphone combination modules around a measurement field domain for transmitting and receiving acoustic signals, and forming an annular acoustic sensor array; step 2, multi-frequency sound transit time and amplitude attenuation measurement; 3, constructing a tomography inversion problem model of sound velocity and sound attenuation distribution; 4, respectively realizing tomography inversion of sound velocity and sound attenuation distribution through an improved joint iterative reconstruction method SIRT, and solving pulverized coal accumulation slowness and sound attenuation coefficient distribution under each frequency; 5, calculating a corresponding multi-frequency complex wave number through the multi-frequency slowness and the sound attenuation coefficient in each grid; step 6, constructing a temperature decoupling model of a complex wave number; and 7, regarding the temperature decoupling model of the complex wave number as a nonlinear objective function optimization problem, establishing a nonlinear objective function, and performing parameter decoupling and temperature calculation.

The invention discloses a road cavity volume calculation method and device based on a three-dimensional FK migration algorithm, equipment and a medium, and relates to the technical field of road cavity detection, and the method comprises the steps: obtaining three-dimensional ground penetrating radar data, and carrying out the preprocessing of the data, and obtaining target radar data; converting the target radar data to a frequency wavenumber domain to obtain frequency wavenumber domain radar data; establishing a road hierarchical velocity model, and performing hierarchical wave field extrapolation and double weight correction based on the frequency wave number domain radar data and the road hierarchical velocity model to obtain corrected wave field data; performing sparse constraint optimization imaging processing on the corrected wave field data to obtain a target imaging result; and carrying out boundary extraction processing on the target imaging result to obtain a three-dimensional cavity boundary point set, and carrying out voxelization processing based on the three-dimensional cavity boundary point set to obtain a road cavity volume. According to the invention, high-precision restoration of the shape of the road cavity and automatic extraction of the volume can be realized in three-dimensional GPR imaging.

The invention relates to a wavenumber domain compressed sensing damage detection method suitable for a plate-shaped structure. The method comprises the following steps: S1, constructing a compressed sensing equation; s2, solving the compressed sensing equation, and reconstructing to obtain a wave number domain signal of a full wave field; and S3, selecting a threshold value for the reconstructed wavenumber domain data for filtering, returning to a spatial domain through inverse Fourier transform, and carrying out merging processing on results by adopting a root mean square method to obtain a damage morphology graph. According to the method, the sparse basis matrix meeting the compressed sensing requirement is constructed by optimizing the sparse representation of the signal, so that the precision of full-wave field signal reconstruction is remarkably improved; the number of sampling points is greatly reduced, the full-wave-field data acquisition efficiency is improved, and the potential of breaking through the limitation of the Nyquist sampling law is achieved; the position and the shape of the damage are detected, the detection effect is good, and the detection efficiency is high.

The invention relates to the technical field of pipeline detection, and provides a non-contact pipeline fluid flow velocity measurement method based on a frequency wavenumber domain spatial spectrogram, and the method comprises the steps: collecting a turbulence signal in a pipeline through a piezoelectric film sensor; carrying out snapshot number segmentation on the turbulence signal to obtain a plurality of segmented turbulence time domain signals; fourier transform and narrowband signal processing are carried out on the segmented turbulence time domain signals, and a spatial spectrum corresponding to each narrowband signal component is obtained; constructing a three-dimensional frequency-wave number domain spectrogram according to the spatial spectrum function; and according to the three-dimensional frequency-wavenumber domain spectrogram, calculating the flow velocity of the pipeline fluid. According to the method, a three-dimensional spectrogram characteristic space is constructed through conjoint analysis of the frequency and the wave number, the angle, distance and frequency characteristics of the turbulence signals are decoupled through spatial spectrum estimation methods such as the MUSIC algorithm, multipath interference and noise are effectively restrained, the flow speed of fluid in a pipeline can be accurately estimated, and limitation in a traditional method is overcome.

The invention relates to the technical field of spectrum detection and optical imaging, and provides a time domain stimulated Raman scattering method and system.The method comprises the steps that pumping pulses or Stokes pulses are modulated, two stimulated Raman excitation pulse pairs with adjustable relative delay are introduced, a to-be-detected sample is excited to generate transient stimulated Raman scattering, and the time domain stimulated Raman scattering is obtained; transient stimulated Raman scattering signals can be obtained, acquisition of hyperspectral data covering hundreds of wave numbers can be achieved only through single excitation, and the method has the advantage of being high in time resolution. Meanwhile, the spectral resolution of the vibration level life broadening limit can be obtained through the method, and the main technical bottleneck faced by the traditional stimulated Raman scattering spectrum technology can be solved at a time. The method can be applied to the fields of Raman scattering imaging, spectrum detection, related Raman probes and the like.

The invention relates to a near-surface defect positioning method based on a laser ultrasonic surface wave, which comprises the following steps of: extracting a defect tiny inverse signal and calculating reflection time delay by utilizing a high-frequency surface wave signal excited by laser through a wavenumber domain method to determine a defect position, and determining the burial depth of the defect through a reflection wave amplitude value of the defect; comprising the steps of S1, laser ultrasonic finite element simulation; s2, acquiring high-frequency laser ultrasonic signals; s3, direct wave and reflected wave extraction and time delay calculation; s4, determining defect space coordinates; and S5, determining the burial depth of the defect. The method makes full use of the advantages of high frequency and long propagation distance of surface waves, can effectively solve the problem of defect detection and positioning of the tested piece, can quickly detect the submillimeter defects, is high in positioning precision, and is suitable for positioning and detecting the submillimeter defects on the near surface of a structure.

The invention discloses a near-fault broadband strong vibration simulation method, which comprises the following steps: S1, data collection and preprocessing: collecting a small vibration record in a research area as a source of an empirical Green function, and simultaneously obtaining information such as a geologic structure and stratum parameters of the area for model construction of a discrete wave number method; s2, correcting an empirical Green function, and selecting records similar to a target earthquake focus mechanism from the collected small earthquake records; and S3, calculating based on a geological structure and stratum parameters of a research area by using a discrete wavenumber method. According to the near-fault broadband strong vibration simulation method, the advantages of the two methods are fused, the description of the discrete wave number method on the medium structure at the low frequency band is accurate, the corrected empirical Green function can consider more actual factors at the high frequency band, accurate simulation of the broadband is achieved, and small vibration records are used as the empirical Green function, so that the simulation accuracy of the near-fault broadband strong vibration simulation method is improved. And information such as a real field effect is included, so that a simulation result is closer to an actual seismic oscillation condition.

The invention discloses a linear array data processing method for surface pulse pressure wave number-frequency spectrum measurement. A one-dimensional linear array is adopted to carry out a pulsating pressure measurement test, a model surface time domain signal is obtained, an array signal cross-spectrum matrix is calculated, and a flow direction wave number-frequency spectrum result is obtained through calculation by taking the cross-spectrum matrix as input. The method is directly obtained through a turbulent boundary layer pulsating pressure flow direction wave number-frequency spectrum physical definition formula and has direct relevance with the physical definition formula, the cross spectrum matrix method is applied, the calculation and analysis capacity is greatly improved, parallel calculation can be adopted in cross spectrum matrix calculation, a direct discrete algorithm and an average discrete algorithm, and the calculation speed is greatly improved.

The invention provides synthetic aperture sonar equipment towed by an unmanned aerial vehicle and a method, and belongs to the technical field of underwater sonar detection. The equipment comprises an unmanned aerial vehicle, a sonar towed body, a towing cable and a water surface control unit, the method comprises: establishing a measurement reference; based on the measurement basis, acquiring state data of the synthetic aperture sonar under the dragging of the unmanned aerial vehicle, and resolving in real time to obtain pose data; extracting a track direction velocity component from the pose data, performing displacement integration on the track direction velocity component, and generating a pulse trigger signal when the displacement integration quantity reaches a sound wave half-wavelength threshold value; synchronously latching current pose data and collecting acoustic echoes based on the pulse trigger signal, and generating a timestamp binding data packet; geometric phase compensation is applied to the timestamp binding data packet, an aliasing-free image is reconstructed through wave number domain resampling, and spatial spectrum aliasing caused by carrier motion disturbance is eliminated.

The invention relates to an unsteady-state sea wave direction spectrum estimation method. The method comprises the steps of S1, signal decomposition; s2, direction distribution calculation; s3, adaptively adjusting the time domain analysis step length; s4, performing nonparametric kernel estimation on the transient wavelet distribution probability; and S5, direction spectrum time domain adaptive reconstruction is carried out. In the time domain analysis process of wavelet analysis, a two-dimensional plug-in non-uniform bandwidth algorithm is adopted to determine a data-adaptive two-dimensional non-uniform and non-parametric kernel function, the probability of different wavelets in a wave number-direction space is analyzed, the lateral propagation sea wave energy in the non-main wave direction is estimated, and the wave energy in the non-main wave direction is obtained. Therefore, the precision and the anti-noise capability of the direction spectrum estimation method are improved. In the time domain analysis process of wavelet analysis, a sea wave signal fluctuation intensity detection algorithm based on the frequency-direction-wave number three-dimensional change rate of wave energy distribution is introduced, the time duration of sea wave signal sections is divided according to the fluctuation intensity, the time domain precision is adjusted, the calculation efficiency is improved, and meanwhile, the sea wave signals which rapidly change in the time domain are accurately captured. Measurement is avoided.

The invention provides an aero-engine blade damage detection method based on an improved ultrasonic sparse reconstruction algorithm, and relates to the technical field of aero-engine damage detection positioning, and the method comprises the steps: S1, setting a sensor array and a detection region of a to-be-detected blade, and obtaining an actual measurement damage reflection signal dictionary through a sensor; s2, wave number distribution of a detection area is determined according to a variable thickness structure Lamb wave propagation model, signal reconstruction is carried out, and a theoretical damage reflection signal dictionary is obtained; s3, an aero-engine blade damage detection model is constructed, an improved sparse reconstruction algorithm is adopted to solve and optimize, and damage judgment is carried out; and S4, according to the pixel value distribution condition of the aero-engine blade detection area, generating an aero-engine blade damage detection image, and positioning a damage position. According to the method, the reflected wave amplitude information can be obtained from the array sensor signals, and the reflected wave amplitude information is compared with the reconstructed theoretical reflected signals, so that the damage condition of the aero-engine blade is identified and positioned.

The invention provides a method for testing flow-induced noise of a pipeline. The method comprises the following steps: S1, determining measurement positions at an upstream flow channel and a downstream flow channel of a measured element, carrying out normalization processing on measured signals, and calculating the pressure power spectrum density of the inner wall surface of a pipe; s2, calculating a wall surface pressure fluctuation time domain cross correlation coefficient, and analyzing correlation characteristics between wall surface pressure pulsation; s3, spectrum intrinsic orthogonal decomposition is carried out on the pressure pulsation field, and space-time decoupling is achieved through time domain-frequency domain conversion; s4, extracting a dominant flow-induced noise mode, constructing an empirical frequency-wave number spectrum, and decomposing fluctuating energy distribution into contributions of different frequency and wave number combinations; s5, analyzing the decoupled flow-induced noise characteristic result, and identifying the amplitude-frequency characteristic of the flow-induced noise; according to the invention, real-time measurement of pressure pulsation in the pipeline can be realized, decoupling analysis of flow field pulsation and sound pressure pulsation can be carried out on the measured signal, and characteristics and a generation mechanism of pipeline flow-induced noise can be identified.

The invention discloses a sparse Bayesian learning shallow sea target depth and horizontal distance passive estimation method, and the method comprises the steps: calculating a modulus depth function hypothesis space through employing a finite difference method based on a prior sound velocity profile and a horizontal wavenumber hypothesis interval; preprocessing the vertical array data to a frequency domain, constructing a linear underdetermined equation set about the frequency domain vertical array data and a modulus depth function of a hypothesis space under sparse constraint, and solving by using a sparse Bayesian learning algorithm to obtain estimation of a horizontal wave number spectrum; performing spectrum peak search on the horizontal wave number spectrum, estimating the order and the horizontal wave number of an actual excitation mode, and estimating a mode depth function of each order of mode by using a finite difference method; vertical array copy sound fields under different depths and horizontal distances are calculated according to a normal wave sound field formula; and respectively performing conjugate transpose multiplication on the vertical array copy sound field and actual data to obtain a depth-horizontal distance ambiguity surface, and searching the position of the maximum surface value as the estimation of the depth and horizontal distance of the source.

The invention relates to the technical field of seismic exploration, in particular to a two-way wave depth migration method for a VTI medium. Respectively forming a detector wave field and a seismic source wave field; respectively reading seismic data and a seismic source wave field of a single shot, and obtaining anisotropic parameters; fourier transform is carried out on the seismic source wave field data and the detector wave field data to obtain frequency-space domain wave field data; calculating parameters required by offset; calculating an anisotropic vertical wave number; solving a wave field of a second depth grid; according to the double-boundary condition, wave field propagation operators are solved layer by layer; eliminating evanescent waves and pseudo S waves at each depth position; migration imaging is carried out; continuation calculation is repeated until the maximum depth is calculated, and an offset result is output; by popularizing the two-way wave migration from an isotropic medium to an anisotropic medium, the application range of the two-way wave depth migration is expanded, and an underground structure can be better imaged.

The invention discloses a multi-scale sea surface electromagnetic scattering calculation method and device based on marine environment elements and a storage medium, and belongs to the technical field of marine detection and physical remote sensing. The method comprises the following steps: calculating a tension wave sea spectrum according to an Elfouhail sea spectrum in combination with a Longuet-Higgins directional spectrum; a JONSWAP sea spectrum is combined with a Longuet-Higgins directional spectrum to calculate a gravity wave sea spectrum; superposing the tension wave on the gravity wave according to the slope by using a linear superposition method, and calculating a self-adaptive cut-off wave number; the contribution of small-scale waves to scattering is calculated by adopting an integral equation method, mirror scattering of large-scale waves is counted by adopting a Kirchhoff approximation method, and statistics and fusion are carried out on each scale wave according to the self-adaptive cut-off wave number, so that an electromagnetic scattering value is obtained. According to the method, a sea surface geometric calculation method based on real marine environment element driving is adopted, and accurate calculation is performed aiming at sea surface scattering characteristics under different sea conditions in combination with an integral equation method-Kirchhoff approximation electromagnetic scattering calculation model of an adaptive cut-off wave number, so that the calculation precision of scattering intensity is improved.

The invention discloses a noise reduction method of a bionic fish scale blade centrifugal pump, which comprises the following steps that a long and short blade centrifugal pump is selected for research, and the main design parameters of the long and short blade centrifugal pump are as follows: the specific speed ns is equal to 40.79, the rated flow Qd is equal to 10.6 m < 3 > / h, the lift H is equal to 34m, and the rotating speed n is equal to 2900r / min; a computational domain three-dimensional model of the long and short vane centrifugal pump is established through Creo8.0; in order to prevent backflow of inlet and outlet sections, ensure that fluid stably flows into an impeller and ensure the accuracy of a simulation result, extension pipelines are additionally arranged at an impeller inlet and a volute outlet of the centrifugal pump, the pipe diameter of the inlet extension section is 4 times that of the outlet extension section, and the pipe diameter of the outlet extension section is 8 times that of the inlet extension section; in order to verify the accuracy of a calculation result, a centrifugal pump experiment platform is built for external characteristic experiments; the bionic fish scale pit has the advantages that the dynamic and static interference between the impeller and the volute is reduced, the pressure fluctuation in the pump is stabilized, the amplitude fluctuation of a wave number-frequency spectrum and the energy region of a low wave number region are reduced, the migration of a wave number-frequency spectrum peak ridge is reduced, the curve change is gentle, and then the noise reduction effect is achieved.

The invention discloses a high-order random wave representation method for a non-Gaussian wind field, and the method comprises the following steps: S1, deducing a frequency-wave number spectrum of a random wave based on a target power spectrum density and a theoretical coherence function model gamma (xi, omega); s2, defining frequency-wave number three spectrums according to the frequency-wave number spectrums, and obtaining a biased three-coherence function of high-order random waves; s3, calculating a pure frequency-wave number spectrum through a partial three-coherence function; s4, performing Fourier-Stetius integral expansion on the random waves, decomposing an orthogonal increment process, and establishing a new orthogonal increment process according to a pure frequency-wave number spectrum and a pure frequency-wave number three-spectrum; and S5, substituting the established orthogonal increment into the Fourier-Stetius integral expansion of the random wave in the step S4 to obtain the high-order spectrum representation of the non-Gaussian random wave. The invention further discloses a non-Gaussian wind field simulation method and system and a storage medium.

The invention provides a method, a medium and a system for calculating ocean wave spectrum parameters of ocean digital twinning, and belongs to the technical field of electrical digital data processing. The method comprises the following steps: acquiring and preprocessing original sea wave spectrum data, selecting a Jonswap frequency sea wave spectrum as a theoretical model, establishing a multi-index fitting criterion, optimizing and solving a peak rise factor parameter, performing statistical analysis according to a significant wave high period level, and distinguishing sea wave spectrum characteristics of typhoon and non-typhoon working conditions. And calculating the wave amplitude wave number and the random phase to construct a wave surface height model, and finally forming a sea wave spectrum parameter database. The technical problem that sea wave spectrum parameters in the ocean digital twinning system are difficult to accurately describe and calculate is solved, the extreme sea condition simulation capability is particularly optimized, and high-precision sea wave field simulation technical support is provided for the ocean digital twinning system.

The invention provides a high-frequency electromagnetic excitation cable detection image reconstruction method and system, and belongs to the technical field of underground facility multi-field coupling detection and imaging, and the method comprises the steps: applying high-frequency electromagnetic excitation to an underground medium; receiving an electromagnetic signal transmitted through an underground medium, synchronously collecting a sound signal generated by electromagnetic excitation, and preprocessing the electromagnetic signal and the sound signal; performing structure imaging on the electromagnetic signal by adopting a frequency wave number migration imaging algorithm to obtain medium structure distribution; performing electrical parameter imaging on the electromagnetic signal by adopting a Born approximate inverse scattering imaging algorithm to obtain medium electromagnetic characteristic parameter distribution; sound source distribution reconstruction is carried out on the sound signals through a time reversal algorithm, and medium sound source response characteristic distribution is obtained; and inputting imaging results of the three algorithms into a pre-constructed deep reinforcement learning model for dynamic fusion to generate an underground medium comprehensive parameter reconstruction image.

The present invention relates to a magnetic susceptibility inversion method and system, the method comprising: S1, obtaining a spatial domain residual based on observed magnetic field intensity and theoretical magnetic field intensity; S2, expressing the spatial domain residual by two-dimensional discrete Fourier inverse transform based on wave number domain residual; S3, obtaining a wave number domain model perturbation based on the spatial domain residual, and further obtaining a spatial domain model perturbation based on Fourier inverse transform; S4, obtaining a magnetic susceptibility inversion model based on the spatial domain model perturbation. The present invention reduces the required memory space, improves the calculation efficiency, and further improves the inversion efficiency.

The invention belongs to the field of optical coherence layer imaging, and particularly discloses a spectrum calibration method for a linear field spectral domain OCT (optical coherence tomography) system, which comprises the following steps: acquiring interference signal instantaneous phases corresponding to two optical path differences by adjusting a lifting platform at a sample mirror in the linear field spectral domain OCT system, and subtracting the interference signal instantaneous phases to obtain an instantaneous phase difference; based on the calibration light source, discrete spectral lines are generated on a detector through a sample arm in the shielding line view field spectral domain OCT system, and spectral intensity distribution corresponding to the calibration light source is obtained; calculating the sub-pixel position of the peak value of each spectral line and the instantaneous phase difference corresponding to the sub-pixel position based on the instantaneous phase difference and the spectral intensity distribution corresponding to the calibration light source; and based on the wave number and the instantaneous phase difference corresponding to each spectral line, determining a linear relationship between the instantaneous phase difference and the wave number through linear fitting, and calibrating the corresponding wave number based on the linear relationship. According to the method and the device, the linear field OCT spectrum calibration result with higher precision can be realized.

The invention relates to the technical field of spectrum detection, in particular to a Raman spectrum-based rice oil detection method, which comprises the following steps: extracting main and auxiliary peaks to construct a peak group structure, performing sliding comparison to output a conformity region, identifying jump to remove an abnormal signal, calibrating a peak position to generate a correction structure, and performing homing numbering to complete a detection scheme. According to the method, the feature matching rule is constructed by combining the peak position spacing, the intensity ratio and the wavenumber index, the recognition precision of the multi-peak structure in the spectrum is enhanced, non-continuous response signals are recognized and eliminated through derivative change trend and proportion deviation constraint, and the accuracy of stable response extraction is improved. Dynamic calibration and position adjustment of peak positions are realized by utilizing an offset trend sequence, the consistency of a peak group structure in an overall direction and the continuity of a local structure are ensured, a numbering system and structure attribution division is completed after homing sorting, and the integrity of feature expression in a complex sample and the orderliness and reliability of spectral data processing are improved.

The invention discloses an ultra-deep reflection imaging method based on passive source seismic data and a related device, and relates to the technical field of earth depth structure imaging, and the method comprises the following steps: firstly, for original continuous seismic waveform data, intercepting seismic waveform data with a first arrival of a P wave and a non-arrival of an S wave, and carrying out the preprocessing of the data, so as to obtain a far earthquake event data segment; the method comprises the following steps: firstly, analyzing and identifying a far earthquake event with a high signal-to-noise ratio by using a time-frequency-wavenumber domain, removing data which are greatly influenced by noise, and then, reconstructing a reflection wave field of a high-quality far earthquake data fragment by using a multi-dimensional deconvolution technology, and carrying out convolution on a seismic source wavelet to obtain ultra-deep reflection wave field data. And finally, performing reflection imaging processing in combination with velocity analysis and superposition imaging technologies, optimizing the analysis capability of a deep reflection interface, and finally obtaining an earth deep reflection imaging profile, so that a continuous structure interface extending from the shallow earth crust to the deep mantle can be clearly revealed, and coupling characteristics between shallow and deep earth structures can be effectively captured.

The invention discloses a method for generating a reflection angle domain point spread function based on an LSTM network, and belongs to the technical field of geophysical exploration, and the method comprises the following steps: S1, selecting a reference point, and constructing a label through ray tracing and wave equation forward modeling; s2, training an LSTM (Long Short Term Memory) network as a ray tracing agent model of the whole space of an imaging domain; s3, predicting a travel time gradient and an amplitude factor through a ray tracing agent model; s4, combining the seismic source and the geophone in pairs, and calculating a reflection angle through the travel time gradient; s5, solving an illumination wave number vector, and solving a wave number domain point spread function in combination with the illumination intensity; s6, performing inverse Fourier transform on the wavenumber domain point spread function to obtain a spatial domain point spread function; and S7, solving a deconvolution problem to optimize the angle gather. According to the method, the calculation efficiency of solving the Green function through seismic wave ray tracing can be improved, the imaging response described by the point spread function is eliminated from the imaging result, and the balance of the resolution of the imaging result and the illumination intensity can be improved.

The present invention relates to the technical field of seismic exploration. Provided are a three-component seismic data vector anti-aliasing interpolation reconstruction method and an apparatus. The method comprises: representing three-component seismic data in a time domain as pure quaternions for vector joint; transforming the pure quaternion data in a time-space domain to a frequency-space domain to obtain frequency-space domain data; transforming the frequency-space domain data into a frequency-wave number domain; in the frequency-wave number domain, using an inclination angle scanning strategy to construct a mask matrix; performing vector anti-aliasing interpolation reconstruction on the frequency-space domain data and the mask matrix, so as to obtain a frequency-space domain reconstruction result; and inversely transforming the frequency-space domain reconstruction result to the time-space domain to obtain a time domain three-component three-dimensional data anti-aliasing reconstruction result. The present invention can effectively solve the problem of vector anti-aliasing reconstruction for three-component seismic data that contain both irregular and regular missing traces, thus improving the accuracy of effective wave event identification.

The invention discloses a real-time wave number space resampling method and system for SS-OCT, and relates to the technical field of signal processing, and the method comprises the steps: synchronously collecting a K clock signal and an OCT interference signal, and carrying out the filtering preprocessing of the signals; carrying out Hilbert transformation on the preprocessed K clock signal, constructing an analysis signal, extracting an instantaneous phase, and forming an instantaneous phase sequence; using a preset ideal equal wave number interval phase sequence as a target phase sequence, matching the target phase sequence with the instantaneous phase sequence, calculating a corresponding position index of each target phase in the instantaneous phase sequence, and obtaining a remapping sampling point index; and processing the preprocessed OCT interference signal by using an SINC interpolation method according to the remapping sampling point index to obtain an OCT signal which is subjected to resampling and has equal wave number intervals. According to the invention, the SINC interpolation algorithm is utilized to complete resampling in real time and with high precision in the data acquisition stage, and the requirement of high-speed real-time imaging is met.

The invention discloses a mobile underwater acoustic communication multipath channel modeling method based on sound field calculation, and relates to underwater sound waves. Channel modeling is carried out by inputting environmental parameters including an ocean sound velocity profile, submarine topography and the like, and modeling parameters including initial coordinates, motion trail equations and the like of a signal transmitting end and a signal receiving end. Performing ocean sound field calculation by adopting a ray tracing model, a normal wave model or a wave number integral model and the like, and calculating a sound field at each track sampling moment according to modeling parameters. And converting a frequency domain sound field result into a channel impulse response through inverse Fourier transform, and recovering the time-varying impulse response of the whole multipath channel by using local spline interpolation. And carrying out convolution on the transmitted signal and the time-varying impulse response of the multipath channel, and adding common propagation time delay to obtain a received signal simulation result. Through accurate sound field calculation and multipath channel modeling, the signal transmission condition in a specific marine environment can be simulated without actual marine environment measurement, the marine actual measurement cost is reduced, and the research and development efficiency is improved.

The invention provides a water surface / underwater target distinguishing method and system based on normal wave limited energy under an ice region half-sound channel, and the method comprises the steps: converting a time domain signal received by a vertical array in an ice region sea area into a frequency domain, and obtaining the frequency domain data of an array receiving signal; simulating and forecasting a normal wave modal depth function, a horizontal wave number and a phase velocity of the test sea area by using the normal wave model; extracting a normal wave modal amplitude which can be distinguished by test measured data; comparing the sound velocity of the phase velocity with the sound velocity of the critical resolution depth to obtain an order index of the limited mode; calculating the energy of a data extraction result limiting mode as a binary detection statistic; and setting a detection threshold according to the ROC curve, and if the limited modal energy is smaller than the detection threshold, determining that the target is an underwater target. Otherwise, determining the target as a water surface target. Compared with a traditional matching field and matching mode method, the method is small in calculation amount, high in environment tolerance and easy to popularize in an actual sonar platform.

The invention provides a method, medium and system for estimating the orientation of a horizontal array in an uncertain shallow sea environment, and belongs to the technical field of geophysics, and the method comprises the steps: firstly, receiving a far-field narrow-band signal through a horizontal array, obtaining a time domain sampling snapshot, then calculating an array sampling covariance matrix, and constructing a signal subspace through the first J feature vectors; meanwhile, according to a standard shallow sea sound velocity profile and a parameter uncertainty range, a horizontal wave number range of a propagation mode is obtained, sampling is carried out in the range to construct a mode matrix, and feature vectors corresponding to the first K maximum feature values are extracted to construct a mode subspace. Next, constructing an expansion matrix by using base vectors in the signal subspace and the modal subspace, performing QR decomposition to obtain an upper triangular matrix, and extracting diagonal elements of the upper triangular matrix to construct a sound source orientation estimation function; and finally, the angle corresponding to the peak value of the sound source orientation estimation function is the sound source orientation estimation value.