35 results about "Bessel function" patented technology

The invention provides a phase generated carrier demodulation method based on two-parameter estimation. The method comprises the following steps of: firstly, realizing estimation of a C value by extracting an amplitude ratio of a carrier frequency multiplication component in an interference signal frequency spectrum and combining a Bessel function characteristic; then, carrying out operation processing by utilizing a first frequency multiplication signal and a second frequency multiplication signal to extract carrier phase delay theta c; and finally, eliminating the influence through dynamic phase compensation. According to the demodulation method provided by the invention, the nonlinear drift of the modulation depth parameter C and the synchronous estimation and compensation of the carrier phase disturbance signal theta c can be realized, the cross interference error in the traditional demodulation system is effectively inhibited, the phase demodulation robustness is remarkably improved, and the harmonic nonlinear distortion is effectively inhibited.

The invention discloses a sine frequency modulation harmonic error compensation method based on Bessel function modulation characteristics, and relates to the technical field of sine frequency modulation nonlinear parameter compensation and high-precision measurement of target absolute distance based on multi-wavelength interference. The problem that the constant phase demodulation precision is insufficient due to the fact that carrier harmonic waves, accompanying amplitude modulation and modulation depth errors cannot be effectively suppressed due to laser modulation nonlinearity in existing sine frequency modulation distance measurement is solved. According to the method, after sine frequency modulation interference signals are collected, high-pass filtering is carried out to remove low-frequency components, Hilbert transform is carried out to remove envelope suppression associated amplitude modulation, window function frequency mixing is carried out to extract modulation depth, a matched carrier containing second harmonics is constructed, and frequency modulation harmonic parameter optimization is achieved through baseband amplitude maximization. And finally, carrying out orthogonal demodulation on the compensated interference signal to obtain a high-precision constant phase. The method is suitable for high-precision multi-wavelength interference ranging, precision displacement measurement and optical metering systems.

This invention discloses a method and system for measuring the frequency response of an electro-optic modulator based on a spectral analyzer. The system includes a laser, a 95:5 single-mode fiber splitter, an optical power meter, an electro-optic modulator, a control power supply, an RF signal generator, and a spectral analyzer. The measurement method includes the following steps: calculating a calibration factor to correct optical power drift through optical path splitting, driving the modulator by applying several discrete single-frequency RF points, and measuring the carrier power P0(f) using the spectral analyzer. i ) and the first-order sideband power P1(f i ), and determined by the power ratio R(f) i )=10 (P1‑P0) / 10 Combining the Jacobi–Anger expansion / Bessel function relation, the phase modulation index β(f) is solved by inverse solution. i Then, the phase transfer function H(f) at each frequency point is calculated and fitted. i The continuous frequency response S of the modulator is obtained. 21 This invention simplifies the testing process, reduces reliance on high-speed photodetectors and vector network analyzers, and enables high-precision electro-optic response testing of electro-optic modulators at high frequencies.

According to the method, a surface gravity wave and inner interface wave velocity potential field is analyzed into a product form of a vertical characteristic item and a horizontal diffraction field, and dimensionality reduction solution of a three-dimensional wave problem is realized; a vertical characteristic function orthogonal system meeting free liquid level, seabed boundary and interface coupling conditions is constructed, horizontal diffraction potential is decomposed into an incident field, a local column and an adjacent column scattering effect, and a Bessel function addition theorem is applied to construct a closed linear equation set of a multi-column interference field; establishing an extended dispersion equation containing density spring layer parameters, and revealing a modulation rule of layered strength on complex wave propagation characteristics; and finally, through tensor product coupling of a vertical mode and a horizontal diffraction field, constructing a full-three-dimensional diffraction velocity potential analytical solution, synthesizing a cross-layer pressure field integral by adopting a domain decomposition technology, and establishing a multi-column array generalized wave excitation force forecasting model. According to the method, efficient analytical calculation of the wave load of the multi-column system in the layered fluid is innovatively achieved.

The application discloses a non-circular sparse array DOA estimation performance evaluation method based on information theory, which comprises the following steps: firstly, a multi-dimensional probability density function (PDF) of a received signal is constructed, a joint PDF of the received signal and DOA and a DOA posteriori PDF are derived through information theory, and a system DOA information quantity is obtained by simplifying the DOA posteriori PDF with a Bessel function; then, a DOA posteriori PDF of given noise is derived, and a DOA information approximate upper limit is obtained by simplifying the DOA posteriori PDF with a Taylor expansion; finally, a DOA estimation performance index entropy error is obtained by using the posteriori differential entropy; the application builds a non-circular sparse array system DOA information theory framework based on information theory, in actual signal processing, the entropy error of a parameter can be calculated only by estimating the posteriori PDF of the parameter, and a performance limit independent of an algorithm is provided; in addition, it is found through simulation that the DOA information quantity approaches the DOA information upper limit under a high signal-to-noise ratio, and the entropy error approaches the Cramer-Rao limit, thereby verifying the rationality of the index.

The invention discloses a single-stage gear parametrically-excited vibration response frequency spectrum prediction method based on modulation feedback, and the method comprises the steps: building a single-degree-of-freedom kinetic equation containing time-varying meshing stiffness based on the Newton second law, and solving a main oscillation frequency through employing a modulation feedback theory and a trigonometric series approximation method; constructing a quantitative prediction model of excitation parameters and dominant side frequency orders on the basis of a first-class Bessel function, and solving system response by applying a four-order Runge-Kutta method in combination with fast Fourier transformation to verify the prediction model; a parameter sensitivity analysis framework is established through a harmonic coefficient method, and the influence of excitation parameters on the main oscillation frequency and the side frequency amplitude ratio is quantified.

The invention belongs to the technical field of ink-jet printing, provides a Bezier curve-based pulse generation method, device and equipment and a medium, and solves the technical problem that the acceleration of an acceleration section and a deceleration section of a printing trolley in the prior art is a fixed value, so that the printing trolley cannot run stably. The method comprises the following steps: acquiring an actual printing trolley speed curve according to a preset printing trolley speed curve, and acquiring acceleration and deceleration coefficients at each moment according to a third-order Bezier function corresponding to the actual printing trolley speed curve; according to the acceleration and deceleration coefficient at each moment, acquiring the frequency of a driving pulse at each moment; and generating the driving pulse according to the frequency of the driving pulse at each moment. The printing trolley can run stably, and the problems of liquid level warning misinformation, overall printer displacement and the like caused by machine shaking are solved.

The invention discloses a pile body torsion dynamic response analysis method and system based on a virtual soil pile in the technical field of pile foundation dynamic response analysis. The method comprises the following steps: constructing a total energy function of the pile-soil coupling system; variation is conducted on the total energy function of the pile-soil coupling system on the basis of the variation principle, and a pile periphery soil body radial control function in a general solution form and a Bessel function form of a pile body control differential equation is obtained; according to the general solution form of the pile body control differential equation and the built-in boundary conditions of the bottom of the virtual soil pile section, an equivalent complex impedance function of the top of the virtual soil pile section is constructed; according to the continuity condition and the boundary condition of the pile body and the equivalent complex impedance function of the top of the virtual soil pile section, a system response linear algebraic equation set is constructed; and iteratively updating the control variable of each layer of soil body to obtain the converged control variable of each layer of soil body, further calculating a torsion power response value at the pile top, and obtaining a torsion impedance function of the pile body based on the torsion power response value at the pile top.

The invention discloses a novel modal shape function construction method suitable for a variable cross-section flexible structure, and the method comprises the following steps: firstly, building a linear kinetic equation of a flexible wing based on an Euler-Bernoulli beam theory and a Lagrange principle; secondly, geometric and physical parameters in the equation are normalized through dimensionless variable substitution, and the coupling relation between different variables is simplified; secondly, formalizing the solution into a linear combination of a first type of Bessel function and a second type of Bessel function by adopting a special function expansion method so as to ensure that boundary conditions are met and effectively reflect the influence of section change on modal characteristics; and finally, solving the inherent frequency and a modal shape function by constructing a novel modal shape function. According to the method, iteration and approximate truncation are not needed, the inherent frequency and the modal shape can be rapidly and accurately calculated, the calculation efficiency is remarkably improved, high calculation precision is guaranteed, and the method is suitable for efficient dynamic analysis of a complex structure.

The invention discloses a micro-motion exploration frequency dispersion curve extraction method, device and medium, and relates to the technical field of micro-motion exploration, and the method comprises the following steps: S1, calculating the average cross-correlation spectrum of each station pair according to the recorded data of a plurality of stations in an observation array; s2, aiming at the average cross-correlation spectrum of each station pair, identifying a first intersection point frequency of the average cross-correlation spectrum and a frequency axis; s3, based on the first intersection point frequency, performing piecewise fitting on the average cross-correlation spectrum to extract a phase velocity: for a frequency band not greater than the first intersection point frequency, performing fitting by adopting a zero-order Bessel function J0; and for a frequency band of which the frequency is greater than the frequency of the first intersection point, performing fitting by adopting a first-class zero-order Hankel function H0 (1). In a single method, the contradiction of artifact interference in a high frequency band and insufficient detection depth in a low frequency band in the prior art is overcome at the same time, the accuracy and reliability of frequency dispersion curve extraction under the sparse array condition are improved, and the detection depth of micro-motion exploration is effectively expanded.

The invention relates to the technical field of geodetic parameter inversion, and discloses a horizontal hierarchical geodetic parameter inversion method, system, equipment and medium, and the method comprises the steps: constructing an inversion objective function according to an apparent resistivity prediction function and an apparent resistivity measurement value which are constructed through employing a kernel function and a first-class zero-order Bessel function; and representing the kernel function by adopting an odd polynomial and an even polynomial, representing the odd polynomial and the even polynomial on the basis of an adjacent layered soil resistivity change rate set and a layered thickness sum set obtained by adopting a permutation and combination algorithm, and then performing partial derivative calculation, according to the obtained first partial derivative relational expression and the second partial derivative relational expression, a complex mirror image method is adopted to obtain a third partial derivative relational expression of the inversion target function to the soil resistivity and a fourth partial derivative relational expression of the inversion target function to the layering thickness; and carrying out horizontal layered geodetic parameter inversion by adopting a trust region algorithm to obtain an inversion result. According to the method provided by the invention, the accuracy and reliability of a horizontal layered geodetic parameter inversion result are improved.

The application discloses a nonlinear phase error compensation method of an optical fiber demodulator, and relates to the technical field of phase error compensation.The method comprises the following steps: synchronously collecting interference signals output by a photoelectric detector by an analog-to-digital converter, performing low-pass filtering and amplitude extraction on the digitized interference sequence by using a carrier wave and a twice-frequency orthogonal reference signal, calculating a current amplitude ratio according to a fundamental wave amplitude and a second harmonic amplitude, estimating a current actual phase modulation depth by using a Bessel function relationship, and generating an error signal by a digital PID.The application adjusts a driving voltage amplitude in real time by using an analog voltage control programmable high-voltage driving source, thereby guaranteeing the safety monitoring of a high-pressure gas pipeline and solving the problem of system failure caused by phase modulation depth drift in long-distance oil and gas pipeline network distributed optical fiber sensing in the prior art.

The invention discloses a carbonate rock gas reservoir productivity prediction method and device. The method comprises the following steps: establishing a first seepage equation of a quasi-steady-state seepage mode and a second seepage equation of an unsteady-state seepage mode in the carbonate rock gas reservoir in combination with a non-darcy effect and a stress sensitive effect; performing dimensionless processing and Laplace transformation to obtain a first comprehensive control differential equation and a second comprehensive control differential equation; performing Fourier cosine transform, combining a Bessel function general solution, and obtaining a pressure distribution solution in the pull-type space through point source and line source functions; a real space quasi-pressure solution is obtained through inversion according to the pressure distribution solution in the pull-type space, and then a yield prediction model under the fixed flowing bottomhole pressure condition is obtained; and based on the established productivity prediction model, predicting the yield of the target well according to the fixed bottomhole flowing pressure of the target well in the carbonate rock gas reservoir. According to the method, complex relevance among triple media of the carbonate gas reservoir is fully considered, and reasonable yield prediction is achieved.

The invention provides a double-interferometer system demodulation method and device, and the method comprises the steps: generating laser through a modulation laser, and inputting the laser into a signal interferometer and a reference interferometer, determining a first orthogonal signal and a second orthogonal signal corresponding to the signal interferometer, and a third orthogonal signal and a fourth orthogonal signal corresponding to the reference interferometer; the arm length differences of the signal interferometer and the reference interferometer are equal; solving a first modulation depth of the signal interferometer and a second modulation depth of the reference interferometer through an ellipse fitting algorithm, calculating Bessel function values of the first modulation depth and the second modulation depth, and obtaining an orthogonal signal without the modulation depth according to the orthogonal signal; and constructing a product sum-difference formula according to the orthogonal signal of which the modulation depth is eliminated, and obtaining an optical phase signal of which the noise is eliminated through arc tangent operation.

The invention relates to the technical field of satellite system engineering and micro-vibration suppression, in particular to a method for quickly constructing and optimizing an optical remote sensing satellite micro-vibration suppression index system, which comprises the following steps of: acquiring an imaging quality index threshold, a system static modulation transfer function value, a structure transfer function data set and a flywheel disturbance source characteristic data set; constructing a system-level frequency domain response model based on a frequency domain transfer function chain, calculating a modulation transfer function attenuation factor caused by high-frequency micro-vibration by utilizing a zero-order Bessel function and a time delay integral smoothing function, and establishing a reverse design inequality to solve an allowable maximum disturbance force upper limit curve of the flywheel under different rotating speed frequencies; and vibration isolation device design parameters or flywheel model selection strategies are output through comparison with a flywheel actual measurement disturbance spectrum. According to the method, the problems that existing time domain simulation calculation is low in efficiency and lacks parameter inverse solution capacity are solved, millisecond-level rapid evaluation and direct design guidance are achieved, and the overall satellite scheme argumentation period is remarkably shortened.

The invention discloses a demodulation stability improving method and system based on primary ellipse fitting, and belongs to the technical field of phase demodulation. According to the invention, based on a one-time ellipse fitting technology, an additional low-frequency sinusoidal signal is superposed in laser modulation, so that an interference signal generates a complete elliptical orbit, and an ellipse fitting algorithm is utilized to extract a direct current amount, an alternating current amplitude and a phase difference of fitting parameters in advance. A modulation depth value C is inquired and matched through fitting parameters in combination with a Bessel function table, accurate calibration of key parameters is achieved, then the phase difference of non-strict orthogonal signals is corrected, direct current is removed, and demodulation distortion caused by phase deviation is eliminated. Meanwhile, an intensity disturbance signal B containing light intensity fluctuation and circuit noise is calculated in real time based on the modulation depth value C and the fitting parameter, and is dynamically removed in the demodulation process, so that the intensity disturbance signal B is prevented from being superposed to a sensing signal, and the stability of a demodulation amplitude value is remarkably improved.

The invention relates to the technical field of transformer station grounding grid evaluation, and provides a dynamic evaluation method and system for power station grounding short circuit transfer potential, and the method comprises the steps: employing a Wenner quadrupole method to measure the apparent resistivity of soil, and constructing an apparent resistivity-depth change curve; reversely deducing the resistivity and the thickness of each layer of soil by adopting an inversion algorithm and combining an analysis kernel function and Bessel function expansion, so as to obtain soil layering parameters of a multi-layer and non-uniform structure; and obtaining a short-circuit current, a grounding grid structure and the obtained soil layering parameters, and obtaining a grounding short-circuit transfer potential of the grounding grid of the to-be-evaluated power station by adopting coupling simulation analysis mixed with a finite element method and a moment method. According to the method, the real soil model is constructed by combining the Wenner quadrupole method measurement and the inversion algorithm, and the finite element and moment method coupling simulation is introduced, so that the grounding short circuit transfer potential is accurately evaluated, and the prediction precision and the system safety under the complex geological condition are remarkably improved.

The present invention discloses a data processing method and system for the layout of a microseismic detection array, which relates to the technical field of microseismic detection. The method includes: setting the measuring point distance, the radius and quantity of the circumscribed circles of nested triangles according to the measured starting point coordinates and the azimuth angle of the survey line, and batch generating the coordinates of field stations; collecting the waveform data of the stations and preprocessing it; calculating the spatial autocorrelation coefficient between any central station and non-central stations; performing DBSCAN clustering based on the station distances; performing spatial averaging on the spatial autocorrelation coefficient according to the clustering results; fitting the averaged spatial autocorrelation coefficient with the zero-order Bessel function to obtain the dispersion spectrum, picking up the surface wave phase velocity in the dispersion spectrum, and calculating the apparent shear wave velocity; adopting an adaptive gradient inversion strategy to invert the shear wave velocity of the waveform data; generating a velocity comprehensive profile in the depth domain and frequency domain according to the surface wave phase velocity, the apparent shear wave velocity, and the inverted shear wave velocity. The present invention provides a new method for microseismic detection from station layout to data processing.

This invention discloses a method and related equipment for detecting geological structures during ocean drilling. The method includes: acquiring preliminary data of a candidate area and collecting seabed background noise in the candidate area using a concentric ring array to obtain raw noise data; preprocessing and segmenting the raw noise data to obtain segmented noise data; performing azimuth-weighted averaging and Bessel function fitting on the segmented noise data based on the principal noise direction of the concentric ring array to obtain a target dispersion curve; using the target dispersion curve as a constraint, performing nonlinear inversion using a genetic algorithm to obtain a shear wave velocity structure model; jointly verifying the shear wave velocity structure model with the preliminary data, and selecting drilling sites based on the joint verification results. This invention provides independent, quantitative, and traceable geophysical evidence for pre-drilling site selection, effectively reducing drilling risks, and can be widely applied in the field of seismic structure detection technology.

The invention discloses a rapid calculation method and system for a grid-connected current THD of a high-voltage direct-hanging energy storage system, and belongs to the technical field of high-voltage direct-hanging energy storage systems. The method aims at solving the problems that existing THD calculation depends on multi-power-frequency period sampling, is slow in updating and cannot adapt to a scene with inconsistent modulation signals of in-phase sub-modules. The method comprises the following steps: sampling to obtain system parameters and modulation signals; calculating modulation signal adjustment quantity amplitude and additional fundamental frequency voltage; a final modulation signal is obtained through alpha-beta axis vector superposition; calculating a harmonic voltage amplitude based on a carrier phase-shifting modulation model and a Bessel function; deriving a harmonic current amplitude and a fundamental current amplitude, and solving THD; and the harmonic calculation frequency is dynamically adjusted through error verification. Only control variables are used for calculation, the calculation speed is high, the precision is high, the relative error is smaller than or equal to 5%, the DSP burden is reduced, control strategy adjustment can be rapidly supported, and stable operation of a power grid is guaranteed.

The invention discloses an electromagnetic sensor signal analysis method and system based on real abrasive particle dynamic trajectory and time-harmonic field intensity coupling, and belongs to the field of mechanical equipment state monitoring and sensor simulation. The method comprises the following steps: firstly, based on a CFD-DEM coupling theory, combining abrasive particle multi-physical field stress, and accurately reconstructing a real three-dimensional dynamic track in a complex flow field of a lubricating oil pipeline; and establishing a three-dimensional time-harmonic magnetic field analysis model based on magnetic vector potential and a modified Bessel function, constructing a dynamic mapping relation between the transient position of the abrasive particles and a non-uniform magnetic field, quantifying the disturbance characteristic of the abrasive particle eddy current to the magnetic field, and solving the transient response of the sensor induced electromotive force. Abrasive particle size, material conductivity and magnetic conductivity parameter correction are fused, simulation distortion caused by a simplified motion model and neglect of radial magnetic field difference in a traditional theory is overcome, and high-precision prediction of dynamic signals of the whole process that the abrasive particles pass through the sensor is achieved. The method is suitable for health monitoring of aero-engines and industrial gearboxes.

This invention discloses a method for spatial domain filtering and spectrum correction of broadband signals in the mode domain of a circular array. The method involves performing a Discrete Fourier Transform (DFT) on the received time-domain data to obtain the array element frequency domain data. The broadband signal frequency band data is then extracted and subjected to a Circular Fourier Transform (CFT) to obtain the sub-band data in the mode domain. The distance between frequency points and the zeros of the Bessel function within the broadband signal frequency band is calculated, and a zero-point correction method is used to replace the frequency points. Sub-band frequency points in the mode domain are selected according to the correction results. The sub-band weighting vectors for mode domain beamforming are obtained according to the minimum variance and distortion-free criterion for data independence, and the beamformer is calculated. The sub-band data in the mode domain are weighted and summed based on the weighting vectors of each sub-band to obtain the beam output frequency domain data. Zeros are padded at all frequency points except for the broadband signal frequency domain data. An Inverse Fourier Transform is then performed on the beam output frequency domain data to obtain the beamforming output data sequence. This invention effectively preserves the original spectral components of the broadband signal while improving the signal-to-noise ratio.

The application provides a behavior recognition method and device based on Doppler frequency shift, an electronic device and a storage medium, comprising: for the frequency domain channel response of a Wi-Fi signal, cleaning the frequency domain channel response to obtain a cleaned power delay spectrum and a cleaned channel impulse response corresponding to the cleaned frequency domain channel response; determining a power sum corresponding to an NLOS channel based on the cleaned power delay spectrum; determining a correlation power value corresponding to the NLOS channel based on the cleaned channel impulse response; determining a multipath channel fading factor according to the power sum and the correlation power value; determining a corresponding Doppler frequency shift in a zero-order Bessel function based on the multipath channel fading factor; and determining a behavior recognition result of a passive target according to Doppler frequency shifts at multiple observation times. The application scheme realizes accurate behavior recognition of a passive target.

The application provides a kind of GIL loss fast calculation method considering electromagnetic space distribution, relate to high voltage gas insulated transmission technology field, solve the problem that existing GIL loss calculation method is difficult to consider physical rigor and calculation efficiency.The method, the coaxial hollow cylindrical structure electromagnetic field model of GIL is established, and then the diffusion equation of time-harmonic electromagnetic field is constructed and the boundary condition is set;Solving time adopts the linear combination of first kind Bessel function and second kind Bessel function as the general solution form, determines the undetermined coefficient, obtains the magnetic field intensity and electric field intensity distribution expression in the conductor region and shell region;Finally, based on the Poynting theorem, the analytical expression of unit length conductor ohmic loss and shell eddy current loss is calculated and derived, and the conductivity parameter is embedded in the form of explicit function of temperature, that is, the GIL loss calculation is completed.The calculation realized by the application can provide a clear, differentiable and traceable heat source coupling interface for the heat network model.

The application provides a quick calculation compensation method for the influence of rainfall on visible light-near infrared band optical communication, and is used for solving the technical problems that the traditional Mie scattering calculation method is difficult to meet the fast and efficient requirements of related research because the programming implementation technology is complex, the calculation speed is slow, numerical instability is prone to occur and the like due to relying on the Bessel function, continued fraction recursion and the like. The quick calculation compensation method for the influence of rainfall on visible light-near infrared band optical communication can realize the quick and stable calculation of Mie scattering optical parameters of three rainfall types and any rainfall intensity in the visible light-near infrared band on the basis of guaranteeing the calculation accuracy, realizes the quick and real-time output of the calculation results, and can be widely applied to the optical communication compensation technology in the fields of laser communication, automatic driving and the like.

The present application relates to the technical field of geodetic parameter inversion, and discloses a horizontal layered geodetic parameter inversion method, system, device and medium, an apparent resistivity prediction function constructed by using a kernel function and a first-order zero Bessel function is used to construct an inversion objective function; the kernel function is represented by using an odd polynomial and an even polynomial, and after the adjacent layered soil resistivity rate of change set and the layered thickness and set obtained by using a permutation and combination algorithm are represented, partial derivative calculation is performed, according to the obtained first partial derivative relationship and second partial derivative relationship, the third partial derivative relationship of the inversion objective function to the soil resistivity and the fourth partial derivative relationship of the inversion objective function to the layered thickness are obtained by using a complex mirror method; the horizontal layered geodetic parameter inversion is performed by using a trust region algorithm, and an inversion result is obtained. The method improves the accuracy and reliability of the horizontal layered geodetic parameter inversion result.

The application provides a directional diagram synthesis method of a multi-mode orbital angular momentum array antenna, relates to the technical field of multi-mode orbital angular momentum array antennas, and comprises the following steps: constructing a concentric circular ring array, determining a first array factor of the concentric circular ring array based on a radiation directional diagram function of a multi-layer nested concentric circular ring; applying a Bessel function theory to derive a second array factor in an integral form; obtaining a third array factor by utilizing the central axis rotational symmetry of the concentric circular ring array and the angular symmetry of the array factor; equivalently converting the third array factor into a straight line array model with a Bessel function factor, and deriving a corresponding equivalent straight line array factor expression; and applying the optimized excitation weight to the first array factor to obtain a final required power directional diagram by setting a constraint condition, using Matlab, and solving by using a convex optimization solving program cvx. The application has good adaptability to large-scale arrays, simplifies two-dimensional array synthesis into a one-dimensional problem, and has the advantages of low calculation complexity and high speed.

The present invention relates to the field of power system grounding measurement and evaluation, and provides a method, system, and storage medium for calculating apparent resistivity of layered earth. The method comprises the following steps: initializing a pole distance a according to the Wenner quadrupole method, inputting a test current I, obtaining a test voltage ΔV corresponding to the pole distance a, and calculating an apparent resistivity measurement value ρ. am ; The Euler transformation method is introduced to obtain the surface potential difference ΔV′(a) between the two voltage poles and the calculated value of apparent resistivity ρ as ; According to whether the relative error ε exceeds the preset threshold, confirm the update of the pole distance a or output the calculated value of the apparent resistivity ρ at this time as The present invention converts the infinite Bessel function integral of the spatial potential into the weighted sum of a finite number of surface characteristic impedances by introducing the Euler transformation method. The processing time is short, it is not limited by the number of soil layers, is applicable to complex earth stratification, and can be quickly applied on engineering sites.

The present application belongs to the technical field of underwater acoustic channel and physical field, and relates to a numerical hybrid calculation method of forward scattering acoustic field of a regular shape object underwater. The method comprises: obtaining an initial projection profile of a three-dimensional geometric model of the regular shape object under a specified projection plane and an incident direction; dividing discrete points of the initial projection profile into upper boundary points and lower boundary points; interpolating and completing the upper boundary points and the lower boundary points, so that the upper boundary points and the lower boundary points have the same number and share the same horizontal coordinates; discretizing the completed projection profile by using trapezoidal plane elements; calculating a forward scattering far field pattern function based on the discretization result; and calculating a forward scattering target strength of the object. The method provided by the present application can cover the shape range of the scattering object applicable to the analytical method, and does not need to calculate complex series expressions composed of special functions such as Bessel functions and spherical harmonic waves, thereby obviously improving the calculation efficiency.

The present application relates to a kind of full-band ESPAC microseismic dispersion imaging method and system based on zero-order stroud function.The method comprises: reading the original microseismic signal recorded by microseismic station array, calculating the average cross-correlation spectrum between station pairs;The Hilbert transform is carried out on the average cross-correlation spectrum, and the complex analytical cross-correlation spectrum is constructed;Unified composite kernel function based on zero-order stroud function is constructed, and the unified composite kernel function is combined by zero-order first kind Bessel function and zero-order stroud function;The complex analytical cross-correlation spectrum is fitted in full band using unified composite kernel function, and the fitting energy residual is calculated;Wherein, the fitting process does not need to be segmented in frequency domain, and does not need to identify the zero-crossing frequency of cross-correlation spectrum;Based on the fitting energy residual of all station pairs, generate dispersion energy spectrum, and extract dispersion curve.Can consider low-frequency band large depth detection capability and high-frequency band artifact suppression capability in full band, improve fitting stability.