CN110134010A - A kind of power attraction repetitive control using equivalent disturbance compensation servo-system - Google Patents
A kind of power attraction repetitive control using equivalent disturbance compensation servo-system Download PDFInfo
- Publication number
- CN110134010A CN110134010A CN201910307598.1A CN201910307598A CN110134010A CN 110134010 A CN110134010 A CN 110134010A CN 201910307598 A CN201910307598 A CN 201910307598A CN 110134010 A CN110134010 A CN 110134010A
- Authority
- CN
- China
- Prior art keywords
- error
- time
- flight
- equivalent disturbance
- sse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000003252 repetitive effect Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000003247 decreasing effect Effects 0.000 claims abstract description 19
- 230000000737 periodic effect Effects 0.000 claims abstract description 18
- 238000013461 design Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 230000014509 gene expression Effects 0.000 claims description 17
- 230000001629 suppression Effects 0.000 claims description 9
- 230000007423 decrease Effects 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 abstract 2
- 230000005764 inhibitory process Effects 0.000 abstract 1
- 230000009897 systematic effect Effects 0.000 abstract 1
- 230000001360 synchronised effect Effects 0.000 description 13
- 238000004088 simulation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
A kind of power attraction repetitive control using equivalent disturbance compensation servo-system, is given module and generates periodic reference signal, construct periodic feedback link, introduced equivalent disturbance compensation in power attracts and restrains, estimated using observer equivalent disturbance;Attract rule building perfect error dynamic based on power, and according to perfect error dynamic design controller, is inputted the signal being calculated as the control of servo-system;Specific attitude conirol can be carried out according to characterization system convergence performance indicator, and give the monotone decreasing region of characterization tracking error convergence process, absolute attractable layer, steady-state error first enter the calculation formula of steady-state error band maximum step number with boundary and tracking error.Power provided by the invention with equivalent disturbance compensation attracts repetitive controller, by the estimation to equivalent disturbance, can be improved systematic tracking accuracy and complete inhibition periodic perturbation.
Description
Technical Field
The invention relates to a power attraction repetitive control method based on equivalent disturbance estimation, which is suitable for a periodic position servo system and other industrial occasions containing periodic operation processes.
Background
When the controller is designed, the internal model principle requires that the closed-loop system comprises an input signal model, namely the input signal model is implanted into the controller to form a feedback control effect, so that the output of the closed-loop system tracks the input reference signal without dead error. Repetitive control provides a controller design method based on an internal model principle, which has 'memory' and 'learning' characteristics to output an error signal to correct the control input of the previous cycle. It can completely suppress periodic disturbance, thereby realizing accurate control. The repetitive control technology is applied to high-precision servo systems such as power electronic circuits, industrial robots, hard disk drives and the like.
The attraction law method directly adopts a tracking error signal, and the controller is more direct and simpler in design. The usual attraction law reflects the desired system error dynamics when perturbations are not considered; for the situation of interference, because of the interference item, the controller designed directly according to the attraction law cannot be realized. The solution is to embed interference suppression measures into the original attraction law, construct an ideal error dynamic with disturbance suppression effect, and design a controller according to the constructed ideal error dynamic equation. Thus, the closed loop system dynamics is determined by the ideal error dynamics and has a desired tracking performance characterized by the ideal error dynamics.
For designing a digital controller by an attraction law method, designing by a discretized continuous time attraction law, and providing performance indexes describing transient and steady-state behaviors of tracking errors by analyzing ideal error dynamic characteristics, the method specifically comprises the following four indexes: an absolute attraction layer, a monotonically decreasing region, a steady state error band, and a maximum number of steps into the steady state error band. The specific values of the indexes depend on the controller parameters and the equivalent interference signal boundaries, so that the controller parameters and the equivalent interference signal boundaries are different, and the values of the three indexes are also different. Once the ideal error dynamic form is given, specific expressions of various indexes can be given in advance for guiding the parameter setting of the controller.
The Extended State Observer (ESO) is a core unit of an active disturbance rejection control system, and the basic method is to define total disturbance (including internal disturbance and external disturbance) as a new state, construct a state observer of an extended state (including total disturbance action) by using a state observation method. The method can estimate the system state, and can also estimate the real-time action quantity of overall disturbance in a system model, so as to compensate the influence of disturbance signals. Since the overall disturbance encompasses uncertainties in the system model, the system model is greatly simplified, and the control gain can also be considered known, facilitating controller design. The extended observer provides a general and practical method for observing uncertain characteristics.
Disclosure of Invention
In order to overcome the defects that the system tracking precision of the existing power attraction repetitive control method is lower and the periodic disturbance cannot be inhibited, the invention provides the power attraction repetitive control method adopting an equivalent disturbance compensation servo system, in order to enable a closed-loop system to have preset expected error tracking performance, a motor servo repetitive controller is designed according to an ideal error dynamic equation of a power attraction structure, the periodic disturbance component is completely inhibited, meanwhile, a disturbance observer is introduced into the closed-loop system to compensate the aperiodic disturbance and further improve the control performance, so that the motor servo system realizes high-speed and high-precision tracking; the invention expands the disturbance effect influencing the system output into a new variable to construct a disturbance observer, the disturbance observer does not need to directly measure a disturbance signal and know a specific model of the disturbance signal, and the invention specifically provides a specific expression of at most four indexes, namely a steady state error band, an absolute attraction layer, a monotone subtraction area and a step number required for a tracking error to enter the steady state error band for the first time, and can be used for guiding the parameter setting of the controller.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a power attraction repetitive controller design method adopting an equivalent disturbance compensation servo system comprises the following steps:
step 1, giving periodic reference signals to satisfy
rk=rk-N(1)
Where N is the period of the reference signal, rkAnd rk-NReference signals respectively representing time k and time k-N;
step 2, defining tracking error
In the formula
A1(q-1)=a1+a2q-1+...+anq-n+1=q(A(q-1)-1)
A(q-1)=1+a1q-1+...+anq-n
B(q-1)=b0q-1+...+bmq-m
Satisfy the requirement of
A(q-1)yk=q-dB(q-1)uk+wk(3)
Wherein e isk+1Represents the tracking error at time k +1, rk+1Reference signal, y, representing the time instant k +1k+1、yk、yk+1-NAnd yk-NRepresenting the output signals at times k +1, k +1-N and k-N, respectively, ukAnd uk-NRepresenting the input signal at times k and k-N, wk+1And wk+1-NRepresenting the interfering signal at times k and k-N, respectively, d represents the delay, A (q)-1) And B (q)-1) Is q-1Polynomial of (a), q-1Denotes a one-step delay operator, n denotes A (q)-1) M represents B (q)-1) Order of (a)1,...,an,b0,...,bmIs a system parameter and b0Not equal to 0, n is more than or equal to m, d is an integer and is more than or equal to 1;
step 3. constructing equivalent disturbance
dk=wk-wk-N(4)
Where N is the period of the reference signal, dkRepresenting the equivalent disturbance signal at time k, wkAnd wk-NRespectively representing interference signals at the k moment and the k-N moment;
expressing the tracking error as
ek+1=rk+1-yk+1-N+A1(q-1)(yk-yk-N)-q-d+1B(q-1)(uk-uk-N)-dk+1(5)
Wherein d isk+1Representing the equivalent disturbance at the moment k + 1;
step 4, designing an observer and estimating equivalent disturbance
Design observer equivalent disturbance dk+1Observing, and compensating equivalent disturbance by the observed value, wherein two observed variables of the observer areAndare used to estimate e respectivelykAnd dkBased on the error dynamics (equation (5)), an observer of the following form is designed
Wherein,represents the error ek+1Is estimated by the estimation of (a) a,represents the error ekIs estimated by the estimation of (a) a,representing equivalent perturbations, β1Indicating observer gain coefficient with respect to error, β2Representing the observer gain coefficient with respect to the equivalent disturbance,an estimation error representing a tracking error;
estimation error of equivalent disturbanceIs composed of
Estimation error of tracking error is
The expressions (7) and (8) are written as follows
Note the bookThe characteristic equation is
|λI-B|=0 (10)
Namely, it is
λ2+(β1-β2-1)λ-β1=0 (11)
Thus, the characteristic root is
For parameter β1And β2Configured so that all feature roots are within the unit circle, then matrix B is a Schur stable matrix, and the estimation error converges asymptotically, i.e.
Step 5. construct the power law of attraction with disturbance suppression measures
Wherein rho and epsilon are both adjustable parameters,denotes an attraction index, and 0 < rho < 1, epsilon > 0,
step 6, constructing a repetitive controller with equivalent disturbance compensation
Combining equation (5) and equation (12), design a repetitive controller with equivalent disturbance compensation
Note the book
Expressing a repetitive controller as
uk=uk-N+vk(14)
Will ukThe controller input signal as servo object can measure and obtain servo system output signal ykFollows the reference signal rkAnd (4) changing.
Further, an expression of four indexes, such as a steady state error band, an absolute attraction layer, a monotone decreasing area, the maximum number of steps required for the tracking error to enter the steady state error band for the first time and the like is given, and the expression is used for describing the tracking performance of the system and guiding the parameter setting of the controller, wherein the steady state error band, the absolute attraction layer, the monotone decreasing area and the maximum convergence number are defined as follows:
1) monotonous decreasing region deltaMDR: when e iskGreater than this boundary, ekThe same number is decreased, namely the following conditions are met:
2) absolute attraction layer ΔAAL: absolute value of system tracking errorkIf | is greater than this boundary, its | ekI, monotonically decreases, i.e. the condition is satisfied:
3) steady state error band ΔSSE: when the system error once converges into the boundary, the error is stabilized in the region, that is, the following condition is satisfied:
4) maximum number of convergence stepsThe tracking error passes through at mostStep into the steady state error band.
Equivalent disturbance compensation error satisfactionThe expression of each index is as follows
Monotonous decreasing region deltaMDR
ΔMDR=max{ΔMDR1,ΔMDR2} (18)
Wherein, DeltaMDR1And ΔMDR2Are all real and are determined by equation (19).
Absolute attraction layer ΔAAL
ΔAAL=max{ΔAAL1,ΔAAL2} (20)
Wherein, DeltaAAL1And ΔAAL2Are all real numbers and are determined by equation (21).
Steady state error band ΔSSE
ΔSSE=max{ΔSSE1,ΔSSE2} (22)
Wherein, DeltaSSE1And ΔSSE2Are all real, and are determined by equation (23);
in addition, given ΔSSEThen, the tracking error enters the maximum number of steps of the steady state error band
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than.
Further, forTwo cases, according to the given ΔMDR、ΔAAL、ΔSSEDetermining a corresponding calculation formula by the expression and the convergence step number expression;
the situation is as follows:
1) monotonous decreasing region deltaMDR
1.1) whenTime of flight
1.2) whenTime of flight
1.3) whenTime of flight
Wherein
2) Absolute attraction layer ΔAAL
2.1) whenTime of flight
2.2) whenTime of flight
2.3) whenTime of flight
Wherein
3) Steady state error band
3.1) whenOr ΔAAL≥δSSETime of flight
ΔSSE=ΔAAL(31)
3.2) whenTime of flight
Wherein deltaSSEIs an equationThe root of Zhengguo.
4) Number of convergence steps
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than;
the situation is as follows:
1) monotonous decreasing region deltaMDR
1.1) whenTime of flight
1.2) whenTime of flight
1.3) whenTime of flight
Wherein
2) Absolute attraction layer ΔAAL
2.1) whenTime of flight
2.2) whenTime of flight
2.3) whenTime of flight
Wherein
3) Steady state error band
3.1) whenOr ΔAAL≥δSSETime of flight
ΔSSE=ΔAAL(40)
3.2) whenTime of flight
Wherein deltaSSEIs an equationRoot of Zhengguo;
4) number of convergence steps
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than.
The technical conception of the invention is as follows: a power-suction repetitive control method using an equivalent disturbance compensation servo system is provided. According to a given reference signal and the constructed equivalent disturbance, an observer is introduced to estimate the equivalent disturbance, and interference suppression measures are embedded into a power law of attraction to form ideal error dynamics with interference suppression effect, so that a repetitive controller with equivalent disturbance compensation is designed, and rapid and high-precision tracking of the given reference signal is realized.
The invention has the following beneficial effects: the method has equivalent disturbance compensation, complete suppression of periodic disturbance, fast convergence performance and high tracking precision.
Drawings
Fig. 1 is a block diagram of an ac permanent magnet synchronous motor servo system.
FIG. 2 is a block diagram of an equivalent disturbance observer.
Fig. 3 is a block diagram of a power attraction repeat controller.
FIG. 4 is a graph of the time when a disturbance w occursk=5sin(2πfkTs) +0.15sgn (sin (2k pi/150)), the simulation results for the controller parameters with ∈ 0.1, ρ 0.3, and Δ 0.3, where Δ is plottedMDR,ΔAALAnd deltaSSE。
FIG. 5 is a graph of the time when a disturbance w occursk=-10sin(2πfkTs) +0.15sgn (sin (2k pi/150)), the simulation results for the controller parameters with ∈ 0.1, ρ 0.3, and Δ 0.3, where Δ is plottedMDR,ΔAALAnd deltaSSE。
FIG. 6 is a graph of the time when a disturbance w occursk=5sin(2πfkTs) +0.15sgn (sin (2k pi/150)), the simulation results for the controller parameters with ∈ 0.15, ρ 0.5, and Δ 0.3, where Δ is plottedMDR,ΔAALAnd deltaSSE。
FIG. 7 is a graph of the time when a disturbance w occursk=-10sin(2πfkTs) +0.15sgn (sin (2k pi/150)), the simulation results for the controller parameters with ∈ 0.15, ρ 0.5, and Δ 0.3, where Δ is plottedMDR,ΔAALAnd deltaSSE。
Fig. 8 to 11 are experimental results of the permanent magnet synchronous motor control device when the feedback controller parameter is ρ ═ 0.3 and ∈ ═ 0.1, where:
FIG. 8 is a reference position signal and an actual position signal under the action of a feedback controller based on the power law of attraction;
FIG. 9 is a control voltage signal under the action of a feedback controller based on the power law of attraction;
FIG. 10 is a graph of position error under the influence of a feedback controller based on the power law of attraction;
fig. 11 is a histogram of a position error distribution by a feedback controller based on the power law of attraction.
Fig. 12-15 show the feedback controller parameters ρ ═ 0.3, ε ═ 0.1, and the observer parameters β1=0.2,β20.5, the experimental result of the permanent magnet synchronous motor control device, wherein:
FIG. 12 is a reference position signal and an actual position signal under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 13 is a control voltage signal under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 14 is a position error under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
fig. 15 is a histogram of the position error distribution under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation.
Fig. 16 to 19 show experimental results of the permanent magnet synchronous motor control device when the repetitive controller parameter is ρ 0.3 and ∈ 0.1, where:
FIG. 16 is a reference position signal and an actual position signal under the influence of a repetitive controller based on the power law of attraction;
FIG. 17 is a control voltage signal under the action of a repetitive controller based on the power law of attraction;
FIG. 18 is a graph of position error under the influence of a repetitive controller based on the power law of attraction;
fig. 19 is a histogram of a position error distribution under the action of a repetitive controller based on the power law of attraction.
Fig. 20 to 23 show that the repetitive controller parameter ρ is 0.3, the parameter ∈ is 0.1, and the observer parameter β1=0.2,β20.5, the experimental result of the permanent magnet synchronous motor control device, wherein:
FIG. 20 is a graph of a reference position signal and an actual position signal under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 21 is a control voltage signal under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 22 is a position error under repetitive controller action based on the power law of attraction and equivalent disturbance compensation;
fig. 23 is a histogram of the position error distribution under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation.
Fig. 24 to 27 are experimental results of the permanent magnet synchronous motor control device when the feedback controller parameter is ρ is 0.5 and ∈ is 0.15, where:
FIG. 24 is a graph of a reference position signal and an actual position signal under the action of a feedback controller based on the power law of attraction;
FIG. 25 is a control voltage signal based on the power law of attraction under the influence of a feedback controller;
FIG. 26 is a graph showing a position error under the action of a feedback controller based on the power law of attraction;
fig. 27 is a histogram of a position error distribution by a feedback controller based on the power law of attraction.
Fig. 28 to 31 show that ρ is 0.5, and ∈ is 0.15 for the feedback controller parameters, and β for the observer parameters1=0.2,β2Experimental results of the permanent magnet synchronous motor control device when 0.5, which is equivalent toThe method comprises the following steps:
FIG. 28 is a graph of the reference position signal and the actual position signal under the influence of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 29 is a control voltage signal under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 30 is a graph of position error under the influence of a feedback controller based on the power law of attraction and equivalent disturbance compensation;
fig. 31 is a position error distribution histogram under the action of a feedback controller based on the power law of attraction and equivalent disturbance compensation.
Fig. 32 to 35 are experimental results of the permanent magnet synchronous motor control device when the repetitive controller parameter is ρ 0.5 and ∈ 0.15, where:
FIG. 32 is a graph of a reference position signal and an actual position signal under the influence of a repetitive controller based on the power law of attraction;
FIG. 33 is a control voltage signal under the influence of a repetitive controller based on the power law of attraction;
FIG. 34 is a graph of position error under the influence of a repetitive controller based on the power law of attraction;
fig. 35 is a histogram of the distribution of position errors under the action of a repetitive controller based on the power law of attraction.
Fig. 36 to 39 show that the repetitive controller parameter ρ is 0.5, the parameter ∈ is 0.15, and the observer parameter β1=0.2,β20.5, the experimental result of the permanent magnet synchronous motor control device, wherein:
FIG. 36 is a graph of a reference position signal and an actual position signal under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 37 is a control voltage signal under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation;
FIG. 38 is a position error under repetitive controller action based on the power law of attraction and equivalent disturbance compensation;
fig. 39 is a histogram of the position error distribution under the influence of a repetitive controller based on the power law of attraction and equivalent disturbance compensation.
Detailed Description
The following further describes embodiments of the present invention with reference to the accompanying drawings.
Referring to fig. 1-39, a power-suction repetitive control method using an equivalent disturbance compensation servo system. Wherein, fig. 1 is a block diagram of a servo system of an alternating current permanent magnet synchronous motor; FIG. 2 is a block diagram of an equivalent disturbance observer; fig. 3 is a block diagram of a power attraction repeat controller.
A power attraction repetitive control method adopting an equivalent disturbance compensation servo system comprises the following steps:
step 1, giving a periodic reference signal, and satisfying the requirement (1);
step 2, defining a tracking error, wherein a second-order difference model of a motor servo system is (2), and a tracking error of the system is (3);
step 3, constructing an equivalent disturbance (4), and expressing a system tracking error as (5) by using the equivalent disturbance (4);
step 4, designing an observer and estimating equivalent disturbance;
step 5, constructing a power attraction law (12) with disturbance suppression measures;
and 6, constructing a repetitive controller with equivalent disturbance compensation, combining the formula (5) and the formula (12), designing a repetitive controller (13) with equivalent disturbance compensation, and expressing the repetitive controller as (14).
The above repetitive controller design is explained as follows:
1) introduction of d into power attraction lawk+1Reflecting suppression measures for disturbing signals of a given periodic pattern, introducedAn estimate of the equivalent disturbance is reflected, thereby providing equivalent disturbance compensation.
2) In the formula (13), ek、yk、yk-1、yk-1-NAll can be obtained by measurement, uk-1、uk-1-NThe stored value for the control signal may be read from memory.
3) When the reference signal satisfies rk=rk-1The discrete repetitive controller is also suitable for the constant value regulation problem, and the equivalent disturbance is dk=wk-wk-1(ii) a Wherein r isk-1Reference signal representing the time instant k-1, wk-1Representing the interference signal at time k-1; the feedback controller with equivalent disturbance compensation is
4) The discrete time controller is designed for a second-order system, and the design result of a higher-order system can be given according to the same method.
Furthermore, an expression of four indexes, namely a steady-state error band, an absolute attraction layer, a monotone decreasing area and the maximum step number required for the tracking error to enter the steady-state error band for the first time is given, and the expression is used for describing the tracking performance of the system and guiding the parameter setting of the controller.
Further, forTwo cases, according to the given ΔMDR、ΔAAL、ΔSSEAnd determining a corresponding calculation formula by the expression and the convergence step number expression.
In this embodiment, for example, the permanent magnet synchronous motor device executes a repetitive tracking task in a fixed interval, and the position reference signal has a periodic symmetry characteristic. TMS320F2812DSP is used as a controller, a Korean LS AC servo motor APM-SB01AGN is used as a control object, and a permanent magnet synchronous motor servo system is formed by the ELMO AC servo driver and an upper mechanism to control the position of the motor. The servo system adopts three-loop control, the current loop and speed loop controller ELMO driver provide, and the position loop is provided by DSP development board.
Obtaining a mathematical model of the servo object by parameter estimation
yk+1-1.8949yk+0.8949yk-1=1.7908uk-0.5704uk-1+wk+1(43)
Wherein, yk,ukPosition output and control input, w, respectively, of the position servo systemkIs an interference signal.
The effectiveness of the repetitive controller given by the present invention will be illustrated in this example by numerical simulation and experimental results.
Numerical simulation this embodiment takes a sinusoidal signal as a system reference signal, and the corresponding repetitive controller expression can be written as
Given a position reference signal of rk=20(sin(2πfkTs-1/2 pi) +1) in degrees (deg), frequency f 1Hz, sampling period Ts0.001s, and the sampling period N is 1000. Selecting proper disturbance amount w during simulationkIt consists of periodic disturbances and non-periodic random disturbances.
Under the action of a repetitive controller (44), different controller parameters rho and epsilon are selected, and three boundary layers of a servo system are different. For purposes of illustrating the invention patent with respect to the monotonically decreasing region ΔMDRAbsolute attraction layer DeltaAALAnd steady state error band ΔSSETheoretical correctness of (1) inNumerical simulations were performed for the examples.
1) When the controller parameter ∈ is 0.1, ρ is 0.3, and Δ is 0.3, the calculation formula of the three boundary values indicates that
ΔSSE=ΔAAL=ΔMDR=0.7035
2) When the controller parameter ∈ is 0.15, ρ is 0.5, and Δ is 0.3, the calculation formula of the three boundary values indicates that
ΔSSE=ΔAAL=0.3823,ΔMDR=0.8884
The simulation results are shown in FIGS. 4-7, where FIGS. 4 and 6 are the disturbance values wk=5sin(2πfkTs) Simulation result of +0.15sgn (2k pi/150),FIG. 5,7 is the disturbance quantity wk=-10sin(2πfkTs) +0.15sgn (2k pi/150).
The numerical results verify the monotonous reduction area delta of the tracking error of the system under the action of the repetitive controller given by the patent under the condition of a given system model, a reference signal and an interference signalMDRAbsolute attraction layer DeltaAALAnd steady state error band ΔSSEThe accuracy of (2).
The block diagram of the permanent magnet synchronous motor control system for the experimental verification experiment is shown in figure 1. And verifying the tracking performance of discrete repetitive control based on the power attraction law by setting different controller parameters. Given position signal rk=A(sin(2πfkTs) +1), where the amplitude a is 135deg, the sampling period Ts5ms, frequency f 1Hz, toExperimental verification was performed for the examples.
The feedback controller adopted has the following form
The feedback controller based on disturbance compensation is adopted and has the following form
The adopted repetitive controller has the following form
The repetitive controller based on disturbance compensation is adopted and has the following form
1) The controller (45) is adopted, the controller parameters are rho is 0.3, epsilon is 0.1, and the experimental results are shown in figures 8-11, wherein delta is shown in the figuresSSE=0.15deg。
2) Adopting a controller (46), wherein the parameters of the controller are rho equal to 0.3, epsilon equal to 0.1, and the parameters of the equivalent disturbance observer are β1=0.2,β2The results are shown in FIGS. 12-15, where Δ is 0.5SSE=0.1deg。
3) The controller (47) is adopted, the controller parameters are rho is 0.3, epsilon is 0.1, and the experimental results are shown in figures 16-19, wherein delta is shown in the figuresSSE=0.1deg。
4) Adopting a controller (48), wherein the parameters of the controller are rho equal to 0.3, epsilon equal to 0.1, and the parameters of the equivalent disturbance observer are β1=0.2,β2The results are shown in FIGS. 20-23, where Δ is 0.5SSE=0.08deg。
5) By means of a controller(45) The controller parameters ρ is 0.5 and ∈ is 0.15, and the experimental results are shown in fig. 24 to 27, where Δ is shown in the graphsSSE=0.15deg。
6) Adopting a controller (46), wherein the parameters of the controller are rho equal to 0.5, epsilon equal to 0.15, and the parameters of the equivalent disturbance observer are β1=0.2,β2The results are shown in FIGS. 28-31, where Δ is 0.5SSE=0.11deg。
7) The controller (47) is adopted, the controller parameters are rho is 0.5, epsilon is 0.15, and the experimental results are shown in figures 32-35, wherein delta is shown in the figuresSSE=0.1deg。
8) Adopting a controller (48), wherein the parameters of the controller are rho 0.5, epsilon 0.15, and the parameters of the equivalent disturbance observer are β1=0.2,β2The results are shown in FIGS. 36-39, where Δ is 0.5SSE=0.07deg。
The experiment result shows that the equivalent disturbance is introduced and is estimated by the equivalent disturbance observer, the compensation for the unmodeled characteristic and the external unknown disturbance of the system is provided, and the influence of the unknown disturbance on the tracking performance can be effectively inhibited; and the periodic disturbance is completely inhibited by adopting repeated control, so that the control performance of the system is further improved.
Claims (3)
1. A power attraction repetitive control method adopting an equivalent disturbance compensation servo system, wherein a controlled object is a periodic servo system, and the method is characterized by comprising the following steps:
step 1, giving periodic reference signals to satisfy
rk=rk-N(1)
Where N is the period of the reference signal, rkAnd rk-NReference signals respectively representing time k and time k-N;
step 2, defining tracking error
In the formula
A1(q-1)=a1+a2q-1+…+anq-n+1=q(A(q-1)-1)
A(q-1)=1+a1q-1+…+anq-n
B(q-1)=b0q-1+…+bmq-m
Satisfy the requirement of
A(q-1)yk=q-dB(q-1)uk+wk(3)
Wherein e isk+1Represents the tracking error at time k +1, rk+1Reference signal, y, representing the time instant k +1k+1、yk、yk+1-NAnd yk-NRepresenting the output signals at times k +1, k +1-N and k-N, respectively, ukAnd uk-NRepresenting the input signal at times k and k-N, wk+1And wk+1-NRepresenting the interfering signal at times k and k-N, respectively, d represents the delay, A (q)-1) And B (q)-1) Is q-1Polynomial of (a), q-1Denotes a one-step delay operator, n denotes A (q)-1) M represents B (q)-1) Order of (a)1,…,an,b0,…,bmIs a system parameter and b0Not equal to 0, n is more than or equal to m, d is an integer and is more than or equal to 1;
step 3. constructing equivalent disturbance
dk=wk-wk-N(4)
Where N is the period of the reference signal, dkRepresenting the equivalent disturbance signal at time k, wkAnd wk-NRespectively representing interference signals at the k moment and the k-N moment;
expressing the tracking error as
ek+1=rk+1-yk+1-N+A1(q-1)(yk-yk-N)-q-d+1B(q-1)(uk-uk-N)-dk+1(5)
Wherein d isk+1Representing the equivalent disturbance at the moment k + 1;
step 4, designing an observer and estimating equivalent disturbance
Design observer equivalent disturbance dk+1Observing and compensating equivalent disturbance by the observed value; two observed variables of the observer areAndare used to estimate e respectivelykAnd dk(ii) a From the error dynamics (equation (5)), an observer of the following form is designed
Wherein,represents the error ek+1Is estimated by the estimation of (a) a,represents the error ekIs estimated by the estimation of (a) a,representing equivalent perturbations, β1Indicating observer gain coefficient with respect to error, β2Representing observer gain coefficients with respect to equivalent disturbances;an estimation error representing a tracking error;
estimation error of equivalent disturbanceIs composed of
Estimation error of tracking error is
The expressions (7) and (8) are written as follows
Note the bookThe characteristic equation is
|λI-B|=0 (10)
Namely, it is
λ2+(β1-β2-1)λ-β1=0 (11)
Thus, the characteristic root isPair β1And β2Is configured such that all feature roots are within the unit circle;
step 5. construct the power law of attraction with disturbance suppression measures
Wherein rho and epsilon are both adjustable parameters,denotes an attraction index, and 0 < rho < 1, epsilon > 0,
step 6, constructing a repetitive controller with equivalent disturbance compensation
Combining equation (5) and equation (12) to obtain a repetitive controller with equivalent disturbance compensation
Note the book
Expressing a repetitive controller as
uk=uk-N+vk(14)
Will ukAs input signal of controller of servo object, measuring to obtain output signal y of servo systemkFollows the reference signal rkAnd (4) changing.
2. The method as claimed in claim 1, wherein the expressions of four indexes, namely a steady state error band, an absolute attraction layer, a monotone decreasing region and the maximum number of steps required for a tracking error to enter the steady state error band for the first time, are given for describing the tracking performance of the system and guiding the parameter setting of the controller, wherein the steady state error band, the absolute attraction layer, the monotone decreasing region and the convergence step are defined as follows:
1) monotonous decreasing region deltaMDR: when e iskGreater than this boundary, ekThe same number is decreased, namely the following conditions are met:
2) absolute attraction layer ΔAAL: absolute value of system tracking errorkIf | is greater than this boundary, its | ekI, monotonically decreases, i.e. the condition is satisfied:
3) steady state error band ΔSSE: when the system error once converges into the boundary, the error is stabilized in the region, that is, the following condition is satisfied:
4) maximum number of convergence stepsThe tracking error passes through at mostEntering a steady state error band;
equivalent disturbance compensation error satisfactionThe expression of each index is as follows
Monotonous decreasing region deltaMDR
ΔMDR=max{ΔMDR1,ΔMDR2} (18)
Wherein, DeltaMDR1And ΔMDR2Are all real and are determined by equation (19);
absolute attraction layer ΔAAL
ΔAAL=max{ΔAAL1,ΔAAL2} (20)
Wherein, DeltaAAL1And ΔAAL2Are all real numbers and are determined by equation (21);
steady state error band ΔSSE
ΔSSE=max{ΔSSE1,ΔSSE2} (22)
Wherein, DeltaSSE1And ΔSSE2Are all real, and are determined by equation (23);
in addition, given ΔSSEThen, the tracking error enters the maximum number of steps of the steady state error band
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than.
3. The method of claim 2, wherein for the power-law pull-in repetitive control using an equivalent perturbation compensation servo systemAndtwo cases, according to the monotone decreasing zone Δ givenMDRAbsolute attraction layer DeltaAALSteady state error band deltaSSEAnd maximum number of convergence stepsThe expression determines a corresponding calculation formula;
the situation is as follows:
1) monotonous decreasing region deltaMDR
1.1) whenTime of flight
1.2) whenTime of flight
1.3) whenTime of flight
Wherein
2) Absolute attraction layer ΔAAL
2.1) whenTime of flight
2.2) whenTime of flight
2.3) whenTime of flight
Wherein
3) Steady state error band
3.1) whenOr ΔAAL≥δSSETime of flight
ΔSSE=ΔAAL(31)
3.2) whenTime of flight
Wherein deltaSSEIs an equationRoot of Zhengguo;
4) number of convergence steps
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than;
the situation is as follows:
1) monotonous decreasing region deltaMDR
1.1) whenTime of flight
1.2) whenTime of flight
1.3) whenTime of flight
Wherein
2) Absolute attraction layer ΔAAL
2.1) whenTime of flight
2.2) whenTime of flight
2.3) whenTime of flight
Wherein
3) Steady state error band
3.1) whenOr ΔAAL≥δSSETime of flight
ΔSSE=ΔAAL(40)
3.2) whenTime of flight
Wherein deltaSSEIs an equationRoot of Zhengguo;
4) number of convergence steps
Wherein e is0In order to be the initial value of the tracking error,represents the smallest integer no less than.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910307598.1A CN110134010B (en) | 2019-04-17 | 2019-04-17 | Power attraction repetitive control method adopting equivalent disturbance compensation servo system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910307598.1A CN110134010B (en) | 2019-04-17 | 2019-04-17 | Power attraction repetitive control method adopting equivalent disturbance compensation servo system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110134010A true CN110134010A (en) | 2019-08-16 |
CN110134010B CN110134010B (en) | 2022-05-03 |
Family
ID=67570057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910307598.1A Active CN110134010B (en) | 2019-04-17 | 2019-04-17 | Power attraction repetitive control method adopting equivalent disturbance compensation servo system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110134010B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456645A (en) * | 2019-08-18 | 2019-11-15 | 台州学院 | A kind of discrete repetitive control for inverter |
CN110579964A (en) * | 2019-09-06 | 2019-12-17 | 浙江大学城市学院 | Robust adaptive repetitive controller design method |
CN110658719A (en) * | 2019-09-27 | 2020-01-07 | 浙江工业大学 | Servo system non-switching attraction repetitive control method adopting equivalent disturbance compensation |
CN110673477A (en) * | 2019-09-27 | 2020-01-10 | 浙江工业大学 | Design method of discrete time power attraction controller adopting interference difference suppression strategy |
CN110716430A (en) * | 2019-09-27 | 2020-01-21 | 浙江工业大学 | Servo system rapid suction repetitive control method adopting equivalent disturbance compensation |
CN110995102A (en) * | 2019-12-31 | 2020-04-10 | 南京工业大学 | Direct torque control method and system for permanent magnet synchronous motor |
CN111752149A (en) * | 2020-06-08 | 2020-10-09 | 浙江工业大学 | Design method of digital repetitive controller for designated time attraction of servo motor |
CN111835254A (en) * | 2020-06-04 | 2020-10-27 | 浙江工业大学 | Permanent magnet synchronous motor speed control method based on finite time attraction law |
CN111835250A (en) * | 2020-06-05 | 2020-10-27 | 浙江工业大学 | 1/2 power finite value attraction repetitive control method for servo motor driving system |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11285283A (en) * | 1998-03-27 | 1999-10-15 | Mitsubishi Heavy Ind Ltd | Disturbance compensation control device |
US20070067051A1 (en) * | 2005-09-22 | 2007-03-22 | Ramos Jesus L | Repetitive controller to compensate for odd harmonics |
CN105446140A (en) * | 2015-12-21 | 2016-03-30 | 内蒙古科技大学 | Finite time convergence-based composite control method of stabilization platform servo system |
WO2017171816A1 (en) * | 2016-03-31 | 2017-10-05 | Etagen, Inc. | Control of piston trajectory in a free-piston combustion engine |
CN107544244A (en) * | 2017-08-25 | 2018-01-05 | 浙江工业大学 | Discrete repetitive control method for motor servo system based on elliptic attraction law and equivalent disturbance expansion state compensation |
CN107544245A (en) * | 2017-08-25 | 2018-01-05 | 浙江工业大学 | Discrete repetitive control method for motor servo system adopting disturbance expansion compensation first-order inertia attraction law |
CN107797448A (en) * | 2017-08-25 | 2018-03-13 | 浙江工业大学 | Motor position discrete repetitive control method adopting disturbance expansion compensation |
CN108646574A (en) * | 2018-07-26 | 2018-10-12 | 台州学院 | A kind of discrete repetitive controller and its control method attracting rule based on power |
WO2018198694A1 (en) * | 2017-04-27 | 2018-11-01 | 日本電産株式会社 | Control device for outputting voltage command value |
CN108983610A (en) * | 2018-08-03 | 2018-12-11 | 华侨大学 | A kind of robust adaptive disturbance rejection control method |
CN108983615A (en) * | 2018-08-17 | 2018-12-11 | 台州学院 | Attract the discrete binary cycle repetitive controller of rule based on asinh |
CN109450307A (en) * | 2018-10-10 | 2019-03-08 | 浙江工业大学 | Attract rule based on class normal distribution and using the discrete repetitive control for motor servo system of disturbance expansion state compensation |
-
2019
- 2019-04-17 CN CN201910307598.1A patent/CN110134010B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11285283A (en) * | 1998-03-27 | 1999-10-15 | Mitsubishi Heavy Ind Ltd | Disturbance compensation control device |
US20070067051A1 (en) * | 2005-09-22 | 2007-03-22 | Ramos Jesus L | Repetitive controller to compensate for odd harmonics |
CN105446140A (en) * | 2015-12-21 | 2016-03-30 | 内蒙古科技大学 | Finite time convergence-based composite control method of stabilization platform servo system |
WO2017171816A1 (en) * | 2016-03-31 | 2017-10-05 | Etagen, Inc. | Control of piston trajectory in a free-piston combustion engine |
WO2018198694A1 (en) * | 2017-04-27 | 2018-11-01 | 日本電産株式会社 | Control device for outputting voltage command value |
CN107544244A (en) * | 2017-08-25 | 2018-01-05 | 浙江工业大学 | Discrete repetitive control method for motor servo system based on elliptic attraction law and equivalent disturbance expansion state compensation |
CN107544245A (en) * | 2017-08-25 | 2018-01-05 | 浙江工业大学 | Discrete repetitive control method for motor servo system adopting disturbance expansion compensation first-order inertia attraction law |
CN107797448A (en) * | 2017-08-25 | 2018-03-13 | 浙江工业大学 | Motor position discrete repetitive control method adopting disturbance expansion compensation |
CN108646574A (en) * | 2018-07-26 | 2018-10-12 | 台州学院 | A kind of discrete repetitive controller and its control method attracting rule based on power |
CN108983610A (en) * | 2018-08-03 | 2018-12-11 | 华侨大学 | A kind of robust adaptive disturbance rejection control method |
CN108983615A (en) * | 2018-08-17 | 2018-12-11 | 台州学院 | Attract the discrete binary cycle repetitive controller of rule based on asinh |
CN109450307A (en) * | 2018-10-10 | 2019-03-08 | 浙江工业大学 | Attract rule based on class normal distribution and using the discrete repetitive control for motor servo system of disturbance expansion state compensation |
Non-Patent Citations (7)
Title |
---|
AHM SAYEM,等: "Discrete extended state observer based repetitive control system for improved disturbance rejection performance", 《2016 IEEE 11TH CONFERENCE ON INDUSTRIAL ELECTRONICS AND APPLICATIONS (ICIEA)》 * |
HITOSHI HAGA,等: "Current control system based on repetitive control and disturbance observer for single-phase five-level inverter", 《 IECON 2015 - 41ST ANNUAL CONFERENCE OF THE IEEE INDUSTRIAL ELECTRONICS SOCIETY》 * |
司刚: "基于混合控制器的飞轮储能系统鲁棒控制方法", 《电器与能效管理技术》 * |
张博等: "基于双幂次滑模趋近律方法的PMLSM精密位置控制", 《微电机》 * |
李捷等: "自抗扰控制在快刀伺服系统控制中的应用", 《机床与液压》 * |
李鹤等: "一类不确定非线性系统的重复学习控制", 《自动化学报》 * |
邹权等: "基于扰动观测和补偿的PMSM伺服系统位置跟踪控制", 《电机与控制学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456645A (en) * | 2019-08-18 | 2019-11-15 | 台州学院 | A kind of discrete repetitive control for inverter |
CN110456645B (en) * | 2019-08-18 | 2022-02-15 | 台州学院 | Discrete repetitive control method for inverter |
CN110579964A (en) * | 2019-09-06 | 2019-12-17 | 浙江大学城市学院 | Robust adaptive repetitive controller design method |
CN110579964B (en) * | 2019-09-06 | 2022-03-29 | 浙江大学城市学院 | Robust adaptive repetitive controller design method |
CN110658719A (en) * | 2019-09-27 | 2020-01-07 | 浙江工业大学 | Servo system non-switching attraction repetitive control method adopting equivalent disturbance compensation |
CN110673477A (en) * | 2019-09-27 | 2020-01-10 | 浙江工业大学 | Design method of discrete time power attraction controller adopting interference difference suppression strategy |
CN110716430A (en) * | 2019-09-27 | 2020-01-21 | 浙江工业大学 | Servo system rapid suction repetitive control method adopting equivalent disturbance compensation |
CN110716430B (en) * | 2019-09-27 | 2022-05-03 | 浙江工业大学 | Servo system rapid suction repetitive control method adopting equivalent disturbance compensation |
CN110995102A (en) * | 2019-12-31 | 2020-04-10 | 南京工业大学 | Direct torque control method and system for permanent magnet synchronous motor |
CN111835254A (en) * | 2020-06-04 | 2020-10-27 | 浙江工业大学 | Permanent magnet synchronous motor speed control method based on finite time attraction law |
CN111835250A (en) * | 2020-06-05 | 2020-10-27 | 浙江工业大学 | 1/2 power finite value attraction repetitive control method for servo motor driving system |
CN111752149A (en) * | 2020-06-08 | 2020-10-09 | 浙江工业大学 | Design method of digital repetitive controller for designated time attraction of servo motor |
Also Published As
Publication number | Publication date |
---|---|
CN110134010B (en) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110134010B (en) | Power attraction repetitive control method adopting equivalent disturbance compensation servo system | |
CN110716430B (en) | Servo system rapid suction repetitive control method adopting equivalent disturbance compensation | |
CN110134014B (en) | Equivalent disturbance compensation method for periodic servo system power attraction repetitive control | |
CN107797448B (en) | Motor position discrete repetitive control method adopting disturbance expansion compensation | |
CN109450307B (en) | Discrete repetitive control method for motor servo system based on similar normal distribution attraction law and adopting disturbance expansion state compensation | |
CN110032073B (en) | 1/2 power attraction repetitive control method with equivalent disturbance compensation | |
CN107544245B (en) | Discrete repetitive control method for motor servo system adopting disturbance expansion compensation first-order inertia attraction law | |
CN107544244B (en) | Discrete repetitive control method for motor servo system based on elliptic attraction law and equivalent disturbance expansion state compensation | |
CN110658719A (en) | Servo system non-switching attraction repetitive control method adopting equivalent disturbance compensation | |
CN105867110A (en) | Discrete repetition control method for motor servo system | |
CN109709792A (en) | Aero-engine stable state circuit pi controller and its design method and device | |
CN109743002B (en) | Servo system controller, feedforward control signal determination method and inertia identification method | |
CN111752149B (en) | Design method of digital repetitive controller for designated time attraction of servo motor | |
CN112068423B (en) | Design method of finite power attraction repetitive controller | |
CN111835250B (en) | 1/2 power finite value attraction repetitive control method for servo motor driving system | |
Kilicarslan et al. | Modeling and hysteresis compensation in a thin SMA wire using ANFIS methods | |
CN112068424B (en) | Discrete repeated control method of elliptic approach law adopting disturbance compensation | |
CN112117945B (en) | Discrete repeated control method of parabolic approach law adopting disturbance compensation | |
Hammouche et al. | Robust and optimal output-feedback control for interval state-space model: application to a two-degrees-of-freedom piezoelectric tube actuator | |
Nassef et al. | Simulation-based parameter identification of a reduced model using neural networks | |
Nyulászi et al. | Comparison of experimental identification methods using measured data from a turbojet engine | |
Coman et al. | Combining the MIT and Lyapunov stability adaptive methods for second order systems | |
CN112099346B (en) | Design method of specified time repetitive controller based on arctangent attraction law | |
CN111340647A (en) | Power distribution network data adjustment method and system | |
Mozaryn et al. | Experimental Verification of Discrete Linear-Quadratic-Gaussian Control System of Electro-Hydraulic Servodrive |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |