31 results about "Chirp rate" patented technology

The invention relates to the technical field of power equipment state monitoring, and discloses a partial discharge detection method and system based on synchronous compression wavelet-synchronous extrusion. The method comprises the following steps: preprocessing a partial discharge signal; calculating a coefficient through continuous wavelet transform and estimating an instantaneous frequency and a chirp rate; synchronous compression is carried out to compensate the phase, and scale remarking is carried out; energy is reallocated to a frequency axis through synchronous extrusion, and a high-aggregation-degree time-frequency diagram is generated; extracting a peak domain based on the saliency spectrum, demodulating along a ridge line, and estimating event parameters; and finally, robust judgment and auditing output are completed by using the joint statistics. According to the method, the sensitivity and robustness of partial discharge signal detection under the conditions of low signal-to-noise ratio and strong interference are remarkably improved.

The application discloses a chirp rate measurement method of a chirped volume grating, and comprises the following steps: constructing a measurement light path; starting a laser and an energy detector; adjusting the included angle and relative position of a to-be-measured chirped volume grating sample and a laser output light beam by adjusting an angle adjusting table and a translation table, recording the angle and position when the chirped volume grating diffracts, and calculating an angle average value; and calculating the chirp rate of the to-be-measured chirped volume grating sample according to a formula. The measurement principle is simple, the calculation method is clear, and the chirp rate of a chirped volume grating of any specification can be measured.

An elastic wave device 1 according to the present invention comprises a support member 3, a piezoelectric film including a piezoelectric layer 6, and at least one IDT electrode 7 having a plurality of electrode fingers. In the at least one IDT electrode 7, the electrode finger pitches and / or the widths of the electrode fingers vary. When one electrode finger pitch is denoted as P1 and the other electrode finger pitch is denoted as P2 in three continuous electrode fingers, and Mmax represents the widest width and Mmin represents the narrowest width among the widths of the three continuous electrode fingers, an adjacent chirp rate represented by at least one of formulae (1) and (2) is larger than 0° in at least a part of the at least one IDT electrode 7. In the entirety of the IDT electrode 7, an adjacent chirp rate represented by each of formulae (1) and (2) is 0.6° or less.　(1): Arctan [|P1-P2| / (P1+P2)] (2): Arctan [|Mmax-Mmin| / (P1+P2)]

This signal parameter estimation device is configured to comprise: a candidate value setting unit (1) that sets a plurality of candidate values for a chirp rate representing the ratio of chirp frequency width to chirp time width of a reception signal, which is an observation signal of a chirp signal that is the target for parameter estimation; and a candidate value selection unit (2) that multiplies the reception signal by opposite-phase components of each of the candidate values set by the candidate value setting unit (1), converts the reception signals after the opposite-phase component multiplication to signals of a frequency region, and selects, from among the plurality of candidate values, the candidate value for which an electrical power peak value of the signal of the frequency region is relatively larger. Moreover, the signal parameter estimation device comprises: a frequency width estimation unit (3) that multiplies the reception signal, which is the observation signal, by the opposite-phase component of the candidate value selected by the candidate value selection unit (2), converts the reception signal after the opposite phase component multiplication to a signal of the frequency region, and, on the basis of frequency intervals included in the post-conversion signal of the frequency region which correspond to a plurality of electrical power peak values, estimates the chirp frequency width; and a time width estimation unit (4) that estimates the chirp time width from the chirp frequency width estimated by the frequency width estimation unit (3) and the candidate value selected by the candidate value selection unit (2).

The invention discloses a three-dimensional Chirplet time-frequency analysis method based on bidirectional sparse search, and belongs to the technical field of signal processing and acoustic target recognition. Comprising the following steps: acquiring a to-be-analyzed non-stationary multi-component signal, and performing three-dimensional Chirplet transformation on the to-be-analyzed non-stationary multi-component signal to obtain three-dimensional time-frequency-chirp rate combined feature mapping fusing time, frequency and chirp rate features; selecting a plurality of initialized candidate points from the three-dimensional time-frequency-chirp rate joint feature mapping, and performing forward and backward path search on the initialized candidate points to obtain a sparse ridge track; and sparse ridge line components are extracted based on the sparse ridge line trajectory, residual signal energy updating is performed on the three-dimensional time-frequency-chirp rate joint features based on the sparse ridge line components, and time domain representation fusing all the sparse ridge line components is output. The method is suitable for sparse feature extraction and component separation of a non-stationary multi-component signal, cross interference and background noise can be effectively suppressed, and the identification capability and reconstruction precision of a signal structure are improved.

The FMCW LiDAR system includes an optical drive electronic circuit to receive a reference frequency signal and a beat frequency signal to generate a drive signal. The optical drive electronic circuit includes a TDC to calculate a phase difference between the reference frequency signal and the beat frequency signal and a digital ramp control to: provided the phase difference is a positive value, produce a ramp down control signal to increase a current chirp rate to an increased chirp rate; provided the phase difference is a negative value, produce a ramp up control signal to decrease the current chirp rate to a decreased chirp rate. The optical drive electronic circuit includes a digital integrator to generate a digital output based on at least one of the ramp down control signal or the ramp up control signal and a DAC to convert the digital output to an analog output.

The application discloses a method for estimating the sweep period of linear frequency modulation interference, and is mainly used for estimating the sweep period of linear frequency modulation interference in the field of signal processing and the like. The method can preset a transform order and a calculation parameter, generates a coefficient required for fractional order Fourier transform calculation, and does not need to normalize the dimension of input data in advance. According to different transform orders, fractional order Fourier transform calculation is performed on input data, until the average peak value interval is calculated after multiple peak values appear in the transform domain. Finally, the sweep period is estimated according to the current searched transform order and the average peak value interval. The method greatly reduces the calculation complexity of the sweep period estimation, multiplexes the anti-interference hardware architecture, saves the hardware resources, does not need to normalize the dimension of input data in advance, increases the compatibility of the method, and does not need to search the period with a fixed length, thereby reducing the calculation amount and increasing the estimation accuracy.

The invention discloses a chirp modulation frequency offset and timing joint estimation signal synchronization method and a chirp modulation frequency offset and timing joint estimation signal synchronization system, which relate to wireless communication and signal processing technologies, and comprise the following steps: completing coarse synchronization by using a demodulation result of an upper chirp in a lead code; the method comprises the following steps: accurately estimating frequency offset and time offset of signals based on a forward chirp signal and a reverse chirp signal in a chirp signal pilot frequency sequence synchronization head, intercepting a section of forward sweep frequency signal at the synchronization head, and performing refined spectrum estimation on an upper chirp spectrum peak position of the forward sweep frequency signal; intercepting a section of reverse frequency sweep signal at the synchronization head, and performing refined frequency spectrum estimation on the lower chirp spectrum peak position of the reverse frequency sweep signal; and resolving a peak value generated by a signal frequency spectrum according to refined frequency spectrum estimation results of the forward frequency sweep signal and the reverse frequency sweep signal to obtain a frequency offset compensation result. According to the method, errors generated by the system can be estimated with fewer symbols, and the method is applicable to different chirp rate conditions.

The application discloses a kind of laser detection device and method, belong to laser detection technical field, the device includes beam output module, light splitting module, optical path multiplexing module, two coherent mixing modules and signal processing module;Beam output module outputs the positive and negative double-chirp laser of frequency-modulated time sequence synchronization of different wavelength, chirp rate absolute value equal, after light splitting, optical path multiplexing, get local light and echo light, two coherent mixing modules generate four-way quadrature beat frequency signal, signal processing module does product subtraction and addition operation to four-way beat frequency signal and constructs complex signal, offset Doppler broadening interference introduced by rotating polygon mirror, again Fourier transform, peak searching obtains intermediate frequency parameter and calculates target vector velocity, the application can completely eliminate the influence of Doppler broadening on speed measurement, significantly improve speed measurement precision.

The application provides a real-time correction method and device for detecting equipment operation parameters, equipment and a medium, and belongs to the technical field of ore body exploration. The method comprises the following steps: acquiring point cloud data of a target area based on a detecting equipment; calling a preset dynamic local quality model to perform spatial division and physical property identification processing on the point cloud data to obtain a quality voxel; performing multipole moment expansion and gravity projection processing on the quality voxel to obtain a terrain gravity correction number; performing frequency shift conversion and modulation processing on a frequency shift parameter of a radar scanner based on the terrain gravity correction number to obtain a compensation frequency shift parameter, which is used for adjusting a laser chirp rate to offset a phase shift of the radar scanner; and correcting an operation parameter of a gravimeter in the detecting equipment based on the compensation frequency shift parameter. The above scheme can effectively improve the stability of the detecting equipment in performing an exploration task.

The invention discloses a multi-focal-length laser focusing system and method, and belongs to the technical field of photoelectrons and lasers. The multi-focal-length laser focusing system comprises a laser beam splitting unit, a light path conversion unit and a laser multi-focal-length focusing unit. The laser beam splitting unit is used for splitting an incident laser beam; the light path conversion unit collimates the plurality of beam splitting lasers and stretches the collimated laser beams; the laser multi-focal-length focusing unit comprises a first acousto-optic driving source and a first acousto-optic device, the first acousto-optic driving source sequentially generates a plurality of chirp signals with adjustable energy, and the first acousto-optic device forms a plurality of acoustic sub-gratings with different grating intervals; and the acoustic sub-grating performs focusing with different focal lengths on the split light beams which are shaped by the light path conversion unit. By adjusting the angle and the intensity of the beam splitting laser and the number and the chirp rate of the chirp signals, the focusing number and the focusing focal length of the beam splitting laser are controlled. The adjustable multi-beam multi-focal-length focusing is realized, and the size of the device is reduced.

To provide a distance and velocity measurement apparatus that can be adopted preferably in a LiDAR or a sensor for a robot, wherein the apparatus can prevent deterioration of SN ratio even in a case where an object in an external environment vibrates. A LiDAR 20 according to the present embodiment includes a first laser apparatus 1a, a second laser apparatus 1b, a polarization-maintaining type optical fiber 2, a WDM filter 6, an optical fiber coupler 3a, an optical amplifier 11, an input / output unit 4, an optical scanner 5, a second optical fiber coupler 3b, a balanced photodetector 7, and a square-law detector 9. Further, a delay line 10 composed of a polarization-maintaining optical fiber is provided on the local port 2b. The first laser apparatus 1a includes a device for generating a first laser light having a first wavelength and a first chirp rate in an interior thereof, and the second laser apparatus 1b includes a device for generating a second laser light having a second wavelength that differs from the first wavelength and a second chirp rate that differs from the first chirp rate.

The invention relates to a system and a method for simultaneous range and velocity measurement in an FMCW LiDAR system. A first light source produces first light having a first frequency that varies according to a first chirp rate. A second light source produces second light having a second frequency that is constant or that varies according to a second chirp rate. A splitter separates the measuring light into reference light and output light, and a scanning unit directs the output light towards an object and receives input light that is obtained by reflection of the output light at the object. A detector detects a superposition of the reference light and the input light. A computing unit computes unambiguously the range and relative velocity by analyzing beat frequencies resulting from the superposition. Ambiguities due to Doppler frequency shifts are removed by performing a decision tree analysis.

The application discloses a time-frequency analysis method and system for bearing fault signal classification processing, electronic equipment and a computer readable storage medium, and belongs to the technical field of time-frequency analysis of non-stationary rotating machinery fault signals. A set of signal classification processing strategies are introduced, and the initial short-time Fourier transform time-frequency representation of the signal is divided into two categories. The time direction compression is adopted for one category with a chirp rate greater than a threshold value, and the frequency direction compression is adopted for the other category with a chirp rate less than the threshold value, so as to reduce the serious error of the frequency direction compression in the case of facing the signal frame with a large chirp rate and the time direction compression in the case of facing the signal frame with a small chirp rate. In the noisy working environment of the rotating machinery, the method can more accurately extract the fault characteristic frequency of the bearing, and effectively diagnose the health condition of the machinery. The problems that the fault characteristic frequency of the bearing cannot be accurately extracted and the health condition of the machinery cannot be accurately diagnosed in the prior art are solved.

The application discloses a three-dimensional extraction method of complex superimposed dense channel reservoir based on statistical linear frequency modulation, and comprises the following steps: S1, inputting a seismic signal x(t) to be analyzed, wherein t is time; S2, constructing a linear chirp transform in a time-frequency-chirp rate (TFC) three-dimensional space, wherein omega is frequency, and beta is a chirp rate; S3, combining a matrix calculation method, and introducing a group of bases to represent an instantaneous attribute value of a single-component signal; S4, introducing a sliding real-value frequency filter H(omega) to realize local representation of an instantaneous attribute value of a multi-component cross signal; S5, according to an inner product formula, constructing a statistical linear frequency modulation filter estimation and a time-frequency coefficient used for extracting an original time-frequency spectrum instantaneous frequency and an instantaneous chirp rate curve in the TFC space, and obtaining a statistical linear frequency modulation three-dimensional extraction transform value T x (t, omega, beta); S6, performing modulus operation on T x (t, omega, beta), and obtaining a time-frequency spectrum of the statistical linear frequency modulation three-dimensional extraction transform, so that the purpose of effectively identifying a complex superimposed dense channel sandstone reservoir is achieved.

This user equipment may acquire a first identifier and a second identifier on the basis of a primary synchronization signal and a secondary synchronization signal received on a cell, and acquire a physical cell identifier of the cell on the basis of the first identifier and the second identifier. The first identifier is one of P first-identifier candidates respectively associated with P primary synchronization signal candidates, where P is a positive integer. The P primary synchronization signal candidates may be respectively based on P orthogonal chirp signals, wherein the P orthogonal chirp signals may include at least a first orthogonal chirp signal and a last orthogonal chirp signal among N orthogonal chirp signals of the same chirp rate.

The invention discloses a gravimeter measuring range automatic matching method and computer equipment, and the method comprises the steps: roughly calculating a local gravity value through employing local latitude and longitude and altitude, and carrying out the rough adjustment of the measuring range of a gravimeter; changing the free evolution time of the gravimeter by controlling a time sequence to obtain atomic interference fringes at different evolution time; the chirp rate is used as an independent variable, the atomic transition probability is used as a dependent variable, data fitting is carried out on the atomic interference fringes, and a chirp rate value corresponding to the sine curve valley bottom is obtained; and obtaining a scatter diagram by taking the free evolution time as an independent variable and the chirp rate value corresponding to the valley bottom of the sine curve as a dependent variable, fitting to obtain a public valley bottom through a shortest path algorithm, obtaining a local gravity value, and finishing the fine adjustment of the range of the gravimeter. The problems that an existing method is complex in operation, low in automation level and long in consumed time can be solved.

The application relates to the field of biometric identification technology and discloses a multi-modal handwriting dynamic feature extraction and identification system. The system synchronously acquires macro kinematic sequences and micro acoustic emission sequences, performs time axis regularization based on a bottom-layer hardware interrupt signal, extracts instantaneous geometric singular points of the macro kinematic sequences, combines stress checking to intercept effective micro windows, performs phase space line integral calculation in the effective micro windows to extract macro dissipation work parameters, applies phase compensation to the micro acoustic emission sequences based on corresponding physical coordinates to realize beam forming, acquires micro friction frequency response chirp rates, fuses the dissipation work parameters and the chirp rates to generate high-dimensional identification feature vectors, and inputs the high-dimensional identification feature vectors into a classifier to output a judgment result. The application quantizes a real dynamic physical work process of writing, effectively resists high-precision mechanical arm impersonation attacks, overcomes the defects that a single-channel acoustic acquisition is susceptible to bottom noise and dispersion interference, and improves identification stability in a complex environment.

The present invention discloses a radar device and a signal processing method for the radar device. The radar device includes a transmitter module, which is configured to generate a transmission wave, including: generating a first chirp chain at a first chirp rate for the transmission wave to be output, including: generating a first transmission signal including at least one modulated signal to be output at a first angle; generating a second transmission signal to be output at a second angle different from the first angle; and generating a second chirp chain at a second chirp rate for the transmission wave to be output, including: generating a third transmission signal including at least one modulated signal to be output at the first angle; and generating a fourth transmission signal including at least one modulated signal to be output at the second angle, wherein the first chirp rate is different from the second chirp rate.

The invention provides a chirp fiber grating preparation device, a preparation method and an optical fiber filter, a femtosecond laser direct writing technology is utilized to prepare a chirp localization grating of an asymmetric structure in a fiber core of a single-mode fiber, and therefore the ultra-wide-band optical fiber filter is achieved. On one hand, chirp is introduced on the basis of the localized grating, on the other hand, inscribing of the point-by-point chirp grating is controlled, preparation of the chirp fiber grating with the controllable initial period and chirp rate is achieved, and the filter with the working bandwidth capable of being flexibly controlled is obtained. The optical fiber filter provided by the invention has larger bandwidth and more flexible bandwidth regulation and control capability, can be prepared conveniently, and is suitable for various broadband filtering application scenes such as a super-continuous light source, a Raman spectrum system, optical coherence tomography and the like.

The present application relates to the technical field of fiber grating preparation, and particularly relates to a system and a method for automatically preparing a chirped fiber grating, the system comprising a laser, an electrically-controlled optical switch, an optical power attenuator, a dichroic mirror, a microscopic objective, and a three-dimensional displacement platform provided with a rotating clamp along an optical path; the system further comprises a control device electrically connected to the laser, the electrically-controlled optical switch, the optical power attenuator, and the three-dimensional displacement platform, respectively; and a CCD camera arranged on a side of the dichroic mirror away from the microscopic objective. The system is used to prepare the chirped fiber grating, and the system's optical pulse timing control and gray-scale lithography are realized through the cooperation of the system components, so that the fiber grating with different apodization modulation functions and chirp rates can be prepared at low cost and high efficiency; and the control of the laser, the electrically-controlled optical switch, the optical power attenuator, and the three-dimensional displacement platform by the controller improves the stability of the laser pulse output power and realizes the customized preparation of complex spectrum shapes.

The application discloses a broadband laser pulse chirp rate measuring method and device. The method is as follows: a broadband chirp pulse is divided into two beams, a reflected beam keeps original beam characteristics as a measured light, a transmitted beam becomes narrow pulse width and narrow spectrum light as a probe light after filtering, the two beams are incident in a non-collinear way in a nonlinear crystal to generate sum frequency through space-time coupling, a crystal is optimally adjusted to obtain a maximum intensity sum frequency signal, a spectrometer is used to measure the sum frequency spectrum and record the data; a probe light path is adjusted to change a relative synchronization relationship with the measured light, so that the probe light scans from a front edge to a rear edge of the measured light, and a delay path amount and a corresponding sum frequency spectrum are recorded; and the chirp rate is calculated according to a corresponding relationship between the sum frequency spectrum and the path delay amount. The device for realizing the method comprises a beam splitter, three mirrors, a narrow-band filter, a path delay device, a nonlinear crystal, an electric angle adjusting frame, an absorber group, a lens and a spectrometer.

The invention relates to the technical field of optical detection, and particularly provides a coherent anti-Stokes Raman scattering imaging method and device, and the method comprises the steps: splitting initial light emitted by a driving light source into driving light and Stokes light; processing the driving light into signal light based on an optical parameter amplifier; performing pulse broadening processing on the signal light on a first light path of the signal light, or performing pulse broadening processing on the Stokes light on a second light path of the Stokes light; combining the signal light and the Stokes light to a third light path, and performing pulse broadening processing on the signal light and the Stokes light on the third light path to obtain the signal light and the Stokes light with matched chirp rates; and performing coherent anti-Stokes Raman scattering imaging on the target sample based on the signal light and the Stokes light which are matched with the chirp rate. The problem that the hyperspectral imaging quality and the hyperspectral imaging speed cannot be considered at the same time in the prior art is solved.

This signal identification device (1) is provided with: a sweep cycle estimation unit (12a) that estimates a sweep cycle in which a sweep of a linear frequency modulation (LFM) wave in a digital reception signal obtained by converting an incoming wave comprising an LFM wave into a digital signal is repeated; a chirp rate estimation unit (12b) that estimates a chirp rate, which is the frequency change rate of the LFM wave in the digital reception signal; and an initial frequency estimation unit (12c) that estimates the initial frequency of the LFM wave in the digital reception signal.

The application discloses a wireless mobile communication method based on generalized orthogonal chirp division multiplexing, and relates to the field of wireless mobile communication. The method comprises the following steps: GOCDM modulation, double-selected channel modeling, received signal modeling, GOCDM demodulation under the double-selected channel, and channel equalization. The GOCDM communication system is applicable to a double-selected channel (i.e. a time-varying multipath channel). Each orthogonal base wave of the GOCDM is a Zadoff-Chu (ZC) sequence. The ZC sequence is determined by a sequence length N and a root coefficient l. The waveform of the GOCDM can be regarded as a digital chirp signal. When the N is fixed, the root coefficient l determines the chirp rate, namely the digital frequency modulation rate. The root coefficient l increases the flexibility of system design. By changing the l, the digital frequency modulation rate can be changed. When the l of the GOCDM is 1, the GOCDM is equivalent to a conventional OCDM system. Since the orthogonal base wave of the GOCDM is a digital chirp signal, the GOCDM has the characteristics of time-frequency domain double diffusion. The application improves the reliability of high-speed wireless communication on the double-selected channel.

To provide a radar device capable of reducing distance calculation errors.SOLUTION: A radar device includes: a light source 20 that outputs an optical signal including a chirp signal; a transmitting / receiving unit 30 that transmits an output signal of the light source 20 and receives the optical signal to generate a received signal; a beat frequency measurement unit 40 that generates a first beat signal by causing interference between the output signal and the received signal and measures a beat frequency of the first beat signal; a chirp rate measurement unit 50 that measures a chirp rate of the output signal; and a ranging calculation unit 60 that calculates a distance to a target by using the beat frequency measured by the beat frequency measurement unit 40 and the chirp rate measured by the chirp rate measurement unit 50.SELECTED DRAWING: Figure 1

The examples herein provide enhanced optical architectures employing techniques and systems for scaling up a signal bandwidth and correcting errors in optical signals to achieve a corrected waveform that can be adjusted for the desired application with broad tuning of bandwidth, chirp rate, chirp duration and symmetry. An optical signal is provided for circulation in an optical loop, and incremental modulations are applied to the optical signal over a quantity of circulations to establish circulating waveforms that form a composite waveform. Detected errors in the circulating waveforms can be corrected during the incremental modulations. In some instances, a portion of the optical signal comprising the composite waveform can then be directed for emission toward a remote target and remote sensing of properties of the remote target.

The invention discloses a matching extraction transformation method which is used for deep carbonate rock sedimentary cycle division. Comprising the following steps: inputting a to-be-analyzed seismic signal; calculating instantaneous frequency (IF) and group delay in a short-time Fourier transform domain; defining a chirp rate effective estimation region (CVD), and solving partial derivatives of time by the instantaneous frequency and the group delay to obtain a matched chirp rate estimator; chirplet transformation is calculated in CVD, and IF is estimated by using a local maximum criterion; extracting the time-frequency coefficient at the IF to obtain matching extraction transformation (MET); performing time-frequency transformation on the seismic data according to the MET to obtain a time-frequency result; and in a time-frequency result obtained by MET, analyzing the increase and decrease trend of the tuning frequency so as to characterize the deposition cycle of the deep carbonate rock.

The invention relates to a linear frequency modulation pulse detection method, which realizes linear frequency modulation pulse detection according to a linear frequency modulation pulse frequency modulation slope invariant characteristic and a phase, frequency and time correlation characteristic. The specific processing process is as follows: firstly, obtaining a data instantaneous phase of node acquisition data through Hilbert transform, and performing phase compensation from front to back; secondly, the instantaneous phase forward difference is obtained through a phase difference method, and frequency data related to the linear frequency modulation pulse starting frequency, the frequency modulation slope, the pulse starting moment and the pulse current moment can be obtained through phase difference data according to the linear frequency modulation pulse frequency modulation slope invariant characteristic and the phase, frequency and time correlation characteristic; performing forward difference on the frequency data to obtain frequency difference data only related to the frequency modulation slope; and finally, the reciprocal of the variance of the frequency difference is used as a detection quantity to realize detection and parameter estimation of the linear frequency modulation pulse, so that the stability of linear frequency modulation pulse detection is improved.

This invention discloses a method for classifying deep carbonate sedimentary cycles based on matched Chirplet compression transform (MCST). The method includes: inputting the seismic signal to be analyzed; calculating the instantaneous frequency and group delay in the short-time Fourier transform domain to obtain the chirp rate function; deriving the IF band center and constructing the time-frequency energy focusing domain in the chirplet transform domain; accurately obtaining the local rate of change of the instantaneous frequency trajectory in the chirplet transform domain by constructing a matched Chirplet rate estimator; redistributing the time-frequency coefficients to the true instantaneous frequency trajectory in the time-frequency energy focusing domain through compression transform, thereby obtaining the time spectrum of matched chirplet compression transform (MCST); applying this method to perform time-frequency characterization of seismic trace data, and classifying deep carbonate sedimentary cycles by analyzing the changing trends of the tuned frequencies in the time spectrum.