88 results about "Time frequency resolution" patented technology

A system and method for monitoring a patient includes a sensor configured to acquire physiological data from a patient and a processor configured to receive the physiological data from the at least one sensor. The processor is also configured to apply a spectral estimation framework that utilizes structured time-frequency representations defined by imposing, to the physiological data, a prior distributions on a time-frequency plane that enforces spectral estimates that are smooth in time and sparse in a frequency domain. The processor is further configured to perform an iteratively re-weighted least squares algorithm to perform yield a denoised time-varying spectral decomposition of the physiological data and generate a report indicating a physiological state of the patient.

An audio encoder for encoding an audio signal has a first coding branch, the first coding branch comprising a first converter for converting a signal from a time domain into a frequency domain. Furthermore, the audio encoder has a second coding branch comprising a second time / frequency converter. Additionally, a signal analyzer for analyzing the audio signal is provided. The signal analyzer, on the hand, determines whether an audio portion is effective in the encoder output signal as a first encoded signal from the first encoding branch or as a second encoded signal from a second encoding branch. On the other hand, the signal analyzer determines a time / frequency resolution to be applied by the converters when generating the encoded signals. An output interface includes, in addition to the first encoded signal and the second encoded signal, a resolution information identifying the resolution used by the first time / frequency converter and used by the second time / frequency converter.

The invention provides a high-precision time-frequency analysis method, and belongs to the field of digital signal processing. The method comprises the following steps of: (1)inputting original digital signals x(t), then carrying out empirical mode decomposition (EMD)to the original digital signals x(t), decomposing the original digital signals x(t) into sets of intrinsic mode function (IMF)components, namely ci(t) is the IMF component; (2) carrying out matching pursuit (MP)decomposition to each IMF component ci(t) respectively to obtain a time-frequency spectrum of each IMF component; and (3) carrying out mutual superposition to the time-frequency spectrum of each IMF component obtained in the step (2) to obtain the time-frequency spectrum of a high time-frequency resolution of the original digital signals x(t). The high-precision time-frequency analysis method disclosed by the invention avoids inaccuracy of time-frequency atoms when MP decomposition matching is implemented for composite signals directly, and MP decomposition is implemented for each IMF component, thereby being capable of obtaining the time-frequency spectrum with higher accuracy and time-frequency resolution and richer information.

The invention relates to a power transmission line short circuit fault diagnosis method based on empirical wavelet transformation and local energy. The method is characterized by comprising the following steps of: collecting short-circuit fault voltage signals; using an empirical wavelet transformation method to process the short-circuit fault voltage signals to obtain an empirical mode; performing short-circuit fault detection on the basis of empirical wavelet transformation; and performing short-circuit fault classification feature extraction on the basis of block local energy. The method has the advantages of scientificity, rationality, high time-frequency resolution, high noise resistance, high adaptability, good application effect and the like.

In an automotive system (300), a knock detection scheme is provided for detecting knock events in an internal combustion engine that are sensed by a transducer structure, such as a non-intrusive acoustic accelerometer sensor (310) to generate sensor signal information which is processed by a signal processing structure (312, 316, 318, 326) which extracts digital signal parameters from the sensor signal information to identify a predetermined pitch frequency (330) and any short-term energy increase (328) in the digital signal information which in combination are used to provide a positive indication of engine knock behavior. When a short term Fourier transform (324) is used to extract the digital signal parameters, time frequency resolution may be improved by appropriately windowing the digital signal being transformed.

The invention discloses a rolling bearing fault intelligent diagnosis method based on improved S-transformation and deep learning. According to the method, a window width adjustment factor is added toS-transformation, so that the window width of a Gaussian window function can be changed, and therefore, the time-frequency resolution of the S-transformation can be improved, and the S-transformationcan accurately detect an impact component in vibration signals, and thus, the fault characteristics of the vibration signals of a rolling bearing can be better extracted. Improved S-transformation isperformed on the vibration signals of the fault of the bearing, so that the feature matrix of the signals can be obtained; the feature matrix is subjected to column-based expansion so as be convertedinto a feature vector, and the feature vector is inputted into a sparse autoencoder model; on the basis of the characteristics of the encoder, the deep features of the data are further extracted, sothat some important hidden information that is not recognized by humans can be mined; the extracted features are accurately classified. With the method of the invention adopted, the accuracy of the fault diagnosis of the rolling bearing can be effectively improved.

The acoustic impedance of the respiratory system can be inferred from oscillations that are generated in an airway of a subject. The impedance describes the frequency-dependent relation between the resulting oscillations in flow and pressure. When the impedance varies from inspiration to expiration, it has to be estimated with a high time resolution. A method is provided that reliably estimates the impedance in time intervals that are short enough for physiological purposes. A simple version of the uncertainty principle has been derived for discrete time and frequency. A discrete time-frequency transform has been developed that gives an optimal time-frequency resolution according to this principle. The transform is orthonormal, which permits an analysis of variance in the discrete time-frequency domain. The impedance follows from bivariate least-squares analysis in the time-frequency domain, under the assumption that noise is present in both flow and pressure.

The invention discloses a microseismic signal noise reduction filtering method based on a VMD and a wavelet packet and belongs to the signal processing technology field. In the invention, a mode of combining the VMD and the wavelet packet is adopted; adaptivity of a VMD decomposition method and characteristics of a high mathematical theory basis, high frequency noise suppression and the like in the algorithm are used; and the wavelet packet has the characteristics of carrying out multilevel division on a frequency band, further decomposing a high frequency portion which is not subdivided in multiresolution analysis, and adaptively selecting a frequency band according to a characteristic of an analyzed signal so that a time frequency resolution is increased. In the invention, based on a condition of maintaining randomness, non-stationary performance and a burst transient characteristic of a microseismic signal, the microseismic signal is filtered. The algorithm is simple and easy and an effect is ideal. Effective noise reduction filtering can be performed on a mine noisy microseismic signal. And a good technical value and an application prospect are possessed.

The invention provides a spectrum analysis method of a seismic signal and an apparatus thereof. The method comprises the following steps: acquiring the seismic signal; carrying out Hilbert transformation operation on the seismic signal to obtain an analytic signal of the seismic signal; carrying out WVD processing on the analytic signal; according to a result of the WVD processing, determining a Gauss window function and an attenuation coefficient; according to the determined Gauss window function and the attenuation coefficient, carrying out smooth pseudo-WVD distribution time frequency analysis processing on the result of the WVD processing. Through the method and the apparatus in the invention, a high time frequency resolution of seismic signal spectrum analysis can be obtained.

The invention discloses an estimation method of a multi-component non-stationary signal instantaneous frequency based on focusing S-transform, and mainly solves the problems that the multi-component non-stationary signal instantaneous frequency with high precision cannot be realized because the energy aggregation is not high and the time frequency resolution ratio cannot be flexible to adjust and cannot realize high precision. The estimation method comprises the realizing steps: 1, inputting a to-be-analyzed non-stationary signal; 2, optimizing a control parameter of a window function in the focusing S-transform according to an aggregation criterion; 3, calculating the optimized focusing S-transform and acquiring a signal time frequency distribution energy diagram; 4, carrying out binarization processing on the time frequency distribution energy diagram of the non-stationary signal; and 5, obtaining the instantaneous frequency estimation on a signal according to a binaryzation time frequency diagram after the binarization processing. According to the method provided by the invention, the self-adaption adjustment of the control parameter of the window function can be realized by improving the S-transform window function, and the time frequency distribution energy aggregation and the instantaneous frequency estimation accuracy can be improved, therefore the method can be used for radar interception, communication countermeasure, voice identification, and medical brain wave signal analysis.

The invention discloses a nonlinear signal time-frequency analysis method based on short-time fractional order Fourier transform, and the nonlinear signal time-frequency analysis method belongs to the technical field of target signal analysis. The nonlinear signal time-frequency analysis method utilizes short-time fractional order Fourier transform, subjects a signal to segmentation processing and combined processing of different orders of fractional order Fourier transform through short-time fractional order Fourier transform, adopts an angle constraint method for fractional-order order searching, improves the time-frequency resolution of the target signal and the precision of time-frequency analysis, completes the time-frequency analysis of single-component and multi-component nonlinear signals, solves the problems of low time-frequency resolution and difficult separation of the target nonlinear echo signal, is conductive to parameter estimation of the weak target signal, improve the precision of signal parameter estimation, has the advantages of high time-frequency resolution and precise separation of different signals in the time-frequency domain, significantly reduces the amount of computing, provides reliable data support for a follow-up system to complete tasks, and is worthy of adoption and popularization.

The invention discloses a seismic spectrum imaging method based on deconvolution generalized S transform, which comprises the steps that a generalized S transform spectrum is acquired by performing two-dimensional convolution on Wigner distribution of original signals and a Gaussian window, a transform spectrum is acquired by performing generalized S transform on seismic data, and time-frequency distribution of the original signal can be acquired through deconvolution when the generalized S transform spectrum and the Wigner distribution of a window function are known. The seismic spectrum imaging method combines advantages of generalized S transform and Wigner-Ville distribution, generation of a cross term of the Wigner-Ville distribution is effectively suppressed through a generalized S transform window, and the generalized S transform spectrum is enabled to acquire high time-frequency aggregation at the same time; and deconvolution generalized S transform can adaptively adjust an analysis window along with variations of a frequency component, is applicable to time-frequency analysis for unstable seismic data, and can acquire high time-frequency resolution; and the seismic spectrum imaging method is applied to detecting the oil-gas possibility of a reservoir, thereby being conducive to improving the reservoir prediction accuracy.

A multirate spectral analyzer receives a data sequence through a tapped delay line (510) onto a plurality of filter channels, each of which includes a first decimator (520), a polyphase filter (530), and a second decimator (540). Each filter channel is coupled to a transform processor (550) which modulates the frequency response of the polyphase filters to a plurality of identical bandpass filters spaced evenly at 2πk / M intervals, where k=1, 2, . . . , M−1, and M is the number of filter channels. The transform processor (550) outputs a plurality of time series representative of the frequency component for each of the corresponding filter channels. By choosing a prototype filter order to be much larger than the number of filter channels, the effects of aliasing in the output are greatly reduced. The reduction in aliasing further allows, by prudent selection of the two decimators, for fractional temporal overlap to be implemented in the spectral estimates.

The invention discloses a method for carrying out separation and time-frequency analysis on respiration and heartbeat signals in non-contact life detection on the basis of HHT (Hilbert Huang Transform), which comprises the following steps of: by adopting a band-pass filter, separating the respiration and heartbeat signals and filtering out noise and harmonic waves to the greatest extent; by adopting EMD (Empirical Mode Decomposition), decomposing the filtered respiration and heartbeat signals into IMF (intrinsic mode functions) of single-frequency signals; filtering out the IMF comprising the respiration harmonic waves by adopting a maximum correlation criterion, i.e. extracting the IMF with the maximum correlation with an original signal and using the IMF as the respiration or heartbeat signal; and directly carrying out Hilbert transformation on the filtered output signal. Compared with a conventional time-frequency analysis technology, the method not only can effectively remove noise and harmonic wave interference of the respiration and heartbeat signals, but also has higher time-frequency resolution so as to effectively describe variation of life characteristics of a human body.

In an automotive system (300), a knock detection scheme is provided for detecting knock events in an internal combustion engine that are sensed by a transducer structure, such as a non-intrusive acoustic accelerometer sensor (310) to generate sensor signal information which is processed by a signal processing structure (312, 316, 318, 326) which extracts digital signal parameters from the sensor signal information to identify a predetermined pitch frequency (330) and any short-term energy increase (328) in the digital signal information which in combination are used to provide a positive indication of engine knock behavior. When a short term Fourier transform (324) is used to extract the digital signal parameters, time frequency resolution may be improved by appropriately windowing the digital signal being transformed.

The invention provides an ECG signal processing method. ECG signals can be analyzed from a plurality of time-frequency resolutions, and SB-LSTM and TD-CNN are separately utilized to mines the overallfluctuation pattern and local fluctuation pattern of signals from ECG sub-layer signals, and a Sum Rule is finally used for fusing the intermediate classification results of each ECG sub-layer to obtain more accurate and robust classification results. The ECG signal processing method avoids individual differences and intra-individual differences caused by manual classification, and can obtain thestable classification results; and in addition, the ECG signal processing method does not need any expert knowledge, does not need to manually design various features, does not need to perform featureselection processing, does not need to separately construct classifiers, is a typical end-to-end method, and has the characteristics of convenient use, high classification accuracy and the like.

The invention aims to provide a method for oil and gas detection by utilizing seismic data instantaneous frequency attribute, which comprises the following steps of: acquiring and recording seismic initial data; carrying out generalized S transformation on the data in a frequency domain to obtain the time-frequency spectrum of each sampling point thereof, and determining the frequencies respectively corresponding to 15 percent and 35 percent of the total energy in the time-frequency spectrum as the initial frequency and the final frequency for computing low-frequency-band absorption attenuation gradient; after taking a natural logarithm from the time-frequency spectrum p (Tk, f) between the initial frequency and the final frequency and fitting by utilizing a least square method to obtain the low-frequency-band absorption attenuation gradient of the sampling point; and corresponding a position with strong low-frequency-band absorption attenuation gradient to a favorable oil and gas contained area, and determining an oil and gas detection result. The invention has lower requirements for signal-to-noise ratio, strong time-frequency resolution capability, simple algorithm and high computational efficiency.

The invention provides an instantaneous rotation speed estimation method for rotation equipment. According to the method, the instantaneous frequency of rotor vibration signals is precisely estimated,and thus rotation speed estimation of the equipment is achieved. The method comprises the steps that firstly, a parameterization time frequency analysis method is used for roughly estimating the instantaneous frequency of the vibration signals, then, a novel envelope tracking filter is used for precisely correcting an initial estimation result, and the finally obtained signal instantaneous frequency directly reflects the instantaneous rotation speed of the equipment; the envelope tracking filter uses a Fourier series model to depict the signal envelope, and a least square method is adopted for optimizing and solving an envelope parameter. According to the method, the instantaneous rotation speed estimation method which combines parameterization time frequency conversion with the envelopetracking filter is adopted, the defect can be overcome that the estimation precision is limited by the time frequency resolution ratio in the prior art, and the method has the advantages of high estimation precision and easy operation.

The invention discloses a method for inhibiting cross terms in the time-frequency distribution of multi-component linear frequency modulation (LFM) signals based on subspace decomposition, which belongs to the technical field of signal processing. In the method, a time-frequency distribution matrix comprising noises and the cross terms is decomposed into signal subspaces and noise subspaces by utilizing singular value decomposition (SVD)-based subspace decomposition. For the problem that the linear modulation signals occupy relatively more bandwidths to make singular values are reduced at relatively lower rates and cannot be separated from the noises effectively, an angle of inclination of the time-frequency distribution of the signals can be obtained by utilizing Wigner-Hough transform. The time-frequency distribution of the signals is rotated according to the obtained angle so as to be parallel to a time base. The method is characterized in that: the singular values are rapidly decreased to zero, and then the signal sub-spaces can be separated effectively. By the method, the cross terms and the noises in the time-frequency distribution of the multi-component LFM signals are inhibited without reducing time-frequency resolution; therefore, the method is vast in application prospect.

The invention discloses a seismic attenuation qualitative estimation method based on self-adaptive optimal kernel time frequency distribution. The method comprises the following steps of 1) collecting original seismic data; 2) picking a position of a seismic data target layer; 3) determining a high frequency fH and a low frequency fL; 4) calculating self-adaptive optimal kernel time frequency distribution for each signal in the superimposed seismic data body to obtain a high-frequency component and a low-frequency component; 5) determining a correction factor; and 6) estimating the seismic attenuation. According to the method, the self-adaptive optimal kernel time frequency distribution is firstly used for carrying out the time frequency decomposition of the seismic data to obtain the high-frequency component and the low-frequency component of the seismic data, and then the seismic attenuation is qualitatively estimated by utilizing the difference of the low-frequency component and the high-frequency component. Compared with the traditional seismic attenuation qualitative estimation method based on the time frequency distribution, the method has high time-frequency resolution ratio, and the seismic attenuation can be more precisely reflected; and in addition, the technical scheme is easy to realize, and the operability is strong.

An embodiment of the invention discloses an anomaly detection method of network flow, comprising: performing line frequency modulation wavelet transformation on detected flow signals, obtaining a transforming spectrum; dividing the transforming spectrum into at least two frequency spectrum coverages; performing signal reconstruction on each frequency spectrum coverage, operating the reconstructed signal through deviation fraction algorithm, obtaining deviation fraction; comparing the deviation fraction with presupposed threshold, determining whether the network flow is abnormal according to the comparing result. An embodiment of the invention also discloses a device realizing the method. The embodiment of the invention decomposes according to different time frequency resolution, thereby effectively detecting continuous anomaly flow for a long time and breaking anomaly flow for a short time, and detecting high frequency and middle frequency attack flow which are not detected by the network flow anomaly detection method based on wavelet transformation, through a second time detection after signal self-adapting reconstitution, anomaly character is further determined, and detecting accuracy is advanced.

Embodiments are described for a system and method for implementing an adaptive time-frequency resolution in audio and video coding systems. A method of adaptively transforming the time-frequency resolution for a defined spectrum comprises dividing the spectrum of the input signal into a into plurality of bands; determining, for each band of the plurality of bands, a characteristic of the content (e.g., tonal or transient content); modifying the time-frequency resolution value to one or more bands of the plurality of bands to increase either a time resolution of the band or a frequency resolution of the band depending on the characteristic of the content; determining a cost associated with modifying the time-frequency resolution value of the one more bands based on an entropy measure of the bands, and altering the modified time-frequency resolution values in a manner that accounts for the coding cost.

The invention discloses a traveling wave time-frequency analysis method based on variational mode decomposition and Wigner-ville distribution. The method includes the following steps that: fault traveling wave signals are detected, and Karebauer phase-mode transformation is performed on three-phase voltage traveling wave signals, so that a traveling wave aerial mode component is obtained; variational mode decomposition is performed on the traveling wave aerial mode component, so that K intrinsic mode components are generated; Wigner-Ville analysis is performed on each intrinsic mode component;and the Wigner-Ville distribution of each intrinsic mode component is linearly superimposed, so that the time-frequency domain distribution of original traveling wave aerial mode signals is obtained.With the method of the invention adopted, the interference of cross terms in Wigner-Ville distribution can be effectively suppressed; a good noise suppression effect enabling good VMD (Variational Mode Decomposition) is preserved; time-frequency resolution enabling high Wigner-Ville distribution and good time-frequency aggregation are preserved; traveling wave time-frequency domain information characteristics are truly and accurately represented; and fault traveling waves can be completely observed. The method is of important theoretical and practical significance for the practical application of fault traveling wave protection and positioning.

The invention discloses a multichannel HHT (Hilbert-Huang Transform) based UWB (Ultra Wideband) radar human body motion micro-Doppler feature extraction method, which comprises the steps of 1, acquiring UWB radar signals; 2, preprocessing the acquired UWB radar signals; 3, performing effective channel selection; and 4, carrying out multichannel time-frequency spectrum fusion. The micro-Doppler time-frequency analysis method provided by the invention has higher time-frequency resolution, micro-Doppler features corresponding to different components are more obvious in component differentiation, and transient features and detail changes of human body motions can be better reflected. The time-frequency analysis method provided by the invention is high in anti-interference ability, and obvious micro-Doppler features can still be captured under the conditions of long-distance through-wall detection and sharp signal attenuation. The time-frequency analysis method provided by the invention has better denoising capability, the acquired time-frequency spectrum is higher in signal-to-noise ratio, and micro-Doppler feature components are more obvious.

The invention discloses an NMF-based micro-seismic weak signal recognition method. The method comprises the following steps of: firstly carrying out time frequency analysis on micro-seismic signals byadoption of S transformation, and rearranging time frequency spectrums in a frequency direction; factorizing a rearranged time frequency matrix by adoption of non-negative matrix factorization (NMF)so as to obtain a frequency domain base vector and a time domain position vector; extracting characteristic parameters such as sharpness, derivative quadratic sum, information entropy and sparseness to construct characteristic spaces of the micro-seismic signals; and finally classifying the micro-seismic signals by adoption of a least squares support vector machine (LSSVM). The method is capable of strengthening low-frequency weak signals and improving the time frequency resolution ratios, and has good time domain and frequency domain localization abilities.

The invention discloses a description method of reservoir sand bodies on an ultra-deep coal layer. The method comprises the following steps of: firstly, analyzing waveform characteristics of erosional truncation type pinchout points and coal layer characteristics by means of forward modeling of different models, and providing bases for boundary determination of practical data sand bodies and coal layer stripping; utilizing overlying strata low trough attribute extreme points to identify pinchout lines of the sand bodies; utilizing a multichannel dynamic matching tracking method under the control of layer positions to carry out coal layer high-shielding stripping, and highlighting reservoir weak signals; on the basis of coal layer high-shielding stripping, utilizing a relatively-high time frequency resolution in GST to carrying out frequency division and reconstruction on frequency spectrum components of main information of the reservoir in a time frequency domain, wherein the processed data can predict the sand bodies more accurately; and finally, according to GR pseudo-acoustic impedance inversion, obtaining three-dimensional wave impedance data bodies for comprehensive analysis, and obtaining a reservoir prediction result. By adopting the method provided by the invention, the pinchout lines of the sand bodies are effectively described, the range of the sand bodies of the ultra-deep reservoir is precisely predicted, and the precision of seismic reservoir prediction is greatly improved.

The invention discloses a fault diagnosis method for an AC / DC charging equipment power device based on wavelet packet analysis, and the method comprises the steps: analyzing a charging equipment powermodule, and determining a fault class; obtaining output signals of the charging equipment power module under the normal condition and various fault conditions; carrying out the high-low frequency decomposition of the output signals, carrying out the multilayer wavelet packet decomposition of the high-low frequency signals obtained through decomposition, and carrying out the reconstruction of a wavelet packet decomposition coefficient; calculating the extracted signal energy at all frequency bands, carrying out the normalized calculation of the signal energy, and finally determining a centralized frequency band range of the signal energy; carrying out the power spectrum analysis of the reconstructed wavelet packet decomposition coefficient, determining the feature frequency and power spectrum value of a frequency band signal, comparing the comparison with the fault analysis, and determining the type of a fault; and carrying out the fault type coding of the type of the fault. The methodimproves the time-frequency resolution of a signal, achieves the precise fault positioning through the spectrum features, and improves the precision of the fault diagnosis of the AC / DC charging equipment power device.

The invention relates to a method for detecting voltage sag based on generalized hyperbolic S-transformation, which comprises the following steps of (1) sampling detected voltage signals; (2) carrying out generalized hyperbolic S-transformation on the sampled signals; (3) calculating a fundamental frequency voltage amplitude value according to an S-model matrix; (4) comparing the fundamental frequency voltage amplitude value with a set voltage sag threshold value, and judging whether voltage sag occurs or not; and (5) extracting a characteristic value if voltage sag occurs, and carrying out detection and analysis on start and stop moments and phase jumps of the voltage sag. Compared with an existing detection method, the method provided by the invention adopts the generalized hyperbolic S-transformation, improvement is carried out on an original hyperbolic window function, and a window width factor and a window width change rate adjustment factor are introduced, thereby being capable of flexibly adjusting the time-frequency resolution, reducing the sensitivity of detection signals for noises, and being capable of accurately detecting the amplitude value, the start and stop moments and the phase jumps of the voltage sag. The method provided by the invention has wide application prospects in fault detection of an electric power system.

The invention discloses a seismic data time-frequency analysis method based on second-order extrusion wavelet transformation. Firstly, a second-order synchronous extrusion operator is introduced on the basis of continuous wavelet transformation, the transformation is rearranged in a time-scale domain according to wavelet transformation coefficients, accordingly energy diffusion is pressed, and analysis results higher in time-frequency resolution ratio can be obtained by adopting the method. Compared with existing synchronous extrusion transformation, the method is more suitable for analysis of seismic signals having stronger frequency modulation, thus the signals are applied to time-frequency analysis of actual seismic data, the boundary of an underground carbon and hydrogen reservoir can be more accurately defined, riverways, faults and other geological information can be more elaborately depicted, and accordingly further seismic data explanation and well position determination are promoted.

The invention discloses a method for estimating a stratum equivalent quality factor based on synchronous extrusion wavelet transform, and the method is a quality factor extraction method based on a Qcompensation principle. By using synchronous extrusion wavelet transform with extremely high time-frequency resolution and energy aggregation to transform seismic data of adjacent horizons (post-stackseismic data) or adjacent seismic records (vertical seismic sections) to a time-frequency domain first, then inverse Q filtering is carried out in a preset certain Q value range, and by calculating amean square error between a result after time-frequency domain Q compensation and a reference waveform frequency spectrum, a corresponding quality factor when the minimum mean square error is satisfied is searched as the best equivalent quality factor, thereby realizing optimal estimation of the quality factor. The method for estimating a stratum equivalent quality factor based on synchronous extrusion wavelet transform avoids the difficulty of frequency band selection for estimating a quality factor in a conventional spectral ratio method, the hypothesis that a seismic source frequency spectrum is a Gaussian spectrum in a centroid frequency method or the like does not exist, methods for quality factor estimation are expanded, and the estimation accuracy of the quality factor is improved.