CN102353119A - Control method of VAV (variable air volume) air-conditioning system - Google Patents

Abstract

The invention provides a control scheme of a VAV (variable air volume) air-conditioning system. A neural network predictive control method is used for a tail end VAV-BOX and an air-conditioning unit, thus the hysteresis characteristic of a VAV system can be overcome, the control accuracy is improved, the resonance phenomena of an actuating mechanism can be greatly reduced, the energy-saving effect is improved by above 13%, and the control parameter tuning problem in the project is solved. A pressure independent cascade stage predictive control method is used for the tail end VAV-BOX, thus the control accuracy can also be improved. The air-conditioning unit is provided with four control loops and can automatically select a static pressure control or total air volume control policy with adjustable setting static pressure by using an all-condition integrated control technology, thus the predictive control on a fan can be realized; the primary air volume of all the tail end VAV-BOXes can be detected and the air supply temperature can be adjusted according to the operation condition, thus the problem of too low temperature in partial air-conditioning area in the lowest fresh air operation is solved; the cascade stage predictive control method is adopted for a fresh air ratio control loop, thus the accurate control on the fresh air ratio can be realized and the energy-saving level can be further improved.

Description

A kind of VAV air conditioning system with variable control method

Technical field

The present invention relates to a kind of VAV air conditioning system with variable control method, belong to civil buildings VAV air quantity variable air conditioner control technology field.

Background technology

The system of VAV air quantity variable air conditioner control at present mainly is based on the independent design in a plurality of loops; And all be to regulate basically with the PID control method; Its robot control system(RCS) needs the commissioning engineer according to the on-the-spot pid parameter of setting of the experience of self, does not possess self-learning capability.Because the time constant of each robot control system(RCS) is different, very easily causes the robot control system(RCS) adaptivity poor, make temperature fluctuation big or produce " resonance " phenomenon of actuator.Simultaneously; Owing to have certain coupling relation between a plurality of loops in the VAV air conditioning system with variable, even debugging of single loop and operation are all out of question, in all loops during co-ordination; Whole system also is difficult for realizing stable control, is prone to " resonance " phenomenon of system.

In addition; The setting of ratio, integration, differential parameter has very big influence to the performance of control system in the PID control method; Because air conditioning system with variable is formed complicacy, equipment is numerous; Various application scenario is also different to the requirement of parameter; Need the commissioning engineer to set according to the experience scene of self, debugging has brought very big difficulty to engineering site.Simultaneously since exist summer, winter and transition season three kinds of different working conditions conditions, the VAV air-conditioning system often needs could satisfy basically more than 1 year client's control performance requirement at least in practical engineering application.

There is minority Design of Variable Air Volume System company to adopt fuzzy PID control method at present,, improves the adaptive ability of control system so that can in control procedure, adjust pid parameter automatically according to environmental change.But, because the PID control method can't realize optimum control in essence, promptly can't make certain performance indications reach optimum, therefore do not realize the purpose of further saving energy consumption.

In sum, at present the VAV air conditioning system with variable exists that control performance is relatively poor, debugging work load big and the energy-saving effect problem of good aspect not, and these problems affect are to the application of VAV air conditioning system with variable.

Summary of the invention

The present invention is directed to the existing this difficult point of multiple-input and multiple-output nonlinear system aspect control that has the large time delay characteristic of VAV air quantity variable air conditioner, the advantage of comprehensive neutral net, optimum control and PREDICTIVE CONTROL has proposed a kind of intelligence control method.This method synthesis Hamilton-Jacobi-Bellman (HJB) and Eular-Lagrange (EL) optimized Algorithm; Utilize predicted roll optimizing idea training multilayer feedforward neural network; Become the optimization feedback solution of multiple-input and multiple-output nonlinear system when then it being found the solution as the optimization feedback controller; Can utilize the multi-step prediction rolling optimization to overcome various uncertainties and the complicated influence that changes simultaneously in amount of calculation and the optimal control problem that takies solution nonlinear system under the moderate situation of memory block capacity.

A kind of VAV air conditioning system with variable control method comprises following step:

The first step: utilize the BP neutral net to set up air conditioning area temperature neural network prediction model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model

1) confirms air conditioning area temperature prediction model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model structure

The input signal of air conditioning area temperature prediction model is outdoor intensity of solar radiation, outdoor temperature, CO2 concentration, indoor temperature, air quantity and terminal valve area, is output as next indoor temperature constantly;

The input signal of terminal air-valve forecast model is terminal valve area and duct static pressure (air-conditioning unit main air duct place), is output as next terminal air quantity constantly;

The input signal of main air duct duct static pressure forecast model is rotation speed of fan, duct static pressure, return air CO2 concentration, outdoor temperature and intensity of solar radiation, is output as next duct static pressure constantly;

The input signal of main air duct air quantity forecast model is rotation speed of fan and duct static pressure, is output as next main air duct air quantity constantly;

The input signal of wind pushing temperature forecast model is wind pushing temperature and water valve aperture, is output as next wind pushing temperature constantly;

The input signal of new wind air-valve forecast model is new wind valve area and duct static pressure (place, fresh wind tube road), is output as next resh air requirement constantly;

The input signal of air quality forecast model is new wind valve area and CO2 concentration, is output as next CO2 concentration constantly;

2) gather sample data;

3) sample data is by formula carried out normalization in (1), (2):

$x_i=\frac{x_{\mathrm{di}}-x_{d\min }}{x_{d\max }-x_{d\min }}---\left(1\right)$

$y_{\mathrm{tl}}=\frac{y_{\mathrm{dl}}-y_{d\min }}{y_{d\max }-y_{d\min }}---\left(2\right)$

Wherein, x _iInput value for neutral net after the normalization; x _DiBe former input value; x _DminBe the minimum of a value in the former input value; x _DmaxBe the maximum in the former input value; y _TlDesired value for neutral net after the normalization; y _DlRepresent former desired value; y _DminRepresent the minimum of a value in the former desired value; y _DmaxBe the maximum in the former desired value;

4) each neural network prediction model is carried out off-line training.

Second step: VAV air conditioning terminal tandem PREDICTIVE CONTROL

1) confirms terminal outer shroud PREDICTIVE CONTROL object function

Terminal outer shroud PREDICTIVE CONTROL object function is: $J_o\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_c-1}{\mathrm{\&Sigma;}}}{L}_o\left(T_o\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_c-1}{\mathrm{\&Sigma;}}}{\left(T_o\right[k]-T_{o,\mathrm{set}}[k\left]\right)}^2$

M wherein _cBe prediction time domain, t ₁Be the initial time in the prediction time domain, T _o[k] is the air conditioning area temperature in k sampling period, T _{O, set}[k] is the air conditioning area desired temperature in k sampling period, L _oIt is terminal outer shroud object function of k sampling period;

2) confirm terminal interior ring PREDICTIVE CONTROL object function

Ring PREDICTIVE CONTROL object function is in terminal: $J_i\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_i-1}{\mathrm{\&Sigma;}}}{L}_i\left(V_i\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_i-1}{\mathrm{\&Sigma;}}}{\left(V_i\right[k]-V_{i,\mathrm{set}}[k\left]\right)}^2$

V wherein _i[k] is the air quantity in k sampling period, V _{I, set}[k] is the air quantity setting value in k sampling period; M _iBe prediction time domain, L _iIt is ring object function in k the sampling period end;

3) confirm terminal interior ring neural network prediction controller and terminal outer shroud neural network prediction controller architecture

As input, terminal valve area is output to ring neural network prediction controller with air quantity setting value, duct static pressure (air-conditioning unit main air duct place) and pipeline air quantity in terminal;

The input parameter of terminal outer shroud neural network prediction controller comprises outdoor temperature, intensity of solar radiation, indoor temperature and air conditioning area desired temperature; Output parameter is the air quantity setting value;

4) ring neural network prediction controller in terminal and terminal outer shroud neural network prediction controller are carried out online optimizing training.

The 3rd step: VAV air-conditioning unit PREDICTIVE CONTROL

1) confirms static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and the VAV air-conditioning unit new wind ratio control loop PREDICTIVE CONTROL object function of VAV air-conditioning unit adjustable settings static pressure

Static pressure control loop PREDICTIVE CONTROL object function is:

$J_s\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_s-1}{\mathrm{\&Sigma;}}}{L}_s\left(P_s\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_s-1}{\mathrm{\&Sigma;}}}\{{\left(P_s\right[k]-P_{s,\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

P wherein _s[k] is the duct static pressure in k sampling period, P _{S, set}[k] is that the air conditioning area in k sampling period is set static pressure, M _sBe prediction time domain, U _Fan[k] is k sampling period blower voltage controlled quentity controlled variable, L _sBe k sampling period static pressure control loop object function;

Total blast volume control loop PREDICTIVE CONTROL object function is:

$J_f\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_f-1}{\mathrm{\&Sigma;}}}{L}_f\left(V_{\mathrm{fan}}\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_f-1}{\mathrm{\&Sigma;}}}\{{\left(V_{\mathrm{fan}}\right[k]-V_{\mathrm{fan}.\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

V wherein _Fan[k] is the pipeline air quantity in k sampling period, M _fBe prediction time domain, V _Fan.set[k] is each terminal prediction air quantity sum in k sampling period, L _fBe k sampling period total blast volume control loop object function;

Wind pushing temperature control loop PREDICTIVE CONTROL object function is:

$J_{\mathrm{st}}\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{st}}\left(T_{\mathrm{st}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{\left(T_{\mathrm{st}}\right[k]-T_{\mathrm{st},\mathrm{set}}[k\left]\right)}^2$

T wherein _St[k] is the temperature in k sampling period, T _{St, set}[k] is the design temperature in k sampling period, M _StBe prediction time domain, L _StBe k sampling period wind pushing temperature control loop object function;

New wind ratio control loop outer shroud PREDICTIVE CONTROL object function is:

$J_q\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_q-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{no}}\left(Q\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_q-1}{\mathrm{\&Sigma;}}}{\left(Q\right[k]-Q_{\mathrm{set}}[k\left]\right)}^2$

Q[k wherein] be the air quality in k sampling period, Q _Set[k] is the air quality setting value in k sampling period, M _qBe prediction time domain, L _NoBe k sampling period new wind ratio control loop outer shroud object function;

Ring PREDICTIVE CONTROL object function is in the new wind ratio control loop:

$J_{\mathrm{ni}}\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{ni}}\left(V_{\mathrm{ni}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{\left(V_{\mathrm{ni}}\right[k]-S[k\left]\right)}^2$

V wherein _i[k] is the prediction air-valve discharge quantity of fan in k sampling period, S[k] be the setting air quantity in k sampling period, M _iBe prediction time domain, L _NiIt is ring object function in k the sampling period new wind ratio control loop;

2) confirm static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and the VAV air-conditioning unit new wind ratio control loop neural network prediction controller architecture of VAV air-conditioning unit adjustable settings static pressure

The input signal of static pressure neural network prediction controller is output as rotation speed of fan for setting duct static pressure and duct static pressure;

Total blast volume neural network prediction controller input signal is pipeline air quantity and total blast volume, is output as rotation speed of fan;

The input signal of wind pushing temperature neural network prediction controller is output as the water valve aperture for setting wind pushing temperature and wind pushing temperature;

New wind ratio outer shroud neural network prediction controller is input as sets CO2 concentration and return air CO2 concentration, is output as the setting resh air requirement; Ring neural network prediction controller will be set resh air requirement, duct static pressure (place, fresh wind tube road) and pipeline air quantity as input in the new wind ratio, and new air valve aperture is output;

3) realize the online optimizing of static pressure neural network prediction controller in the static pressure control loop;

4) realize the online optimizing of total blast volume neural network prediction controller in the total blast volume control loop;

5) realize the online optimizing of wind pushing temperature neural network prediction controller in the wind pushing temperature control loop;

6) realize ring and the online optimizing of outer shroud neural network prediction controller in the new wind ratio in the new wind ratio control loop.

The invention has the advantages that:

(1) the present invention utilizes neutral net to set up forecast model, and model can onlinely be revised, and disturbance is had the good adaptive ability; Simultaneously, through the Neural Network Self-learning function, can solve the engineering middle controller on-site parameters difficult problem of adjusting;

(2) the present invention adopts the PREDICTIVE CONTROL scheme, can solve the bigger problem of actual environment temperature fluctuation that air-conditioning system produces owing to hysteresis characteristic, and can realize optimization control, can further reduce the energy consumption of VAV air-conditioning system;

(3) because the control method that the present invention proposes is taken all factors into consideration comfort index and energy consumption index is optimized PREDICTIVE CONTROL as optimizing performance indications; Compare with the PID control method; Under the suitable situation of selected performance indications parameter, can make VAV air quantity variable air conditioner blower fan system energy efficient more than 13%;

(4) the pressure independent type that the present invention is directed to VAV air conditioning system with variable end equipment is controlled process characteristic, adopts the tandem forecast Control Algorithm based on neutral net, has improved control accuracy;

(5) the present invention is according to the operating mode characteristics in winter, summer and transition season, and system selects the static pressure control strategy or the total blast volume control strategy of adjustable settings static pressure automatically, realizes the accurate control to the blower fan total blast volume;

(6) the present invention has adopted new wind pushing temperature control strategy, and through detecting all terminal VAV-BOX primary air flows, the operating condition adjustment wind pushing temperature based on VAV-BOX has solved because minimum new wind moves, and caused the low excessively problem of part air conditioning area temperature;

(7) new wind ratio control loop of the present invention adopts the tandem forecast Control Algorithm, realizes the accurate control to new wind ratio, thereby has further improved the energy-saving effect of system;

(8) the neural network prediction system optimizing control that adopts of the present invention can the non-linear and probabilistic influence of resolution system, and the algorithm real-time is good, and it is few to take the memory block, is easy to Project Realization.When realizing on the controller at the scene; It is 1.1MB that whole procedure takies space, field controller memory block; Process committed memory size is 1.3MB only during operation, is applicable to networked control technologys such as the extensive fieldbus that adopts in present building building automatic control field, EPA.

Description of drawings

Fig. 1 is air conditioning area temperature prediction model structure figure

Fig. 2 is terminal air-valve forecast model structure chart

Fig. 3 is a main air duct duct static pressure forecast model structure chart

Fig. 4 is a main air duct air quantity forecast model structure chart

Fig. 5 is a wind pushing temperature forecast model structure chart

Fig. 6 is new wind air-valve forecast model structure chart

Fig. 7 is an air quality forecast model structure chart

Fig. 8 is a VAV air conditioning terminal schematic diagram

Fig. 9 is terminal tandem predictive control loop figure

Figure 10 is terminal interior ring neural network prediction controller architecture figure

Figure 11 is terminal outer shroud neural network prediction controller architecture figure

Figure 12 is terminal outer shroud and interior ring neural network prediction controller searching process figure

Figure 13 is the static pressure predictive control loop figure of adjustable settings static pressure

Figure 14 is total blast volume predictive control loop figure

Figure 15 is wind pushing temperature predictive control loop figure

Figure 16 is new wind ratio predictive control loop figure

Figure 17 is static pressure neural network prediction controller architecture figure

Figure 18 is total blast volume neural network prediction controller architecture figure

Figure 19 is a wind pushing temperature predictive controller structure chart

Figure 20 is new wind ratio outer shroud neural network prediction controller architecture figure

Figure 21 is ring neural network prediction controller architecture figure in the new wind ratio

Figure 22 is static pressure neural network prediction controller searching process figure

Figure 23 is total blast volume neural network prediction controller searching process figure

Figure 24 is wind pushing temperature neural network prediction controller searching process figure

Figure 25 is new wind ratio outer shroud and inside and outside neural network prediction controller searching process figure

The specific embodiment

To combine accompanying drawing and embodiment that the present invention is done further detailed description below.

The present invention is a kind of VAV air conditioning system with variable control method, specifically comprises following step:

The first step: utilize the BP neutral net to set up air conditioning terminal regional temperature forecast model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model

1) confirms air conditioning terminal regional temperature forecast model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model structure

As shown in Figure 1, the input signal of air conditioning area temperature prediction model is outdoor intensity of solar radiation, outdoor temperature, CO ₂Concentration, indoor temperature, air quantity and terminal valve opening are output as next indoor temperature constantly;

As shown in Figure 2, the input signal of terminal air-valve forecast model is terminal valve opening and duct static pressure (air-conditioning unit main air duct place), is output as next terminal air quantity constantly;

As shown in Figure 3, the input signal of main air duct duct static pressure forecast model is VAV air-conditioning unit rotation speed of fan, duct static pressure, return air CO2 concentration, outdoor temperature and intensity of solar radiation, is output as next duct static pressure constantly;

As shown in Figure 4, the input signal of main air duct air quantity forecast model is rotation speed of fan and duct static pressure, is output as next main air duct air quantity constantly;

As shown in Figure 5, the input signal of wind pushing temperature forecast model is wind pushing temperature and water valve aperture, is output as next wind pushing temperature constantly;

As shown in Figure 6, the input signal of new wind air-valve forecast model is new air valve aperture and duct static pressure (place, fresh wind tube road), is output as next resh air requirement constantly;

As shown in Figure 7, the input signal of air quality forecast model is new wind air-valve valve opening and CO2 concentration, is output as next CO2 concentration constantly.

2) gather sample data

The sampling time scope be at 8 in the morning to 6 pm, sampling time interval 150 seconds, each forecast model are gathered about 2000 groups of data; Getting the time interval when setting up forecast model is 5 minutes;

Air conditioning area temperature prediction model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and also 1V is divided into ten grades to terminal valve area to the 10V interval by 0V, gathers outdoor intensity of solar radiation, outdoor temperature, CO ₂Concentration, indoor temperature and air quantity;

Terminal air-valve forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and also 1V is divided into ten grades to terminal valve area to the 10V interval by 0V, gathers duct static pressure (air-conditioning unit main air duct place) and air quantity;

Main air duct duct static pressure forecast model: after each terminal debugging was accomplished, all 1V was divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, gathers VAV air-conditioning unit rotation speed of fan, duct static pressure, return air CO2 concentration, outdoor temperature, intensity of solar radiation;

Main air duct air quantity forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, gathers rotation speed of fan, duct static pressure, main air duct air quantity;

The wind pushing temperature forecast model: 1V is divided into ten grades to the 10V interval by 0V with the water valve aperture, gathers wind pushing temperature and water valve aperture;

New wind air-valve forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and also 1V is divided into ten grades to new wind valve area to the 10V interval by 0V, gathers duct static pressure (place, fresh wind tube road) and air quantity;

The air quality forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and also 1V is divided into ten grades to new wind valve area to the 10V interval by 0V, gathers CO2 concentration;

3) sample data is by formula carried out normalization in (1), (2):

$x_i=\frac{x_{\mathrm{di}}-x_{d\min }}{x_{d\max }-x_{d\min }}---\left(1\right)$

$y_{\mathrm{tl}}=\frac{y_{\mathrm{dl}}-y_{d\min }}{y_{d\max }-y_{d\min }}---\left(2\right)$

X wherein _iBe the input value of neutral net after the normalization, x _DiBe former input value, x _DminBe the minimum of a value in the former input value, x _DmaxBe the maximum in the former input value; y _TlBe the desired value of neutral net after the normalization, y _DlRepresent former desired value, y _DminRepresent the minimum of a value in the former desired value, y _DmaxBe the maximum in the former desired value;

4) by table 1 parameter neutral net is carried out off-line training:

Table 1 Neural Network Training Parameter table

The neutral net type	Single latent layer forward direction BP network
		The input layer number	A *
Output layer node number	1
		The number of hidden nodes	B *
The neuron excitation function	Latent layer ' tansig ', output layer ' purelin '
		Learning function	′learngdm′
Performance function	' msereg ' (weighted mean square is poor)
		The network training function	' trainbr ' (Bayes's normalization method)
Power (threshold) value initialization method	' initnw ' (Nguyen-Widrow method)
		The maximum training time	2000Epochs
Target error	0

A ^*: air conditioning area temperature prediction model value 6, terminal air-valve forecast model, wind pushing temperature forecast model, main air duct air quantity forecast model, air quality forecast model and new wind air-valve forecast model value 2, main air duct duct static pressure forecast model value 5;

B ^*：

Second step: VAV air conditioning terminal tandem PREDICTIVE CONTROL

The present invention adopts the tandem forecast Control Algorithm for the VAV air conditioning terminal.Fig. 8 shows VAV air conditioning terminal schematic diagram, and air-valve, pressure sensor and air velocity transducer are installed in the air delivery duct.Temperature sensor and CO ₂Concentration sensor is separately positioned in air conditioning area and the return airway; At outdoor mounting temperature sensor and the intensity of solar radiation sensor of also needing.Terminal tandem control loop as shown in Figure 9; Cascade control system inner and outer ring controller is all selected the neural network prediction controller for use; According to the air conditioning area temperature of gathering; Input outer shroud neural network prediction controller; Calculate the setting air quantity, ring neural network prediction controller is adjusted terminal valve area in utilizing again.Predictive controller not only can improve the trace performance of air quantity setting value, eliminates the influence of static pressure to air quantity simultaneously, has improved control accuracy.

1) confirms terminal outer shroud PREDICTIVE CONTROL object function

Terminal outer shroud PREDICTIVE CONTROL object function is: $J_o\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_c-1}{\mathrm{\&Sigma;}}}{L}_o\left(T_o\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_c-1}{\mathrm{\&Sigma;}}}{\left(T_o\right[k]-T_{o,\mathrm{set}}[k\left]\right)}^2$

M wherein _cBe prediction time domain, t ₁Be the initial time in the prediction time domain, T _o[k] is the air conditioning area temperature in k sampling period, T _{O, set}[k] is the air conditioning area desired temperature in k sampling period, L _oIt is terminal outer shroud object function of k sampling period;

2) confirm terminal interior ring PREDICTIVE CONTROL object function

Ring PREDICTIVE CONTROL object function is in terminal: $J_i\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_i-1}{\mathrm{\&Sigma;}}}{L}_i\left(V_i\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_i-1}{\mathrm{\&Sigma;}}}{\left(V_i\right[k]-V_{i,\mathrm{set}}[k\left]\right)}^2$

V wherein _i[k] is the air quantity in k sampling period, V _{I, set}[k] is the air quantity setting value in k sampling period; M _iBe prediction time domain, L _iIt is ring object function in k the sampling period end;

3) confirm terminal interior ring neural network prediction controller and terminal outer shroud neural network prediction controller architecture

Encircle the neural network prediction controller architecture as shown in figure 10 in terminal, air quantity setting value, duct static pressure (air-conditioning unit main air duct place) are input with the pipeline air quantity, and valve opening is for exporting.The number of hidden nodes of controller neutral net is 5;

The structure of terminal outer shroud neural network prediction controller as shown in figure 11, input parameter comprises outdoor temperature, intensity of solar radiation, indoor temperature and desired temperature, output parameter is the air quantity setting value.The number of hidden nodes of controller neutral net is 8;

4) ring neural network prediction controller in terminal and terminal outer shroud neural network prediction controller are carried out online optimizing training

Each of interior ring of initialization end and outer shroud neural network prediction controller is connected weights, and assignment is a small random number in [1,1] scope, computing controller output then.Prediction time domain M _iAnd M _oBe respectively 3 and 6, predetermined period was got 5 minutes;

Searching process as shown in figure 12, x[k wherein] be k the relevant state variables parameter of outer shroud controlled device (being air conditioning area), i.e. indoor temperature constantly;

Be k+1 forecast model output constantly, i.e. k+1 indoor predicted temperature constantly; x ^*It is desired temperature; S[k] be k air quantity setting value constantly; Y[k] be the interior constantly ring relevant state variables of k parameter, i.e. air quantity; C[k] for encircling the k terminal valve area of the moment that the neural network prediction controller calculates in terminal; U[k] be that the terminal valve area that i.e. optimizing obtains is exported in terminal interior ring neural network prediction controller optimizing end back k control constantly;

The k+1 that exports for terminal air-valve forecast model predicts air quantity constantly.With x[t ₁], x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount S[t ₁], with S[t ₁], Y[t ₁] and-1 act on terminal in ring neural network prediction controller, obtain C[t ₁].Then with C[t ₁] act on terminal air-valve forecast model, obtain predicting air quantity

Keep terminal interior ring neural network prediction controller weights constant, will

S[t ₁] and-1 act on terminal in ring neural network prediction controller, obtain C ' [t ₁+ 1].With C ' [t ₁+ 1] the terminal air-valve forecast model of input obtains Will

S[t ₁] and-1 act on terminal in ring neural network prediction controller obtain C ' [t ₁+ 2]; With C ' [t ₁+ 2] the terminal air-valve forecast model of input obtains

The data that calculated are preserved.Make λ _i[k+M _i]=0, from after the λ the calculating formula (3) respectively forward _i[k] and q _i[k]:

$q_i\left[k\right]=\frac{\&PartialD;f_v{\left(k\right)}^T}{\&PartialD;C\left[k\right]}{\mathrm{\&lambda;}}_i[k+1]+\frac{\&PartialD;L_i{\left(k\right)}^T}{\&PartialD;C\left[k\right]}$

${\mathrm{\&lambda;}}_i\left[k\right]=\frac{\&PartialD;f_v{\left(k\right)}^T}{\&PartialD;y\left[k\right]}{\mathrm{\&lambda;}}_i[k+1]+\frac{\&PartialD;L_i{\left(k\right)}^T}{\&PartialD;y\left[k\right]}+\frac{\&PartialD;g_i{(k;W_i)}^T}{\&PartialD;y\left[k\right]}{q}_i\left[k\right]---\left(3\right)$

F in the formula _v(k) the terminal air-valve forecast model of setting up before the representative; g _i(k; W _i) be terminal interior ring neural network prediction controller equation.According to the q that calculates _i[k], through type (4) and formula (5) are revised the weights of terminal interior ring neural network prediction controller:

$\mathrm{\&Delta;}{W}_i=-{\mathrm{\&mu;}}_i\frac{\&PartialD;g_i{(k,W_i)}^T}{\&PartialD;W_i}{q}_i\left[k\right]---\left(4\right)$

W _i＝W _i+ΔW _i (5)

W wherein _iBe the weights battle array of terminal interior ring neural network prediction controller, μ _iBe the right value update rate, μ _iSelect 0.05.Constantly revise the weights of terminal interior ring neural network prediction controller, until Δ W _i＜0.001;

With S[t ₁], Y[t ₁] and-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁].With u[t] act on terminal air-valve forecast model, obtain

Afterwards with x[t ₁], Act on the air conditioning area forecast model, obtain

Keep terminal outer shroud neural network prediction controller weights constant, will

x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 1]; Ring optimizing in carrying out again is with s ' [t ₁+ 1],

With-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁+ 1].Again with u ' [t ₁+ 1] acts on terminal air-valve forecast model, obtain

Afterwards with x[t ₁+ 1], Act on the air conditioning area forecast model, obtain

Will

x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 2]; Ring optimizing in carrying out again is with s ' [t ₁+ 2],

With-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁+ 2].Again with u ' [t ₁+ 2] act on terminal air-valve forecast model, obtain Afterwards, will

Act on the air conditioning area forecast model, obtain

And the rest may be inferred, utilizes controller neutral net and object forecast model to extrapolate following u ' [t ₁+ i] and

Value, i=3 wherein ..., 6, and the data that calculated are preserved.Make λ _o[k+M _o]=0, from after the λ the calculating formula (6) respectively forward _o[k] and q _o[k]:

$q_o\left[k\right]=\frac{\&PartialD;f_z{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_o[k+1]+\frac{\&PartialD;L_o{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_o\left[k\right]=\frac{\&PartialD;f_z{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_o[k+1]+\frac{\&PartialD;L_o{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_o{(k;W_o)}^T}{\&PartialD;x\left[k\right]}{q}_o\left[k\right]---\left(6\right)$

F in the formula _z(k) the air conditioning area forecast model of setting up before the representative; g _o(k; W _o) be terminal outer shroud neural network prediction controller equation.According to the q that calculates _o[k] is for k=t ₁+ M _o-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (7) and formula (8) are revised the weights of terminal outer shroud neural network prediction controller:

$\mathrm{\&Delta;}{W}_o=-{\mathrm{\&mu;}}_o\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_o-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_o{(k,W_o)}^T}{\&PartialD;W_o}{q}_o\left[k\right]---\left(7\right)$

W _o＝W _o+ΔW _o (8)

W wherein _oBe the weights battle array of terminal outer shroud neural network prediction controller, μ _oBe the right value update rate, μ _oSelect 0.05.Constantly revise the weights of terminal outer shroud neural network prediction controller, until Δ W _o＜0.001; After the outer shroud optimizing finishes, with x[t ₁] and x ^*Import terminal outer shroud neural network prediction controller, again with the S[t that obtains ₁] and Y[t ₁] the terminal interior ring neural network prediction controller of input, will export u[t ₁] directly act on terminal air-valve;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

The 3rd step: VAV air-conditioning unit PREDICTIVE CONTROL

The present invention is respectively static pressure control loop (Figure 13), VAV air-conditioning unit total blast volume control loop (Figure 14), VAV air-conditioning unit wind pushing temperature control loop (Figure 15) and the VAV air-conditioning unit new wind ratio control loop (Figure 16) of VAV air-conditioning unit adjustable settings static pressure through the control of 4 control loops completion air-conditioning units.

The present invention realizes the control to air-conditioning unit air quantity according to different working conditions through two kinds of control strategies: when the static pressure decline of system static pressure monitoring point reaches setting value, select the static pressure PREDICTIVE CONTROL strategy (Figure 13) of adjustable settings static pressure; When the static pressure of system static pressure monitoring point is higher than setting value, then move total blast volume PREDICTIVE CONTROL strategy (Figure 14).

The present invention is different from conventional VAV air-conditioning system and all adopts the control strategy of deciding wind pushing temperature, variable air rate; But adopt the Different control strategy according to the different working conditions condition; When moving owing to minimum resh air requirement operating mode, cause the low excessively problem of temperature with solution transition season part air conditioning area.Concrete control scheme is: the control system detects all terminal VAV-BOX primary air flows, when a certain VAV-BOX primary air flow of appearance is lower than nominal air delivery 30%, reduces the water valve aperture, improves 0.5 ℃ of wind pushing temperature; When a certain terminal VAV-BOX primary air flow greater than 70% the time, increase the water valve aperture, reduce by 0.5 ℃ of wind pushing temperature (Figure 15).

The conventional control strategy of VAV air-conditioning system resh air requirement is the aperture of new wind air-valve of control and return air air-valve; And for the recuperation of heat unit is set on the roof; Adopt the vertical system that concentrates new wind; Because the dynamic and static pressure of vertical VMC relation, the resh air requirement adjustment of adjacent floor can influence the quantity delivered of the new wind of this layer.The present invention adopts the tandem PREDICTIVE CONTROL of pressure independent type in the new wind ratio control loop; According to the return air CO2 concentration of gathering; Input outer shroud neural network prediction controller calculates the setting resh air requirement, and ring neural network prediction controller is adjusted new wind valve area in utilizing again.Outer shroud is to carry out PREDICTIVE CONTROL according to air conditioning area CO2 concentration; And interior ring is the PREDICTIVE CONTROL of resh air requirement being carried out the pressure independent type; Thereby reach accurate control, under the prerequisite that satisfies the air conditioning area environmental quality, realize maximum energy-saving control (Figure 16) resh air requirement.

1) confirms VAV air-conditioning unit static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and VAV air-conditioning unit new wind ratio control loop PREDICTIVE CONTROL object function

Static pressure control loop PREDICTIVE CONTROL object function is:

$J_s\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_s-1}{\mathrm{\&Sigma;}}}{L}_s\left(P_s\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_s-1}{\mathrm{\&Sigma;}}}\{{\left(P_s\right[k]-P_{s,\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

P wherein _s[k] is the duct static pressure in k sampling period, P _{S, set}[k] is that the air conditioning area in k sampling period is set static pressure, M _sBe prediction time domain, U _Fan[k] is k sampling period blower voltage controlled quentity controlled variable, L _sBe k sampling period static pressure control loop object function;

Total blast volume control loop PREDICTIVE CONTROL object function is:

$J_f\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_f-1}{\mathrm{\&Sigma;}}}{L}_f\left(V_{\mathrm{fan}}\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_f-1}{\mathrm{\&Sigma;}}}\{{\left(V_{\mathrm{fan}}\right[k]-V_{\mathrm{fan}.\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

V wherein _Fan[k] is the pipeline air quantity in k sampling period, M _fBe prediction time domain, V _Fan.set[k] is each terminal prediction air quantity sum in k sampling period, L _fBe k sampling period total blast volume control loop object function;

Wind pushing temperature control loop PREDICTIVE CONTROL object function is:

$J_{\mathrm{st}}\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{st}}\left(T_{\mathrm{st}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{\left(T_{\mathrm{st}}\right[k]-T_{\mathrm{st},\mathrm{set}}[k\left]\right)}^2$

T wherein _St[k] is the temperature in k sampling period, T _{St, set}[k] is the design temperature in k sampling period, M _StBe prediction time domain, L _StBe k sampling period wind pushing temperature control loop object function;

New wind ratio control loop outer shroud PREDICTIVE CONTROL object function is:

$J_q\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_q-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{no}}\left(Q\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_q-1}{\mathrm{\&Sigma;}}}{\left(Q\right[k]-Q_{\mathrm{set}}[k\left]\right)}^2$

Q[k wherein] be the air quality in k sampling period, Q _Set[k] is the air quality setting value in k sampling period, M _qBe prediction time domain, L _NoBe k sampling period new wind ratio control loop outer shroud object function;

Ring PREDICTIVE CONTROL object function is in the new wind ratio control loop:

$J_{\mathrm{ni}}\left[k\right]=\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{ni}}\left(V_{\mathrm{ni}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{\left(V_{\mathrm{ni}}\right[k]-S[k\left]\right)}^2$

V wherein _i[k] is the prediction air-valve discharge quantity of fan in k sampling period, S[k] be the setting air quantity in k sampling period, M _iBe prediction time domain, L _NiIt is ring object function in k the sampling period new wind ratio control loop;

2) confirm VAV air-conditioning unit static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and VAV air-conditioning unit new wind ratio control loop neural network prediction controller architecture

The input signal of static pressure neural network prediction controller is output as rotation speed of fan (control signal of motor frequency conversion device is 0～10V voltage) for setting duct static pressure and duct static pressure; Hidden node is 3 (Figure 17);

Total blast volume neural network prediction controller input signal is pipeline air quantity and total blast volume, is output as rotation speed of fan (control signal of motor frequency conversion device is 0～10V voltage); Hidden node is 3 (Figure 18);

The input signal of wind pushing temperature neural network prediction controller is output as water valve aperture (control signal of electric valve executing mechanism is 0～10V voltage) for setting wind pushing temperature and wind pushing temperature; Hidden node is 3 (Figure 19);

New wind ratio outer shroud neural network prediction controller input signal is output as the setting resh air requirement for setting CO2 concentration and return air CO2 concentration; Hidden node is 4 (Figure 20); Ring neural network prediction controller will be set resh air requirement, duct static pressure (place, fresh wind tube road) and pipeline air quantity as input in the new wind ratio, and new air valve aperture is output (control signal of electric valve executing mechanism is 0～10V voltage); Hidden node is 5 (Figure 21);

3) the online optimizing of static pressure neural network prediction controller in the static pressure predictive control loop of realization adjustable settings static pressure

Each of initialization static pressure neural network prediction controller connects weights, and assignment is a small random number in [1,1] scope, computing controller output then.Prediction time domain M _sIt was 2 steps, for the right value update rate μ of static pressure neural network prediction controller _sSelect 0.05, predetermined period was got 5 minutes.

Searching process as shown in figure 22, x[k wherein] parameter be duct static pressure, u _Fan[k] is the blower fan controlled quentity controlled variable, x ^*Be to set duct static pressure.With x[t ₁], x ^*With-1 act on static pressure neural network prediction controller, controlled amount u _Fan[t ₁], then with u _Fan[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u _Fan[t ₁], x[t ₁] and-1 act on main air duct duct static pressure forecast model, obtain Keep static pressure neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain

Will

X[t ₁+ 1] and-1 act on main air duct duct static pressure forecast model, obtains

The data that calculated are preserved.Make λ _s[k+M _s]=0, from after the λ the calculating formula (9) respectively forward _s[k] and q _s[k]:

$q_s\left[k\right]=\frac{\&PartialD;f_s{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}{\mathrm{\&lambda;}}_s[k+1]+\frac{\&PartialD;L_s{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}$

${\mathrm{\&lambda;}}_s\left[k\right]=\frac{\&PartialD;f_s{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_s[k+1]+\frac{\&PartialD;L_s{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_s{(k;W_s)}^T}{\&PartialD;x\left[k\right]}{q}_s\left[k\right]---\left(9\right)$

F in the formula _s(k) the main air duct duct static pressure forecast model of setting up before the representative, g _s(k; W _s) be static pressure neural network prediction controller equation.According to the q that calculates _s[k] is for k=t ₁+ M _s-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (10) and formula (11) are revised the weights of static pressure neural network prediction controller:

$\mathrm{\&Delta;}{W}_s=-{\mathrm{\&mu;}}_s\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_s-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_s{(k,W_s)}^T}{\&PartialD;W_s}{q}_s\left[k\right]---\left(10\right)$

W _s＝W _s+ΔW _s (11)

W wherein _sBe the weights battle array of static pressure neural network prediction controller, μ _sBe the right value update rate, μ _sSelect 0.05.Constantly revise the weights of static pressure neural network prediction controller, until Δ W _s＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

4) realize the online optimizing of total blast volume neural network prediction controller in the total blast volume control loop

Each of initialization total blast volume neural network prediction controller connects weights, and assignment is a small random number in [1,1] scope, computing controller output then.Prediction time domain M _fIt was 2 steps, for the right value update rate μ of total blast volume neural network prediction controller _fSelect 0.05, predetermined period was got 5 minutes.

Searching process as shown in figure 23, x[k wherein] parameter be the main air duct air quantity, u _Fan[k] is the blower fan controlled quentity controlled variable, x ^*It is each regional air quantity setting value sum.With x[t ₁], x ^*With-1 act on total blast volume neural network prediction controller, controlled amount u _Fan[t ₁], then with u _Fan[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u _Fan[t ₁], x[t ₁] and-1 act on main air duct air quantity forecast model, obtain

Keep total blast volume neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain

Will

X[t ₁+ 1] and-1 act on main air duct air quantity forecast model, obtains

The data that calculated are preserved.Make λ _f[k+M _f]=0, from after the λ the calculating formula (12) respectively forward _f[k] and q _f[k]:

$q_f\left[k\right]=\frac{\&PartialD;f_f{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}{\mathrm{\&lambda;}}_f[k+1]+\frac{\&PartialD;L_f{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}$

${\mathrm{\&lambda;}}_f\left[k\right]=\frac{\&PartialD;f_f{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_f[k+1]+\frac{\&PartialD;L_f{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_f{(k;W_f)}^T}{\&PartialD;x\left[k\right]}{q}_f\left[k\right]---\left(12\right)$

F in the formula _f(k) the main air duct air quantity forecast model of setting up before the representative, g _f(k; W _f) be total blast volume neural network prediction controller equation.According to the q that calculates _f[k] is for k=t ₁+ M _f-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (13) and formula (14) are revised the weights of total blast volume neural network prediction controller:

$\mathrm{\&Delta;}{W}_f=-{\mathrm{\&mu;}}_f\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_f-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_f{(k,W_f)}^T}{\&PartialD;W_f}{q}_f\left[k\right]---\left(13\right)$

W _f＝W _f+ΔW _f (14)

W wherein _fBe the weights battle array of total blast volume neural network prediction controller, μ _fBe the right value update rate, μ _fSelect 0.05.Constantly revise the weights of total blast volume neural network prediction controller, until Δ W _f＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

5) realize the online optimizing of wind pushing temperature neural network prediction controller in the wind pushing temperature control loop

Each of initialization wind pushing temperature neural network prediction controller connects weights, and assignment is a small random number in [1,1] scope, computing controller output then.Prediction time domain M _StIt was 2 steps, for the right value update rate μ of wind pushing temperature neural network prediction controller _StSelect 0.05, predetermined period was got 5 minutes.

Searching process as shown in figure 24, x[k wherein] parameter be wind pushing temperature, u[k] be the water valve aperture, x ^*Be to set wind pushing temperature.With x[t ₁], x ^*With-1 act on wind pushing temperature neural network prediction controller, controlled amount u[t ₁], then with u[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u[t ₁], x[t ₁] and-1 act on the wind pushing temperature forecast model, obtain

Keep wind pushing temperature neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain u ' [t ₁+ 1], with u ' [t ₁+ 1], x[t ₁+ 1] and-1 act on the wind pushing temperature forecast model, obtains

The data that calculated are preserved.Make λ _St[k+M _St]=0, from after the λ the calculating formula (15) respectively forward _St[k] and q _St[k]:

$q_{\mathrm{st}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{st}}[k+1]+\frac{\&PartialD;L_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{st}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{st}}[k+1]+\frac{\&PartialD;L_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_{\mathrm{st}}{(k;W_{\mathrm{st}})}^T}{\&PartialD;x\left[k\right]}{q}_{\mathrm{st}}\left[k\right]---\left(15\right)$

F in the formula _St(k) the wind pushing temperature forecast model of setting up before the representative, g _St(k; W _St) be wind pushing temperature neural network prediction controller equation.According to the q that calculates _St[k] is for k=t ₁+ M _St-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (16) and formula (17) are revised the weights of wind pushing temperature neural network prediction controller:

$\mathrm{\&Delta;}{W}_{\mathrm{st}}=-{\mathrm{\&mu;}}_{\mathrm{st}}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_{\mathrm{st}}{(k,W_{\mathrm{st}})}^T}{\&PartialD;W_{\mathrm{st}}}{q}_{\mathrm{st}}\left[k\right]---\left(16\right)$

W _st＝W _st+ΔW _st (17)

W wherein _StBe the weights battle array of wind pushing temperature neural network prediction controller, μ _StBe the right value update rate, μ _StSelect 0.05.Constantly revise the weights of wind pushing temperature neural network prediction controller, until Δ W _St＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

6) realize the online optimizing of neural network prediction controller in the new wind ratio control loop

Each of interior ring of initialization new wind ratio and outer shroud neural network prediction controller is connected weights, and assignment is a small random number in [1,1] scope; Computing controller output then is if output valve is all in [0.7,1] scope; Then use this group controller initial weight, otherwise random initializtion again.Prediction time domain M _NiAnd M _qBe respectively 3 and 6, for the right value update rate μ selection 0.05 of ring and outer shroud neural network prediction controller in the new wind ratio, predetermined period was got 5 minutes;

Searching process as shown in figure 25, x[k wherein] be k CO2 concentration constantly;

Be k forecast model output constantly, i.e. k+1 return air CO2 concentration constantly;

Be t ₁Return air CO2 concentration constantly; x ^*It is the CO2 concentration set point; S[k] set resh air requirement constantly for k; Y[k] be that k encircles the relevant state variables parameter constantly, comprise pipeline air quantity and duct static pressure; O[k] the constantly new wind valve area of the k that calculates for interior ring neural network prediction controller; U[k] be that interior ring neural network prediction controller optimizing finishes back k control output constantly, the i.e. new wind valve area that optimizing obtains;

The k+1 that exports for new wind air-valve forecast model predicts air quantity constantly.With x[t ₁], x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount S[t ₁], with S[t ₁], Y[t ₁] and-1 act on ring neural network prediction controller in the new wind ratio, obtain O[t ₁].Then with O[t ₁] act on new wind air-valve forecast model respectively, obtain predicting air quantity

It is constant to keep encircling neural network prediction controller weights in the new wind ratio, will S[t ₁] and-1 act on ring neural network prediction controller in the new wind ratio, obtain O ' [t ₁+ 1].With O ' [t ₁+ 1] the new wind air-valve forecast model of input obtains

Will S[t ₁] and-1 act on that ring neural network prediction controller obtains O ' [2] in the new wind ratio.The O '[2] Enter the new air damper prediction model, get

to the calculated data saved.Make λ _Ni[k+2]=0, from after the λ the calculating formula (18) respectively forward _Ni[k] and q _Ni[k]:

$q_{\mathrm{ni}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{nv}}{\left(k\right)}^T}{\&PartialD;O\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{ni}}[k+1]+\frac{\&PartialD;L_{\mathrm{ni}}{\left(k\right)}^T}{\&PartialD;O\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{ni}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{nv}}{\left(k\right)}^T}{\&PartialD;y\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{ni}}[k+1]+\frac{\&PartialD;L_{\mathrm{ni}}{\left(k\right)}^T}{\&PartialD;y\left[k\right]}+\frac{\&PartialD;g_{\mathrm{ni}}{(k;W_{\mathrm{ni}})}^T}{\&PartialD;y\left[k\right]}{q}_{\mathrm{ni}}\left[k\right]---\left(18\right)$

F in the formula _Nv(k) the new wind air-valve forecast model of setting up before the representative; g _Ni(k; W _Ni) be ring neural network prediction controller equation in the new wind ratio; According to the q that calculates _Ni[k], through type (19) and formula (20) are revised the weights of ring neural network prediction controller in the new wind ratio:

$\mathrm{\&Delta;}{W}_{\mathrm{ni}}=-{\mathrm{\&mu;}}_{\mathrm{ni}}\frac{\&PartialD;g_{\mathrm{ni}}{(k,W_{\mathrm{ni}})}^T}{\&PartialD;W_{\mathrm{ni}}}{q}_{\mathrm{ni}}\left[k\right]---\left(19\right)$

W _ni＝W _ni+ΔW _ni (20)

W wherein _NiBe the weights battle array of ring neural network prediction controller in the new wind ratio, μ _NiBe the right value update rate, μ _NiSelect 0.05.Constantly revise the weights of ring neural network prediction controller in the new wind ratio, until Δ W _Ni＜0.001;

With S[t ₁], Y[t ₁] and-1 new wind ratio that acts on after the optimizing in ring neural network prediction controller, obtain u[t ₁].Then with u[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u[t ₁] act on new wind air-valve forecast model, obtain

Afterwards with x[t ₁],

Act on the air conditioning area forecast model, obtain

Keep new wind ratio outer shroud neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 1]; Ring optimizing in carrying out again is with s ' [t ₁+ 1],

With ring neural network prediction controller in-1 new wind ratio that acts on after the optimizing, obtain u ' [t ₁+ 1].Again with u ' [t ₁+ 1] acts on new wind air-valve forecast model, obtain

Afterwards, with x[t ₁+ 1],

Act on the air conditioning area forecast model, obtain

Will

x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 2]; Ring optimizing in carrying out again is with s ' [t ₁+ 2],

With ring neural network prediction controller in-1 new wind ratio that acts on after the optimizing, obtain u ' [t ₁+ 2].Again with u ' [t ₁+ 2] act on new wind air-valve forecast model, obtain Afterwards, will

Act on the air conditioning area forecast model, obtain

And the rest may be inferred, utilizes controller neutral net and object forecast model to extrapolate following u ' [t ₁+ i] and

Value, i=3 wherein ..., 6, and the data that calculated are preserved.Make λ _No[k+M _q]=0, from after the λ the calculating formula (21) respectively forward _No[k] and q _No[k]:

$q_{\mathrm{no}}\left[k\right]=\frac{\&PartialD;f_n{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{no}}[k+1]+\frac{\&PartialD;L_{\mathrm{no}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{no}}\left[k\right]=\frac{\&PartialD;f_n{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{no}}[k+1]+\frac{\&PartialD;L_{\mathrm{no}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_{\mathrm{no}}{(k;W_{\mathrm{no}})}^T}{\&PartialD;x\left[k\right]}{q}_{\mathrm{no}}\left[k\right]---\left(21\right)$

F in the formula _n(k) the air quality forecast model of setting up before the representative; g _No(k; W _No) be new wind ratio outer shroud neural network prediction controller equation.According to the q that calculates _No[k] is for k=t ₁+ M _q-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (22) and formula (23) are revised the weights of new wind ratio outer shroud neural network prediction controller:

$\mathrm{\&Delta;}{W}_{\mathrm{no}}=-{\mathrm{\&mu;}}_{\mathrm{no}}\underset{k=t_1}{\overset{{\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000082058780000023.png,HDA0000082058780000033.png"\; state="\{\{state\}\}">t</figure-callout>}}_1+M_q-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_{\mathrm{no}}{(k,W_{\mathrm{no}})}^T}{\&PartialD;W_{\mathrm{no}}}{q}_{\mathrm{no}}\left[k\right]---\left(22\right)$

W _no＝W _no+ΔW _no (23)

W wherein _NoBe the weights battle array of new wind ratio outer shroud neural network prediction controller, μ _NoBe the right value update rate, μ _NoSelect 0.05.Constantly revise the weights of new wind ratio outer shroud neural network prediction controller, until Δ W _No＜0.001; After the outer shroud optimizing finishes, with x[t ₁] and x ^*Input new wind ratio outer shroud neural network prediction controller is again with the S[t that obtains ₁] and Y[t ₁] the interior ring of input new wind ratio neural network prediction controller, will export u[t ₁] directly act on new wind air-valve;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

Claims (8)

1. a VAV air conditioning system with variable control method is characterized in that, comprises following step:

The first step: utilize the BP neutral net to set up air conditioning area temperature prediction model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model

1) confirms air conditioning area temperature prediction model, terminal air-valve forecast model, air-conditioning unit main air duct duct static pressure forecast model, main air duct air quantity forecast model, wind pushing temperature forecast model, new wind air-valve forecast model and air quality forecast model structure

The input signal of air conditioning area temperature prediction model is outdoor intensity of solar radiation, outdoor temperature, CO2 concentration, indoor temperature, air quantity and terminal valve opening, is output as next indoor temperature constantly;

The input signal of terminal air-valve forecast model is terminal valve area and air-conditioning unit main air duct duct static pressure, is output as next terminal air quantity constantly;

The input signal of main air duct duct static pressure forecast model is VAV air-conditioning unit rotation speed of fan, duct static pressure, return air CO2 concentration, outdoor temperature and intensity of solar radiation, is output as next duct static pressure constantly;

The input signal of main air duct air quantity forecast model is rotation speed of fan and duct static pressure, is output as next main air duct air quantity constantly;

The input signal of wind pushing temperature forecast model is wind pushing temperature and water valve aperture, is output as next wind pushing temperature constantly;

The input signal of new wind air-valve forecast model is new wind valve area and fresh wind tube road pipeline static pressure, is output as next resh air requirement constantly;

The input signal of air quality forecast model is new wind valve area and CO2 concentration, is output as next CO2 concentration constantly;

2) gather sample data;

3) sample data is by formula carried out normalization in (1), (2):

$x_i=\frac{x_{\mathrm{di}}-x_{d\min }}{x_{d\max }-x_{d\min }}---\left(1\right)$

$y_{\mathrm{tl}}=\frac{y_{\mathrm{dl}}-y_{d\min }}{y_{d\max }-y_{d\min }}---\left(2\right)$

X wherein _iBe the input value of neutral net after the normalization, x _DiBe former input value, x _DminBe the minimum of a value in the former input value, x _DmaxBe the maximum in the former input value; y _TlBe the desired value of neutral net after the normalization, y _DlRepresent former desired value; y _DminRepresent the minimum of a value in the former desired value; y _DmaxBe the maximum in the former desired value;

4) neutral net is carried out off-line training;

Second step: VAV air conditioning terminal tandem PREDICTIVE CONTROL

1) confirms terminal outer shroud PREDICTIVE CONTROL object function

Terminal outer shroud PREDICTIVE CONTROL object function is: $J_o\left[k\right]=\underset{k=t_1}{\overset{t_1+M_c-1}{\mathrm{\&Sigma;}}}{L}_o\left(T_o\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_c-1}{\mathrm{\&Sigma;}}}{\left(T_o\right[k]-T_{o,\mathrm{set}}[k\left]\right)}^2$

M wherein _cBe prediction time domain, t ₁Be the initial time in the prediction time domain, T _o[k] is the air conditioning area temperature in k sampling period, T _{O, set}[k] is the air conditioning area desired temperature in k sampling period, L _oIt is terminal outer shroud object function of k sampling period;

2) confirm terminal interior ring PREDICTIVE CONTROL object function

Ring PREDICTIVE CONTROL object function is in terminal: $J_i\left[k\right]=\underset{k=t_1}{\overset{t_1+M_i-1}{\mathrm{\&Sigma;}}}{L}_i\left(V_i\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_i-1}{\mathrm{\&Sigma;}}}{\left(V_i\right[k]-V_{i,\mathrm{set}}[k\left]\right)}^2$

V wherein _i[k] is the air quantity in k sampling period, V _{I, set}[k] is the air quantity setting value in k sampling period; M _iBe prediction time domain, L _iIt is ring object function in k the sampling period end;

3) confirm terminal interior ring neural network prediction controller and terminal outer shroud neural network prediction controller architecture

Terminal tandem control loop inner and outer ring controller is all selected the neural network prediction controller for use.As input, valve opening is output to ring neural network prediction controller with air quantity setting value, air-conditioning unit main air duct duct static pressure and pipeline air quantity in terminal; The input parameter of terminal outer shroud neural network prediction controller comprises outdoor temperature, intensity of solar radiation, indoor temperature and air conditioning area desired temperature; Output parameter is the air quantity setting value;

4) ring neural network prediction controller in terminal and terminal outer shroud neural network prediction controller are carried out online optimizing training; Based on the air conditioning area temperature of gathering; Input outer shroud neural network prediction controller; Calculate the setting air quantity, ring neural network prediction controller is adjusted terminal valve area in utilizing again;

The 3rd step: VAV air-conditioning unit PREDICTIVE CONTROL

1) confirms static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and the VAV air-conditioning unit new wind ratio tandem control loop PREDICTIVE CONTROL object function of VAV air-conditioning unit adjustable settings static pressure

Static pressure control loop PREDICTIVE CONTROL object function is:

$J_s\left[k\right]=\underset{k=t_1}{\overset{t_1+M_s-1}{\mathrm{\&Sigma;}}}{L}_s\left(P_s\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_s-1}{\mathrm{\&Sigma;}}}\{{\left(P_s\right[k]-P_{s,\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

P wherein _s[k] is the duct static pressure in k sampling period, P _{S, set}[k] is that the air conditioning area in k sampling period is set static pressure, M _sBe prediction time domain, U _Fan[k] is k sampling period blower voltage controlled quentity controlled variable, L _sBe k sampling period static pressure control loop object function;

Total blast volume control loop PREDICTIVE CONTROL object function is:

$J_f\left[k\right]=\underset{k=t_1}{\overset{t_1+M_f-1}{\mathrm{\&Sigma;}}}{L}_f\left(V_{\mathrm{fan}}\right[k],U_{\mathrm{fan}}[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_f-1}{\mathrm{\&Sigma;}}}\{{\left(V_{\mathrm{fan}}\right[k]-V_{\mathrm{fan}.\mathrm{set}}[k\left]\right)}^2+U_{\mathrm{fan}}^2[k\left]\right\}$

V wherein _Fan[k] is the pipeline air quantity in k sampling period, M _fBe prediction time domain, V _Fan.set[k] is each terminal prediction air quantity sum in k sampling period, L _fBe k sampling period total blast volume control loop object function;

Wind pushing temperature control loop PREDICTIVE CONTROL object function is:

$J_{\mathrm{st}}\left[k\right]=\underset{k=t_1}{\overset{t_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{st}}\left(T_{\mathrm{st}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}{\left(T_{\mathrm{st}}\right[k]-T_{\mathrm{st},\mathrm{set}}[k\left]\right)}^2$

T wherein _St[k] is the temperature in k sampling period, T _{St, set}[k] is the design temperature in k sampling period, M _StBe prediction time domain, L _StBe k sampling period wind pushing temperature control loop object function;

New wind ratio control loop outer shroud PREDICTIVE CONTROL object function is:

$J_q\left[k\right]=\underset{k=t_1}{\overset{t_1+M_q-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{no}}\left(Q\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_q-1}{\mathrm{\&Sigma;}}}{\left(Q\right[k]-Q_{\mathrm{set}}[k\left]\right)}^2$

Q[k wherein] be the air quality in k sampling period, Q _Set[k] is the air quality setting value in k sampling period, M _qBe prediction time domain, L _NoBe k sampling period new wind ratio control loop outer shroud object function;

Ring PREDICTIVE CONTROL object function is in the new wind ratio control loop:

$J_{\mathrm{ni}}\left[k\right]=\underset{k=t_1}{\overset{t_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{L}_{\mathrm{ni}}\left(V_{\mathrm{ni}}\right[k],k)=\frac{1}{2}\underset{k=t_1}{\overset{t_1+M_{\mathrm{ni}}-1}{\mathrm{\&Sigma;}}}{\left(V_{\mathrm{ni}}\right[k]-S[k\left]\right)}^2$

V wherein _i[k] is the prediction air-valve discharge quantity of fan in k sampling period, S[k] be the setting air quantity in k sampling period, M _iBe prediction time domain, L _NiIt is ring object function in k the sampling period new wind ratio control loop;

2) confirm static pressure control loop, VAV air-conditioning unit total blast volume control loop, VAV air-conditioning unit wind pushing temperature control loop and the VAV air-conditioning unit new wind ratio tandem control loop neural network prediction controller architecture of VAV air-conditioning unit adjustable settings static pressure

The input signal of static pressure neural network prediction controller is output as rotation speed of fan for setting duct static pressure and duct static pressure;

Total blast volume neural network prediction controller input signal is pipeline air quantity and total blast volume, is output as rotation speed of fan;

The input signal of wind pushing temperature neural network prediction controller is output as the water valve aperture for setting wind pushing temperature and wind pushing temperature;

New wind ratio tandem control loop inner and outer ring controller is all selected the neural network prediction controller for use.New wind ratio outer shroud neural network prediction controller is input as sets CO2 concentration and return air CO2 concentration, is output as the setting resh air requirement; Ring neural network prediction controller will be set air quantity, fresh wind tube road static pressure and pipeline air quantity as input in the new wind ratio, and new air valve aperture is output;

3) when the static pressure decline of system static pressure monitoring point reaches setting value, select the static pressure PREDICTIVE CONTROL strategy of adjustable settings static pressure, realize the online optimizing of static pressure neural network prediction controller in the static pressure control loop;

4) when the static pressure of system static pressure monitoring point is higher than setting value, operation total blast volume PREDICTIVE CONTROL strategy is realized the online optimizing of total blast volume neural network prediction controller in the total blast volume control loop;

5) realize the online optimizing of wind pushing temperature neural network prediction controller in the wind pushing temperature control loop;

6) realize ring and the online optimizing of outer shroud neural network prediction controller in the new wind ratio in the new wind ratio tandem control loop; According to the return air CO2 concentration of gathering; Input outer shroud neural network prediction controller; Calculate the setting resh air requirement, ring neural network prediction controller is adjusted new wind valve area in utilizing again.

2. a kind of VAV air conditioning system with variable control method according to claim 1 is characterized in that, the first step 2) gather sample data and be specially:

The sampling time scope be at 8 in the morning to 6 pm, sampling time interval 150 seconds, each forecast model are gathered about 2000 groups of data; Getting the time interval when setting up forecast model is 5 minutes;

Air conditioning area temperature prediction model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V; Simultaneously; Also 1V is divided into ten grades to terminal valve area to the 10V interval by 0V, gathers outdoor intensity of solar radiation, outdoor temperature, CO2 concentration, indoor temperature and air quantity;

Terminal air-valve forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and simultaneously, also 1V is divided into ten grades to terminal valve area to the 10V interval by 0V, gathers air-conditioning unit main air duct duct static pressure and air quantity;

Main air duct duct static pressure forecast model: after each terminal debugging was accomplished, all 1V was divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, gathers VAV air-conditioning unit rotation speed of fan, duct static pressure, return air CO2 concentration, outdoor temperature, intensity of solar radiation;

Main air duct air quantity forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, gathers rotation speed of fan, duct static pressure, main air duct air quantity;

The wind pushing temperature forecast model: 1V is divided into ten grades to the 10V interval by 0V with the water valve aperture, gathers wind pushing temperature and water valve aperture;

New wind air-valve forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and simultaneously, also 1V is divided into ten grades to new wind valve area to the 10V interval by 0V, gathers fresh wind tube road pipeline static pressure and air quantity;

The air quality forecast model: all 1V is divided into ten grades to blower fan controlled quentity controlled variable signal to the 10V interval by 0V, and is same, and also 1V is divided into ten grades to new wind valve area to the 10V interval by 0V, gathers CO2 concentration.

3. a kind of VAV air conditioning system with variable control method according to claim 1 is characterized in that, the first step 4) carry out off-line training according to table 1 pair neutral net:

Table 1 Neural Network Training Parameter table

The neutral net type Single latent layer forward direction BP network The input layer number A ^* Output layer node number 1 The number of hidden nodes B ^* The neuron excitation function Latent layer ' tansig ', output layer ' purelin ' Learning function ′learngdm′ Performance function ′msereg′ The network training function ′trainbr′ Power (threshold) value initialization method ′initnw′ The maximum training time 2000Epochs Target error 0

A ^*: air conditioning area temperature prediction model value 6, terminal air-valve forecast model, air quality forecast model, wind pushing temperature forecast model, main air duct air quantity forecast model and new wind air-valve forecast model value 2, main air duct duct static pressure forecast model value 5;

B ^*：

4. a kind of VAV air conditioning system with variable control method according to claim 1; It is characterized in that; Second step 4) be specially: each of ring and outer shroud nerve network controller is connected weights in the initialization end; Assignment is [1; 1] small random number in the scope, computing controller output then; Prediction time domain M _iAnd M _oBe respectively 3 and 6, predetermined period was got 5 minutes;

Searching process is: establish x[k] be the k relevant state variables parameter of outer shroud controlled device (being air conditioning area), i.e. indoor temperature constantly;

Be k+1 forecast model output constantly, i.e. k+1 indoor predicted temperature constantly; x ^*It is desired temperature; S[k] be k air quantity setting value constantly; Y[k] be the interior constantly ring relevant state variables of k parameter, i.e. air quantity; C[k] for encircling the k terminal air-valve valve opening of the moment that the neural network prediction controller calculates in terminal; U[k] be that the terminal air-valve valve opening that i.e. optimizing obtains is exported in terminal interior ring neural network prediction controller optimizing end back k control constantly;

The k+1 that exports for valve end air-valve forecast model predicts air quantity constantly.With x[t ₁], x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount S[t ₁], with S[t ₁], Y[t ₁] and-1 act on terminal in ring neural network prediction controller, obtain C[t ₁].Then with C[t ₁] act on terminal air-valve forecast model, obtain predicting air quantity

Keep terminal interior ring neural network prediction controller weights constant, will

S[t ₁] and-1 act on terminal in ring neural network prediction controller, obtain C ' [t ₁+ 1].With C ' [t ₁+ 1] the terminal air-valve forecast model of input obtains

Will

S[t ₁] and-1 act on terminal in ring neural network prediction controller, obtain C ' [t ₁+ 2]; With C ' [t ₁+ 2] the terminal air-valve forecast model of input obtains

The data that calculated are preserved.Make λ _i[k+M _i]=0, from after the λ the calculating formula (3) respectively forward _i[k] and q _i[k]:

$q_i\left[k\right]=\frac{\&PartialD;f_v{\left(k\right)}^T}{\&PartialD;C\left[k\right]}{\mathrm{\&lambda;}}_i[k+1]+\frac{\&PartialD;L_i{\left(k\right)}^T}{\&PartialD;C\left[k\right]}$

${\mathrm{\&lambda;}}_i\left[k\right]=\frac{\&PartialD;f_v{\left(k\right)}^T}{\&PartialD;y\left[k\right]}{\mathrm{\&lambda;}}_i[k+1]+\frac{\&PartialD;L_i{\left(k\right)}^T}{\&PartialD;y\left[k\right]}+\frac{\&PartialD;g_i{(k;W_i)}^T}{\&PartialD;y\left[k\right]}{q}_i\left[k\right]---\left(3\right)$

F in the formula _v(k) the terminal air-valve forecast model of setting up before the representative; g _i(k; W _i) be terminal interior ring neural network prediction controller equation; According to the q that calculates _i[k], through type (4) and formula (5) are revised the weights of terminal interior ring neural network prediction controller:

$\mathrm{\&Delta;}{W}_i=-{\mathrm{\&mu;}}_i\frac{\&PartialD;g_i{(k,W_i)}^T}{\&PartialD;W_i}{q}_i\left[k\right]---\left(4\right)$

W _i＝W _i+ΔW _i (5)

W wherein _iBe the weights battle array of terminal interior ring neural network prediction controller, μ _iBe the right value update rate, μ _iSelect 0.05.Constantly revise the weights of terminal interior ring neural network prediction controller, until Δ W _i＜0.001;

With S[t ₁], Y[t ₁] and-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁].With u[t] act on terminal air-valve forecast model, obtain

Afterwards with x[t ₁],

Act on the air conditioning area forecast model, obtain

Keep terminal outer shroud neural network prediction controller weights constant, will x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 1]; Ring optimizing in carrying out again is with s ' [t ₁+ 1], With-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁+ 1].Again with u ' [t ₁+ 1] acts on terminal air-valve forecast model, obtain

Afterwards with x[t ₁+ 1] x[1],

Act on the air conditioning area forecast model, obtain

Will

x ^*With-1 act on terminal outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 2]; Ring optimizing in carrying out again is with s ' [t ₁+ 2],

With-1 act on ring neural network prediction controller in terminal after the optimizing, obtain u ' [t ₁+ 2].Again with u ' [t ₁+ 2] act on terminal air-valve forecast model, obtain

Afterwards, will

Act on the air conditioning area forecast model, obtain And the rest may be inferred, utilizes controller neutral net and object forecast model to extrapolate following u ' [t ₁+ i] and

Value, i=3 wherein ..., 6, and the data that calculated are preserved.Make λ _o[k+M _o]=0, from after the λ the calculating formula (6) respectively forward _o[k] and q _o[k]:

$q_o\left[k\right]=\frac{\&PartialD;f_z{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_o[k+1]+\frac{\&PartialD;L_o{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_o\left[k\right]=\frac{\&PartialD;f_z{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_o[k+1]+\frac{\&PartialD;L_o{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_o{(k;W_o)}^T}{\&PartialD;x\left[k\right]}{q}_o\left[k\right]---\left(6\right)$

F in the formula _z(k) the air conditioning area forecast model of setting up before the representative, g _o(k; W _o) be terminal outer shroud neural network prediction controller equation.According to the q that calculates _o[k] is for k=t ₁+ M _o-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (7) and formula (8) are revised the weights of terminal outer shroud neural network prediction controller:

$\mathrm{\&Delta;}{W}_o=-{\mathrm{\&mu;}}_o\underset{k=t_1}{\overset{t_1+M_o-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_o{(k,W_o)}^T}{\&PartialD;W_o}{q}_o\left[k\right]---\left(7\right)$

W _o＝W _o+ΔW _o (8)

W wherein _oBe the weights battle array of terminal outer shroud neural network prediction controller, μ _oBe the right value update rate, μ _oSelect 0.05.Constantly revise the weights of terminal outer shroud neural network prediction controller, until Δ W _o＜0.001; After the outer shroud optimizing finishes, with x[t ₁] and x ^*Import terminal outer shroud neural network prediction controller, again with the S[t that obtains ₁] and Y[t ₁] the terminal interior ring neural network prediction controller of input, will export u[t ₁] directly act on terminal air-valve;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

5. a kind of VAV air conditioning system with variable control method according to claim 1; It is characterized in that, the 3rd step 3) be specially: each of initialization static pressure neural network prediction controller connects weights, and assignment is [1; 1] small random number in the scope, computing controller output then; Prediction time domain M _sIt was 2 steps, for the right value update rate μ of static pressure neural network prediction controller _sSelect 0.05, predetermined period was got 5 minutes;

Searching process: establish x[k] parameter be duct static pressure, u _Fan[k] is the blower fan controlled quentity controlled variable, x ^*Be to set duct static pressure.With x[t ₁], x ^*With-1 act on static pressure neural network prediction controller, controlled amount u _Fan[t ₁], then with u _Fan[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u _Fan[t ₁], x[t ₁] and-1 act on main air duct duct static pressure forecast model, obtain

Keep static pressure neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain

Will

X[t ₁+ 1] and-1 act on main air duct duct static pressure forecast model, obtains

The data that calculated are preserved.Make λ _s[k+M _s]=0, from after the λ the calculating formula (9) respectively forward _s[k] and q _s[k]:

$q_s\left[k\right]=\frac{\&PartialD;f_s{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}{\mathrm{\&lambda;}}_s[k+1]+\frac{\&PartialD;L_s{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}$

${\mathrm{\&lambda;}}_s\left[k\right]=\frac{\&PartialD;f_s{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_s[k+1]+\frac{\&PartialD;L_s{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_s{(k;W_s)}^T}{\&PartialD;x\left[k\right]}{q}_s\left[k\right]---\left(9\right)$

F in its Chinese style _s(k) the main air duct duct static pressure forecast model of setting up before the representative, g _s(k; W _s) be static pressure neural network prediction controller equation.According to the q that calculates _s[k] is for k=t ₁+ M _s-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (10) and formula (11) are revised the weights of static pressure neural network prediction controller:

$\mathrm{\&Delta;}{W}_s=-{\mathrm{\&mu;}}_s\underset{k=t_1}{\overset{t_1+M_s-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_s{(k,W_s)}^T}{\&PartialD;W_s}{q}_s\left[k\right]---\left(10\right)$

W _s＝W _s+ΔW _s (11)

W wherein _sBe the weights battle array of static pressure neural network prediction controller, μ _sBe the right value update rate, μ _sSelect 0.05.Constantly revise the weights of static pressure neural network prediction controller, until Δ W _s＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

6. a kind of VAV air conditioning system with variable control method according to claim 1; It is characterized in that, the 3rd step 4) be specially: each of initialization total blast volume neural network prediction controller connects weights, and assignment is [1; 1] small random number in the scope, computing controller output then; Prediction time domain M _fIt was 2 steps, for the right value update rate μ of total blast volume neural network prediction controller _fSelect 0.05, predetermined period was got 5 minutes;

Searching process: establish x[k] parameter be the main air duct air quantity, u _Fan[k] is the blower fan controlled quentity controlled variable, x ^*It is each regional air quantity setting value sum.With x[t ₁], x ^*With-1 act on total blast volume neural network prediction controller, controlled amount u _Fan[t ₁], then with u _Fan[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u _Fan[t ₁], x[t ₁] and-1 act on main air duct air quantity forecast model, obtain

Keep total blast volume neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain

Will

X[t ₁+ 1] and-1 act on main air duct air quantity forecast model, obtains The data that calculated are preserved.Make λ _f[k+M _f]=0, from after the λ the calculating formula (12) respectively forward _f[k] and q _f[k]:

$q_f\left[k\right]=\frac{\&PartialD;f_f{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}{\mathrm{\&lambda;}}_f[k+1]+\frac{\&PartialD;L_f{\left(k\right)}^T}{\&PartialD;u_{\mathrm{fan}}\left[k\right]}$

${\mathrm{\&lambda;}}_f\left[k\right]=\frac{\&PartialD;f_f{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_f[k+1]+\frac{\&PartialD;L_f{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_f{(k;W_f)}^T}{\&PartialD;x\left[k\right]}{q}_f\left[k\right]---\left(12\right)$

F wherein _f(k) the main air duct air quantity forecast model of setting up before the representative, g _f(k; W _f) be total blast volume neural network prediction controller equation.According to the q that calculates _f[k] is for k=t ₁+ M _f-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (13) and formula (14) are revised the weights of total blast volume neural network prediction controller:

$\mathrm{\&Delta;}{W}_f=-{\mathrm{\&mu;}}_f\underset{k=t_1}{\overset{t_1+M_f-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_f{(k,W_f)}^T}{\&PartialD;W_f}{q}_f\left[k\right]---\left(13\right)$

W _f＝W _f+ΔW _f (14)

W wherein _fBe the weights battle array of total blast volume neural network prediction controller, μ _fBe the right value update rate, μ _fSelect 0.05.Constantly revise the weights of total blast volume neural network prediction controller, until Δ W _f＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

7. a kind of VAV air conditioning system with variable control method according to claim 1; It is characterized in that, the 3rd step 5) be specially: each of initialization wind pushing temperature neural network prediction controller connects weights, and assignment is [1; 1] small random number in the scope, computing controller output then; Prediction time domain M _StIt was 2 steps, for the right value update rate μ of wind pushing temperature neural network prediction controller _StSelect 0.05, predetermined period was got 5 minutes;

Searching process: establish x[k] be wind pushing temperature, u[k] be the water valve aperture, x ^*Be to set wind pushing temperature.The control system detects all terminal VAV-BOX primary air flows, when a certain VAV-BOX primary air flow of appearance is lower than nominal air delivery 30%, reduces the water valve aperture, will set wind pushing temperature x ^*Improve 0.5 ℃; When a certain terminal VAV-BOX primary air flow greater than 70% the time, increase the water valve aperture, will set wind pushing temperature x ^*Reduce by 0.5 ℃.

With x[t ₁], x ^*With-1 act on wind pushing temperature neural network prediction controller, controlled amount u[t ₁], then with u[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u[t ₁], x[t ₁] and-1 act on the wind pushing temperature forecast model, obtain

Keep wind pushing temperature neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on controller, obtain u ' [t ₁+ 1], with u ' [t ₁+ 1], x[t ₁+ 1] and-1 act on the wind pushing temperature forecast model, obtains

The data that calculated are preserved.Make λ _St[k+M _St]=0, from after the λ the calculating formula (15) respectively forward _St[k] and q _St[k]:

$q_{\mathrm{st}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{st}}[k+1]+\frac{\&PartialD;L_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{st}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{st}}[k+1]+\frac{\&PartialD;L_{\mathrm{st}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_{\mathrm{st}}{(k;W_{\mathrm{st}})}^T}{\&PartialD;x\left[k\right]}{q}_{\mathrm{st}}\left[k\right]---\left(15\right)$

F wherein _St(k) the wind pushing temperature forecast model of setting up before the representative, g _St(k; W _St) be wind pushing temperature neural network prediction controller equation.According to the q that calculates _St[k] is for k=t ₁+ M _St-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (16) and formula (17) are revised the weights of wind pushing temperature neural network prediction controller:

$\mathrm{\&Delta;}{W}_{\mathrm{st}}=-{\mathrm{\&mu;}}_{\mathrm{st}}\underset{k=t_1}{\overset{t_1+M_{\mathrm{st}}-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_{\mathrm{st}}{(k,W_{\mathrm{st}})}^T}{\&PartialD;W_{\mathrm{st}}}{q}_{\mathrm{st}}\left[k\right]---\left(16\right)$

W _st＝W _st+ΔW _st (17)

W wherein _StBe the weights battle array of wind pushing temperature neural network prediction controller, μ _StBe the right value update rate, μ _StSelect 0.05.Constantly revise the weights of wind pushing temperature neural network prediction controller, until Δ W _St＜0.001;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

8. a kind of VAV air conditioning system with variable control method according to claim 1; It is characterized in that; The 3rd step 6) be specially ring and outer shroud neural network prediction controller in the initialization new wind ratio each be connected weights; Assignment is a small random number in [1,1] scope, computing controller output then; If output valve is all [0.7; 1] in the scope, then use this group controller initial weight, otherwise random initializtion again; Prediction time domain M _NiAnd M _qBe respectively 3 and 6, for the right value update rate μ selection 0.05 of ring and outer shroud neural network prediction controller in the new wind ratio, predetermined period was got 5 minutes;

Searching process: establish x[k] be k CO2 concentration constantly;

Be k forecast model output constantly, i.e. k+1 return air CO2 concentration constantly;

Be t ₁Return air CO2 concentration constantly; x ^*It is the CO2 concentration set point; S[k] set resh air requirement constantly for k; Y[k] be that k encircles the relevant state variables parameter constantly, comprise pipeline air quantity and duct static pressure; O[k] the constantly new wind valve area of the k that calculates for interior ring neural network prediction controller; U[k] be that interior ring neural network prediction controller optimizing finishes back k control output constantly, the i.e. new wind valve area that optimizing obtains;

The k+1 that exports for new wind air-valve forecast model predicts air quantity constantly.With x[t ₁], x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount S[t ₁], with S[t ₁], Y[t ₁] and-1 act on ring neural network prediction controller in the new wind ratio, obtain O[t ₁].Then with O[t ₁] act on new wind air-valve forecast model respectively, obtain predicting air quantity

It is constant to keep encircling neural network prediction controller weights in the new wind ratio, will

S[t ₁] and-1 act on ring neural network prediction controller in the new wind ratio, obtain O ' [t ₁+ 1].With O ' [t ₁+ 1] the new wind air-valve forecast model of input obtains

Will

S[t ₁] and-1 act on that ring neural network prediction controller obtains O ' [2] in the new wind ratio.The O '[2] Enter the new air damper prediction model, get

to the calculated data saved.Make λ _Ni[k+2]=0, from after the λ the calculating formula (18) respectively forward _Ni[k] and q _Ni[k]:

$q_{\mathrm{ni}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{nv}}{\left(k\right)}^T}{\&PartialD;O\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{ni}}[k+1]+\frac{\&PartialD;L_{\mathrm{ni}}{\left(k\right)}^T}{\&PartialD;O\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{ni}}\left[k\right]=\frac{\&PartialD;f_{\mathrm{nv}}{\left(k\right)}^T}{\&PartialD;y\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{ni}}[k+1]+\frac{\&PartialD;L_{\mathrm{ni}}{\left(k\right)}^T}{\&PartialD;y\left[k\right]}+\frac{\&PartialD;g_{\mathrm{ni}}{(k;W_{\mathrm{ni}})}^T}{\&PartialD;y\left[k\right]}{q}_{\mathrm{ni}}\left[k\right]---\left(18\right)$

F in the formula _Nv(k) the new wind air-valve forecast model of setting up before the representative; g _Ni(k; W _Ni) be ring neural network prediction controller equation in the new wind ratio; According to the q that calculates _Ni[k], through type (19) and formula (20) are revised the weights of ring neural network prediction controller in the new wind ratio:

$\mathrm{\&Delta;}{W}_{\mathrm{ni}}=-{\mathrm{\&mu;}}_{\mathrm{ni}}\frac{\&PartialD;g_{\mathrm{ni}}{(k,W_{\mathrm{ni}})}^T}{\&PartialD;W_{\mathrm{ni}}}{q}_{\mathrm{ni}}\left[k\right]---\left(19\right)$

W _ni＝W _ni+ΔW _ni (20)

W wherein _NiBe the weights battle array of ring neural network prediction controller in the new wind ratio, μ _NiBe the right value update rate, μ _NiSelect 0.05.Constantly revise the weights of ring neural network prediction controller in the new wind ratio, until Δ W _Ni＜0.001;

With S[t ₁], Y[t ₁] and-1 new wind ratio that acts on after the optimizing in ring neural network prediction controller, obtain u[t ₁].Then with u[t ₁] act on controlled device, obtain x[t ₁+ 1], again with u[t ₁] act on new wind air-valve forecast model, obtain Afterwards with x[t ₁],

Act on the air conditioning area forecast model, obtain

Keep new wind ratio outer shroud neural network prediction controller weights constant, with x[t ₁+ 1], x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 1]; Ring optimizing in carrying out again is with s ' [t ₁+ 1],

With ring neural network prediction controller in-1 new wind ratio that acts on after the optimizing, obtain u ' [t ₁+ 1].Again with u ' [t ₁+ 1] acts on new wind air-valve forecast model, obtain Afterwards, with x[t ₁+ 1], Act on the air conditioning area forecast model, obtain

Will

x ^*With-1 act on new wind ratio outer shroud neural network prediction controller, obtain set amount s ' [t ₁+ 2]; Ring optimizing in carrying out again is with s ' [t ₁+ 2],

With ring neural network prediction controller in-1 new wind ratio that acts on after the optimizing, obtain u ' [t ₁+ 2].Again with u ' [t ₁+ 2] act on new wind air-valve forecast model, obtain

Afterwards, will

Act on the air conditioning area forecast model, obtain

And the rest may be inferred, utilizes controller neutral net and object forecast model to extrapolate following u ' [t ₁+ i] and

Value, i=3 wherein ..., 6, and the data that calculated are preserved.Make λ _No[k+M _q]=0, from after the λ the calculating formula (21) respectively forward _No[k] and q _No[k]:

$q_{\mathrm{no}}\left[k\right]=\frac{\&PartialD;f_n{\left(k\right)}^T}{\&PartialD;u\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{no}}[k+1]+\frac{\&PartialD;L_{\mathrm{no}}{\left(k\right)}^T}{\&PartialD;u\left[k\right]}$

${\mathrm{\&lambda;}}_{\mathrm{no}}\left[k\right]=\frac{\&PartialD;f_n{\left(k\right)}^T}{\&PartialD;x\left[k\right]}{\mathrm{\&lambda;}}_{\mathrm{no}}[k+1]+\frac{\&PartialD;L_{\mathrm{no}}{\left(k\right)}^T}{\&PartialD;x\left[k\right]}+\frac{\&PartialD;g_{\mathrm{no}}{(k;W_{\mathrm{no}})}^T}{\&PartialD;x\left[k\right]}{q}_{\mathrm{no}}\left[k\right]---\left(21\right)$

F in the formula _n(k) the air quality forecast model of setting up before the representative, g _No(k; W _No) be new wind ratio outer shroud neural network prediction controller equation.According to the q that calculates _No[k] is for k=t ₁+ M _q-1 ..., t ₁+ 2, t ₁+ 1, t ₁, through type (22) and formula (23) are revised the weights of new wind ratio outer shroud neural network prediction controller:

$\mathrm{\&Delta;}{W}_{\mathrm{no}}=-{\mathrm{\&mu;}}_{\mathrm{no}}\underset{k=t_1}{\overset{t_1+M_q-1}{\mathrm{\&Sigma;}}}\frac{\&PartialD;g_{\mathrm{no}}{(k,W_{\mathrm{no}})}^T}{\&PartialD;W_{\mathrm{no}}}{q}_{\mathrm{no}}\left[k\right]---\left(22\right)$

W _no＝W _no+ΔW _no (23)

W wherein _NoBe the weights battle array of new wind ratio outer shroud neural network prediction controller, μ _NoBe the right value update rate, μ _NoSelect 0.05.Constantly revise the weights of new wind ratio outer shroud neural network prediction controller, until Δ W _No＜0.001; After the outer shroud optimizing finishes, with x[t ₁] and x ^*Input new wind ratio outer shroud neural network prediction controller is again with the S[t that obtains ₁] and Y[t ₁] the interior ring of input new wind ratio neural network prediction controller, will export u[t ₁] directly act on new wind air-valve;

The next sampling period then repeats aforesaid operations, and the value of each moment controlled quentity controlled variable finishes until control procedure after calculating respectively.

Images