Application of Simulated Annealing Particle Swarm Optimization Based on Correlation in Parameter Identification of Induction Motor

The strengths and weaknesses of correlation algorithm, simulated annealing algorithm, and particle swarm optimization algorithm are studied in this paper. A hybrid optimization algorithm is proposed by drawing upon the three algorithms, and the specific application processes are given. To extract the current fundamental signal, the correlation algorithm is used. To identify the motor dynamic parameter, the filtered stator current signal is simulated using simulated annealing particle swarm algorithm. The simulated annealing particle swarm optimization algorithm effectively incorporates the global optimization ability of simulated annealing algorithm with the fast convergence of particle swarm optimization by comparing the identification results of asynchronous motor with constant torque load and step load.


Introduction
A certain motor parameter identification model is used, and the residual square sum of the predicted output value and the actual motor output value at the discrete time of the motor model under the parameter identification value serves as the minimization objective function in the identification of motor parameter follows for iterative optimization [1][2][3][4][5].In recent years, plenty of studies using generalized Kalman filtering, least square method, and genetic algorithm have been performed on the identification of asynchronous motor parameters [6][7][8][9].The weakness of the generalized Kalman filter method in these methods is that it requires considerable computation in every step of vector or matrix operation.The objective function of the motor parameters of the derivative is used by least square method in the optimization, and speed fluctuations are very susceptible.Genetic algorithm is easier to reach local optimum, which makes it difficult to have higher recognition accuracy in terms of convergence.Besides, the fundamental information will be impacted by the noise signal component in the stator current.Given these, the more effective algorithm should be studied to identify asynchronous motor parameters.
The correlation algorithm is used in this paper to extract the amplitude and phase of the stator current fundamental signal, and the harmonic components are filtered out to emphasize the fundamental signal given the characteristics of the stator voltage and current fundamental frequency equal to construct the same frequency with the stator voltage reference signal [10][11][12].Then, to identify the parameters of the dynamic mathematical model of asynchronous motor, a novel optimization algorithm incorporating simulated annealing with the advantages of particle swarm optimization is used [13][14][15].Some particles in the particle swarm get a new target position by flight in the identification iteration.To obtain a new target position, some particles are randomly sampled following simulated annealing [16][17][18][19][20].By calculating the acceptance probability of a new position, each particle can decide whether to reach a new position to ultimately update the position status of the entire particle swarm.The ability of simulated annealing algorithm to jump out of the local optimal value effectively incorporates the rapid optimization of particle swarm optimization, as suggested from the simulation results, which greatly increases the accuracy of parameter identification of the dynamic mathematical model of asynchronous motor. Hypothesis where  denotes the initial phase angle of () at time t = 0;  presents the phase difference between () and ().
Given that the period T is a finite value, the estimated value of the cross-correlation function   () for () and () is expressed as Cross-correlation analysis of two periodic signals with the same frequency maintains both the same frequency component and the phase information, as suggested from (5).As the nonfrequency periodic signals are irrelevant, the same frequency of the reference signal and the measured signal can be processed using cross-correlation.Given that the interference signal differs from the reference signal frequency, the interference can be eliminated after the measured signal and the amplitude and phase information can be extracted.

Related Algorithms in the Application of Noise Signal
Cancellation.The fundamental components of the stator voltage and the stator current signal are considered to have the same frequency when the grid has sufficient capacity.A signal with the frequency identical to that of the stator voltage can serve as a reference signal.The noise component and the reference signal have different frequencies by crosscorrelating the grid voltage signal and the stator current signal of the motor.Accordingly, the amplitude and phase information of the stator current fundamental component can be obtained.The implementation is detailed as follows: the stator current signal of the measured motor is denoted as (), the phase of the fundamental component expressed as , the amplitude is represented by A, the angular frequency refers to , and   () denotes the harmonic component of the stator signal and the other components of the noise.The following is available: Two-reference signals with the same frequency as the stator voltage are created, and the rated current of the motor to be measured is the effective value and B is the reference signal, as expressed below The stator current signal of the measured motor () goes through the cross-correlation with the reference signals () and  1 (), respectively.Given that   () is not associated with the two-reference signals () and  1 (),  = 0, the estimated value of the correlation function between the measured signal and the reference signal () and  1 () is expressed as Solving (9) to (11) results in (12) to ( 14) Assume that the number of sampling points in the period T is N, the discrete time series of the stator current signal s (t) expressed as (  ), and the discrete time series of the two-reference signals (),  1 () are denoted as (  ),  1 (  ), respectively   = 0, 1, 2 ⋅ ⋅ ⋅ ,  − 1.Then, the following is formed: ) ,   1 , and   are calculated through the discrete time series (), (), and  1 ().By introducing them into (12) to ( 14), the amplitude A and the phase  of the stator current fundamental wave signal can be obtained.And then according to (18) (  ) denotes the fundamental signal sequence, and the   (  ) represents the harmonic and noise signal sequence.
SVD filtering is an algorithm on the basis of the singular value decomposition.Keep The information of induction motor can be known through the parameter identification of   (  ).

SA-PSO
3.1.SA Algorithm.SA (short for Simulated annealing algorithm) refers to a heuristic random search method following the strategy of Monte-Carlo iterative solution.Using the Metropolis sampling criterion characterized by probabilistic sudden jump, the SA is first searched randomly in the solution space at a comparatively high initial temperature.The sampling process is repeated with the continuous decrease of temperature, and eventually the global optimal solution of the problem is obtained.
The general steps of the SA algorithm are as follows: (1) optional initial solution x 0 ; x i = x 0 ; k = 0; t 0 = t max (2) selecting a value x j randomly from the field N(x i ) and calculating (3) calculating; terminating the calculation if the termination conditions are satisfied; otherwise, back to step two.

PSO Algorithm.
Each latent solution in D-dimensional space is deemed as a particle in PSO (short for particle swarm optimization) algorithm.m particles are first generated in the solution space in the optimization calculation [21][22][23].The position vector of the ith particle in the particle swarm is x i = (x i1 , x i2 , . . ., x iD ), and the velocity vector is vi = (v i1 , v i2 , . . ., v iD ).The optimal position searched by the particle is denoted by p i = (p i1 , p i2 , . . ., p iD ).The optimal position searched by the entire particle group is denoted by p g = (p g1 , p g2 , . . ., p gD ).Particle swarm optimization algorithm can be expressed below: x where t denotes an iterative algebra; c 1 and c 2 represent a learning factor, generally c 1 = c 2 = 2, and rand( ) refers to a random number between 0 and 1 independent of each other.There is a maximum speed of particle flight [24][25][26].The maximum speed is used instead when the calculated speed surmounts this maximum.
An inertia weight term w is added to the basic evolution equation to obtain wPSO given that the flight speed of particles affects the global convergence performance of the algorithm.The improved equation is as follows: In the basic evolution equation to increase the contraction factor r CPSO algorithm, the iterative equation is expressed as 3.3.SA-PSO Hybrid Algorithm.This paper proposes an optimization algorithm following the particle swarm optimization algorithm and simulated annealing algorithm given the characteristics of SA and PSO.The two algorithms are mixed and iterated in the calculation.Using the Metropolis criterion of simulated annealing, the acceptance of particles for each of their leapfrog results is modeled.The optimal solution in a particle swarm serves as the initial solution to the next temperature sampling process of the simulated annealing algorithm.The number of particles reaching a certain level of similarity to the optimal particle is considered as the length of the randomly sampled Markov chain at that temperature, and they are not flown.Yet a new solution is generated by the simulated annealing random disturbance to reach a new location, which increases the diversity of particle populations to avoid falling into the local optimum.
As for SA-PSO particle swarm optimization, Fast Convergence, and global convergence of simulated annealing algorithms, the "premature" phenomenon of PSO is effectively overcome, and the convergence speed is fast by the powerful combination of the two-algorithm search.

SA-PSO Used in Asynchronous Motor
Parameter Identification

Asynchronous Motor Dynamic Mathematical Model.
Motor voltage equation is Electromagnetic torque equation is Equation of motion is In ( 27)-( 29), u ds and u qs denote the stator voltages of d and q axes, u dr and u qr denote d and q-axis rotor voltage on the stator conversion value, i ds and i qs denote d and q-axis stator current, and i dr and i qr , denote d and q-axis rotor current stator converted value.R s refers to the stator's resistance per phase.R r represents the rotor's commutation value per stator resistance.Ls denotes the inductance of each phase of the stator.L r represents the rotor's inductance per phase of the stator conversion value.L m refers to the magnetizing inductance of each phase of the stator.X r denotes the rotor electrical angular velocity.T e refers to the electromagnetic torque.p represents the asynchronous motor pole pairs.T r denotes the load torque.B represents the damping coefficient.J refers to the moment of inertia.
To be identified parameter vector is Input vector is State vector is where X denotes the state variable; V refers to the input variable; H represents the identification parameter variable; D denotes the differential operator,  = /.

The identification parameter variable matrix (N particles) is denoted as
where  i = [ i1  i2 . . . im ] suggests that the ith particle is encompassed by m identification parameters.
where Y(t k ) denotes the measured value.C represents the coefficient matrix associated with selecting those measurements.k refers to the order of each measurement during the measurement.t k denotes the instantaneous time for the kth measurement.The entire measurement starts from t k = 0 and ends at t k = TN.
Two different sets of measurements are selected, i.e., where i as denotes the A-phase winding current; i bs represents the stator B-phase winding current; i cs refers to the stator C-phase winding current.The corresponding T , respectively, is Assume that iter is an iterative algebra in SA-PSO algorithm and the steps of applying the algorithm to perform the foregoing parameter identification of dynamic mathematical model are as follows.
Step 1. Randomly generate the initial generation particle group all individuals P(0) = [ θ(0) , in which each particle identification parameter is under the following conditions: where θ()  denotes the jth identification parameter of the ith particle in the first iteration. min  and  max  refer to the upper and lower limits of the jth identification parameter of all particles, respectively.The upper and lower limits should comply with the experience given.The parameter is from iter = 0.
Step 2. Introduce each particle θ()  into (29), find the state variable X(  ), and then introduce (34) into Ŷ(  ).Then find the error between the measurement result and the calculation result: Obtain the objective function value of the particle θ() where ∧ denotes the identity matrix; t k refers to the discrete computation time exactly identical to the measurement time.
Step 3. Calculate the objective function value of all individuals in a particle swarm.Calculate the initial temperature T emp0 if t iter = 0, so that the initial acceptance rate of various particles in the initial iteration is large (∑   [exp(−Δ  / 0 ) > rand( )] > (0.95).At least 95% of the particles pass ( 24) and ( 25) can be accepted as a new location after iteration).
Step 4. Determine the optimal position of each particle  pbest and the global optimum position  gbest .The sampling initial value of the simulated annealing algorithm is established at  gbest .
Step 7. Calculate the respective objective function value for N new particles.And determine whether to accept the new state by following the Metropolis criteria.Number of iterations increases by 1. Perform the annealing operation.Annealing temperature selection equation Step 8. Determine whether to reach the maximum number of iterations; otherwise turn Step 2; output the optimal value.

Experimental Equipment.
The sampling frequency of the stator current signal is 8kHz in this paper.Table 1 lists the induction motor parameters.

Induction Motor Parameter Identification under Load
Torque Change.In most cases the load is changing given the actual asynchronous motor operation.Thus, the load torque is merely changed, and the three classic loads are described, such as the constant torque load, the step load, and the constant power load to study the performance of the CSA-PSO algorithm under different loads.

Induction Motor Parameter Identification under Load
Torque Change.Constant torque load refers to the load torque T L at any speed being always constant or basically constant, as well as the speed-independent load.Constant torque load characteristics primarily fall into resistive torque load characteristics, the potential of the torque load characteristics, and potential energy load characteristics of three.
The constant torque load is mostly resistance under normal circumstances; i.e., the load torque T L polarity changes with the direction of the speed variation.Such load torque depends merely on the weight of the load, irrespective of the speed, and the load characteristics in the one or three quadrants.
The characteristic taken on by the potential load is opposite to that of the reactive load; i.e., the polarity of the load torque T L does not change with the change of the rotational speed.The potential energy load characteristic is similar to that of the potential load.The load torque T L takes on a fixed direction, and its polarity is not varied with the variation of the rotational speed direction.Its load characteristic is in a quadrant or quadrant.In the simulation of asynchronous motor speed of 1000rpm, there is a constant torque load of 20Nm.
Asynchronous motor parameter identification results under constant torque load curve, which are stator resistance, rotor resistance, stator inductance, rotor inductance, and mutual inductance identification curve are presented in Figures 1-5, the stator resistance identification results are presented in Figure 1, the rotor resistance identification results are presented in Figure 2, the stator inductance identification results are presented in Figure 3, the rotor inductance identification results are presented in Figure 4, and the result of mutual inductance is identified in Figure 5.Because the stator current contains the noise signal, the value of the objective function in the induction motor parameter identification is always in a state of fluctuation.At this point, the identification result of the induction motor parameters has been fluctuating.After the stator current is filtered, the recognition results of the induction motor parameters are quickly converged to the identification value after a brief fluctuation and the error between the identification value and the actual value is very small.When there is no stator current filtering, the induction motor parameter identification result is always in a fluctuation state and the error between the identification value and the actual value is large.After 60 iterations of identification, the identification parameters are stable.The accuracy of the induction motor parameter identification is notably greater than when the stator current is not filtered as the stator current is filtered.The identification of the motor parameters stability error is small, as suggested from the identification curve and the actual curve comparison.

Parameter Identification for
Step Load Torque Variation.The step input is usually the most unfavorable input situation under the actual situation.Practically, many inputs are similar to step inputs.As a result, the system performance index  of the control system is usually obtained by analyzing the step response of the system as a stepping number is applied.Similarly, fast tracking and accuracy can be provided using the load torque as an input to the control system if the system's step response is high.It is suggested that CSA-PSO algorithm identification method is applicable under the step load.
The speed of asynchronous motor was 1000rpm, and the load torque suddenly increased from 10Nm to 20Nm at 0.5s in the simulation.The curve of asynchronous motor parameter identification results is presented in Figures 6-10 under step load, i.e., stator resistance, rotor resistance, stator inductance, rotor inductance, and mutual inductance identification curve.The stator resistance identification results are presented in Figure 6, the rotor resistance identification results are presented in Figure 7, the stator inductance identification results are presented in Figure 8, the rotor inductance identification results are presented in Figure 9, and the result of mutual inductance is identified in Figure 10.
At 0.5 s, the load torque changes from 10 Nm to 20 Nm.Correspondingly, the output torque increases and the output current also increases.In the armature circuit, the voltage drop across the inductor and resistor will increase, and the effect of the skin effect will increase, causing the induction motor parameters to change.However, it can be seen from Figures 6-10 that the difference between the identification result of the induction motor parameter and its actual value is not significant.In the case of abrupt changes in load torque, the CSA-PSO algorithm can accurately search the true value of the estimated parameter after 60 iterations.After 100 iterations, relatively stable output results can be obtained.Therefore, it can be considered that the step load disturbance has almost no effect on the convergence speed and accuracy of the CSA-PSO algorithm.After the stator current is filtered, the accuracy of the induction motor parameter identification is significantly higher than when the stator current is not filtered.According to the identification curve and the actual curve comparison, the identification of the motor parameters stability error is small.Accordingly, the disturbance of step load exerts no effect on the convergence speed and precision of identification.
.1.Correlation Function Theory Overview.The dependency of the instantaneous value () of a time t and the instantaneous value ( + ) is expressed by the correlation function at another time.For periodic signals, the integral average time T is the signal period.The relevant autocorrelation function   () and cross-correlation function   () are defined as 2 large singular values (complying with energy-focused signals or large energy signals such as harmonic signals).Discard small singular values (complying with energy dispersive signals or small energy signals, such as noise signals).Singular value decomposition is performed on   (  ), and large singular values are reserved.To obtain a pure harmonic component   (  ), harmonic restoration is performed based on these large singular values.The induction motor current signal is   (  ) as expressed below: