CN109188908B - Digital controller design method based on exponential type non-switching attraction law - Google Patents

Digital controller design method based on exponential type non-switching attraction law Download PDF

Info

Publication number
CN109188908B
CN109188908B CN201811112071.5A CN201811112071A CN109188908B CN 109188908 B CN109188908 B CN 109188908B CN 201811112071 A CN201811112071 A CN 201811112071A CN 109188908 B CN109188908 B CN 109188908B
Authority
CN
China
Prior art keywords
interference
error
signal
reference signal
controller
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.)
Active
Application number
CN201811112071.5A
Other languages
Chinese (zh)
Other versions
CN109188908A (en
Inventor
孙明轩
胡志云
张钰
李鹤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201811112071.5A priority Critical patent/CN109188908B/en
Publication of CN109188908A publication Critical patent/CN109188908A/en
Application granted granted Critical
Publication of CN109188908B publication Critical patent/CN109188908B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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 digital controller non-switching attraction law design method adopting an interference differential suppression strategy is characterized in that a given module generates a reference signal; constructing a corresponding interference difference compensation feedback link according to the specific form of a given reference signal, wherein the output signal of the interference difference compensation feedback link is used for correcting a digital controller; and constructing an ideal error dynamic state based on a non-switching attraction law, designing a digital controller according to the ideal error dynamic state, and taking a signal obtained by calculation of the current controller as the input of a servo object. The specific controller parameter setting can be carried out according to the indexes of the convergence performance of the representation system, and a steady-state error band, an absolute attraction layer and a monotone decreasing area of the representation tracking error convergence process are provided. And then the anti-interference capability and the tracking precision of the system are improved, so that the servo system can follow the change of the reference signal. The invention provides a method for designing a non-switching suction law of a digital controller, which can eliminate the phenomenon of system buffeting and has good control precision.

Description

Digital controller design method based on exponential type non-switching attraction law
Technical Field
The invention relates to an exponential type non-switching attraction law design method adopting an interference difference suppression strategy digital controller, which is suitable for a position servo system and other industrial application occasions.
Background
The approach law method is an effective tool for designing a sliding mode controller of a servo system, and due to the adoption of the approach law, the dynamic process of a closed-loop system is represented by an approach process and a sliding mode, and the stability and the convergence of the closed-loop system are determined by the specific approach law and a switching function form. The actual controller design needs to consider the influence of various disturbances, the disturbance suppression measures are 'embedded' in the original approach law, and the modified approach law forms ideal switching dynamics. In this way, a controller dynamically designed according to ideal switching can effectively suppress disturbances.
The attraction law method directly adopts the tracking error signal without defining a switching function, and the design of the controller is more direct and simpler. The attraction law reflects the expected dynamic characteristics of the system error when disturbance is not considered; in the presence of disturbances, a controller directly based on the attraction law cannot be implemented. The interference suppression measures can be 'embedded' into the attraction law, and ideal error dynamics with the disturbance suppression effect are constructed. The digital controller is designed according to the constructed ideal error dynamic equation, and the dynamic process of the closed-loop system is determined by the ideal error dynamic equation and has the expected tracking performance represented by the ideal error dynamic equation.
The attraction law method is different from an approach law method of discrete sliding mode control. The main differences between the two are as follows: replacing a switching function by the tracking error and replacing a switching surface by the original point by the attraction law method; the approach law method requires a finite time to reach the switching surface, while the attraction law method requires a finite time to reach the origin; the closed-loop system designed by the attraction law method still has robust performance related to parameter drift and external interference, only the sliding mode control focuses on invariance of sliding mode motion, and the attraction law method pursues invariance of system steady state.
When the discrete controller is designed by an attraction law method, indexes describing transient and steady-state behaviors of the tracking error can be dynamically given by an ideal error, and the indexes specifically comprise the following four indexes: a steady state error band, an absolute attraction layer, a monotonically decreasing area, and a maximum number of steps required for the tracking error to first enter the steady state error band. In fact, the specific values of the four indexes depend on the controller parameters, the controller parameters are different, and the values of the four indexes are also different. Once the ideal error dynamic form is given, specific expressions of four indexes can be given in advance for parameter setting of the controller. In the currently published attraction law method, four indexes all depend on the boundary of the equivalent interference signal. The boundary that effectively inhibits interference and reduces equivalent interference signals is a difficult problem to be solved urgently by an attraction law method.
Disclosure of Invention
In order to overcome the defects that the influence of interference signals on the performance of a servo system cannot be inhibited and the tracking control precision is low in the existing design method suitable for the digital controller of the position servo system, the invention provides an attraction law method suitable for designing the digital controller of the position servo system. In order to inhibit the influence of interference signals on the performance of a servo system and improve the tracking control precision, an interference differential compensation technology is adopted and embedded into an exponential attraction law so as to construct an ideal error dynamic state with disturbance inhibition capability. The digital controller is dynamically designed according to the ideal error, so that the closed-loop system has the characteristic of being dynamically drawn by the ideal error, and the anti-interference capability and the tracking performance of the position servo system are improved. The controller is designed by adopting a non-switching suction law, so that the buffeting phenomenon can be eliminated. The invention specifically provides a specific expression of four indexes, namely a steady-state error band, an absolute attraction layer, a monotone decreasing area, the maximum step number required for a tracking error to enter the steady-state error band for the first time, and the like, 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:
an exponential attraction law design for a digital controller with interference differential compensation, comprising the following steps:
1) given a reference signal rkFor a polynomial of the time variable k, M representing the highest power of the polynomial, three specific reference signals are as follows:
square wave signal, M ═ 0:
Figure BDA0001809381180000021
triangular wave signal, M ═ 1:
Figure BDA0001809381180000022
s-curve, M ═ 3:
Figure BDA0001809381180000031
wherein, A is amplitude, and N is sampling times of the reference signal in one period;
2) structure ideal error dynamics
For exponential type no-switching attraction law
ek+1=(1-ρ(ek))ek (4)
Figure BDA0001809381180000032
Wherein e isk=rk-ykSystematic tracking error at time k, ykThe actual output signal of the system at the moment k;
embedding interference suppression measures into the attraction law to construct an ideal error dynamic
Figure BDA0001809381180000033
Wherein, delta is more than 0 and less than 1, alpha is more than 0, beta is more than 1, dk+1Is the equivalent interference at the time k +1,
Figure BDA0001809381180000036
compensating the signal for equivalent interference; due to the parameter ρ (e)k) At ekAway from the origin, approaches δ at ekApproaches to the origin and approaches to 1, so that the system error starts from the initial position under the action of an ideal error dynamic state (6), the convergence speed is gradually increased, and delta, alpha and beta are parameters for adjusting the suction speed;
3) interference differential compensation strategy
Take the equivalent interference compensation effect as
Figure BDA0001809381180000035
Defining an equivalent interference for a particular reference signal in the form of
Equivalent interference signal when M is 0
dk+1=wk+1 (7)
Equivalent interference signal when M is 1
dk+1=wk+1-wk (8)
Equivalent interference signal when M is 3
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (9)
Wherein, wk+1Interference at time k + 1;
defining the step number of the interference difference as L, wherein L represents that the interference compensation error comprises L interferences at successive moments, and the two-step interference difference d is shown as formula (7)k+1-dk=wk+1-wkComprising wk+1And wkTwo-time interference, in order to effectively suppress interference, when selecting equivalent interference, the following condition should be satisfied
Figure BDA0001809381180000046
Wherein the content of the first and second substances,
Figure BDA0001809381180000041
is the smallest integer not less than;
4) controller design
According to the ideal error dynamics (6) and the equivalent interference signal dk+1The following controller expression is given:
for a square wave reference signal equation (1),
Figure BDA0001809381180000042
wherein, F (q)-1)=B(q-1)-b0
For the triangular wave reference signal equation (2),
Figure BDA0001809381180000043
wherein the content of the first and second substances,
Figure BDA0001809381180000044
for the S-curve reference signal equation (3),
Figure BDA0001809381180000045
in formulae (11) to (13), A (q)-1)、B(q-1) For serving systems
A(q-1)yk=q-1B(q-1)uk+wk (14)
With respect to q-1The parameter polynomial of (2):
A(q-1)=1+a1q-1+a2q-2+……+anq-n
B(q-1)=b0+b1q-1+b2q-2+……+bmq-m
wherein u iskAnd ykInput and output signals, q, respectively, at time k of the servo system-1For one-step delay operator, m and n are respectively A (q)-1)、B(q-1) Order of (b)0≠0,1≤m≤n。
Further, the method comprises the following steps:
5) analysis of properties
And giving specific expressions of a steady-state error band, an absolute attraction layer, a monotone subtraction area and at most four indexes of steps required by the tracking error entering the steady-state error band for the first time, and being 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 and the monotone subtraction area are defined as follows:
steady state error band ΔSSE
Figure BDA0001809381180000051
Absolute attraction layer ΔAAL
Figure BDA0001809381180000052
Monotonous decreasing region deltaMDR
Figure BDA0001809381180000053
Equivalent interference compensation error satisfaction
Figure BDA0001809381180000054
The expression of each index is as follows:
steady state error band (Δ)SSE)
ΔSSE=Δ (15)
Absolute attraction layer (. DELTA.)AAL)
ΔAAL=max{ΔAAL1AAL2} (16)
In the formula,. DELTA.AAL1,ΔAAL2Is a real number, determined by
Figure BDA0001809381180000055
Monotonous decreasing area (delta)MDR)
ΔMDR=max{ΔMDR1MDR2} (18)
In the formula,. DELTA.MDR1,ΔMDR2Is real, and satisfies
Figure BDA0001809381180000056
Number of convergence steps
Figure BDA0001809381180000057
Figure BDA0001809381180000058
Wherein e is0Is the initial value of the tracking error.
The technical conception of the invention is as follows: an exponential type no-switching attraction law method for the design of a tracking controller of a position servo system is provided. Equivalent interference is defined according to a given reference signal, and interference suppression measures are embedded into the non-switching attraction law to form ideal error dynamics with interference suppression effect. The digital controller is dynamically designed according to the ideal error to realize accurate tracking of a given reference signal.
The control effect of the invention is mainly shown in that: according to the given reference signal, a corresponding interference difference compensation measure is adopted, and the tracking accuracy is improved by suppressing interference. Meanwhile, a non-switching discrete time suction law is adopted to eliminate system buffeting.
Drawings
FIG. 1 is a block diagram of a controller system.
In FIG. 2 to FIG. 3, r isk=10sin(2πfkTs)deg,f=2Hz,TsWhen δ is 0.01, δ is 0.1, α is 5, β is 10, and Δ is 0.6, the numerical simulation result of controller equation (22) is used:
FIG. 2 shows a tracking error signal ek
Fig. 3 shows interference compensation errors
Figure BDA0001809381180000061
FIGS. 4-7 show a reference signal rkAs in formula (1), A ═ 5, TsWhen δ is 0.1, δ is 0.5, α is 10, and β is 10, the numerical simulation result of controller equation (22) is used:
FIG. 4 shows a reference signal rkAnd the output signal yk
FIG. 5 shows a tracking error signal ek
FIG. 6 shows interference compensation errors
Figure BDA0001809381180000062
FIG. 7 shows a control signal uk
FIGS. 8-11 show a reference signal rkWhen a is 20, δ is 0.5, α is 10, and β is 10, the numerical simulation result of the controller equation (22) is used as equation (2):
FIG. 8 shows a reference signal rkAnd the output signal yk
FIG. 9 shows a tracking error signal ek
FIG. 10 shows interference compensation errors
Figure BDA0001809381180000063
FIG. 11 shows a control signal uk
FIGS. 12-15 show a reference signal rkAs in formula (2), A ═ 20, TsWhen δ is 0.1, δ is 0.5, α is 10, and β is 10, the numerical simulation result of controller equation (23) is used:
FIG. 12 shows a reference signal rkAnd the output signal yk
FIG. 13 shows a tracking error signal ek
FIG. 14 shows the interference compensation error dk-dk-1
FIG. 15 shows a control signal uk
FIGS. 16-19 show a reference signal rkAs in formula (3), A ═ 5, TsWhen δ is 0.1, δ is 0.5, α is 10, and β is 10, the numerical simulation result of controller equation (22) is used:
FIG. 16 shows a reference signal rkAnd the output signal yk
FIG. 17 shows a tracking error signal ek
FIG. 18 shows interference compensation errors
Figure BDA0001809381180000071
FIG. 19 shows a control signal uk
FIGS. 20-23 show a reference signal rkWhen a is 5, δ is 0.5, α is 10, β is 10, and δ is 0.5, the numerical simulation result of the controller formula (23) is adopted as formula (3):
FIG. 20 shows a reference signal rkAnd the output signal yk
FIG. 21 shows a tracking error signal ek
FIG. 22 shows interference compensation errors
Figure BDA0001809381180000072
FIG. 23 shows a control signal uk
FIGS. 24-27 show a reference signal rkAs in formula (3), A ═ 5, TsWhen α is 0.1, β is 10, and δ is 0.5, the numerical simulation result of controller equation (24) is used:
FIG. 24 shows a reference signal rkAnd the output signal yk
FIG. 25 shows a tracking error signal ek
FIG. 26 shows interference compensation errors
Figure BDA0001809381180000073
FIG. 27 shows a control signal uk
FIGS. 28-29 show reference signals rkAs in formula (1), A ═ 15deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller equation (22) is used:
FIG. 28 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000074
FIG. 29 shows a tracking error signal ekA histogram of (a).
FIGS. 30-31 show a reference signal rkSuch as formula (2), A ═ 180deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller equation (22) is used:
FIG. 30 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000075
FIG. 31 shows a tracking error signal ekA histogram of (a).
FIGS. 32-33 show a reference signal rkSuch as formula (2), A ═ 180deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller equation (23) is used:
FIG. 32 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000076
FIG. 33 shows a tracking error signal ekA histogram of (a).
FIGS. 34-35 show a reference signal rkSuch as formula (3), A ═ 180deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller equation (22) is used:
FIG. 34 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000081
FIG. 35 shows a tracking error signal ekA histogram of (a).
FIGS. 36-37 show a reference signal rkSuch as formula (3), A ═ 180deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller equation (23) is used:
FIG. 36 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000082
FIG. 37 shows a tracking error signal ekA histogram of (a).
FIGS. 38-39 show the reference signal rkSuch as formula (3), A ═ 180deg, TsWhen α is 0.2ms, β is 10, and δ is 0.2, the experimental result of the controller formula (24) is used:
FIG. 38 shows the reference signal r in turnkAnd the output signal ykControl signal ukTracking error signal ekAnd interference compensation error
Figure BDA0001809381180000083
FIG. 39 shows a tracking error signal ekA histogram of (a).
Detailed Description
The following further describes embodiments of the present invention with reference to the accompanying drawings.
FIG. 1 is a block diagram of a servo system. Referring to fig. 2 to 39, an exponential type non-switching attraction law method for a discrete servo system with interference differential compensation and effective interference suppression is designed by a controller, and comprises the following steps:
1) given reference signal rk
Reference signal rkIs a polynomial of a time variable k, M representing the highest power of the polynomial; three specific reference signals are as follows:
square wave signal, M ═ 0:
Figure BDA0001809381180000084
triangular wave signal, M ═ 1:
Figure BDA0001809381180000091
s-curve, M ═ 3:
Figure BDA0001809381180000092
wherein, A is amplitude, and N is sampling times of the reference signal in one period;
2) structure ideal error dynamics
For exponential type no-switching attraction law
ek+1=(1-ρ(ek))ek (4)
Figure BDA0001809381180000093
Wherein e isk=rk-ykSystematic tracking error at time k, ykThe actual output signal of the system at the moment k;
embedding interference suppression measures into the attraction law to construct an ideal error dynamic
Figure BDA0001809381180000094
Wherein, delta is more than 0 and less than 1, alpha is more than 0, beta is more than 1, dk+1Is the equivalent interference at the time k +1,
Figure BDA0001809381180000097
the signal is compensated for equivalent interference. Due to the parameter ρ (e)k) At ekAway from the origin, approaches δ at ekApproaches to the origin and approaches to 1, so under the action of an ideal error dynamic state (6), the system error starts from an initial position, the convergence speed gradually increases, and delta, alpha and beta are parameters for adjusting the suction speed;
3) differential compensation measures for interference
The equivalent interference compensation measures are provided as
Figure BDA0001809381180000096
Defining an equivalent interference for a given reference signal in the form of
Equivalent interference signal when M is 0
dk+1=wk+1 (7)
Equivalent interference signal when M is 1
dk+1=wk+1-wk (8)
Equivalent interference signal when M is 3
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (9)
Wherein, wk+1Interference at time k + 1;
and defining the interference difference step number as L, wherein L represents that the interference compensation error comprises L interferences at successive time instants. As in equation (7), a two-step interference differential dk+1-dk=wk+1-wkComprising wk+1And wkTwo moments interfere; in order to effectively suppress interference, the following condition should be satisfied when selecting the equivalent interference
Figure BDA0001809381180000101
Wherein the content of the first and second substances,
Figure BDA0001809381180000102
is the smallest integer not less than.
4) Controller design
According to the ideal error dynamics (6) and the equivalent interference signal dk+1The following controller is easily given:
for the square wave reference signal formula (1), the controller is
Figure BDA0001809381180000103
For the triangular wave reference signal formula (2), the controller is
Figure BDA0001809381180000104
For the S-curve reference signal formula (3), the controller is
Figure BDA0001809381180000105
Formula (11) and formula(12) And in formula (13), a1,a2,b0,b1For serving systems
yk+1+a1yk+a2yk-1=b0uk+b1uk-1+wk+1 (14)
The parameters of (1);
5) analysis of properties
Giving specific expressions of four indexes such as 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, and the like, and being 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 and the monotone decreasing area are defined as follows
Steady state error band ΔSSE
Figure BDA0001809381180000111
Absolute attraction ΔAAL
Figure BDA0001809381180000112
Monotonous decreasing region deltaMDR
Figure BDA0001809381180000113
Equivalent interference compensation error satisfaction
Figure BDA0001809381180000114
Specific expressions of the indexes are as follows:
steady state error band (Δ)SSE)
ΔSSE=Δ (15)
Absolute attraction layer (. DELTA.)AAL)
ΔAAL=max{ΔAAL1AAL2} (16)
Wherein, DeltaAAL1,ΔAAL2Is a real number, determined by
Figure BDA0001809381180000115
Monotonous decreasing area (delta)MDR)
ΔMDR=max{ΔMDR1MDR2} (18)
Wherein, DeltaMDR1,ΔMDR2Is real, and satisfies
Figure BDA0001809381180000116
Number of convergence steps
Figure BDA0001809381180000117
Figure BDA0001809381180000118
Wherein e is0Is the initial value of the tracking error.
As is clear from equation (15), since the steady-state error band of the tracking error decreases with a decrease in Δ, the control accuracy can be improved by adopting the disturbance difference compensation measure.
In the embodiment, the permanent magnet synchronous motor device executes a position accurate tracking task, and a digital controller is designed for position loop control, wherein a current loop and a speed loop controller are provided by an ELMO driver; the position loop controller is provided by a DSP development board TMS320F 2812.
The mathematical model of the servo system is as follows
yk+1-1.5001yk+0.4989yk-1=2.1589uk-0.5113uk-1+wk+1 (21)
When M is 0, the compound is represented by formula (10)
Figure BDA0001809381180000121
When M is 1, the compound is represented by formula (11)
Figure BDA0001809381180000122
When M is 3, the compound is represented by formula (12)
Figure BDA0001809381180000123
The effectiveness of interference difference compensation measures in a discrete servo system is verified through numerical simulation and experimental results.
The simulation is divided into two parts, the first part verifies the specific expression of the performance index given by the formula (15), and the second part verifies the interference suppression effect of the interference difference compensation measure.
(1) Given reference letter rk=10sin(2πfkTs) deg, frequency f 2Hz, sampling period Ts0.01 interference wk0.28| mod (k,20) -10| +0.32| mod (k +7,20) -10 |. Under the action of the controller formula (22), when the controller parameters Δ is 0.6, δ is 0.1, α is 5, and β is 10 (see fig. 2 and 3), the performance index is
ΔSSE=0.6
Figure BDA0001809381180000124
Through simulation, the result shows that the steady-state error band deltaSSEThe formula (15) is satisfied.
(2) The reference signals are respectively a square wave signal formula (1), a triangular wave signal formula (2) and an S curve formula (3), and the amplitude A is respectively 5, 10 and 5. To verify that the interference difference step number L is in the full equation (10), the controller is designed to accurately track the corresponding reference signal, α is 10, β is 10, and δ is 0.5, and the disturbance signal is selected as wk=0.1rkSampling period Ts=0.1。
1) Reference signal rkFor equation (1), controller equation (22) is used, and the simulation results are shown in FIGS. 4-7, where ΔSSE=0deg。
2) Reference signal rkFor equation (2), the controller is adopted as equation (22), and the simulation results are shown in FIGS. 8-11, in which ΔSSE=0.1deg。
3) Reference signal rkFor equation (2), the controller is adopted as equation (23), and the simulation results are shown in FIGS. 12-15, in which ΔSSE=0deg。
4) Reference signal rkFor equation (3), the controller is adopted as equation (22), and the simulation results are shown in FIGS. 16-19, in which ΔSSE=0.2deg。
5) Reference signal rkFor equation (3), the controller is adopted as equation (23), and the simulation results are shown in fig. 20-23, wherein Δ isSSE=0.07deg。
6) Reference signal rkFor equation (3), the controller is adopted as equation (24), and the simulation results are shown in FIGS. 24-27, in which ΔSSE=0deg。
Through simulation (2), it is shown that a controller adopting the interference difference compensation design can realize accurate tracking of a given reference signal when M and L satisfy the formula (10), and the closer L is to each other when M and L do not satisfy the formula (10)
Figure BDA0001809381180000131
The better the tracking effect and the control process has no buffeting.
The control method provided by the invention is verified on a position servo device, and FIG. 1 is a block diagram of a position servo system. The reference signals are respectively a square wave signal formula (1), a triangular wave signal formula (2) and an S curve formula (3), and the maximum powers of the 3 reference signals with respect to the variable k are respectively 0, 1 and 3. The effects of the interference difference compensation technique were verified using controller equations (20), (21), and (22), respectively. Sampling period T in the experimentsThe controller parameter α is 10, β is 10, δ is 0.2ms, and the amplitudes a of the reference signals equation (1), equation (2), and equation (3) are 15deg, 180deg, and 180deg, respectively. The experimental results are as follows:
1) reference signal rkAs shown in formula (1), the controller (22) is adopted to test the knotThe results are shown in fig. 28-29. In the figure, ΔSSE=0.05deg。
2) Reference signal rkThe experimental results are shown in fig. 30-fig. 31 using the controller (22) as in formula (2). In the figure, ΔSSE=0.2deg。
3) Reference signal rkThe experimental results are shown in fig. 32-33 using the controller (23) as in formula (2). In the figure, ΔSSE=0.05deg。
4) Reference signal rkThe experimental results are shown in fig. 34-35 using the controller (22) as in equation (3). In the figure, ΔSSE=0.1deg。
5) Reference signal rkThe results of the experiment are shown in fig. 36-37 using the controller (23) as in equation (3). In the figure, ΔSSE=0.07deg。
6) Reference signal rkThe experimental results are shown in fig. 38-39 using the controller (24) as in equation (3). In the figure, ΔSSE=0.05deg。
Experimental results show that when M and L satisfy the formula (10), the controller adopting the interference difference compensation design can realize accurate tracking of a given reference signal, and when M and L do not satisfy the formula (10), the closer L is
Figure BDA0001809381180000141
The better the system tracking performance and the control process is buffeting free.

Claims (2)

1. An exponential type non-switching attraction law design method of a discrete time controller adopting an interference difference suppression strategy is characterized by comprising the following steps:
1) given a reference signal rkFor a polynomial of the time variable k, M representing the highest power of the polynomial, three specific reference signals are as follows:
square wave signal, M ═ 0:
Figure FDA0002771779250000011
triangular wave signal, M ═ 1:
Figure FDA0002771779250000012
s-curve, M ═ 3:
Figure FDA0002771779250000013
wherein, A is amplitude, and N is sampling times of the reference signal in one period;
2) structure ideal error dynamics
For exponential type no-switching attraction law
ek+1=(1-ρ(ek))ek (4)
Figure FDA0002771779250000021
Wherein e isk=rk-ykSystematic tracking error at time k, ykThe actual output signal of the system at the moment k;
embedding interference suppression measures into the attraction law to construct an ideal error dynamic
Figure FDA0002771779250000022
Wherein, delta is more than 0 and less than 1, alpha is more than 0, beta is more than 1, dk+1Is the equivalent interference at the time k +1,
Figure FDA0002771779250000023
compensating the signal for equivalent interference; due to the parameter ρ (e)k) At ekAway from the origin, approaches δ at ekApproaches to the origin, and thus the system error starts from the initial position under the action of the ideal error dynamics (6), the convergence rate gradually increases, and delta, alpha and beta are used for adjusting the suction speedA parameter of degree;
3) interference differential compensation strategy
Take the equivalent interference compensation effect as
Figure FDA0002771779250000024
Defining an equivalent interference for a particular reference signal in the form of
Equivalent interference signal when M is 0
dk+1=wk+1 (7)
Equivalent interference signal when M is 1
dk+1=wk+1-wk (8)
Equivalent interference signal when M is 3
dk+1=[(wk+1-wk)-(wk-wk-1)]-[(wk-wk-1)-(wk-1-wk-2)] (9)
Wherein, wk+1Interference at time k + 1;
defining the step number of the interference difference as L, wherein L represents that the interference compensation error comprises L interferences at successive moments, and the two-step interference difference d is shown as formula (7)k+1-dk=wk+1-wkComprising wk+1And wkTwo-time interference, in order to effectively suppress interference, when selecting equivalent interference, the following condition should be satisfied
Figure FDA0002771779250000025
Wherein the content of the first and second substances,
Figure FDA0002771779250000026
is the smallest integer not less than;
4) controller design
According to the ideal error dynamics (6) and the equivalent interference signal dk+1The following controller expression is given:
for a square wave reference signal equation (1),
Figure FDA0002771779250000027
wherein, F (q)-1)=B(q-1)-b0
For the triangular wave reference signal equation (2),
Figure FDA0002771779250000031
wherein the content of the first and second substances,
Figure FDA0002771779250000032
for the S-curve reference signal equation (3),
Figure FDA0002771779250000033
in formulae (11) to (13), A (q)-1)、B(q-1) For serving systems
A(q-1)yk=q-1B(q-1)uk+wk (14)
With respect to q-1The parameter polynomial of (2):
A(q-1)=1+a1q-1+a2q-2+……+anq-n
B(q-1)=b0+b1q-1+b2q-2+……+bmq-m
wherein u iskAnd ykInput and output signals, q, respectively, at time k of the servo system-1For one-step delay operator, m and n are respectively A (q)-1)、B(q-1) Order of (b)0≠0,1≤m≤n。
2. The method of claim 1, wherein the method further comprises the following steps:
5) analysis of properties
And giving specific expressions of a steady-state error band, an absolute attraction layer, a monotone subtraction area and at most four indexes of steps required by the tracking error entering the steady-state error band for the first time, and being 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 and the monotone subtraction area are defined as follows:
steady state error band ΔSSE
Figure FDA0002771779250000034
Absolute attraction layer ΔAAL
Figure FDA0002771779250000035
Monotonous decreasing region deltaMDR
Figure FDA0002771779250000036
Equivalent interference compensation error satisfaction
Figure FDA0002771779250000041
The expression of each index is as follows:
steady state error band ΔSSE
ΔSSE=Δ (15)
Absolute attraction layer ΔAAL
ΔAAL=max{ΔAAL1AAL2} (16)
In the formula,. DELTA.AAL1,ΔAAL2Is a real number, determined by
Figure FDA0002771779250000042
Monotonous decreasing region deltaMDR
ΔMDR=max{ΔMDR1MDR2} (18)
In the formula,. DELTA.MDR1,ΔMDR2Is real, and satisfies
Figure FDA0002771779250000043
Number of convergence steps
Figure FDA0002771779250000044
Figure FDA0002771779250000045
Wherein e is0Is the initial value of the tracking error.
CN201811112071.5A 2018-09-25 2018-09-25 Digital controller design method based on exponential type non-switching attraction law Active CN109188908B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811112071.5A CN109188908B (en) 2018-09-25 2018-09-25 Digital controller design method based on exponential type non-switching attraction law

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811112071.5A CN109188908B (en) 2018-09-25 2018-09-25 Digital controller design method based on exponential type non-switching attraction law

Publications (2)

Publication Number Publication Date
CN109188908A CN109188908A (en) 2019-01-11
CN109188908B true CN109188908B (en) 2021-02-26

Family

ID=64909439

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811112071.5A Active CN109188908B (en) 2018-09-25 2018-09-25 Digital controller design method based on exponential type non-switching attraction law

Country Status (1)

Country Link
CN (1) CN109188908B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110297426B (en) * 2019-08-14 2022-02-25 台州学院 Parallel type half-cycle repetitive control method adopting periodic interference differential suppression strategy
CN110456645B (en) * 2019-08-18 2022-02-15 台州学院 Discrete repetitive control method for inverter
CN110412876B (en) * 2019-08-18 2022-02-08 台州学院 Inverter repetitive control method based on attraction law
CN111752149B (en) * 2020-06-08 2022-12-20 浙江工业大学 Design method of digital repetitive controller for designated time attraction of servo motor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3292951B2 (en) * 1998-11-13 2002-06-17 住友重機械工業株式会社 Position control device
WO2004059400A1 (en) * 2002-12-27 2004-07-15 Inea D.O.O. Self-tuning controller for non-linear processes described by set of local linear models
CN103048921A (en) * 2012-11-12 2013-04-17 浙江工业大学 Half-cycle repetitive controller for position servo system
CN104898427A (en) * 2015-05-19 2015-09-09 浙江工业大学 Segmentation filtering iterative learning control method of motor servo system
CN107490966A (en) * 2017-08-22 2017-12-19 浙江工业大学 A kind of aircraft finite time Adaptive Attitude control method based on improvement power Reaching Law
CN107544244A (en) * 2017-08-25 2018-01-05 浙江工业大学 Based on the oval discrete repetitive control for motor servo system for attracting rule and the compensation of equivalent disturbance expansion state

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3292951B2 (en) * 1998-11-13 2002-06-17 住友重機械工業株式会社 Position control device
WO2004059400A1 (en) * 2002-12-27 2004-07-15 Inea D.O.O. Self-tuning controller for non-linear processes described by set of local linear models
CN103048921A (en) * 2012-11-12 2013-04-17 浙江工业大学 Half-cycle repetitive controller for position servo system
CN104898427A (en) * 2015-05-19 2015-09-09 浙江工业大学 Segmentation filtering iterative learning control method of motor servo system
CN107490966A (en) * 2017-08-22 2017-12-19 浙江工业大学 A kind of aircraft finite time Adaptive Attitude control method based on improvement power Reaching Law
CN107544244A (en) * 2017-08-25 2018-01-05 浙江工业大学 Based on the oval discrete repetitive control for motor servo system for attracting rule and the compensation of equivalent disturbance expansion state

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Convergence Performance of Discrete Power Attracting Law;lingwei wu等;《2018 IEEE 7th Data Driven Control and Learning Systems Conference》;20180531;第974-978页 *
带扰动补偿的无抖振离散重复控制器设计;周文委等;《控制与决策》;20170331;第32卷(第3期);第467-473页 *

Also Published As

Publication number Publication date
CN109188908A (en) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109188908B (en) Digital controller design method based on exponential type non-switching attraction law
CN103048921B (en) For the semiperiod repetitive controller of positional servosystem
CN107425548B (en) Interpolation H∞Dynamic state estimation method for extended Kalman filter generator
CN107797448B (en) Motor position discrete repetitive control method adopting disturbance expansion 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
CN110134010B (en) Power attraction repetitive control method adopting equivalent disturbance compensation servo system
CN108828958B (en) Non-switching attraction law design method of discrete time controller adopting interference differential suppression strategy
CN108155833B (en) Motor servo system asymptotic stable control method considering electrical characteristics
CN108958041B (en) Discrete double-period repetitive control method based on hyperbolic secant attraction law
CN103399485B (en) For the partial periodicity repetitive controller of positional servosystem
CN103197556B (en) Based on 1/2nd cycle repetitive control attracting rule
US20120019217A1 (en) Control concept for a digitally controlled magnetic supply device
CN111546346B (en) Flexible joint disturbance observation method, moment control method and equipment
CN110658719A (en) Servo system non-switching attraction repetitive control method adopting equivalent disturbance compensation
CN109100938B (en) Discrete time controller attraction law design method adopting interference differential compensation
Takemura et al. Proposal of novel rolling friction compensation with data-based friction model for ball screw driven stage
CN107544250B (en) Less-conservative permanent magnet spherical motor trajectory tracking control method
CN111752149B (en) Design method of digital repetitive controller for designated time attraction of servo motor
CN110673477A (en) Design method of discrete time power attraction controller adopting interference difference suppression strategy
CN115313931A (en) Sensor-free vector control method of permanent magnet synchronous motor based on AEKF
CN111835250B (en) 1/2 power finite value attraction repetitive control method for servo motor driving system
CN112117945A (en) Discrete repetitive control method of parabola approximation law adopting disturbance compensation
CN111856941B (en) Adaptive terminal dynamic sliding mode control method based on active disturbance rejection
CN113411026A (en) Design method of self-adaptive sliding film disturbance observer for permanent magnet synchronous linear motor

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