PARTIAL-FORM MODEL-FREE ADAPTIVE DISTURBANCE COMPENSATION CONTROL IN THE PRESENCE OF MEASURABLE DISTURBANCES

Abstract

A method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, includes establishing a dynamic data model of a controlled plant subject to measurable disturbances, wherein the dynamic data model is described by a pseudo Jacobian input matrix and a pseudo Jacobian disturbance matrix; constructing cost functions and solving their optimization problems to find optimal values of the pseudo Jacobian input matrix and the pseudo Jacobian disturbance matrix; designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances; constructing an energy function and solving it by using a momentum gradient descent method to find optimal values of the partial-form adaptive input matrix and the partial-form adaptive disturbance matrix; controlling the controlled plant by using the control law. The control method of the present invention provides significant improvements in disturbance compensation control performance and achieves effective tracking of desired system outputs.

Claims

1. A method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, executed on a hardware platform for controlling a controlled plant subject to measurable disturbances, said controlled plant being a multi-input multi-output (MIMO) system with a predetermined number of control inputs and a predetermined number of system outputs, said method comprising: step 1: obtaining measurable disturbances at time k, establishing a dynamic data model of said controlled plant subject to measurable disturbances, wherein said dynamic data model is described by a pseudo Jacobian input matrix ?(k) and a pseudo Jacobian disturbance matrix ?(k); step 2: constructing cost functions and solving optimization problems for said cost functions to find an optimal value of said pseudo Jacobian input matrix ?(k) in said step 1 and an optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 1; step 3: utilizing said measurable disturbances at time k, employing said dynamic data model described by said optimal value of said pseudo Jacobian input matrix ?(k) and said optimal value of said pseudo Jacobian disturbance matrix x(k) in said step 2, designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances, wherein said control law comprising a partial-form adaptive input matrix ?.sub.p(k) and a partial-form adaptive disturbance matrix ?.sub.p(k); step 4: constructing an energy function and solving said energy function by using a momentum gradient descent method to find an optimal value of said partial-form adaptive input matrix ?.sub.p(k) in said step 3 and an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 3; step 5: controlling said controlled plant by using said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of partial-form adaptive input matrix ?.sub.p(k) and said optimal value of partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4, weakening the effect of measurable disturbances on actual system outputs of said controlled plant, achieving effective tracking of desired system outputs of said controlled plant.

2. The method as claimed in claim 1 wherein said step 1, obtaining measurable disturbances at time k, establishing a dynamic data model of said controlled plant subject to measurable disturbances as
?y(k+1)=?(k)?u(k)+?(k)?d(k) where k is a sampling time, k is a positive integer; y (k+1) is an actual system output vector of said controlled plant at time k+1, y(k+1)=[y.sub.1(k+1), . . . , y.sub.n(k +1)].sup.T, ?y(k+1)=y(k+1)?y(k); n is a total number of system outputs in said controlled plant, n is a positive integer greater than 1; u(k) is a control input vector of said controlled plant at time k, u(k)=[u.sub.1(k), . . . , ?.sub.m(k)].sup.T, ?u(k)=u(k)?u(k?1); m is a total number of control inputs in said controlled plant, m is a positive integer greater than 1; d(k) is a measurable disturbance vector in said controlled plant at time k, d(k)=[d.sub.1(k), . . . , d.sub.q(k)].sup.T, ?d(k)=d(k)?d(k?1); q is a total number of measurable disturbances in said controlled plant, q is a positive integer; ?k) is said pseudo Jacobian input matrix at time k and ?(k) is said pseudo Jacobian disturbance matrix at time k.

3. The method as claimed in claim 1 wherein said step 2, constructing cost functions and solving optimization problems for said cost functions to find an optimal value of said pseudo Jacobian input matrix ?(k) in said step 1 and an optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 1, comprising: step 2.1: constructing a cost function for said pseudo Jacobian input matrix ?(k) as
J(?(k))=??y(k)??(k)?u(k?1)??(k?1)?d(k?1)?.sup.2+?.sub.1???(k)?.sup.2 where ?.sub.1 is the first weighting factor; step 2.2: constructing a cost function for said pseudo Jacobian disturbance matrix ?(k) as
J(?(k))=??y(k)??(k?1)?u(k?1)??(k)?d(k?1)?.sup.2+?.sub.2???(k)?.sup.2 where ?.sub.2 is the second weighting factor; step 2.3 solving an optimization problem for said J(?(k)) in said step 2.1, finding an optimal value of said pseudo Jacobian input matrix ?(k) as $?\left(k\right)=?\left(k-1\right)+\frac{{?}_1\left(?y\left(k\right)-?\left(k-1\right)?u\left(k-1\right)-?\left(k-1\right)?d\left(k-1\right)\right)?{u\left(k-1\right)}^T}{{?}_1+{.Math.?u\left(k-1\right).Math.}^2}$ where ?.sub.1 is the first step size factor; step 2.4: solving an optimization problem for said J(?(k)) in said step 2.2, finding an optimal value of said pseudo Jacobian disturbance matrix ?(k) as $?\left(k\right)=?\left(k-1\right)+\frac{{?}_2\left(?y\left(k\right)-?\left(k-1\right)?u\left(k-1\right)-?\left(k-1\right)?d\left(k-1\right)\right)?{d\left(k-1\right)}^T}{{?}_2+{.Math.?d\left(k-1\right).Math.}^2}$ where ?.sub.2 is the second step size factor.

4. The method as claimed in claim 1 wherein said step 3, utilizing said measurable disturbances at time k, employing said dynamic data model described by said optimal value of said pseudo Jacobian input matrix ?(k) and said optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 2, designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances as
u(k)=u(k?1)+?.sub.p(k)?E(k)+?.sub.p(k)?D(k) where ?E(k)=[?e(k).sup.T, ?e(k).sup.T, . . . , ?e(k?L+2).sup.2].sup.T, ?D(k)=[?d(k).sup.T, . . . , ?d(k?L+1).sup.T].sup.T; e(k) is a system error vector of said controlled plant at time k, e(k)=y*(k)?y(k), e(k)=[e.sub.1(k), . . . , e.sub.n(k)].sup.T, ?e(k)=e(k)?e(k?1); L is linearized length constant and is positive integer; ?.sub.p(k) is said partial-form adaptive input matrix at time k and ?.sub.p(k) is said partial-form adaptive disturbance matrix at time k.

5. The method as claimed in claim 1 wherein said step 4, constructing an energy function and solving said energy function by using a momentum gradient descent method to find an optimal value of said partial-form adaptive input matrix ?.sub.p(k) in said step 3 and an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 3, comprising: step 4.1: constructing an energy function as $W=\frac{1}{2}{.Math.y^{*}\left(k+1\right)-y\left(k+1\right).Math.}^2+\frac{1}{2}?{.Math.?u\left(k\right).Math.}^2$ where y*(k+1) is a desired system output vector of said controlled plant at time k+1; y*(k+1)=[y.sub.1*(k+1), . . . , y.sub.n*(k+1)].sup.T; ? is a penalty factor; step 4.2: solving said energy function in said step 4.1 by using a momentum gradient descent method, finding an optimal value of said partial-form adaptive input matrix ?.sub.p(k) as ${?}_p\left(k\right)={?}_p\left(k-1\right)-{?}_1\left(1-{?}_1\right)\frac{?W}{?{?}_p\left(k-1\right)}+{?}_1?{?}_p\left(k-1\right)$ where ?.sub.1 is the first learning rate, ?.sub.1 is the first momentum factor; $?{?}_p\left(k-1\right)={?}_p\left(k-1\right)-{?}_p\left(k-2\right);\frac{?W}{?{?}_p\left(k-1\right)}$ is a partial derivative of said energy function W to ?.sub.p(k?1); solving said energy function by using a momentum gradient descent method, finding an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) as ${?}_p\left(k\right)={?}_p\left(k-1\right)-{?}_2\left(1-{?}_2\right)\frac{?W}{?{?}_p\left(k-1\right)}+{?}_2{\mathrm{??}}_p\left(k-1\right)$ where ?.sub.2 is the second learning rate, ?.sub.2 is the second momentum factor; ${\mathrm{??}}_p\left(k-1\right)={?}_p\left(k-1\right)-{?}_p\left(k-2\right);\frac{?W}{?{?}_p\left(k-1\right)}$ is a partial derivative of said energy function W to ?.sub.p(k?1).

6. The method as claimed in claim 5 wherein said partial derivative of said energy function W to ?.sub.p(k?1) in said step 4.2 is calculated as $\frac{?W}{?{?}_p\left(k-1\right)}=-\frac{?y\left(k\right)}{?u\left(k-1\right)}e\left(k\right)?{E\left(k-1\right)}^T+??u\left(k-1\right)?{E\left(k-1\right)}^T;$ said partial derivative of said energy function W to ?.sub.p(k?1) in said step 4.3 is calculated as $\frac{?W}{?{?}_p\left(k-1\right)}=-\frac{?y\left(k\right)}{?u\left(k-1\right)}e\left(k\right)?{D\left(k-1\right)}^T+??u\left(k-1\right)?{D\left(k-1\right)}^T$

7. The method as claimed in claim 6 wherein said $\frac{?y\left(k\right)}{?u\left(k-1\right)}$ is calculated as $\frac{?y\left(k\right)}{?u\left(k-1\right)}={?\left(k\right)}^T.$

8. The method as claimed in claim 1 wherein said step 5, controlling said controlled plant by using said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of partial-form adaptive input matrix ?.sub.p(k) and said optimal value of partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4, comprising: step 5.1: obtaining a measurable disturbance vector d(k); step 5.2: obtaining a desired system output vector y*(k) and an actual system output vector y(k), calculating a system error vector e(k); step 5.3: based on said steps 5.1 and 5.2, calculating a control input vector u(k) according to said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of said partial-form adaptive input matrix ?.sub.p(k) and said optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4; step 5.4: generating an actual system output vector of said controlled plant based on application of said control input vector u(k).

9. A non-transitory computer-readable storage medium having a computer program stored thereon, wherein when said computer program is executed by a processor, causing said processor to carry out said method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances in claim 1.

10. An electronic device comprising a memory, a processor and a computer program stored on said memory and runnable on said processor, wherein when said processor executes said computer program, causing said processor to carry out said method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances in claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0034] FIG. 1 shows a schematic diagram according to the embodiments of the present invention;

[0035] FIG. 2 shows a block diagram of an engineering application according to the embodiments of the present invention;

[0036] FIG. 3 shows a schematic diagram of a hardware platform for running the embodiments of the present invention;

[0037] FIG. 4 shows the tracking performance of the first system output when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0038] FIG. 5 shows the tracking performance of the second system output when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0039] FIG. 6 shows the first control input when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0040] FIG. 7 shows the second control input when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0041] FIG. 8 shows a refrigeration cycle of a vapor compression refrigeration system;

[0042] FIG. 9 shows the two measurable disturbances in the vapor compression refrigeration system;

[0043] FIG. 10 shows the tracking performance of the first system output when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0044] FIG. 11 shows the tracking performance of the second system output when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method;

[0045] FIG. 12 shows the first control input when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; and

[0046] FIG. 13 shows the second control input when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method.

DETAILED DESCRIPTION OF THE INVENTION

[0047] The present invention is hereinafter described in detail with reference to the embodiments and accompanying drawings. It is to be understood that other embodiments may be utilized and structural changes may be made without departing from the scope of the present invention.

[0048] FIG. 1 shows a schematic diagram according to the embodiments of the present invention. A method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, includes establishing a dynamic data model of a controlled plant subject to measurable disturbances, wherein the dynamic data model is described by a pseudo Jacobian input matrix and a pseudo Jacobian disturbance matrix; constructing cost functions and solving their optimization problems to find optimal values of the pseudo Jacobian input matrix and the pseudo Jacobian disturbance matrix; designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances; constructing an energy function and solving it by using a momentum gradient descent method to find optimal values of the partial-form adaptive input matrix and the partial-form adaptive disturbance matrix; controlling the controlled plant by using the control law.

[0049] In the following, the implementation steps of the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances are further explained: [0050] the method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, executed on a hardware platform for controlling a controlled plant subject to measurable disturbances, said controlled plant being a multi-input multi-output (MIMO) system with a predetermined number of control inputs and a predetermined number of system outputs, said method comprising:

[0051] step 1: obtaining measurable disturbances at time k, establishing a dynamic data model of said controlled plant subject to measurable disturbances, wherein said dynamic data model is described by a pseudo Jacobian input matrix ?(k) and a pseudo Jacobian disturbance matrix ?(k);

[0052] step 2: constructing cost functions and solving optimization problems for said cost functions to find an optimal value of said pseudo Jacobian input matrix ?(k) in said step 1 and an optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 1;

[0053] step 3: utilizing said measurable disturbances at time k, employing said dynamic data model described by said optimal value of said pseudo Jacobian input matrix ?(k) and said optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 2, designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances, wherein said control law comprising a partial-form adaptive input matrix ?.sub.p(k) and a partial-form adaptive disturbance matrix ?.sub.p(k);

[0054] step 4: constructing an energy function and solving said energy function by using a momentum gradient descent method to find an optimal value of said partial-form adaptive input matrix ?.sub.p(k) in said step 3 and an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 3;

[0055] step 5: controlling said controlled plant by using said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of partial-form adaptive input matrix ?.sub.p(k) and said optimal value of partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4, weakening the effect of measurable disturbances on actual system outputs of said controlled plant, achieving effective tracking of desired system outputs of said controlled plant.

[0056] Said step 1, obtaining measurable disturbances at time k, establishing a dynamic data model of said controlled plant subject to measurable disturbances as

?y(k+1)??(k)?u(k)+?(k)?d(k)

where k is a sampling time, k is a positive integer; y(k+1) is an actual system output vector of said controlled plant at time k+1, y(k+1)=[y.sub.1(k+1), . . . , y.sub.n(k+1)].sup.T, ?y(k+1)=y(k+1)?y(k); n is a total number of system outputs in said controlled plant, n is a positive integer greater than 1; u(k) is a control input vector of said controlled plant at time k, u(k)=[u.sub.1(k), . . . , u.sub.m(k)].sup.T, ?u(k)=u(k)?u(k?1); m is a total number of control inputs in said controlled plant, m is a positive integer greater than 1; d(k) is a measurable disturbance vector in said controlled plant at time k, d(k)=d (k)=[d.sub.1(k), . . . , d.sub.q(k)].sup.T, ?d(k)=d(k)?d(k?1); q is a total number of measurable disturbances in said controlled plant, q is a positive integer; ?(k) is said pseudo Jacobian input matrix at time k and ?(k) is said pseudo Jacobian disturbance matrix at time k.

[0057] Said step 2, constructing cost functions and solving optimization problems for said cost functions to find an optimal value of said pseudo Jacobian input matrix ?(k) in said step 1 and an optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 1, comprising:

[0058] step 2.1: constructing a cost function for said pseudo Jacobian input matrix ?(k) as

J(?(k)=??y(k)??(k)?u(k?1)??(k?1)?d(k?1)?.sup.2+?.sub.1???(k)?.sup.2

where ?.sub.1 is the first weighting factor; ??(k)=?(k)??(k?1); ??? is a Euclidean norm;

[0059] step 2.2: constructing a cost function for said pseudo Jacobian disturbance matrix ?(k) as

J(?(k))=??y(k)??(k?1)?u(k?1)??(k)?d(k?1)?.sup.2+?.sub.2???(k)?.sup.2

where ?.sub.2 is the second weighting factor; ??(k)=?(k)??(k?1);

[0060] step 2.3: solving an optimization problem for said J(?(k)) in said step 2.1, finding an optimal value of said pseudo Jacobian input matrix ?(k) as

[00012]$?\left(k\right)=?\left(k-1\right)+\frac{{?}_1\left(?y\left(k\right)-?\left(k-1\right)?u\left(k-1\right)-?\left(k-1\right)?d\left(k-1\right)\right)?{u\left(k-1\right)}^T}{{?}_1+{.Math.?u\left(k-1\right).Math.}^2}$

where ?.sub.1 is the first step size factor;

[0061] step 2.4: solving an optimization problem for said J(?(k)) in said step 2.2, finding an optimal value of said pseudo Jacobian disturbance matrix ?(k) as

[00013]$?\left(k\right)=?\left(k-1\right)+\frac{{?}_2\left(?y\left(k\right)-?\left(k-1\right)?u\left(k-1\right)-?\left(k-1\right)?d\left(k-1\right)\right)?{d\left(k-1\right)}^T}{{?}_2+{.Math.?d\left(k-1\right).Math.}^2}$

where ?.sub.2 is the second step size factor.

[0062] Said step 3, utilizing said measurable disturbances at time k, employing said dynamic data model described by said optimal value of said pseudo Jacobian input matrix ?(k) and said optimal value of said pseudo Jacobian disturbance matrix ?(k) in said step 2, designing a partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances as

u(k)=u(k?1)+?.sub.p(k)?E(k)+?.sub.p(k)?D(k)

where ?E(k)=[?e(k).sup.T, ?e(k).sup.T, . . . , ?e(k?L+2).sup.T].sup.T, ?D(k)=[?d(k).sup.T, . . . , ?d(k?L+1).sup.T].sup.T; e(k) is a system error vector of said controlled plant at time k, e(k)=y*(k)?y(k), y*(k) is a desired system output vector of said controlled plant at time k, e(k)=[e.sub.1(k), . . . , e.sub.n(k)].sup.T, ?e(k)=e(k)?e(k?1); L is linearized length constant and is positive integer; ?.sub.p(k) is said partial-form adaptive input matrix at time k and ?.sub.p(k) is said partial-form adaptive disturbance matrix at time k.

[0063] Said step 4, constructing an energy function and solving said energy function by using a momentum gradient descent method to find an optimal value of said partial-form adaptive input matrix ?.sub.p(k) in said step 3 and an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 3, comprising:

[0064] step 4.1: constructing an energy function as

[00014]$W=\frac{1}{2}{.Math.y^{*}\left(k+1\right)-y\left(k+1\right).Math.}^2+\frac{1}{2}?{.Math.?u\left(k\right).Math.}^2$

where y*(k+1) is a desired system output vector of said controlled plant at time k+1; y*(k+1)=[y.sub.1*(k+1), . . . , y.sub.n*(k+1)].sup.T; ? is a penalty factor;

[0065] step 4.2: solving said energy function in said step 4.1 by using a momentum gradient descent method, finding an optimal value of said partial-form adaptive input matrix ?.sub.p(k) as

[00015]${?}_p\left(k\right)={?}_p\left(k-1\right)-{?}_1\left(1-{?}_1\right)\frac{?W}{?{?}_p\left(k-1\right)}+{?}_1?{?}_p\left(k-1\right)$

where ?.sub.1 is the first learning rate, n is the first momentum factor;

[00016]$?{?}_p\left(k-1\right)={?}_p\left(k-1\right)-{?}_p\left(k-2\right);\frac{?W}{?{?}_p\left(k-1\right)}$

is a partial derivative of said energy function W to ?.sub.p(k?1);

[0066] solving said energy function by using a momentum gradient descent method, finding an optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) as

[00017]${?}_p\left(k\right)={?}_p\left(k-1\right)-{?}_2\left(1-{?}_2\right)\frac{?W}{?{?}_p\left(k-1\right)}+{?}_2{\mathrm{??}}_p\left(k-1\right)$

where ?.sub.2 is the second learning rate, ?.sub.2 is the second momentum factor;

[00018]$?{?}_p\left(k-1\right)={?}_p\left(k-1\right)-{?}_p\left(k-2\right);\frac{?W}{?{?}_p\left(k-1\right)}$

is a partial derivative of said energy function W to ?.sub.p(k?1).

[0067] Said partial derivative of said energy function W to ?.sub.p(k?1) in said step 4.2 is calculated as

[00019]$\frac{?W}{?{?}_p\left(k-1\right)}=-\frac{?y\left(k\right)}{?u\left(k-1\right)}e\left(k\right)?{E\left(k-1\right)}^T+??u\left(k-1\right)?{E\left(k-1\right)}^T;$

[0068] said partial derivative of said energy function W to ?.sub.p(k?1) in said step 4.3 is calculated as

[00020]$\frac{?W}{?{?}_p\left(k-1\right)}=-\frac{?y\left(k\right)}{?u\left(k-1\right)}e\left(k\right)?{D\left(k-1\right)}^T+??u\left(k-1\right)?{D\left(k-1\right)}^T.$

said

[00021]$\frac{?y\left(k\right)}{?u\left(k-1\right)}$

is calculated as

[00022]$\frac{?y\left(k\right)}{?u\left(k-1\right)}={?\left(k\right)}^T.$

[0069] Said step 5, controlling said controlled plant by using said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of partial-form adaptive input matrix ?.sub.p(k) and said optimal value of partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4, comprising:

[0070] step 5.1: obtaining a measurable disturbance vector d(k);

[0071] step 5.2: obtaining a desired system output vector y*(k) and an actual system output vector y(k), calculating a system error vector e(k);

[0072] step 5.3: based on said steps 5.1 and 5.2, calculating a control input vector u(k) according to said partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with said optimal value of said partial-form adaptive input matrix ?.sub.p(k) and said optimal value of said partial-form adaptive disturbance matrix ?.sub.p(k) in said step 4;

[0073] step 5.4: generating an actual system output vector of said controlled plant based on application of said control input vector u(k).

[0074] FIG. 2 shows a block diagram of an engineering application according to the embodiments of the present invention. For the hardware platform in FIG. 2, a schematic diagram of the hardware platform for running the embodiments of the present invention is shown in FIG. 3; specifically, a non-transitory computer-readable storage medium having a computer program stored thereon, wherein when said computer program is executed by a processor, causing said processor to carry out the method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances; an electronic device comprising a memory, a processor, and a computer program stored on said memory and runnable on said processor, wherein when said processor executes said computer program, causing said processor to carry out the method of partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances.

[0075] Two exemplary embodiments of the present invention are given for further explanation.

[0076] The first exemplary embodiment:

[0077] a two-input two-output controlled plant:

[00023]$x_{11}\left(k+1\right)=\frac{x_{11}^2\left(k\right)}{1+x_{11}^2\left(k\right)}+0.3x_{12}\left(k\right)$$x_{12}\left(k+1\right)=\frac{x_{11}^2\left(k\right)}{1+x_{12}^2\left(k\right)+x_{21}^2\left(k\right)+x_{22}^2\left(k\right)}+a\left(k\right)u_1\left(k\right)+d_1\left(k\right)$$x_{21}\left(k+1\right)=\frac{x_{21}^2\left(k\right)}{1+x_{21}^2\left(k\right)}+0.2x_{22}\left(k\right)$$x_{22}\left(k+1\right)=\frac{x_{21}^2\left(k\right)}{1+x_{11}^2\left(k\right)+x_{12}^2\left(k\right)+x_{22}^2\left(k\right)}+b\left(k\right)u_2\left(k\right)+d_2\left(k\right)$$y_1\left(k+1\right)=x_{11}\left(k+1\right)$$y_2\left(k+1\right)=x_{21}\left(k+1\right)$

where a(k)=1+0.1 sin(2?k/1500), b(k)=1+0.1 cos(2?k/1500) are two time-varying parameters; d.sub.1(k)=0.15 sin(k/10), d.sub.2(k)=0.15 sin(k/10) are two measurable disturbances. The two-input two-output controlled plant is subject to measurable disturbances.

[0078] The desired system outputs are as follows:

[00024]$y_1^{*}\left(k\right)=\{\begin{array}{cc}2,& k?550\\ -0.15,& k>550\end{array}$$y_2^{*}\left(k\right)=\{\begin{array}{c}0,k?700\\ 1,k>700\end{array}$

[0079] In this embodiment, m=n=q=2.

[0080] As is known to all, the existing PID is a well-established and widely used control method in the field of control theory and engineering, which is used for comparison in the embodiments of the present invention. To quantitatively compare the control performance of the control method of the present invention with the existing PID control method, the integral time-weighted absolute error (ITAE) is used as the control performance index for evaluation:

[00025]$\mathrm{ITAE}\left(e_j\right)=\underset{k=1}{\overset{N}{.Math.}}k.Math.e_j\left(k\right).Math.$

where e.sub.j(k)=y.sub.j*(k)?y.sub.j(k), y.sub.j*(k) is the j-th desired system output at time k, y.sub.j(k) is the j-th actual system output at time k, j=1, . . . , n . The smaller the value of ITAE(e.sub.j), the smaller the error between the j-th actual system output and the j-th desired system output, the higher control accuracy and response speed, and the better the control performance.

[0081] The hardware platform for running the embodiment of the present invention is the industrial control computer.

[0082] The embodiment of the present invention adopts the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances; set the parameters of the control method of the present invention as L=3, ?(1)=[0.6, ?0.05; 0.1, 0.6], ?(1)=[0.2, 0; 0, 0.2], ?.sub.p(1)=[?0.35, 0, 0.1, 0, 0.08, 0; 0, ?0.27, 0, 0.1, 0, 0.05], ?.sub.p(1)=[?0.35, 0, 0.01, 0, 0.01, 0; 0, ?0.9, 0, 0.05, 0, 0.01], ?.sub.1=0.5, ?.sub.2=0.5 , ?.sub.1=1, ?.sub.2=0.9, ?.sub.1=0.7, ?.sub.2=0.5, ?.sub.1=0.25, ?.sub.2=0.25, ?=2.

[0083] When controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, the following steps are included at each time k: a) obtaining the measurable disturbance vector d(k); b) obtaining the desired system output vector y*(k) and the actual system output vector y (k), calculating the system error vector e(k); c) based on steps a) and b), calculating the control input vector u(k) according to the partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with the optimal value of the partial-form adaptive input matrix ?.sub.p(k) and the optimal value of the partial-form adaptive disturbance matrix ?.sub.p(k); d) controlling the controlled plant by applying the control input vector u(k), generating the actual system output vector of the controlled plant at time k+1 based on the application of the control input vector; e) repeating steps a) to d) until the end of the control time.

[0084] The control performance of the control method of the present invention and the PID control method are given in FIGS. 4-7. FIG. 4 shows the tracking performance of the first system output when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; FIG. 5 shows the tracking performance of the second system output when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; FIG. 6 shows the first control input when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; and FIG. 7 shows the second control input when controlling the two-input two-output controlled plant in the first exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; comparing the control performance indexes, the ITAE(e.sub.1) of the first system output controlled by the control method of the present invention is 12090, the ITAE(e.sub.2) of the second system output controlled by the control method of the present invention is 10029, the ITAE(e.sub.1) of the first system output controlled by the PID control method is 24067, the ITAE(e.sub.2) of the second system output controlled by the PID control method is 25998, the comparison results of the control performance indexes are listed in Table 1; comparing the system output curves, the control method of the present invention can effectively suppress the effect of measurable disturbances on actual system outputs of the two-input two-output controlled plant, and the control performance using the control method of the present invention is better than that using the PID control method. From the above comparison results, it can be seen that the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances provides a significant improvement in disturbance compensation control performance and achieves effective tracking of the desired system outputs.

TABLE-US-00001 TABLE 1 Comparison Results of The Control Performance Index ITAE(e.sub.1) ITAE(e.sub.2) PID 24067 25998 The present invention 12090 10029

[0085] The second exemplary embodiment:

[0086] Vapor Compression Refrigeration Systems (VCRS) are the most common refrigeration cycle equipment used in homes (e.g., home refrigerators, air conditioners), commercial (e.g., building and automotive air conditioning, refrigerated warehouses) and industrial (e.g., petrochemical plants, natural gas processing plants), and the refrigeration cycle of the vapor compression refrigeration system is shown in FIG. 8. Two disturbances of the vapor compression refrigeration system are the inlet temperature of the condenser secondary flux and the inlet temperature of the evaporator secondary flux. Nowadays, high energy-consuming refrigeration equipment is used in large quantities, and the realization of disturbance compensation control of VCRS is of great significance to the promotion of energy conservation and consumption reduction in the world.

[0087] The vapor compression refrigeration system is a two-input two-output controlled plant. Two control inputs u.sub.1 and u.sub.2 of the vapor compression refrigeration system are the compressor frequency (Hz) and valving opening (%), respectively. Two system outputs y.sub.1 and y.sub.2 of the vapor compression refrigeration system are the degree of superheating (? C.) and the outlet temperature of the evaporator secondary flux (? C.), respectively. Two disturbances d.sub.1 and d.sub.2 of the vapor compression refrigeration system are the inlet temperature of the condenser secondary flux (? C.) and the inlet temperature of the evaporator secondary flux (? C.), respectively, where d.sub.1, d.sub.2 are measured online via temperature sensors and are therefore measurable disturbances. FIG. 9 shows the two measurable disturbances in the vapor compression refrigeration system. The vapor compression refrigeration system is a two-input two-output controlled plant subject to measurable disturbances. In this embodiment, m=n=q=2. The hardware platform for running the embodiment of the present invention is the industrial control computer.

[0088] The initial conditions of the vapor compression refrigeration system are given as: u.sub.1 (0)=36.45 Hz, u.sub.2(0)=48.79%, y.sub.1(0)=14.65? C., y.sub.2(0)=?22.15? C. To meet the cooling demand of the evaporator secondary flux, the desired system output y.sub.1* is adjusted from 14.65? C. to 7.2? ? C. at the 2nd min, then adjusted from 7.2? C. to 22.2? C. at the 9th min, and finally adjusted from 22.2? C. to 11.65? C. at the 16th min, the desired system output y.sub.2 is adjusted from ?22.15? C. to ?22.65? C. at the 2nd min.

[0089] The embodiment of the present invention adopts the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances; set the parameters of the control method of the present invention as L=2, ?(1)=[2, 0; 0, 0.1], ?(1)=[0.2, 0; 0, 0.2], ?.sub.p(1)=[?2, 0, 0.01, 0; 0, ?1, 0, 0], ?.sub.p(1)=[?1.3, 0, 0.2, 0.3; 0, ?0.05, 0.03, 0], ?.sub.1=0.5, ?.sub.2=0.5, ?.sub.1=1, ?.sub.2=1, ?.sub.1=0.5, ?.sub.2=0.9, ?.sub.1=0.2, ?.sub.2=0.2, ?=0.1.

[0090] When controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances, the following steps are included at each time k: a) obtaining the measurable disturbance vector d(k); b) obtaining the desired system output vector y*(k) and the actual system output vector y (k), calculating the system error vector e(k); c) based on steps a) and b), calculating the control input vector u(k) according to the partial-form model-free adaptive disturbance compensation control law in the presence of measurable disturbances with the optimal value of the partial-form adaptive input matrix ?.sub.p(k) and the optimal value of the partial-form adaptive disturbance matrix ?.sub.p(k); d) controlling the vapor compression refrigeration system by applying the control input vector u(k), generating the actual system output vector of the vapor compression refrigeration system at time k+1 based on the application of the control input vector; e) repeating steps a) to d) until the end of the control time.

[0091] The control performance of the control method of the present invention and the existing PID control method are given in FIGS. 10-13. FIG. 10 shows the tracking performance of the first system output when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; FIG. 11 shows the tracking performance of the second system output when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; FIG. 12 shows the first control input when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; and FIG. 13 shows the second control input when controlling the vapor compression refrigeration system in the second exemplary embodiment by using the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances and the PID control method; comparing the control performance indexes, the ITAE(e.sub.1) of the first system output controlled by the control method of the present invention is 120374, the ITAE(e.sub.2) of the second system output controlled by the control method of the present invention is 5609, the ITAE(e.sub.1) of the first system output controlled by the PID control method is 471529, the ITAE(e.sub.2) of the second system output controlled by the PID control method is 52531, the comparison results of the control performance indexes are listed in Table 2; comparing the system output curves, the control method of the present invention can effectively suppress the effect of measurable disturbances on actual system outputs of the vapor compression refrigeration system, and the control performance using the control method of the present invention is better than that using the PID control method. From the above comparison results, it can be seen that the partial-form model-free adaptive disturbance compensation control in the presence of measurable disturbances provides a significant improvement in disturbance compensation control performance and achieves effective tracking of the desired system outputs.

TABLE-US-00002 TABLE 2 Comparison Results of The Control Performance Indexes of VCRS ITAE(e.sub.1) ITAE(e.sub.2) PID 471529 52531 The present invention 120374 5609

[0092] Furthermore, the following two points should be noted in particular:

[0093] (1) Disturbances widely exist in actual control systems, such as controlled plants in oil refining, petrochemical, chemical, pharmaceutical, food, paper, water treatment, thermal power, metallurgy, cement, rubber, machinery, electrical, transportation, and robotics industries, including reactors, distillation towers, machines, equipment, devices, production lines, workshops, factories, unmanned vehicles, unmanned ships, unmanned aircraft, and autonomous mobile robots. By way of example, the vapor compression refrigeration system is affected by two persistent complex disturbances, i.e., the inlet temperature of the condenser secondary flux and the inlet temperature of the evaporator secondary flux. From the second exemplary embodiment, the control method of the present invention provides a significant improvement in disturbance compensation control performance and achieves effective tracking of the desired system outputs. To give another example, the unmanned ship in operation is extremely vulnerable to the impact of the wind field on the water, changes in wind speed and wind direction will not only affect the speed and heading of the unmanned ship, but may also lead to capsizing of the unmanned ship in serious cases; based on the wind speed sensor and wind direction sensor, the two measurable disturbances, i.e., wind speed and wind direction, can be monitored online, and the control method of the present invention can be used to compensate for the measurable disturbances and realize the smooth operation of the unmanned ship, which is of great significance to improve the safety and reliability of the unmanned ship.

[0094] (2) In the first and second exemplary embodiments, the hardware platform for running the embodiments of the present invention is the industrial control computer; in practical applications, according to the specific circumstance, a single chip microcomputer controller, a microprocessor controller, a field programmable gate array controller, a digital signal processing controller, an embedded system controller, a programmable logic controller, a distributed control system, a fieldbus control system, an industrial control system based on internet of things, or an industrial internet control system, can also be used as the hardware platform for running the control method of the present invention.

[0095] From the above detailed description of the invention, it is clear to those skilled in the art that the implementation of the present invention can be achieved with the help of software and the necessary hardware platform. Embodiments of the present invention can be implemented by using the existing processor, or by a dedicated processor being used for this or other purposes in an appropriate system, or by a hardwired system. Embodiments of the present invention also include a non-transitory computer-readable storage medium comprising a machine-readable medium for carrying or having machine-executable instructions or data structures stored thereon; the machine-readable medium can be any available medium accessible by a general purpose or the dedicated computer or other machines with a processor. By way of example, the machine-readable medium includes RAM, ROM, EPROM, EEPROM, CD-ROM or other optical disk memory, disk memory or other magnetic storage devices, or any other medium that can carry or store the required computer program code in the form of machine-executable instructions or data structures, and that can be accessed by a general purpose or the dedicated computer or other machines with a processor. When information is transmitted or made available to a machine over a network or other communication connection (hardwired, wireless, or a combination of hardwired and wireless), the connection is also considered a machine-readable medium.

[0096] It should be appreciated that the foregoing is only preferred embodiments of the invention and is not for use in limiting the invention. Any modification, equivalent substitution, and improvement without departing from the spirit and principle of this invention should be covered in the protection scope of the invention.