123results about How to "Shorten calculation cycle" patented technology

The invention discloses a probability calculation-based multiple input multiple output detector, which comprises a matrix QR composer, the matrix QR composer is respectively connected with two random sequence generators, the two random sequence generators are respectively connected with a probability complex multiplier, one of the probability complex multipliers is connected with a Gibbs sampling update unit, and both the other probability complex multiplier and the Gibbs sampling update unit are connected with a log-likelihood ratio calculation unit. Because the multiple input multiple output detector uses probability calculation to carry out the MCMC (Markov chain Monte Carlo) algorithm, the complexity of operation is greatly decreased, the transition probability of the Markov chain in the MCMC algorithm is increased, and the problem of locking under a high signal-to-noise ratio is solved. A sliding window sequence generation method is utilized to carry out Gibbs sampling update, so the length of a probability sequence is reduced. When the multiple input multiple output detector is applied to construct a full-parallel detector, a high throughput rate can be achieved with the full-parallel mode.

In one embodiment, the present invention discloses an apparatus and method for providing efficient implementations of Finite Impulse Response (FIR) digital filters. Specifically, a result from a FIR digital filter can be efficiently computed by using an AVG operation or instruction in conjunction with one or more other operations. The unique use of the AVG operation will allow FIR filters of various types, e.g., Types 1-4, to significantly reduce computational cycles.

The invention relates to an interface method in the field of electric power system dispatching, in particular to a standard data exchange interface method of a parallel cooperative system. Data transmission between a database in the parallel cooperative system and an external data source is carried out through a data exchange interface; transmitted data adopt standard data formats; standard data adopt power grid data in E-language formats. The standard data exchange interface method comprises the following steps of (1) initializing an Oracle database; (2) reading a file in the E-language format to a memory; (3) copying in a list-to-list way; (4) writing the data to the Oracle database. According to the standard data exchange interface method disclosed by the invention, the problem of importing the data in different formats to the database is solved, and the defects such as manual data conversion and repeated operation which are caused by discrete of a daily traditional mainstream algorithm can be eliminated.

The invention discloses a model parameter adaptive method for the inverter model prediction control, and the method comprises the steps: firstly substituting system parameters calculated at a last moment into a prediction model, and calculating prediction current values corresponding to all vectors; secondly building a cost function according to a control requirement, selecting an optimal vector which enables the cost function to be minimum, and applying the optimal vector to a control system; thirdly carrying out the online correction of parameters in the model through voltage and current measurement values at this moment and all parameter values at the last moment after the optimal vector is selected. If inductance and resistance and other parameters need to be corrected on line, the method can correct the parameters at different moments so as to avoid excessive calculation burden caused by formula iteration. After the model prediction control is finished, a parameter recognition algorithm is added, and the method carries out the online correction of all parameters in the model through the prediction model of the system, and improves the operation robustness and control flexibility of an inverter.

The invention discloses a non-linear load calculation method of a flying wing unmanned aerial vehicle, which comprises the steps of solving flight attitude parameters. A multi-degree-of-freedom coupling flight dynamics simulation model is adopted for aerodynamic and dynamic flight calculation of the aircraft based on aerodynamic non-linear coefficients in solving the flight attitude parameters, Atthe same time, an objective function is established to optimize the flight attitude parameters which do not meet the requirements by using the optimization algorithm, and the output flight attitude parameters are used as the input of the nonlinear concentrated load calculation. As for that aerodynamic non-linearity characteristic of the whole aircraft of the unmanned aerial vehicle (UAV), the invention adopt a flight dynamics simulation model of multi-degree-of-freedom coupling to solve the flight attitude parameters of three-axis couple non-linearity on the basis of the aerodynamic non-linearity coefficient; Thus, the nonlinear load calculation of flying-wing UAV is realized, which has good practicability.

An industrial model predictive control method realizing dynamic optimization based on linear programming comprises the following steps: 1) carrying out steady state target optimization: computing the expected target value output in a controlled way in the steady state according to the input quantity of the sampling period; and 2) carrying out dynamic optimization: 2.1) obtaining step response models of the industrial process by a least square method, wherein the step response models are established in N-numbered periods; 2.2) predicting the output quantity at each future moment from the k+1 moment, referring to formula (2); 2.3) ensuring the input and the output in the dynamic process to be in own upper and lower bound ranges, referring to formula (3), setting target constraint, referring to formula (4) and creating the following performance indexes, referring to formula (5); forming the dynamic optimization problem by combining the formulas (2), (3), (4) and (5), solving all the control increments deltau(i) by linear programming, wherein i=0, 1,..., N-m-1, and then substituting the control increments into the formula (2) to compute the predictive output values y(k+1), y(k+2),..., y(N). The method simplifies computation, shortens the computation period and is good in practicability.

The invention discloses a live broadcast platform user interestingness mining system and mining method based on association rules and relates to the technical field of data mining. The system comprises a data acquisition module for acquiring live broadcast platform user behavior information from a server and generating sampling transaction event library and a testing transaction event library, an association rule model establishing module for adopting a spark calculation framework to conduct frequency mode mining on the testing transaction event library to obtain an association rule model, a user interestingness output module for using the testing transaction event library as an input variable to input the association rule model, obtaining an output variable of the association rule model and using the output variable as user interestingness contents. The live broadcast platform user interestingness mining system has the advantages of being short in computing cycle and high in testing result accuracy and practicability.

The invention discloses a method for acquiring discrete state parameters of a power transformer, which solves problems of difficult power transformer state parameter recognition and difficult power transformer multi-source data mining. According to the invention, a discrete gas state parameter is acquired by adopting an improved equal width method according to the type and the content of dissolved gases in tested power transformer oil; a discrete partial discharge state parameter is acquired by adopting an information entropy decision tree method according to a discharge capacity statistical parameter, a discharge phase statistical parameter and a number of discharge statistical parameter; a discrete winding state parameter is acquired by adopting a binning method according to winding deformation and short-circuit current; and a discrete mechanical state parameter is acquired by adopting a chi-square test method according to the transformer temperature, vibration and noises. The method disclosed by the invention is applicable to the fields of power equipment detection and data mining, has the characteristics o f accurate data, high work efficiency, wide fault coverage, easy computer processing and the like, and has high economic values and broad market application prospects.

The invention relates to a method for calculating the static strength of a cabin door structure and belongs to the field of airplane structures strength design. Static strength calculation is performed based on a linear contact analysis method, and a finite element model of the cabin door structure is firstly established; according to the cooperative relation of the cabin door structure and an airplane body structure, boundary conditions are determined; then, the loading conditions applicable to constraint conditions of the cabin door structure are determined according to the boundary conditions, finally constraint of a linear contact analysis model is performed, a constraint control file is established, and finally the result of linear contact constraint reaction is calculated. By comparing a calculation result of linear contact analysis with a calculation result of nonlinear contact analysis, the static strength of the cabin door structure can be more conveniently, quickly and accurately calculated through linear contact analysis.

The invention relates to the technical field of explosion damage assessment, in particular to a Tianhe supercomputer-based city explosion damage effect assessment method. The method comprises the steps of extracting typical analysis characteristics of city buildings through a BIM; establishing and solving a setting rule according to a simulation model; developing modeling and solving display programs on a simulation software platform and a Tianhe supercomputer platform; realizing automatic simulation model building, solving calculation and display function customization; and realizing city explosion damage effect assessment. According to the method, large-scale and high-precision multi-task parallel calculation of the complex simulation model can be performed based on a Tianhe-1 supercomputer, so that the calculation cycle is shortened; by virtue of the BIM, accurate city geometric information can be quickly extracted, and a corresponding geometric model is generated according to analysis needs; and based on this, a simulation solving model is automatically built according to the analysis needs, and corresponding software and hardware resources are matched, so that the time for adjusting the model by personnel is shortened and relatively good calculation configuration is achieved.

The invention discloses an MATLAB based cycloid gear parameter optimization method. The method comprises the following steps: step1. obtaining number ZP of needle teeth and needle teeth central circle radius rp; step2. according to restricting conditions, obtaining an eccentric distance and a needle gear sleeve outer circle radius rrp; step3. determining tooth number range and maximum engagement force involving in engagement simultaneously at the moment when the cycloid gear transmits; and step4. when needle teeth contact stress reaches strength requirement and rotary arm bearing life meets the requirement, outputting the eccentric distance a and the needle gear sleeve outer circle radius rrp, otherwise, executing step 2 until the requirement is met. According to the cycloid gear parameter optimization method, cycloidal gear parameter optimization is completed by analyzing basic principles of cycloid pin gear transmission and a design method of the cycloid gear. In the optimization, main parameters affecting cycloid gear transmission carrying capacity, the eccentic distance and the needle tooth radius rrp, are used as design variables, cycloid pin gear engagement strength and rotary arm bearing life are to-be-solved targets, and certain constraints are met.

The invention discloses a clustering center selection method. The method includes the following steps that: the local densities of sample points in a sample set and high density distances between eachsample point and sample points with higher local densities are obtained on the basis of a plurality of threads set by the OpenMp; on the basis of preset rules, a density threshold is calculated withthe local densities adopted as parameters, and a high density distance threshold is calculated with the high density distances adopted as parameters; a target sample point of which the local density and the density threshold satisfy a first preset relationship and of which the high-density distance and the high-density distance threshold satisfy a second preset relationship is selected from the sample set and is adopted as a clustering center. With the clustering center selection method of the invention adopted, the time complexity of calculation is reduced, a calculation period is correspondingly reduced, and therefore, the execution efficiency of a clustering algorithm is ensured. The present invention also provides a clustering center selection device and a medium. The clustering centerselection device and the medium have the same advantages as the above method.

The present invention discloses a C4.5 decision tree algorithm-based specific user mining system and a method thereof, and relates to the live broadcast website data mining field. The system comprises a sample selection module used for selecting the user samples, wherein the user samples are divided into the paying user samples and the non-paying user samples according to the class labels; a behavior attribute statistics module used for gathering the attribute values of the classification attributes of the user samples; a sample processing module used for carrying out the normalization processing on the attribute values of the classification attributes as the training sample data; and an algorithm platform used for receiving the training sample data, wherein the algorithm platform comprises a C4.5 decision tree algorithm, provides an algorithm interface for the C4.5 decision tree algorithm, and trains a C4.5 decision tree model based on the training sample data and the C4.5 decision tree algorithm. The C4.5 decision tree algorithm-based specific user mining system of the present invention can gather the user interest degree correlated indexes comprehensively on the conditions of many data dimensions and large data size, and is convenient to mine the specific users.

The invention discloses a system and a method of double-scalar multiplication of a streamlined elliptic curve. The system and the method realize the calculation of the scalar multiplication of an elliptic curve through a structure comprising two or three finite domain square operation units, three finite domain multiplication units, a finite domain addition operation unit, a register set, a plurality of multipath gates and a control logic unit. The system and the method of the double-scalar multiplication of a streamlined elliptic curve have simple interface mode, greatly improved operation speed and no hardware idle time, can meet the encryption requirement of data with different widths, shorten a calculation period of an elliptic curve encryption algorithm and in parallel calculate finite domain element inversion and is quite suitable for an elliptic curve encryption system working with a high throughput rate.

The invention is a large-scale Kirchhoff pre-stack time migration parallel method with no speed-up ratio bottleneck, and the method is used for seismic data processing in geophysical prospecting for petroleum. An output space is subdivided into four dimensions such as wire size and geophone offset, the number of the geophone offsets is calculated according to the minimum geophone offset, the maximum geophone offset and the incrementation parameter of the geophone offset, all calculated nodes are grouped, the numbers of nodes in different groups are even, the nodes in the same group are in consecutive order in the node name, and the pre-stack time migration of one geophone offset component is calculated by one group of nodes in one time. According to the parallel method provided by the invention, the minimum network transfer of the seismic data is achieved, the problem of network transfer intensity of the Kirchhoff pre-stack time migration is completely solved, and the scale of computing hardware is expanded, so that the computing interval of the pre-stack time migration is shortened.

The invention provides an overspeed protection method and device for a train. The method comprises the following steps of obtaining an initial speed limit position point of each speed limit area in apreset number of speed limit areas and a first speed limit value corresponding to each initial speed limit position point when an emergency braking of the train is triggered to obtained a plurality offirst speed limit values; obtaining the current driving position of the train and the corresponding second speed limit value, and obtaining the current driving speed of the train; determining the initial speed limit position point corresponding to a deceleration that satisfies a preset condition in the deceleration of a current driving speed relative to each first speed limit value as a target speed limit position point; determining an emergency braking speed according to the deceleration of the train based on the second speed limit value and the first speed limit value corresponding to the target speed limit position point, and conducting overspeed protection on the train according to the emergency braking speed. The method and device can effectively improve the accuracy of the emergencybraking speed calculation, reduce the calculation period, improve the driving safety and efficiency of the train, and improve the user's riding experience.

The invention discloses a compressible wall function calculation method considering a boundary layer combustion heat release effect, and belongs to the field of computational fluid mechanics. A velocity wall law and a temperature wall law under a combustion chemical reaction condition are proposed through theoretical derivation; based on this, a unified wall function calculation method capable of being compatible with compressible turbulent boundary layers with / without a chemical reaction simultaneously is established; and a detailed process coupled with a CFD program for solving is designed. The calculation method is low in technical realization difficulty and is suitable for the condition with or without the boundary layer combustion heat release effect. The numerical simulation efficiency of frictional resistance / heat flow of a flow passage in a scramjet engine can be greatly improved, so that the calculation cycle is greatly shortened.

The invention provides a reverse time migration method and apparatus for a VTI medium, and a medium. The method comprises: a focus coordinate is extracted in shot gather data and a focus wavelet is generated; according to a forward modeling scheme, the focus wavelet is extended based on a qP wave fluctuation equation along the time forward direction, wherein a plurality of wave field parts are selected for storage during the extending process; extension is carried out by using the plurality of wave field parts based on the qP wave fluctuation equation along the reverse time direction, therebyobtaining a focus wave field; seismic data are extracted from the shot gather data and interpolated seismic data are extended along the reverse time direction based on the qP wave fluctuation equationto obtain a demodulator probe wave field; correlation imaging is carried out on the focus wave field and the demodulator probe wave field to obtain a shot migration profile, and the above-mentioned steps are performed to obtain a plurality of shot migration profiles; and the multiple shot migration profiles are superposed to obtain a final imaging profile. Therefore, under the circumstance that the computing load is reduced substantially, the longitudinal wave component of the same kinematic feature can be obtained and thus the accurate reflection event returning of the imaging profile is realized, and the imaging precision of the complex medium is improved.

The invention discloses an evapotranspiration remote sensing inversion application based on a reflectivity-vegetation coverage two-dimensional space, and belongs to the technical field of surface evapotranspiration remote sensing estimation, and the method comprises the steps: S1, constructing an evapotranspiration remote sensing inversion network based on a convolutional neural network; S2, training the evapotranspiration remote sensing inversion network by using the training sample; S3, inputting the meteorological data and the satellite remote sensing image of the target area in the presettime scale into the trained evapotranspiration remote sensing inversion network to obtain evapotranspiration data of the target area in the preset time scale. According to the evapotranspiration remote sensing inversion application based on the reflectivity-vegetation coverage two-dimensional space, the neural network and the method for accurately obtaining the surface evapotranspiration of the target time scale through the remote sensing means are perfectly combined, the evapotranspiration calculation speed is increased, the calculation period is shortened, operation is easy, and intelligenceis good.

The invention provides a rectangular slab continuous casting cooling bed equipment foundation electric calculation combination fine adjustment method. The rectangular slab continuous casting cooling bed equipment foundation electric calculation combination fine adjustment method comprises the steps that pixels to be finely adjusted are defined and divided to obtain cooling bed equipment foundationpixels and slag runner transition section pixels; electric calculation fine adjustment is conducted on slag runner upper ramp pixels through upper ramp slag runner auxiliary data pixels, setting adjustment is conducted on slag runner body side wall horizontal steel bar pixels based on data extracted by the upper ramp slag runner auxiliary data pixels, and setting adjustment is conducted on upperramp slag runner peripheral pixels, outer edge foundation side wall horizontal steel bar pixels and foundation top low wall steel bar pixels; and summary calculation is conducted on steel bars, and the steel bars subjected to summary calculation are locked. According to the rectangular slab continuous casting cooling bed equipment foundation electric calculation combination fine adjustment method,by continuously exploring the combination of experimental thinking and a software calculation mode, existing functions of software are applied flexibly, a modeling method is innovated, and rectangular slab continuous casting cooling bed equipment foundation electric calculation combination modeling and electric calculation pixel combination modeling are completed.

The invention provides a boundary layer grid solving method of the underwater glider lift-drag ratio at different attack angles. The method aims at determining the optimal sliding angle of an underwater glider, obtaining resistance and lifting force at different attack angles and providing theory support for later selection of a control system and a power source of the underwater glider. Finally, the accuracy of boundary layer grid numerical simulation is verified through an experiment, and the obtained matching condition is good. The analysis method can be popularized to other series of underwater motion machines, and compared with traditional complex formula derivation, the calculation period is saved, calculating efficiency can be improved by two times or above, and the efficiency and other problems of the underwater machines can be predicted and solved more effectively. The method comprises the steps that A, an underwater glider physical model is simplified; b, an underwater glider fluid outer area model is created; C, an underwater glider bonding layer grid model is created; D, numerical solution calculation of the underwater glider lift-drag ratio at different attack angles is carried out; E, data result postprocessing is carried out. The method is used for effective solving of the underwater glider lift-drag ratio.

The invention discloses a dual-path radon transformation based test curve inflection point interval calculating method. By means of a LabView software platform, industrial detection data are generated into a test curve, the test curve is subjected to dual-path radon transformation to fit multiple sections of smooth and regular straight lines, all inflection points of the fit straight lines are gotten by means of an intersecting line relationship, and intervals between adjacent inflection points are calculated. The dual-path radon transformation based test curve inflection point interval calculating method is small in the dependency on smoothness of the curve, and suitable for extracting inflection points of noisy industrial test curves with complicated characteristics; the calculation method is capable of simultaneously performing radon projection transformation on two paths and has a halved counting period and improved operating efficiency by comparison to a linear radon method; compared with traditional algorithms, the calculating method has the advantages that calculation errors are reduced by over ten times, identification accuracy is high, robustness is good, and the calculating method is convenient to achieve. When an industrial test system detects complicated data, the dual-path radon transformation based test curve inflection point interval calculating method can meet the requirement of industrial monitoring and fine detecting.

The present invention discloses a hypersonic flight target flow field model and electromagnetic model coupling method. The method comprises: taking an appropriate volume to contain a target and a surrounding flow field assembly model of the volume, carrying out required electromagnetic grid division on the volume, mapping each flow field grid node to the electromagnetic grids, and taking the meanvalue of flow field node data contained in the electromagnetic grid for the grid containing the flow field nodes in the electromagnetic grid to replace the flow field data at the location of the electromagnetic grid; using the existing data electromagnetic grid linear interpolation technology for the remaining electromagnetic grids to obtain flow field data at the locations of the electromagneticgrids; converting the flow field data at locations of all electromagnetic grids into electromagnetic parameters required by the electromagnetic simulation through the formula; and establishing a plasma Drude model, and converting the plasma Drude model into a time domain for simulation and calculation. The coupling method disclosed by the present invention is more reasonable, and can better show the superiority of the time domain algorithm.

The invention provides a rectangular slab continuous casting equipment foundation serialization computerization combination modeling adjustment method which comprises the steps that a rectangular slabcontinuous casting equipment foundation is divided and defined to distinguish a bale rotary table primitive, an outer arc tangent primitive, a slag flushing channel inlet primitive, a slag flushing channel transition section primitive and a continuous casting cooling bed equipment foundation primitive; defining a large ladle turret primitive, an outer arc tangent primitive, a slag sluice inlet primitive, a slag sluice transition section primitive and a continuous casting cooling bed equipment basic primitive respectively, and performing computerized combination fine adjustment; and summarizing and calculating reinforcing steel bars in the ladle turret primitive, the outer arc tangent primitive, the slag sluice inlet primitive, the slag sluice transition section primitive and / or the continuous casting cooling bed equipment basic primitive, and locking the reinforcing steel bars subjected to summarizing and calculating. According to the method, the reinforcing steel bars meet the strokerule design and specification requirements in the foundation, and the edited stressed reinforcing steel bars are drawn in the model and compared with manual drawing sample turning data to form positive values.

The invention discloses a method for determining the position of a bridge closure section of a cooperative system, which belongs to the technical field of bridge design and comprises the following steps that a full-bridge finite element calculation model is established, and a reasonable bridge forming state is determined; the closure section is preliminarily set to be located on the outer side ofthe side sling, and based on the unstressed state quantity construction control theory, the beam end vertical displacement delta s of the suspension end of the closure section is obtained through finite element calculation; if delta s is within the set interval range, it is determined that the preliminarily set closure segment position is basically reasonable; and if delta s is not in the set interval range, the overlapped section of the cable-stayed section and the suspension section serves as an interval, the position of the closure section is reset through a dichotomy, the full-bridge finite element calculation model is updated, and the vertical displacement delta s of the beam ends at the two ends of the closure section is recalculated till delta s is within the set interval range. Theoptimal closure position can be calculated, the position of the closure section and the closure scheme can be quickly calculated and determined, the operability is high, the calculation period is shortened, and the design efficiency is improved.

The invention provides a simulation system and method suitable for an integrated environment test of electronic equipment in a desert area. The simulation system and method aim to perform environment test simulation on electronic equipment and improve simulation efficiency. The system comprises a preprocessing module, a network model construction module, a simulation calculation file generation module, a simulation calculation and postprocessing module and a file management library module, wherein the five modules and software involved in simulation are all integrated in the system. According to the method, an electronic equipment geometric model is constructed, and parameterized solving programs of different environments and parameterized log files in different forms of different environments are written at the same time; flow field meshing is performed on the electronic equipment geometric model; a simulation calculation log file and a simulation calculation self-defining function UDF file are generated; and the simulation calculation log file and a simulation calculation self-defining UDF function are read to perform automatically-driven simulation calculation on a flow field meshing model, and a result cloud chart is output.

The invention discloses a gust response retarding method of a wing based on a variable-camber flexible trailing edge, and relates to the field of gust mitigation control of flexible aircrafts. The method comprises the steps of firstly, designing a flexible trailing edge structure based on a corrugated plate form for a certain type of wing; then, establishing a wing finite element model to obtain corresponding modal information; thirdly, identifying training data obtained through CFD by an observer method, and creating a wing trailing edge deflection generalized unsteady aerodynamic force reduced-order model; utilizing a surface element method to calculate modal and gust aerodynamic force; establishing a wing open-loop gust response model by the aerodynamic model and the modal information; and finally, adopting a flexible trailing edge as a control surface, taking the wingtip acceleration as a feedback input signal, calculating and designing a closed-loop control law by a GPC method, and adding a closed loop based on open-loop response to obtain an overall wing gust response closed-loop active control system. The variable-camber flexible trailing edge is adopted to slow down the gust response of the wing, the gust response can be effectively slowed down, and compared with a conventional wing, the variable-camber flexible trailing edge wing has greater advantages in the aspects of aerodynamic characteristics and structural quality.

The invention provides an ensemble Kalman filter dynamic reservoir history matching method based on hyper-sphere transform. The method comprises the following steps: 1, initializing a set of reservoir models, wherein the set comprises static parameters, dynamic parameters and oil well production data; 2, performing hyper-sphere transformation on permeability in the static parameters, and constructing a novel state vector set; 3, inputting each state vector in the novel state vector set into a reservoir simulator for predicting so as to obtain a state prediction value, wherein each state vector and the state prediction value thereof form a prediction set; 4, calculating a Kalman gain matrix according to the prediction set; 5, updating the prediction set according to the prediction set, the Kalman gain matrix and the measured oil well production data, thereby obtaining the updated static parameters, dynamic parameters and oil well production data. According to the method disclosed by the invention, the accuracy of reservoir history matching precision can be improved, and the blindness of manual history matching can be reduced.

The invention belongs to the field of distributed computing, and provides a multiple application function asynchronous concurrent scheduling method, which is applied to a distributed computing management platform. The distributed computing management platform combines the time consumption characteristic and number of the computing tasks of each application function and the scale and performance information of a computer cluster node, realizes the computing operation of each application function by independently setting a proper computing operation scheduling grain size for each application function, and adds the computing operation into the scheduling sequence of the distributed computer management platform to realize the computing task asynchronous concurrent submission of the plurality of application functions, the uniform scheduling of the computing operation of the plurality of application functions and the computing result asynchronous recovery of the plurality of application functions, thereby fully utilizing the computing performance of the computer cluster and reducing computing time.

A pulse compression method for CHIRP signal and a wireless transceiver thereof; pulse compression processing of a cycle C k has the following steps: analytically obtaining all instantaneous frequency components of instantaneous signals and then isolating out a corresponding amplitude; using storage spaces M n to store the amplitude and phase of each of the instantaneous frequency components; setting T n as a length of time from t n to an end time of the corresponding cycle; calculating T n according to formula T n =T-t n ; using another set of storage spaces M T to rearrange amplitudes of the different instantaneous frequency components; storing the amplitude of each of the instantaneous frequency components in a set of T n value corresponding storage spaces of M T space; and accumulating overlapping amplitudes of instantaneous frequency components into another space in the storage spaces M T .