METHOD FOR TUNING PREDICTIVE CONTROL PARAMETERS OF BUILDING ENERGY CONSUMPTION SYSTEM BASED ON FUZZY LOGIC

Abstract

A method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic: 1) constructing a controlled building energy consumption system, performing generalized predictive control on the building energy consumption system, and initializing an tuned parameter λ of a generalized predictive controller; 2) collecting the output slope y.sub.k (t), the actual output y(t), the set value y.sub.r(t) and the predicted output value ŷ(t+i) of the controlled building energy consumption system in the control process, and then taking them as fuzzy target parameters; 3) constructing a membership function for the fuzzy target parameters in step 2), and then optimally selecting the parameters of the fuzzy membership function by using a particle swarm optimization algorithm to obtain membership function parameters of each fuzzy target parameter; 4) carrying out fuzzy reasoning operation on the membership function parameters, and tuning the adjusted parameter λ by using the results of fuzzy reasoning operation.

Claims

1. A method for tuning one or more predictive control parameters of a building energy consumption system based on fuzzy logic, comprising the steps of: constructing a controlled building energy consumption system, performing generalized predictive control on the controlled building energy consumption system, and initializing a tuned parameter λ of a generalized predictive controller; collecting the output slope y.sub.k(t), the actual output y(t), the set value y.sub.r(t) and the predicted output value ŷ(t+i) of the controlled building energy consumption system as one or more fuzzy target parameters; constructing a membership function for the one or more fuzzy target parameters, and optimally selecting the fuzzy target parameters of the fuzzy membership function by using a particle swarm optimization algorithm to obtain membership function parameters of a given one of the one or more fuzzy target parameters, thereby determining the membership function; carrying out a fuzzy reasoning operation on the membership function parameters of the given one of the one or more fuzzy target parameters; and tuning an adjusted parameter λ by using the results of the fuzzy reasoning operation, thus completing tuning predictive control parameters of the building energy consumption system based on fuzzy logic.

2. The method for tuning one or more predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, wherein the adjusted parameter λ is substituted into a cost function of a next cycle to improve the performance of the controlled building energy consumption system in the next cycle.

3. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, wherein constructing the controlled building energy consumption system comprises constructing a variable air volume air-conditioning system, wherein a dynamic model process transfer function of the controlled building energy consumption system is a first-order time-delay model, namely: $\frac{y\left(s\right)}{u\left(s\right)}=\frac{K}{1+\mathrm{Ts}}{e}^{-\tau s}$ where y is a static pressure of an air supply pipe; u is a fan input frequency; K is a process gain; T is a time constant; and τ is a delay time.

4. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, wherein the adjusted parameter λ is initialized as ${\left[\begin{array}{c}0\\ \mathrm{.Math.}\\ 0\end{array}\right]}_{\left(\mathrm{tend}+N+d\right)\times 1};$ where t.sub.end is a number of cycles of the predictive control process, N is a predictive time domain, and d is a pure delay in the controlled system model.

5. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, comprising: taking the output slope y.sub.k(t) of the controlled building energy consumption system as a first fuzzy target parameter: $y_k\left(t\right)=\{\begin{array}{cc}0,& t=1\\ \frac{y\left(t\right)-y\left(t-1\right)}{T},& t>1\end{array}$ y.sub.k(t) is constrained by y.sub.k.sub.min≤y.sub.k(t)≤y.sub.k.sub.max, where t is a sampling order number, T is a sampling time interval, y(t) is an actual output at a given sampling time in the predictive control process, and y(t−1) is an actual output at a previous time in the predictive control process; on the basis of the collected actual output y(t) and the set value y.sub.r(t), taking the absolute deviation e(t) between the actual output y(t) and the set value y.sub.r(t) a second fuzzy target parameter:
e(t)=y(t)−y.sub.r(t)(0≤e(t)≤e.sub.max) a time ts(t) when an output of the change rate of the absolute deviation e(t) reaches a set value is acquired taking ts(t) is taken as a third fuzzy target parameter: $\begin{array}{cc}\mathrm{ts}\left(t\right)=\{\begin{array}{cc}e\left(t\right)/.Math.\stackrel{.}{e}\left(t\right).Math.,& \stackrel{.}{e}\left(t\right)<0\\ 0,& \stackrel{.}{e}\left(t\right)=0,e\left(t\right)=0\\ M,& \stackrel{.}{e}\left(t\right)\ge 0,e\left(t\right)\ne 0\end{array}& \left(4\right)\end{array}$ where ts.sub.min≤ts(t)≤ts.sub.max, ė(t)=d(e(t))/dt, M is a constant; and on the basis of the collected predicted output value, taking ŷ(t+i), ŷ(t+i) as a fuzzy target parameter, which is constrained by:
ŷ.sub.min≤ŷ(t+i)≤ŷ.sub.max(i=1,2, . . . N) where i=1, 2, . . . , N.

6. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, wherein the membership function parameters of the given one of the one or more fuzzy target parameters are v.sub.min, v.sub.max, p.sub.1 and p.sub.2, where v.sub.min and v.sub.max are a minimum value and a maximum value of given one of the one or more fuzzy target parameters, and p.sub.1 and p.sub.2 are referred to as fuzzy width.

7. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 1, wherein a Mamdani fuzzy reasoning method is used to perform the fuzzy reasoning operation on membership.

8. The method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to claim 2, wherein the cost function is:
J−E{Σ.sub.j-N.sub.1.sup.N.sup.2(y(k+j)−y.sub.r(k+j)).sup.2+Σ.sub.j-1.sup.N.sup.uλ(j)(Δu(k+j−1)).sup.2}  (13) where N.sub.1 is a minimum prediction time domain length, when a system time delay d is known, N.sub.1=d, N.sub.2 is a maximum prediction time domain length, N.sub.u is a control time domain length, and λ(j) is a control weighting coefficient matrix greater than zero, λ(j)=λ.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0027] FIG. 1 is a flow chart of the present disclosure;

[0028] FIG. 2 is a flow chart of a particle swarm optimization algorithm;

[0029] FIG. 3 is a schematic diagram of a fuzzy membership parameter table according to the present disclosure;

[0030] FIG. 4 is a schematic diagram of a fuzzy membership function according to the present disclosure;

[0031] FIG. 5 is a comparison diagram of system simulation output curves according to the present disclosure;

[0032] FIG. 6 is a graph showing the variation of the adjusted parameters according to the present disclosure.

DETAILED DESCRIPTION

[0033] The present disclosure will be described in further detail with reference to the accompanying drawings.

[0034] Referring to FIG. 1 and FIG. 2, the method for tuning predictive control parameters of a building energy consumption system based on fuzzy logic according to the present disclosure comprises the following steps.

[0035] 1) A controlled building energy consumption system is constructed, generalized predictive control is performed on the building energy consumption system, and an adjusted parameter λ of a generalized predictive controller is initialized.

[0036] In step 1), the controlled building energy consumption system is a variable air volume air-conditioning system, and the dynamic model process transfer function of the controlled building energy consumption system is a first-order time-delay model, namely:

[00005]$\begin{array}{cc}\frac{y\left(s\right)}{u\left(s\right)}=\frac{K}{1+\mathrm{Ts}}{e}^{-\tau s}& \left(1\right)\end{array}$

where y is the static pressure of the air supply pipe; u is the fan input frequency; K is the process gain; T is a time constant; τ is the delay time.

[0037] In step 1), the adjusted parameter λ of a generalized predictive controller is initialized as

[00006]${\left[\begin{array}{c}0\\ \mathrm{.Math.}\\ 0\end{array}\right]}_{\left(\mathrm{tend}+N+d\right)\times 1};$

where tend is the number of cycles of the predictive control process, N is the predictive time domain, and d is the pure delay in the controlled system model.

[0038] 2) The output slope y.sub.k(t), the actual output y(t), the set value y.sub.r(t) and the predicted output value ŷ(t+i) of the controlled building energy consumption system in the control process are collected, and then the output slope y.sub.k(t), the actual output y(t), the set value y.sub.r(t) and the predicted output value ŷ(t+1) of the controlled building energy consumption system are taken as fuzzy target parameters.

[0039] In step 2),

[0040] on the basis of the collected output slope y.sub.k(t) of the controlled building energy consumption system, the output slope y.sub.k(t) of the controlled building energy consumption system is taken as a first fuzzy target parameter:

[00007]$\begin{array}{cc}{y}_k\left(t\right)=\{\begin{array}{cc}0,& t=1\\ \frac{y\left(t\right)-y\left(t-1\right)}{T},& t>1\end{array}& \left(2\right)\end{array}$

[0041] y.sub.k(t) is constrained by y.sub.k.sub.min≤y.sub.k(t)≤y.sub.k.sub.max, where t is the sampling order number, T is the sampling time interval, y(t) is the actual output at the sampling time in the predictive control process, and y(t−1) is the actual output at the previous time in the predictive control process.

[0042] On the basis of the collected actual output y(t) and the set value y.sub.r(t), the absolute deviation e(t) between the actual output y(t) and the set value y.sub.r(t) is taken as a second fuzzy target parameter:

e(t)=y(t)−y.sub.r(t)(0≤e(t)≤e.sub.max)  (3)

[0043] The time ts(t) when the output of the change rate of the absolute deviation e(t) reaches the set value is acquired, and ts(t) is taken as a third fuzzy target parameter:

[00008]$\begin{array}{cc}\mathrm{ts}\left(t\right)=\{\begin{array}{cc}e\left(t\right)/.Math.\stackrel{.}{e}\left(t\right).Math.,& \stackrel{.}{e}\left(t\right)<0\\ 0,& \stackrel{.}{e}\left(t\right)=0,e\left(t\right)=0\\ M,& \stackrel{.}{e}\left(t\right)\ge 0,e\left(t\right)\ne 0\end{array}& \left(4\right)\end{array}$

where ts.sub.min≤ts(t)≤ts.sub.max, ė(t)=d(e(t))/dt, M is a constant.

[0044] On the basis of the collected predicted output value ŷ(t+i), ŷ(t+i) is taken as a fuzzy target parameter, which is constrained by:

ŷ.sub.min≤ŷ(t+i)≤ŷ.sub.max(i=1,2, . . . ,N)

where i=1, 2, . . . , N.

[0045] 3) A membership function is constructed for the fuzzy target parameters in step 2), and then the parameters of the fuzzy membership function are optimally selected by using a particle swarm optimization algorithm to obtain membership function parameters of each fuzzy target parameter, thereby determining the membership function.

[0046] In step 3), the membership function parameters of each fuzzy target parameter are v.sub.min, v.sub.max, p.sub.1 and p.sub.2, where v.sub.min and v.sub.max are the minimum value and the maximum value of fuzzy target parameters, and p.sub.1 and p.sub.2 are referred to as fuzzy width.

[0047] Referring to FIG. 3, the specific process of step 3) is as follows.

[0048] 31) The number of particles in the particle swarm is set, the speed and the position of all particles are initialized, and the maximum speed interval is set. The position information of each particle comprises 16 membership function parameters, namely:

[00009]$\hspace{1em}\left[\begin{array}{cccc}{v}_{1\min }& v_{1\max }& p_{11}& p_{21}\\ v_{2\min }& v_{2\max }& p_{12}& p_{22}\\ v_{3\min }& v_{3\max }& p_{13}& p_{23}\\ v_{4\min }& v_{4\max }& p_{14}& p_{24}\end{array}\right]$

[0049] 32) The fitness function of each particle is calculated, the current individual extremum of each particle is found, and a global optimal solution is found from these individual historical optimal solutions and is compared with the historical optimal solution to select the optimal solution as the current historical optimal solution.

[0050] 33) The speed and position information of each particle is updated, and the update formula is:

V.sub.id=ωV.sub.id+C.sub.1random(0,1)(P.sub.id−X.sub.id)+C.sub.2random(0,1)(P.sub.gd−X.sub.id)  (5)

X.sub.id=X.sub.id+V.sub.id  (6)

where ω is an inertia factor, ω is non-negative, V.sub.id is the particle velocity, X.sub.id is the current position of a particle, a four-dimensional matrix consists of four groups of membership function parameters, P.sub.id is the historical optimal position of particles, P.sub.gd is the global optimal position of a swarm, C.sub.1 and C.sub.2 are learning factors which are usually 0-4, and random(0,1) is a random number in the interval [0,1].

[0051] 34) It is detected whether the updated particles meet the condition of ending the cycle, if not, continuing the cycle, if so, outputting the optimal solution as the parameter of the fuzzy membership function, and establishing the membership function to obtain the required membership. The membership functions established by the present disclosure are as follows:

[00010]$\begin{array}{cc}{\mu }_{y_k}=\{\begin{array}{cc}0& y_k<y_{k_{\min }}-p_1\\ 1+\frac{y_k\left(t\right)-y_{k_{\min }}}{p_1},& y_{k_{\min }}-p_1\le y_k\left(t\right)<y_{k_{\min }}\\ 1,& y_{k_{\min }}\le y_k\left(t\right)<y_{k_{\max }}\\ 1-\frac{y_k\left(t\right)-y_{k_{\min }}}{p_2},& y_{k_{\max }}\le y_k\left(t\right)<y_{k_{\max }}+p_2\\ 0,& y_{k_{\max }}+p_2\le y_k\left(t\right)\end{array}& \left(7\right)\\ {\mu }_e={\mu }_{\tilde{e}\left(t\right)}\left(e\left(t\right)\right)=\{\begin{array}{cc}1,& 0\le e\left(t\right)<e_{\max }\\ 1-\frac{e\left(t\right)-e_{\max }}{p_2},& e_{\max }\le e\left(t\right)<e_{\max }+p_2\\ 0,& e\left(t\right)\ge e_{\max }+p_2\end{array}& \left(8\right)\\ {\mu }_{\mathrm{ts}}=\{\begin{array}{cc}\frac{{\mathrm{ts}}_{\min }-\mathrm{ts}\left(t\right)}{p_1},& {\mathrm{ts}}_{\min }-p_1\le \mathrm{ts}\left(t\right)<{\mathrm{ts}}_{\min },\stackrel{.}{e}\left(t\right)<0\\ 0,& {\mathrm{ts}}_{\min }\le \mathrm{ts}\left(t\right)<{\mathrm{ts}}_{\max },\stackrel{.}{e}\left(t\right)<0\\ \frac{\mathrm{ts}\left(t\right)-{\mathrm{ts}}_{\max }}{p_2},& {\mathrm{ts}}_{\max }\le \mathrm{ts}\left(t\right)<{\mathrm{ts}}_{\max }+p_2,\stackrel{.}{e}\left(t\right)<0\\ 1,& \mathrm{ts}\left(s\right)<{\mathrm{ts}}_{\min }-p_1,\stackrel{.}{e}\left(i\right)<0\\ 1,& \mathrm{ts}\left(t\right)\ge {\mathrm{ts}}_{\max }+p_s,\stackrel{.}{e}\left(t\right)<0\\ 0,& \stackrel{.}{e}\left(t\right)=0,e\left(t\right)=0\\ 1,& \stackrel{.}{e}\left(t\right)\ge 0,e\left(t\right)\ne 0\end{array}& \left(9\right)\\ {\mu }_{\hat{y}}\left(i\right)={\mu }_{\hat{y}\left(t+i\right)}\left(\hat{y}\left(t+i\right)\right)=\{\begin{array}{cc}0,& \hat{y}\left(t+i\right)<{\hat{y}}_{\min }-p_1\\ 1+\frac{\hat{y}\left(t+i\right)-{\hat{y}}_{\min }}{p_1},& {\hat{y}}_{\min }-p_1\le \hat{y}\left(t+i\right)<{\hat{y}}_{\min }\\ 1,& {\hat{y}}_{\min }\le \hat{y}\left(t+i\right)<{\hat{y}}_{\max }\\ 1-\frac{\hat{y}\left(t+i\right)-{\hat{y}}_{\max }}{p_2},& {\hat{y}}_{\max }\le \hat{y}\left(t+i\right)<{\hat{y}}_{\max }+p_2\\ 0,& {\hat{y}}_{\max }+p_2\le \hat{y}\left(t+i\right)\end{array}& \left(10\right)\end{array}$

where i=1, 2, . . . , N.

[0052] 4) Fuzzy reasoning operation is carried out on the membership, and the adjusted parameter λ is tuned by using the results of fuzzy reasoning operation, thus completing tuning predictive control parameters of the building energy consumption system based on fuzzy logic.

[0053] In step 4), Mamdani fuzzy reasoning method is used to perform fuzzy reasoning operation on membership, which is specifically as follows.

[0054] 41) The most commonly used Mamdani fuzzy reasoning method is used, and the obtained membership is subject to Cartesian product operation, namely:

μ.sub.min=μ.sub.yk∧μ.sub.e∧μ.sub.ts∧(min{μ.sub.ŷ(1),μ.sub.ŷ(2), . . . ,μ.sub.ŷ(N)})  (11)

where the index y.sub.k(t) is the slope of the output curve, which is directly related to the increase and decrease of the output curve at each moment. Taking this index as a fuzzy target parameter to adjust the system output can directly influence the trend of the output curve. The index ts(t) contains the prediction amount of the rise time of the system. The indexes e(t) and ŷ(t+i) describe the difference between the current and predicted system output and the set value.

[0055] 42) According to the μ.sub.min value obtained in step 31), λ is tuned between λ.sub.min and λ.sub.max according to a certain exponential law, and the weight of the control quantity is changed, whose algebraic expression is:

λ=λ.sub.max×exp(μ.sub.min×Ig(λ.sub.min/λ.sub.max))  (12)

where λ.sub.max>λ.sub.min>0.

[0056] 5) The tuned adjusted parameter λ is substituted into the cost function of the next cycle to improve the performance of the controlled building energy consumption system in the next cycle, wherein the cost function is:

J−E{Σ.sub.j-N.sub.1.sup.N.sup.2(y(k+j)−y.sub.r(k+j)).sup.2+Σ.sub.j-1.sup.N.sup.uλ(j)(Δu(k+j−1)).sup.2}  (13)

where N.sub.1 is the minimum prediction time domain length, when the system time delay d is known, N.sub.1=d, N.sub.2 is the maximum prediction time domain length, N.sub.u is the control time domain length, λ(j) is the control weighting coefficient matrix greater than zero, λ(j)=λ.

Embodiment 1

[0057] In the variable air volume air-conditioning system, the static pressure of an air supply pipe in the static pressure control loop of central air-conditioning supply is taken as the controlled object, and the fan frequency-static pressure model is obtained by system identification with the objective that the output of the static pressure prediction model can quickly and accurately follow the set value of static pressure, as shown in the following formula:

[00011]$\begin{array}{cc}G\left(s\right)=\frac{4.193e^{-4.87s}}{6.219s+1}& \left(10\right)\end{array}$

[0058] The method specifically comprises the following steps.

[0059] 1) The building energy consumption system is modeled, generalized predictive control is performed, and an tuned parameter λ is initialized.

[0060] The tuned parameter λ of a generalized predictive controller is initialized as

[00012]${\left[\begin{array}{c}0\\ \mathrm{.Math.}\\ 0\end{array}\right]}_{\left(\mathrm{tend}+N+d\right)\times 1};$

where tend is 200, N is the predictive time domain which is 5, and d is the pure delay in the controlled system model, which is 4.

[0061] 2) In the process of collecting control, the system includes the output slope y.sub.k(t), the actual output y(t), the set value y.sub.r(t) and the predicted output value ŷ(t+i). On the basis of the collected output slope y.sub.k(t) of the controlled system, the output slope is taken as a first fuzzy target parameter and is constrained by 0≤y.sub.k(t)≤0.4:

[00013]$\begin{array}{cc}{y}_k\left(t\right)=\{\begin{array}{cc}0,& t=1\\ \frac{y\left(t\right)-y\left(t-1\right)}{T},& t>1\end{array}& \left(11\right)\end{array}$

where t is the sampling order number which is 200 times in total, T is the sampling time interval which is 1 s, y(t) is the actual output at the sampling time in the predictive control process, and y(t−1) is the actual output at the previous time in the predictive control process.

[0062] On the basis of the collected actual output y(t) and the set value y.sub.r(t), the absolute deviation between the actual output y(t) and the set value y.sub.r(t) is taken as a fuzzy target parameter:

e(t)=y(t)−y.sub.r(t)(0≤e(t)≤0.3)  (12)

[0063] The time ts(t) when the output based on the change rate of the absolute deviation e(t) reaches the set value is acquired:

[00014]$\begin{array}{cc}\mathrm{ts}\left(t\right)=\{\begin{array}{cc}e\left(t\right)/.Math.\stackrel{.}{e}\left(t\right).Math.,& \stackrel{.}{e}\left(t\right)<0\\ 0,& \stackrel{.}{e}\left(t\right)=0,e\left(t\right)=0\\ M,& \stackrel{.}{e}\left(t\right)\ge 0,e\left(t\right)\ne 0\end{array}& \left(13\right)\end{array}$

where 0.15≤ts(t)≤0.75, ė(t)=d(e(t))/dt, M is a big constant.

[0064] On the basis of the collected predicted output value ŷ(t+i), ŷ(t+i) is taken as a fuzzy target parameter, which is constrained by:

0.3≤y(t+i)≤1(i=1,2, . . . ,N)

where i=1, 2, . . . , N.

[0065] 3) A membership function is constructed for the four collected fuzzy target parameters. In the present disclosure, the parameters of the fuzzy membership function are optimally selected by using a particle swarm optimization algorithm, so that each fuzzy target parameter can obtain four membership function parameters, namely, v.sub.min, v.sub.max, p.sub.1, p.sub.2, where v.sub.min and v.sub.max are the minimum value and the maximum value of fuzzy target parameters, and p.sub.1 and p.sub.2 are referred to as fuzzy width.

[0066] The specific operation process of step 3) is as follows.

[0067] 31) The number of particles in the particle swarm is set as 50, the speed and the position of all particles are initialized, and the maximum speed interval is set. The position information of each particle comprises 16 membership function parameters, namely:

[00015]$\hspace{1em}\left[\begin{array}{cccc}{v}_{1\min }& v_{1\max }& p_{11}& p_{21}\\ v_{2\min }& v_{2\max }& p_{12}& p_{22}\\ v_{3\min }& v_{3\max }& p_{13}& p_{23}\\ v_{4\min }& v_{4\max }& p_{14}& p_{24}\end{array}\right]$

[0068] 32) The fitness function of each particle is calculated, the current individual extremum of each particle is found, and a global optimal solution is found from these individual historical optimal solutions and is compared with the historical optimal solution to select the optimal solution as the current historical optimal solution.

[0069] 33) The speed and position information of each particle is updated, and the update formula is:

V.sub.id=ωV.sub.id+C.sub.1random(0,1)(P.sub.id−X.sub.id)+C.sub.2random(0,1)(P.sub.gd−X.sub.id)  (14)

X.sub.id=X.sub.id+V.sub.id  (15)

where ω is an inertia factor, ω has a value of 0.85, V.sub.id is the particle velocity, X.sub.id is the current position of a particle, a four-dimensional matrix consists of four groups of membership function parameters, P.sub.id is the historical optimal position of particles, P.sub.gd is the global optimal position of a swarm, C.sub.1 and C.sub.2 are learning factors which are usually 0.5, and random(0,1) is a random number in the interval [0,1].

[0070] 34) It is detected whether the updated particles meet the condition of ending the cycle, if not, continuing the cycle, if so, outputting the optimal solution as the parameter of the fuzzy membership function to obtain the parameter as shown in FIG. 3, and establishing the membership function as shown in FIG. 4 to obtain the required membership. The membership functions established by the present disclosure are as follows:

[00016]$\begin{array}{cc}{\mu }_{y_k}=\{\begin{array}{cc}0,& y_k<-0.6\\ 1+\frac{y_k\left(t\right)-0.4}{1},& -0.6\le y_k\left(t\right)<0.4\\ 0,& 0.4<y_k\left(t\right)\end{array}& \left(16\right)\\ {\mu }_e={\mu }_{e\left(t\right)}\left(e\left(t\right)\right)=\{\begin{array}{cc}1,& 0\le e\left(t\right)<0.1\\ 1-\frac{e\left(t\right)-0.3}{0.75},& 0.3\le e\left(t\right)<1.05\\ 0,& e\left(t\right)\ge 1.05\end{array}& \left(17\right)\\ {\mu }_{\mathrm{ts}}=\{\begin{array}{cc}\frac{0.15-\mathrm{ts}\left(t\right)}{0.1},& 0.05\le \mathrm{ts}\left(t\right)<0.15,\stackrel{.}{e}\left(t\right)<0\\ 0,& 0.15\le \mathrm{ts}\left(t\right)<0.75,\stackrel{.}{e}\left(t\right)<0\\ \frac{\mathrm{ts}\left(t\right)-0.75}{0.1},& 0.75\le \mathrm{ts}\left(t\right)<0.85,\stackrel{.}{e}\left(t\right)<0\\ 1,& \mathrm{ts}\left(s\right)<0.05,\stackrel{.}{e}\left(t\right)<0\\ 1,& \mathrm{ts}\left(t\right)\ge 0.85,\stackrel{.}{e}\left(t\right)<0\\ 0,& \stackrel{.}{e}\left(t\right)=0,e\left(t\right)=0\\ 1,& \stackrel{.}{e}\left(t\right)\ge 0,e\left(t\right)\ne 0\end{array}& \left(18\right)\\ u_{\hat{y}}\left(i\right)={\mu }_{\hat{y}\left(t+i\right)}\left(\hat{y}\left(t+i\right)\right)=\{\begin{array}{cc}0,& \hat{y}\left(t+i\right)<0.03\\ 1,& 0.3\le \hat{y}\left(t+i\right)<1\\ 1-\frac{\hat{y}\left(t+i\right)-1}{0.3},& 1\le \hat{y}\left(t+i\right)<1.3\\ 0,& 1.3\le \hat{y}\left(t+i\right)\end{array}& \left(19\right)\end{array}$

where i−1, 2, . . . , N.

[0071] 4) Fuzzy reasoning operation is carried out on the membership of fuzzy target parameters, and the adjusted parameter λ is tuned by using the results of fuzzy reasoning operation.

[0072] The specific process of step 4) is as follows.

[0073] 41) The Mamdani fuzzy reasoning method is used, and the obtained membership is subject to Cartesian product operation, namely:

μ.sub.min=μ.sub.yk∧μ.sub.e∧μ.sub.ts∧(min{μ.sub.ŷ(1),μ.sub.ŷ(2), . . . ,μ.sub.ŷ(N)})  (20)

where, y.sub.k(t) is the slope of the output curve, which is directly related to the increase and decrease of the output curve at each moment. Taking this index as a fuzzy target parameter to tune the system output can directly influence the trend of the output curve. The index ts(t) contains the prediction amount of the rise time of the system. The indexes e(t): and ŷ(t+i) describe the difference between the current and predicted system output and the set value.

[0074] 42) According to the μ.sub.min value, λ is tuned between λ.sub.min and λ.sub.max according to a certain exponential law, and the weight of the control quantity is changed, whose algebraic expression is:

λ=λ.sub.max×exp(μ.sub.min×Ig(λ.sub.min/λ.sub.max))  (21)

where λ.sub.max=3; λ.sub.min=0.001.

[0075] 5) The tuned weighting coefficient λ is substituted into the cost function of the next cycle again, wherein the cost function is:

J−E{Σ.sub.j-N.sub.1.sup.N.sup.2(y(k+j)−y.sub.r(k+j)).sup.2+Σ.sub.j-1.sup.N.sup.uλ(j)(Δu(k+j−1)).sup.2}  (22)

where N.sub.1 is the minimum prediction time domain length, which is selected as 1, when the system time delay d is known, N.sub.1=d=4, N.sub.2 is the maximum prediction time domain length, which is 5, N.sub.u is the control time domain length, which is 2.

[0076] According to the present disclosure, the variable air volume air-conditioning system is predicted, controlled and simulated, and the result is shown in FIG. 5. It can be seen that the rapidity and stability of the variable air volume air-conditioning system are obviously improved, and the result shows the effectiveness of the present disclosure.

[0077] In addition, as shown in FIG. 6, compared with FIG. 5, it can be seen that with the change of λ value from large to small, the response speed of the corresponding system is accelerated, and the control function is enhanced, so that the system can reach the set stable state faster, which proves reversely that the weighting coefficient λ of the controller has a direct influence on the control process.

[0078] Compared with the prior art, the present disclosure has the following advantages.

[0079] The underutilized output slope, predicted output and other information containing many system characteristics in the generalized predictive control process are applied to the fuzzy logic algorithm to tune the control system parameters, which improves the utilization rate of the control system.

[0080] By constructing the fuzzy membership function, the constraints on performance indexes such as rise time and overshoot in the system are transformed into constraints on four fuzzy control targets, which greatly reduces the calculation amount in the process of tuning parameters.

[0081] Only the weighting coefficient λ of the control quantity is taken as an adjustable parameter, and other control parameters such as a flexibility coefficient and a control time domain are fixed quantities, so that the contradiction between rapidity and stability of the control system can be solved most effectively.

[0082] By applying the system output slope at each moment and the time when the output based on the current absolute error change rate reaches the set value to the parameter tuning, the adjustment time of the system can be effectively shortened.

[0083] The absolute deviation between the actual output and the set output and the predicted output are applied to the fuzzy logic algorithm to tune the system parameters, so that the overshoot of the system is smaller and the robustness is stronger.

[0084] By using the particle swarm optimization algorithm, the parameters of the fuzzy membership function are found more accurately, so that the present disclosure is more universal.