# Study of Error Propagation in the Transformations of Dynamic Thermal Models of Buildings

Authors: Loïc Raillon; Christian Ghiaus
Journal: Journal of Control Science and Engineering
 (2017)

---

## Abstract

Dynamic behaviour of a system may be described by models with different forms: thermal (RC) networks, state-space representations, transfer functions, and ARX models. These models, which describe the same process, are used in the design, simulation, optimal predictive control, parameter identification, fault detection and diagnosis, and so on. Since more forms are available, it is interesting to know which one is the most suitable by estimating the sensitivity of the model to transform into a physical model, which is represented by a thermal network. A procedure for the study of error by Monte Carlo simulation and of factor prioritization is exemplified on a simple, but representative, thermal model of a building. The analysis of the propagation of errors and of the influence of the errors on the parameter estimation shows that the transformation from state-space representation to transfer function is more robust than the other way around. Therefore, if only one model is chosen, the state-space representation is preferable.


---

## Body

## 1. Introduction

Dynamic models for thermal behaviour of buildings are of interest for model predictive control (MPC) and measurement of performance. These applications use different representations (thermal networks, state-space models, autoregressive models, etc.) for the same physical system. However, the study of error propagation in the models due to the model transformations needs more clarifications.According to the literature of the last few years, MPC is a widespread strategy with promising energy saving potential [1–3]. The algorithm minimizes the energy consumption while optimizing the thermal comfort by incorporating weather forecast, occupancy schedule, and other factors. Since MPC is a model-based strategy, it means that the performances rely on the model accuracy; therefore the model must be capable of catching the building dynamics. The information on the model of the building may be used to tackle down the influence on the energy consumption of inputs such as weather and occupancy. This allows distinguishing the intrinsic part of the model from its disturbances. Therefore, it is interesting to keep the physical meaning in the model for two reasons: (1) calibration of model on actual data by using the parameter identification based on in situ measurements and (2) physical interpretation of the parameters of the building such as time constants and envelop thermal conductance. This strategy is attractive because nowadays we are not able to measure the energy efficiency of a building. Indeed actual methods present discrepancies from 50% to 200% between the estimated consumption and the real one [4–6]. In the conception phase, simulations are used to determine if a new or a refurbished building matches the performance specifications. If it is not the case, the building parameters are optimized until the requirements are met. Nevertheless these detailed simulations use conventional data such as weather records, occupancy schedules, and set-point temperatures, which can be far from the reality as pointed out by the magnitude of the gaps between the simulation results and the measurements. This problem in the conception phase cannot be fixed since it is not possible to have in situ measurement when the building is not even built. The same methodology is used to quantify the energy efficiency gained by retrofitting. It follows naturally the reassessment on the relevance of the predicted gains and how the return of investment can be certified. A solution could be to measure the energy consumption in a reference situation, before and after refurbishment and determine the savings [7]. During the baseline period, the model is calibrated on measurements. Thereafter, in the reporting period, the difference between the consumption measured by sensors and that estimated by the model represents the energy savings due to the refurbishment.The deadlock in these applications is the experimental identification of the model. Two model structures are widely used for control analysis and synthesis: state-space and transfer function. When the first is established from a thermal circuit, the relations between the physical parameters (i.e., thermal resistances,R, and capacities, C) are explicit contrary to the transfer function, where the coefficients have to follow some transformations to find the values of the physical parameters. Figure 1 summarises the feasible transformations and the process to follow. Both state-space representations and transfer functions are suitable for control strategy and fault detection and diagnosis. The question is which one is better since their parameters are identified by using different methods. Therefore we are interested in studying the propagation of error through the transformations from a form of the model to another.Figure 1
Transformation flow chart.According to Figure1, we consider three model types. For ease of notation, we refer to the thermal network by RC, the continuous state-space by SS, and the ARX model by ARX. For instance, if a thermal network is put in state-space representation, it is a RC to SS transformation.


## 2. Error Propagation by Monte Carlo Simulation

The transformations in Figure1 can be represented by a function of the form y=f(X), which relates k input parameters, X, from an initial representation, to the output parameter, y, of a new representation. In practice, the input parameter vector X is often estimated from noisy measurements and consequently is uncertain.How are the uncertainties of an estimated parameter propagated from one representation to another? Is it better to identify the parameters of a state-space model and then to have the possibility of converting it in a transfer function or to do the procedure the other way around? Which transformation is the most robust to uncertainties? The error propagation between model transformations is quantified by Monte Carlo simulations where the vector,X, is randomly chosen such that(1)X~Npin∗,δpin∗,where Nμ,σ2 is a normal distribution with mean μ and standard deviation σ, pin∗ is the true input parameters vector, and δ represents the percentage of pin∗ introduced as variance.Then, givenX, the output parameter y is obtained. This process is repeated N times, for N sufficiently large such that the input parameter space is well explored and the output parameter distribution is well defined.The transformations between different model representations are first presented, starting from a dynamic thermal network in which the parameters have a physical meaning.


### 2.1. Transformation between Thermal Network and State-Space Model

To illustrate the methodology presented in this section [8], the second-order thermal network in Figure 2 is considered. The heat transfers through the envelope and the ventilation are modelled by the following heat transfer rates:(i)
q1: infiltration and windows conduction (W)(ii)
q2: outside-air convection (W)(iii)
q3: wall conduction (W)(iv)
q4: wall conduction and inside-air convection (W).The thermal capacities Cw and Ca represent, respectively, the energy accumulated in the envelope and the energy accumulated in the medium and the indoor air.Figure 2
Thermal network for the heat balance method.A temperature source and two heat rate sources are considered as inputs:(i)
To: outdoor temperature (°C)(ii)
Q˙o: solar radiation on the outside-wall surfaces (W)(iii)
Q˙H+Q˙g: heat flow from the HVAC system and internal gains (W).For the thermal network in Figure 2, the temperature differences over each resistance can be written in a matrix form:(2)e=-Aθ+be1e2e3e4=010100-10101-1θsoθaθw+ToTo00,where e is the vector of temperature drops over the thermal resistances, θ is the vector of temperature in the nodes, and b is the vector of temperature sources on the branches.The heat transfer rates in the branches can be expressed in matrix form:(3)q=Geq1q2q3q4=Rv-10000Rco-10000Rw2-10000Rw2+Rci-1e1e2e3e4,where q is the vector of heat rates in the branches and G is the diagonal matrix of thermal conductivities.The heat balance in each temperature node gives(4)Cθ˙=Kθ+Kbb+fwhere K=-ATGA, Kb=ATG, C is the diagonal matrix of thermal capacities, and f is the vector of heat rate sources connected to the temperature nodes.The matrixC partitions the set of equations (4) into nodes with thermal capacities θ˙s and nodes without thermal capacities θ˙0:(5)000Ccθ˙0θ˙s=K11K12K21K22θ0θs+Kb1Kb2b+I1100I22f0fc.The heat balance equation can be transformed in state-space model by eliminating the nodes without thermal capacity and therefore separating differential and algebraic equations.The state-space model of the thermal network is then given by(6)θ˙s=Asθs+Bsu,where(7)As=Cc-1-K21K11-1K12+K22isthestatematrix(8)Bs=Cc-1-K21K11-1Kb1+Kb2-K21K11-1I22istheinputmatrix(9)u=ToTo00Q˙oQ˙H+Q˙g0T.The input matrix (8) can be reduced by taking only the columns corresponding to nonzero elements in the input vector (9) and by adding the first and the second column which correspond to the same input, To, which gives(10)As=-2Rci+2Rv+RwCaRv2Rci+Rw2Ca2Rci+Rw2Cw2Rci+Rw-4Rci+Rco+RwCw2Rci+Rw2Rco+RwBs=1CaRv01Ca2Cw2Rco+Rw2RcoCw2Rco+Rw0.This transformation maps six physical parameters into eight state-space parameters, which means that two parameters are dependent on the others, such as a22=-a21-b21 and a11=-b11-a12, where aij and bij represent, respectively, the parameter of the state matrix As and input matrix Bs in row i and column j. Reciprocally, two parameters which do not act directly on the solution have to be removed to balance the system. Here, a12 and b21 are not used, which gives the following expressions: (11)Rco=-b22a22a21Rci=-b11+a11b22+a22+2a21b13a22+a21b11+a11a22+a11a21Rv=b13b11Rw=2b22b11+a11+2a21b13b11a22+a21+a11a22+a11a21Cw=-b11+a11a21b13Ca=1b13.Finally, the state-space model (10) is discretized with the forward Euler method to conserve the explicit relations between the physical parameters(12)Ad=I+AstsBd=Bsts,where ts is the sampling time; ts=60s is considered for the rest of the paper.


### 2.2. Transformation between State-Space Model and ARX Model

The purpose of this section is to find a relationship between parameters of the state-space model and parameters of the discrete transfer function. Although there is a well-known formula to pass from state-space to transfer function, reciprocally, the latter is not unique since the input-output relationship is invariant to similarity transformation [9]. The state-space model (12) can be transformed in a new representation (13)θ˙c=Acθc+Bcuif a new state vector is chosen according to the following relation:(14)θc=T-1θd,where T is an invertible matrix and θd is the state vector of the discretized state-space.The new representation is linked to the original state-space model by(15)Ac=T-1AdT(16)Bc=T-1Bd(17)Oc=OdT,where O is the observability matrix computed by(18)Os=CsCsAsCsAs2⋯CsAsn-1Twith n, the number of states.One of the easiest ways to get an input-output representation is to describe it as a linear difference equation in which the next output depends on previous observations (outputs and inputs) [10]. Since observed data are always collected by sampling, only the discrete form is considered. The system can be represented by a discrete transfer function, where the numerator and the denominator are polynomial in q-1, the shift operator; this type of model is called autoregressive with exogenous inputs (ARX).A multi-input multioutput (MIMO) system can be described by a difference equation [11]:(19)yk+α1yk-1+α2yk-2+⋯+αnyk-n=β1uk-1+β2uk-2+⋯+βnuk-n,where u∈Rm is the input, y∈Rp is the output, y(k) is the output at time instant k, yk-n is the output delayed of n samples, αi∈Rp×p, and βi∈Rp×m.By using the properties of the shift operatorq-1yk=yk-1, (19) is rewritten:(20)I+α1q-1+α2q-2+⋯+αnq-nyk=β1q-1+β2q-2+⋯+βnq-nuk.By defining(21)αq=I+α1q-1+α2q-2+⋯+αnq-nβq=β1q-1+β2q-2+⋯+βnq-nthe difference equation (20) is transformed in ARX model with(22)ARXq=ykuk=α-1qβq.One solution to convert a state-space model in ARX model and vice versa could be to get the transfer function from the continuous state-space and then discretize it or discretize the state-space before the transformation, which gives the structure to impose in the ARX model. Relations between both transfer functions can be found afterwards.This technique presents some limitations. First the transfer function must be found to know the ARX’s structure; the number of ARX parameters might not be equal to the number of the state-space ones. This is the case for system (10) where the state-space model has eight parameters and the ARX model has seven parameters. An internal relation in the state-space model has to be found to balance the system before solving it; this method is not elegant and not unique as mentioned before. In order to keep the physical meaning of the state-space, the ARX model must be transformed in the same system of coordinates as the state-space model. This is possible if the state-space model is transformed in observability canonical form, where the ARX model coefficients appear explicitly. The observability canonical form has the following structure:(23)Ac=0I0⋯000I⋱⋮⋮⋮⋱⋱000⋯0I-αn-αn-1…-α2-α1∈Rpn×pnBc=γ1γ2⋯γnT∈Rpn×mwith γi∈Rp×m given by(24)γ1γ2⋮γn=Iα1Iα2α1I⋮⋱⋱-αn-1-αn-2…α1I-1β1β2⋮βn.The canonical state-space (23) can be projected in the physical state-space coordinates by using similarity transformation. The observability matrices of the physical state-space and the canonical state-space model are linked by relation (17). Furthermore the observability matrix of the canonical state-space is the identity matrix; therefore the similarity matrix reduces to(25)T=Os-1.By using the knowledge of the physical state-space, a transformation which projects a realization in a given coordinates system and then keeps the physical meaning can be found. Moreover the transformation works for MIMO systems and stays straightforward even if the order of the model increases because only matrix operations are involved. In the next section the error propagation through model conversion is studied in order to test the robustness of the method.


### 2.3. Simulation for Error Propagation

The input parameterspin are normally distributed using (1) with δ=0.1 for the transformations from RC to ARX and δ=0.01 in the other way. The following results are given for N=106 simulations.Three criteria are used to characterize the parameter transformations: the expected value of the output parameter distributionEpout which gives the location of the highest occurrence; the absolute mean error (MAE) between the true output parameter values pout∗; and the simulated output parameters pout such that(26)MAE=1N∑i=1Npouti-pout∗.The Kolmogorov-Smirnov test [12] is used to know if the output parameters are normally distributed. To determine if it is the case, the test finds the maximum distance d between two curves:(27)d=maxx⁡Fx-Sx,where d is the absolute difference between Fx, the empirical cumulative distribution function calculated from x and S(x), and the cumulative distribution function (CDF) of the hypothetical distribution. d is then compared to the critical level dα, the probability of rejecting the hypothesis when it is true (α=0.05 gives dα=1.36/N=1.36e-3). dα<d indicates the rejection of the hypothesis at significance level α and dα>d indicates a failure to reject the hypothesis at significance level α, which means that the data belong to the hypothetical distribution. The Kolmogorov-Smirnov test is illustrated in Figure 3 to check if the parameter a11 is normally distributed. The red line is the empirical cumulative distribution function of the data and the black dotted lines represent the critical level boundaries. If Fx crosses a critical level boundary, it means that the test fails to reject the null hypothesis; that is, the data set is not normally distributed. N=104 has been used to widen the critical level boundaries only for the illustration.Figure 3
Kolmogorov-Smirnov test for parametera11 (RC to SS transformation, α=0.05,N=104).The simulation results for the six transformations are given in Figure4 and Table 1 for the RC to SS transformation, in Figure 5 and Table 2 for the SS to ARX transformation, in Figure 6 and Table 3 for the RC to ARX transformation, in Figure 7 and Table 4 for the ARX to SS transformation, in Figure 8 and Table 5 for the SS to RC transformation, and in Figure 9 and Table 6 for the ARX to RC transformation. If histograms are not centred on x-axis scale, it means that outliers are present but are not visible because of the low occurrence.Table 1
RC to SS transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
a
11
-
5.02
e
-
4
-
5.04
e
-
4
5.66
e
-
5
4.27
e
-
2
a
12
3.46
e
-
5
3.48
e
-
5
3.98
e
-
6
4.29
e
-
2
a
21
7.09
e
-
7
7.13
e
-
7
8.16
e
-
8
4.30
e
-
2
a
22
-
1.43
e
-
6
-
1.44
e
-
6
1.65
e
-
7
4.32
e
-
2
b
11
4.67
e
-
4
4.71
e
-
4
5.44
e
-
5
4.34
e
-
2
b
13
1.22
e
-
5
1.21
e
-
5
9.93
e
-
7
4.07
e
-
2
b
21
7.17
e
-
7
7.23
e
-
7
8.30
e
-
8
4.33
e
-
2
b
22
2.87
e
-
9
2.80
e
-
9
4.04
e
-
10
3.92
e
-
2Table 2
SS to TF transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
n
11
2.80
e
-
2
2.80
e
-
2
2.24
e
-
3
6.75
e
-
4
n
12
-
2.80
e
-
2
-
2.80
e
-
2
2.24
e
-
3
6.69
e
-
4
n
22
3.57
e
-
10
3.49
e
-
10
4.04
e
-
11
1.44
e
-
2
n
31
7.32
e
-
4
7.25
e
-
4
5.84
-
5
5.22
e
-
4
n
32
-
7.32
e
-
4
-
7.25
e
-
4
5.84
-
5
5.19
e
-
4
d
1
-
1.97
-
1.97
2.40
e
-
3
5.34
e
-
4
d
2
9.70
e
-
1
9.70
e
-
1
2.39
e
-
3
5.33
e
-
4Table 3
RC to ARX transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
n
11
2.80
e
-
2
2.72
e
-
2
3.27
e
-
3
4.35
e
-
2
n
12
-
2.80
e
-
2
-
2.72
e
-
2
3.27
e
-
3
4.35
e
-
2
n
22
3.57
e
-
10
3.64
e
-
10
7.87
e
-
11
6.38
e
-
2
n
31
7.32
e
-
4
7.25
e
-
4
5.96
-
5
4.09
e
-
2
n
32
-
7.32
e
-
4
-
7.25
e
-
4
5.96
-
5
4.09
e
-
2
d
1
-
1.97
-
1.97
3.40
e
-
3
4.28
e
-
2
d
2
9.70
e
-
1
9.70
e
-
1
3.39
e
-
3
4.28
e
-
2Table 4
ARX to SS transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
a
21
7.09
e
-
7
-
7.50
e
-
3
1.37
e
-
1
6.09
e
-
4
a
22
-
1.43
e
-
6
2.50
e
-
5
2.62
e
-
4
8.19
e
-
4
b
11
4.67
e
-
4
4.68
e
-
4
3.73
e
-
6
1.11
e
-
3
b
13
1.22
e
-
5
1.22
e
-
5
9.74
e
-
8
4.07
e
-
4
b
21
7.17
e
-
7
-
5.00
e
-
4
4.35
e
-
3
1.18
e
-
3
b
22
2.87
e
-
9
2.87
e
-
9
2.29
e
-
11
8.70
e
-
4Table 5
SS to RC transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
R
c
o
4.00
e
-
3
3.99
e
-
3
7.78
e
-
5
8.02
e
-
3
R
c
i
8.00
e
-
3
9.15
e
-
3
8.15
e
-
3
1.42
e
-
2
R
v
2.61
e
-
2
2.61
e
-
2
2.90
e
-
4
3.37
e
-
3
R
w
6.89
e
-
1
6.33
e
-
1
1.22
e
-
1
8.96
e
-
2
C
w
4.00
e
6
4.04
e
6
6.35
e
5
1.09
e
-
3
C
a
8.20
e
4
8.19
e
4
6.55
e
2
4.21
e
-
3Table 6
ARX to RC transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
R
c
o
4.00
e
-
3
-
1
e
-
6
0
4.00
e
-
3
4.84
e
-
1
R
c
i
8.00
e
-
3
6.85
e
-
1
7.12
e
-
1
3.24
e
-
1
R
v
2.61
e
-
2
2.61
e
-
2
2.99
e
-
4
2.92
e
-
3
R
w
6.89
e
-
1
-
6.75
e
-
1
1.41
4.00
e
-
1
C
w
4.00
e
6
0
300
4.00
e
6
4.66
e
-
1
C
a
8.20
e
4
8.20
e
4
6.54
e
2
4.22
e
-
3Figure 4
Histograms, RC to SS transformation.Figure 5
Histograms, SS to ARX transformation.Figure 6
Histograms, RC to ARX transformation.Figure 7
Histograms, ARX to SS transformation.Figure 8
Histograms, SS to RC transformation.Figure 9
Histograms, ARX to RC transformation.The three transformations from RC to ARX have expected values very close to the true parameter valuespout∗ and all the absolute mean errors are approximately around 10% of pout∗, which corresponds to the variance introduced to the input parameters pin. The majority of the distributions are a little bit skewed but close to the normal distribution. The parameter n22 in the RC to ARX transformation (Figure 6, Table 3) is the worst case with a mean absolute error superior to 22%, which is still a satisfactory result. Indeed, by looking at the analytical expression, n22 depends on fives parameters with 10% of their values as variances. Furthermore, it is more difficult to obtain accurate results with very small values like n22=3.57·10-10.In the other way, from ARX to RC, the transformations are more sensitive to parameter deviations; this is why the variance of (1) has been scaled with δ=0.01 instead of δ=0.1 as previously. The histograms of Rco, Rci, and Rw in Figure 8 are zoomed to see most of the distribution. Otherwise, outliers far from the true parameter values give an improper graphical scale. The same situation occurs for Rco, Rci, Rw, and Cw in Figure 9.From ARX to SS, the parametersa11 and a12 are not present in the results (Figure 7 and Table 4) because they are not a function of the ARX parameters; they are known directly from the similarity matrix T. In this transformation, all the parameters are normally distributed but the mean absolute errors of a21, a22, and b21 are far greater than their respective true values. These transformations depend on the poles d1 and d2 which are several orders of magnitude higher than the output parameters; consequently uncertainties introduced in the input parameters have strong impacts on the output parameters.The transformations from SS to RC involve only operations with values relatively small and with different orders of magnitude which explains the skewed forms ofRco, Rci, and Rw (Figure 8) and the presence of outliers. Moreover, Rci and Rw have complicated transformations involving six input parameters.The problem gets even worse when both transformations are combined in the ARX to RC transformation. The expected values ofRco, Rci, Rw, and Cw are far from their respective true values and the normal distribution shapes of Rco and Cw are completely lost and therefore cannot be shown properly by histograms. Only two parameters out of six are robust against uncertainties and some R and C parameters have negative values which has no physical meaning. These results show that the transformation from ARX to RC is not suitable even if only 1% of the parameter values are introduced as variance.In this section, the most sensitive relations to parameter deviations have been located. Sometimes, functions relating input and output parameters are complicated and it is hard to identify which parameters are most responsible for the output variance.


## 2.1. Transformation between Thermal Network and State-Space Model

To illustrate the methodology presented in this section [8], the second-order thermal network in Figure 2 is considered. The heat transfers through the envelope and the ventilation are modelled by the following heat transfer rates:(i)
q1: infiltration and windows conduction (W)(ii)
q2: outside-air convection (W)(iii)
q3: wall conduction (W)(iv)
q4: wall conduction and inside-air convection (W).The thermal capacities Cw and Ca represent, respectively, the energy accumulated in the envelope and the energy accumulated in the medium and the indoor air.Figure 2
Thermal network for the heat balance method.A temperature source and two heat rate sources are considered as inputs:(i)
To: outdoor temperature (°C)(ii)
Q˙o: solar radiation on the outside-wall surfaces (W)(iii)
Q˙H+Q˙g: heat flow from the HVAC system and internal gains (W).For the thermal network in Figure 2, the temperature differences over each resistance can be written in a matrix form:(2)e=-Aθ+be1e2e3e4=010100-10101-1θsoθaθw+ToTo00,where e is the vector of temperature drops over the thermal resistances, θ is the vector of temperature in the nodes, and b is the vector of temperature sources on the branches.The heat transfer rates in the branches can be expressed in matrix form:(3)q=Geq1q2q3q4=Rv-10000Rco-10000Rw2-10000Rw2+Rci-1e1e2e3e4,where q is the vector of heat rates in the branches and G is the diagonal matrix of thermal conductivities.The heat balance in each temperature node gives(4)Cθ˙=Kθ+Kbb+fwhere K=-ATGA, Kb=ATG, C is the diagonal matrix of thermal capacities, and f is the vector of heat rate sources connected to the temperature nodes.The matrixC partitions the set of equations (4) into nodes with thermal capacities θ˙s and nodes without thermal capacities θ˙0:(5)000Ccθ˙0θ˙s=K11K12K21K22θ0θs+Kb1Kb2b+I1100I22f0fc.The heat balance equation can be transformed in state-space model by eliminating the nodes without thermal capacity and therefore separating differential and algebraic equations.The state-space model of the thermal network is then given by(6)θ˙s=Asθs+Bsu,where(7)As=Cc-1-K21K11-1K12+K22isthestatematrix(8)Bs=Cc-1-K21K11-1Kb1+Kb2-K21K11-1I22istheinputmatrix(9)u=ToTo00Q˙oQ˙H+Q˙g0T.The input matrix (8) can be reduced by taking only the columns corresponding to nonzero elements in the input vector (9) and by adding the first and the second column which correspond to the same input, To, which gives(10)As=-2Rci+2Rv+RwCaRv2Rci+Rw2Ca2Rci+Rw2Cw2Rci+Rw-4Rci+Rco+RwCw2Rci+Rw2Rco+RwBs=1CaRv01Ca2Cw2Rco+Rw2RcoCw2Rco+Rw0.This transformation maps six physical parameters into eight state-space parameters, which means that two parameters are dependent on the others, such as a22=-a21-b21 and a11=-b11-a12, where aij and bij represent, respectively, the parameter of the state matrix As and input matrix Bs in row i and column j. Reciprocally, two parameters which do not act directly on the solution have to be removed to balance the system. Here, a12 and b21 are not used, which gives the following expressions: (11)Rco=-b22a22a21Rci=-b11+a11b22+a22+2a21b13a22+a21b11+a11a22+a11a21Rv=b13b11Rw=2b22b11+a11+2a21b13b11a22+a21+a11a22+a11a21Cw=-b11+a11a21b13Ca=1b13.Finally, the state-space model (10) is discretized with the forward Euler method to conserve the explicit relations between the physical parameters(12)Ad=I+AstsBd=Bsts,where ts is the sampling time; ts=60s is considered for the rest of the paper.


## 2.2. Transformation between State-Space Model and ARX Model

The purpose of this section is to find a relationship between parameters of the state-space model and parameters of the discrete transfer function. Although there is a well-known formula to pass from state-space to transfer function, reciprocally, the latter is not unique since the input-output relationship is invariant to similarity transformation [9]. The state-space model (12) can be transformed in a new representation (13)θ˙c=Acθc+Bcuif a new state vector is chosen according to the following relation:(14)θc=T-1θd,where T is an invertible matrix and θd is the state vector of the discretized state-space.The new representation is linked to the original state-space model by(15)Ac=T-1AdT(16)Bc=T-1Bd(17)Oc=OdT,where O is the observability matrix computed by(18)Os=CsCsAsCsAs2⋯CsAsn-1Twith n, the number of states.One of the easiest ways to get an input-output representation is to describe it as a linear difference equation in which the next output depends on previous observations (outputs and inputs) [10]. Since observed data are always collected by sampling, only the discrete form is considered. The system can be represented by a discrete transfer function, where the numerator and the denominator are polynomial in q-1, the shift operator; this type of model is called autoregressive with exogenous inputs (ARX).A multi-input multioutput (MIMO) system can be described by a difference equation [11]:(19)yk+α1yk-1+α2yk-2+⋯+αnyk-n=β1uk-1+β2uk-2+⋯+βnuk-n,where u∈Rm is the input, y∈Rp is the output, y(k) is the output at time instant k, yk-n is the output delayed of n samples, αi∈Rp×p, and βi∈Rp×m.By using the properties of the shift operatorq-1yk=yk-1, (19) is rewritten:(20)I+α1q-1+α2q-2+⋯+αnq-nyk=β1q-1+β2q-2+⋯+βnq-nuk.By defining(21)αq=I+α1q-1+α2q-2+⋯+αnq-nβq=β1q-1+β2q-2+⋯+βnq-nthe difference equation (20) is transformed in ARX model with(22)ARXq=ykuk=α-1qβq.One solution to convert a state-space model in ARX model and vice versa could be to get the transfer function from the continuous state-space and then discretize it or discretize the state-space before the transformation, which gives the structure to impose in the ARX model. Relations between both transfer functions can be found afterwards.This technique presents some limitations. First the transfer function must be found to know the ARX’s structure; the number of ARX parameters might not be equal to the number of the state-space ones. This is the case for system (10) where the state-space model has eight parameters and the ARX model has seven parameters. An internal relation in the state-space model has to be found to balance the system before solving it; this method is not elegant and not unique as mentioned before. In order to keep the physical meaning of the state-space, the ARX model must be transformed in the same system of coordinates as the state-space model. This is possible if the state-space model is transformed in observability canonical form, where the ARX model coefficients appear explicitly. The observability canonical form has the following structure:(23)Ac=0I0⋯000I⋱⋮⋮⋮⋱⋱000⋯0I-αn-αn-1…-α2-α1∈Rpn×pnBc=γ1γ2⋯γnT∈Rpn×mwith γi∈Rp×m given by(24)γ1γ2⋮γn=Iα1Iα2α1I⋮⋱⋱-αn-1-αn-2…α1I-1β1β2⋮βn.The canonical state-space (23) can be projected in the physical state-space coordinates by using similarity transformation. The observability matrices of the physical state-space and the canonical state-space model are linked by relation (17). Furthermore the observability matrix of the canonical state-space is the identity matrix; therefore the similarity matrix reduces to(25)T=Os-1.By using the knowledge of the physical state-space, a transformation which projects a realization in a given coordinates system and then keeps the physical meaning can be found. Moreover the transformation works for MIMO systems and stays straightforward even if the order of the model increases because only matrix operations are involved. In the next section the error propagation through model conversion is studied in order to test the robustness of the method.


## 2.3. Simulation for Error Propagation

The input parameterspin are normally distributed using (1) with δ=0.1 for the transformations from RC to ARX and δ=0.01 in the other way. The following results are given for N=106 simulations.Three criteria are used to characterize the parameter transformations: the expected value of the output parameter distributionEpout which gives the location of the highest occurrence; the absolute mean error (MAE) between the true output parameter values pout∗; and the simulated output parameters pout such that(26)MAE=1N∑i=1Npouti-pout∗.The Kolmogorov-Smirnov test [12] is used to know if the output parameters are normally distributed. To determine if it is the case, the test finds the maximum distance d between two curves:(27)d=maxx⁡Fx-Sx,where d is the absolute difference between Fx, the empirical cumulative distribution function calculated from x and S(x), and the cumulative distribution function (CDF) of the hypothetical distribution. d is then compared to the critical level dα, the probability of rejecting the hypothesis when it is true (α=0.05 gives dα=1.36/N=1.36e-3). dα<d indicates the rejection of the hypothesis at significance level α and dα>d indicates a failure to reject the hypothesis at significance level α, which means that the data belong to the hypothetical distribution. The Kolmogorov-Smirnov test is illustrated in Figure 3 to check if the parameter a11 is normally distributed. The red line is the empirical cumulative distribution function of the data and the black dotted lines represent the critical level boundaries. If Fx crosses a critical level boundary, it means that the test fails to reject the null hypothesis; that is, the data set is not normally distributed. N=104 has been used to widen the critical level boundaries only for the illustration.Figure 3
Kolmogorov-Smirnov test for parametera11 (RC to SS transformation, α=0.05,N=104).The simulation results for the six transformations are given in Figure4 and Table 1 for the RC to SS transformation, in Figure 5 and Table 2 for the SS to ARX transformation, in Figure 6 and Table 3 for the RC to ARX transformation, in Figure 7 and Table 4 for the ARX to SS transformation, in Figure 8 and Table 5 for the SS to RC transformation, and in Figure 9 and Table 6 for the ARX to RC transformation. If histograms are not centred on x-axis scale, it means that outliers are present but are not visible because of the low occurrence.Table 1
RC to SS transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
a
11
-
5.02
e
-
4
-
5.04
e
-
4
5.66
e
-
5
4.27
e
-
2
a
12
3.46
e
-
5
3.48
e
-
5
3.98
e
-
6
4.29
e
-
2
a
21
7.09
e
-
7
7.13
e
-
7
8.16
e
-
8
4.30
e
-
2
a
22
-
1.43
e
-
6
-
1.44
e
-
6
1.65
e
-
7
4.32
e
-
2
b
11
4.67
e
-
4
4.71
e
-
4
5.44
e
-
5
4.34
e
-
2
b
13
1.22
e
-
5
1.21
e
-
5
9.93
e
-
7
4.07
e
-
2
b
21
7.17
e
-
7
7.23
e
-
7
8.30
e
-
8
4.33
e
-
2
b
22
2.87
e
-
9
2.80
e
-
9
4.04
e
-
10
3.92
e
-
2Table 2
SS to TF transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
n
11
2.80
e
-
2
2.80
e
-
2
2.24
e
-
3
6.75
e
-
4
n
12
-
2.80
e
-
2
-
2.80
e
-
2
2.24
e
-
3
6.69
e
-
4
n
22
3.57
e
-
10
3.49
e
-
10
4.04
e
-
11
1.44
e
-
2
n
31
7.32
e
-
4
7.25
e
-
4
5.84
-
5
5.22
e
-
4
n
32
-
7.32
e
-
4
-
7.25
e
-
4
5.84
-
5
5.19
e
-
4
d
1
-
1.97
-
1.97
2.40
e
-
3
5.34
e
-
4
d
2
9.70
e
-
1
9.70
e
-
1
2.39
e
-
3
5.33
e
-
4Table 3
RC to ARX transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
n
11
2.80
e
-
2
2.72
e
-
2
3.27
e
-
3
4.35
e
-
2
n
12
-
2.80
e
-
2
-
2.72
e
-
2
3.27
e
-
3
4.35
e
-
2
n
22
3.57
e
-
10
3.64
e
-
10
7.87
e
-
11
6.38
e
-
2
n
31
7.32
e
-
4
7.25
e
-
4
5.96
-
5
4.09
e
-
2
n
32
-
7.32
e
-
4
-
7.25
e
-
4
5.96
-
5
4.09
e
-
2
d
1
-
1.97
-
1.97
3.40
e
-
3
4.28
e
-
2
d
2
9.70
e
-
1
9.70
e
-
1
3.39
e
-
3
4.28
e
-
2Table 4
ARX to SS transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
a
21
7.09
e
-
7
-
7.50
e
-
3
1.37
e
-
1
6.09
e
-
4
a
22
-
1.43
e
-
6
2.50
e
-
5
2.62
e
-
4
8.19
e
-
4
b
11
4.67
e
-
4
4.68
e
-
4
3.73
e
-
6
1.11
e
-
3
b
13
1.22
e
-
5
1.22
e
-
5
9.74
e
-
8
4.07
e
-
4
b
21
7.17
e
-
7
-
5.00
e
-
4
4.35
e
-
3
1.18
e
-
3
b
22
2.87
e
-
9
2.87
e
-
9
2.29
e
-
11
8.70
e
-
4Table 5
SS to RC transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
R
c
o
4.00
e
-
3
3.99
e
-
3
7.78
e
-
5
8.02
e
-
3
R
c
i
8.00
e
-
3
9.15
e
-
3
8.15
e
-
3
1.42
e
-
2
R
v
2.61
e
-
2
2.61
e
-
2
2.90
e
-
4
3.37
e
-
3
R
w
6.89
e
-
1
6.33
e
-
1
1.22
e
-
1
8.96
e
-
2
C
w
4.00
e
6
4.04
e
6
6.35
e
5
1.09
e
-
3
C
a
8.20
e
4
8.19
e
4
6.55
e
2
4.21
e
-
3Table 6
ARX to RC transformation characteristics.
p
o
u
t
∗
E
p
o
u
t
MAE
d
R
c
o
4.00
e
-
3
-
1
e
-
6
0
4.00
e
-
3
4.84
e
-
1
R
c
i
8.00
e
-
3
6.85
e
-
1
7.12
e
-
1
3.24
e
-
1
R
v
2.61
e
-
2
2.61
e
-
2
2.99
e
-
4
2.92
e
-
3
R
w
6.89
e
-
1
-
6.75
e
-
1
1.41
4.00
e
-
1
C
w
4.00
e
6
0
300
4.00
e
6
4.66
e
-
1
C
a
8.20
e
4
8.20
e
4
6.54
e
2
4.22
e
-
3Figure 4
Histograms, RC to SS transformation.Figure 5
Histograms, SS to ARX transformation.Figure 6
Histograms, RC to ARX transformation.Figure 7
Histograms, ARX to SS transformation.Figure 8
Histograms, SS to RC transformation.Figure 9
Histograms, ARX to RC transformation.The three transformations from RC to ARX have expected values very close to the true parameter valuespout∗ and all the absolute mean errors are approximately around 10% of pout∗, which corresponds to the variance introduced to the input parameters pin. The majority of the distributions are a little bit skewed but close to the normal distribution. The parameter n22 in the RC to ARX transformation (Figure 6, Table 3) is the worst case with a mean absolute error superior to 22%, which is still a satisfactory result. Indeed, by looking at the analytical expression, n22 depends on fives parameters with 10% of their values as variances. Furthermore, it is more difficult to obtain accurate results with very small values like n22=3.57·10-10.In the other way, from ARX to RC, the transformations are more sensitive to parameter deviations; this is why the variance of (1) has been scaled with δ=0.01 instead of δ=0.1 as previously. The histograms of Rco, Rci, and Rw in Figure 8 are zoomed to see most of the distribution. Otherwise, outliers far from the true parameter values give an improper graphical scale. The same situation occurs for Rco, Rci, Rw, and Cw in Figure 9.From ARX to SS, the parametersa11 and a12 are not present in the results (Figure 7 and Table 4) because they are not a function of the ARX parameters; they are known directly from the similarity matrix T. In this transformation, all the parameters are normally distributed but the mean absolute errors of a21, a22, and b21 are far greater than their respective true values. These transformations depend on the poles d1 and d2 which are several orders of magnitude higher than the output parameters; consequently uncertainties introduced in the input parameters have strong impacts on the output parameters.The transformations from SS to RC involve only operations with values relatively small and with different orders of magnitude which explains the skewed forms ofRco, Rci, and Rw (Figure 8) and the presence of outliers. Moreover, Rci and Rw have complicated transformations involving six input parameters.The problem gets even worse when both transformations are combined in the ARX to RC transformation. The expected values ofRco, Rci, Rw, and Cw are far from their respective true values and the normal distribution shapes of Rco and Cw are completely lost and therefore cannot be shown properly by histograms. Only two parameters out of six are robust against uncertainties and some R and C parameters have negative values which has no physical meaning. These results show that the transformation from ARX to RC is not suitable even if only 1% of the parameter values are introduced as variance.In this section, the most sensitive relations to parameter deviations have been located. Sometimes, functions relating input and output parameters are complicated and it is hard to identify which parameters are most responsible for the output variance.


## 3. Factor Prioritization

As in Section2, the transformations are represented by a function of the form y=f(X), which relates k uncertain input parameters, X, to the output parameter y. This section introduces a variance-based method which ranks the input parameters, X, according to their contribution to the output parameter variance [13]. This information allows fixing parameters with negligible influence and, on the other hand, focusing on the most sensitive parameters to reduce the output variance. Indeed, constraints can be added to these parameters in the identification process, to increase the robustness of the transformations.The variance of the outputV(y) is decomposable in 2k summands (28)Vy=∑iVi+∑i∑j>iVij+V12,…,K,where(29)Vi=VXiEX~iy∣XiVij=VXiXjEX~ijy∣Xi,Xj-VXiEX~iy∣Xi-VXjEX~jy∣Xj,where Xi is the ith input and X~i denotes the vector of all inputs but Xi, EX~iy∣Xi is the expectation of y taken over all possible values of X~i while keeping Xi fixed, and VXi· is the variance over all possible values of Xi.Dividing both sides of (28) by Vy yields(30)∑iSi+∑i∑j>iSij+S12,…,K=1,where S are the sensitivity indices which express the share of variance of y that is due to a given input (Si) or input combination (Sij).The number of sensitivity indices of a model withk inputs grows as 2k since all the input interactions are taken into account. A solution to overcome this cumbersome approach is to only compute the k first-order (Si) and k total sensitivity indices (STi). These two indices are sufficient to describe synthetically the sensitivity pattern of the model.(31)Si=VXiEX~iy∣XiVy(32)STi=EX~iVXiy∣X~iVy,where VXiEX~iy∣Xi is the expected reduction in variance that would be obtained if Xi could be fixed and EX~iVXiy∣X~i is the expected variance that would be left if all factors but Xi could be fixed. STi represents Si plus all the interactions linked to Xi.The indicesSi and STi can be computed efficiently by the following sampling procedure [14]:(a)
Generate two independent sample matricesA and B∈RN×k, with k being the number of inputs and N being the number of simulations. Each sample matrix is defined in the k-dimensional unit hypercube.(b)
Generatek further matrices ABi∈RN×k, for i=1,2,…,k, such that all columns are from A except the ith column which is from B. This resampling strategy is called radial design.(c)
ComputeN(k+2) model evaluation, one for each row of the matrices A, B, and ABi, which gives, respectively, fA, fB, and f(ABi)(d)
Compute the first (31) and total (32) sensitivity indices using(33)VXiEX~iy∣Xi=1N∑j=1NfBjfABij-fAjEX~iVXiy∣X~i=12N∑j=1NfAj-fABij2.The factor prioritization is now applied to the transformations presented in Section 2 to have a better explanation of the results in Section 2.3.


### 3.1. Factor Prioritization Applied to Model Transformation

All the sensitivity indices were estimated withN=214 simulations to converge to a sufficient accuracy of three decimals. If the model transformation y=fX is additive, the sums of the first order and total indices must be equal to one. If the first index is inferior to one, it means there are interactions between inputs and therefore total sensitivity indices are superior to one. The gap between the actual sum and one is an indicator of the input’s interactions. In the following tables, parameters on the left column (outputs) are function of the parameters in the header (inputs); each output has two lines; the first one represents the first-order indices and the bold one the total-order indices.The factor prioritization results for the six transformations are given in Table7 for the RC to SS transformation, in Table 8 for the SS to ARX transformation, in Table 9 for the RC to ARX transformation, in Table 11 for the ARX to SS transformation, in Table 10 for the SS to RC transformation, and in Table 12 for the ARX to RC transformation.Table 7
First-order and total sensitivity indices, RC to SS transformation.
R
c
o
R
c
i
R
v
R
w
C
w
C
a
Sum
a
11
0
0.017
0.224
0.005
0
0.509
0.756
0
0.038
0.444
0.019
0
0.751
1.251
a
12
0
0.137
0
0.042
0
0.558
0.737
0
0.365
0
0.180
0
0.780
1.324
a
21
0
0.137
0
0.042
0.558
0
0.737
0
0.365
0
0.180
0.780
0
1.324
a
22
0.046
0.042
0
0.160
0.678
0
0.820
0.110
0.111
0
0.160
0.839
0
1.220
b
11
0
0
0.337
0
0
0.337
0.674
0
0
0.667
0
0
0.666
1.334
b
13
0
0
0
0
0
0.997
0.997
0
0
0
0
0
1.001
1.001
b
21
0.137
0
0
0.042
0.558
0
0.737
0.365
0
0
0.180
0.780
0
1.324
b
22
0.045
0
0
0.033
0.845
0
0.922
0.091
0
0
0.067
0.921
0
1.080Table 8
First-order and total sensitivity indices, SS to ARX transformation.
a
11
a
12
a
21
a
22
b
11
b
13
b
21
b
22
Sum
n
11
0
0
0
0
1.000
0
0
0
1.000
0
0
0
0
1.000
0
0
0
1.000
n
12
0
0.218
0
0.217
0.231
0
0.219
0
0.885
0
0.276
0
0.275
0.290
0
0.276
0
1.117
n
22
0
0.441
0
0
0
0
0
0.440
0.881
0
0.560
0
0
0
0
0
0.560
1.120
n
31
0
0
0
0
0
0.999
0
0
0.999
0
0
0
0
0
1.000
0
0
1.000
n
32
0
0
0
0.426
0
0.459
0
0
0.885
0
0
0
0.543
0
0.573
0
0
1.116
d
1
0.501
0
0
0.501
0
0
0
0
1.002
0.500
0
0
0.500
0
0
0
0
1.000
d
2
0.228
0.216
0.215
0.229
0
0
0
0
0.888
0.286
0.272
0.272
0.286
0
0
0
0
1.115Table 9
First-order and total sensitivity indices, RC to ARX transformation.
R
c
o
R
c
i
R
v
R
w
C
w
C
a
Sum
n
11
0
0
0.352
0
0
0.340
0.693
0
0
0.667
0
0
0.684
1.351
n
12
0.015
0.017
0.091
0.017
0.152
0.144
0.436
0.082
0.092
0.392
0.107
0.527
0.495
1.696
n
22
0.009
0.056
0
0.064
0.182
0.164
0.475
0.035
0.394
0
0.248
0.526
0.545
1.748
n
31
0
0
0
0
0
1.000
1.000
0
0
0
0
0
1.001
1.001
n
32
0.018
0.022
0
0.027
0.281
0.274
0.622
0.075
0.103
0
0.101
0.613
0.586
1.478
d
1
0.015
0.044
0.171
0.031
0.175
0.373
0.809
0.028
0.082
0.319
0.072
0.208
0.555
1.264
d
2
0.015
0.017
0.091
0.017
0.152
0.144
0.436
0.082
0.092
0.393
0.107
0.527
0.495
1.696Table 10
First-order and total sensitivity indices, SS to RC transformation.
a
11
a
12
a
21
a
22
b
11
b
13
b
21
b
22
Sum
R
c
o
0
0
0.184
0.180
0
0
0
0.362
0.726
0
0
0.405
0.408
0
0
0
0.496
1.309
R
c
i
0.126
0
0.019
0.100
0.130
0.243
0
0.108
0.726
0.276
0
0.091
0.177
0.277
0.335
0
0.146
1.302
R
v
0
0
0
0
0.644
0.182
0
0
0.826
0
0
0
0
0.819
0.354
0
0
1.174
R
w
0.036
0
0.047
0.254
0.037
0.070
0
0.267
0.709
0.090
0
0.226
0.441
0.090
0.106
0
0.363
1.315
C
w
0.018
0
0.287
0
0.020
0.275
0
0
0.560
0.076
0
0.662
0
0.075
0.635
0
0
1.449
C
a
0
0
0
0
0
1.005
0
0
1.005
0
0
0
0
0
1.004
0
0
1.004Table 11
First-order and total sensitivity indices, ARX to SS transformation.
n
11
n
12
n
22
n
31
n
32
d
1
d
2
Sum
a
21
0
0
0
0
0
0.484
0.516
1.000
0
0
0
0
0
0.484
0.516
1.000
a
22
0
0
0
0
0
0.997
0
0.997
0
0
0
0
0
1.000
0
1.000
b
11
1.000
0
0
0
0
0
0
1.000
1.000
0
0
0
0
0
0
1.000
b
13
0
0
0
1.002
0
0
0
1.002
0
0
0
1.001
0
0
0
1.001
b
21
0.621
0.287
0
0
0
0.072
0
0.979
0.641
0.287
0
0
0
0.091
0
1.019
b
22
0
0
0.999
0
0
0
0
0.999
0
0
0.999
0
0
0
0
0.999Table 12
First-order and total sensitivity indices, ARX to RC transformation.
n
11
n
12
n
22
n
31
n
32
d
1
d
2
Sum
R
c
o
0
0
0.836
0
0
0.061
0.063
0.960
0
0
0.870
0
0
0.082
0.087
1.040
R
c
i
0.716
0
0
0.095
0
0
0
0.811
0.902
0
0
0.286
0
0
0
1.189
R
v
0.646
0
0
0.180
0
0
0
0.826
0.819
0
0
0.355
0
0
0
1.173
R
w
0.715
0
0.002
0.094
0
0
0
0.811
0.901
0
0.002
0.286
0
0
0
1.188
C
w
0.180
0
0
0.577
0
0.014
0.012
0.783
0.380
0
0
0.788
0
0.029
0.034
1.231
C
a
0
0
0
1.003
0
0
0
1.003
0
0
0
1.003
0
0
0
1.003It has been highlighted in Section2.3 that the transformation from ARX to RC is highly sensitive to parameter deviations but, thanks to the factor prioritization, the influence of each parameter on the output variance is known.It has been shown in Section2.3 that the transformation from ARX to SS is very sensitive: good results are obtained only for Rv and Ca because they depend, respectively, on two and three parameters (Table 12). The parameters Rco and Cw present the worst results; by looking at Table 12, they are the only input parameters which are function of the poles d1 and d2. Even if the impact of d1 and d2 is very small, the difference of orders of magnitude between the parameters is unfavourable. The parameters Rci and Rw are function of two and three input parameters, respectively (Table 12); however almost 90% of the output variance is due to the ARX parameter n11. If a priori information is available for these parameters, boundaries on the ARX parameters can be determined such that the physical meaning is not violated. Another example is the parameter n22 in the RC to ARX transformation with an absolute mean error equal to 22% of its value. Both thermal capacities in this transformation, Ca and Cw, account for approximately 30% of the output variance; constraints can be added in the identification of the state-space to limit their variations and therefore make the transformation into ARX model more robust.


## 3.1. Factor Prioritization Applied to Model Transformation

All the sensitivity indices were estimated withN=214 simulations to converge to a sufficient accuracy of three decimals. If the model transformation y=fX is additive, the sums of the first order and total indices must be equal to one. If the first index is inferior to one, it means there are interactions between inputs and therefore total sensitivity indices are superior to one. The gap between the actual sum and one is an indicator of the input’s interactions. In the following tables, parameters on the left column (outputs) are function of the parameters in the header (inputs); each output has two lines; the first one represents the first-order indices and the bold one the total-order indices.The factor prioritization results for the six transformations are given in Table7 for the RC to SS transformation, in Table 8 for the SS to ARX transformation, in Table 9 for the RC to ARX transformation, in Table 11 for the ARX to SS transformation, in Table 10 for the SS to RC transformation, and in Table 12 for the ARX to RC transformation.Table 7
First-order and total sensitivity indices, RC to SS transformation.
R
c
o
R
c
i
R
v
R
w
C
w
C
a
Sum
a
11
0
0.017
0.224
0.005
0
0.509
0.756
0
0.038
0.444
0.019
0
0.751
1.251
a
12
0
0.137
0
0.042
0
0.558
0.737
0
0.365
0
0.180
0
0.780
1.324
a
21
0
0.137
0
0.042
0.558
0
0.737
0
0.365
0
0.180
0.780
0
1.324
a
22
0.046
0.042
0
0.160
0.678
0
0.820
0.110
0.111
0
0.160
0.839
0
1.220
b
11
0
0
0.337
0
0
0.337
0.674
0
0
0.667
0
0
0.666
1.334
b
13
0
0
0
0
0
0.997
0.997
0
0
0
0
0
1.001
1.001
b
21
0.137
0
0
0.042
0.558
0
0.737
0.365
0
0
0.180
0.780
0
1.324
b
22
0.045
0
0
0.033
0.845
0
0.922
0.091
0
0
0.067
0.921
0
1.080Table 8
First-order and total sensitivity indices, SS to ARX transformation.
a
11
a
12
a
21
a
22
b
11
b
13
b
21
b
22
Sum
n
11
0
0
0
0
1.000
0
0
0
1.000
0
0
0
0
1.000
0
0
0
1.000
n
12
0
0.218
0
0.217
0.231
0
0.219
0
0.885
0
0.276
0
0.275
0.290
0
0.276
0
1.117
n
22
0
0.441
0
0
0
0
0
0.440
0.881
0
0.560
0
0
0
0
0
0.560
1.120
n
31
0
0
0
0
0
0.999
0
0
0.999
0
0
0
0
0
1.000
0
0
1.000
n
32
0
0
0
0.426
0
0.459
0
0
0.885
0
0
0
0.543
0
0.573
0
0
1.116
d
1
0.501
0
0
0.501
0
0
0
0
1.002
0.500
0
0
0.500
0
0
0
0
1.000
d
2
0.228
0.216
0.215
0.229
0
0
0
0
0.888
0.286
0.272
0.272
0.286
0
0
0
0
1.115Table 9
First-order and total sensitivity indices, RC to ARX transformation.
R
c
o
R
c
i
R
v
R
w
C
w
C
a
Sum
n
11
0
0
0.352
0
0
0.340
0.693
0
0
0.667
0
0
0.684
1.351
n
12
0.015
0.017
0.091
0.017
0.152
0.144
0.436
0.082
0.092
0.392
0.107
0.527
0.495
1.696
n
22
0.009
0.056
0
0.064
0.182
0.164
0.475
0.035
0.394
0
0.248
0.526
0.545
1.748
n
31
0
0
0
0
0
1.000
1.000
0
0
0
0
0
1.001
1.001
n
32
0.018
0.022
0
0.027
0.281
0.274
0.622
0.075
0.103
0
0.101
0.613
0.586
1.478
d
1
0.015
0.044
0.171
0.031
0.175
0.373
0.809
0.028
0.082
0.319
0.072
0.208
0.555
1.264
d
2
0.015
0.017
0.091
0.017
0.152
0.144
0.436
0.082
0.092
0.393
0.107
0.527
0.495
1.696Table 10
First-order and total sensitivity indices, SS to RC transformation.
a
11
a
12
a
21
a
22
b
11
b
13
b
21
b
22
Sum
R
c
o
0
0
0.184
0.180
0
0
0
0.362
0.726
0
0
0.405
0.408
0
0
0
0.496
1.309
R
c
i
0.126
0
0.019
0.100
0.130
0.243
0
0.108
0.726
0.276
0
0.091
0.177
0.277
0.335
0
0.146
1.302
R
v
0
0
0
0
0.644
0.182
0
0
0.826
0
0
0
0
0.819
0.354
0
0
1.174
R
w
0.036
0
0.047
0.254
0.037
0.070
0
0.267
0.709
0.090
0
0.226
0.441
0.090
0.106
0
0.363
1.315
C
w
0.018
0
0.287
0
0.020
0.275
0
0
0.560
0.076
0
0.662
0
0.075
0.635
0
0
1.449
C
a
0
0
0
0
0
1.005
0
0
1.005
0
0
0
0
0
1.004
0
0
1.004Table 11
First-order and total sensitivity indices, ARX to SS transformation.
n
11
n
12
n
22
n
31
n
32
d
1
d
2
Sum
a
21
0
0
0
0
0
0.484
0.516
1.000
0
0
0
0
0
0.484
0.516
1.000
a
22
0
0
0
0
0
0.997
0
0.997
0
0
0
0
0
1.000
0
1.000
b
11
1.000
0
0
0
0
0
0
1.000
1.000
0
0
0
0
0
0
1.000
b
13
0
0
0
1.002
0
0
0
1.002
0
0
0
1.001
0
0
0
1.001
b
21
0.621
0.287
0
0
0
0.072
0
0.979
0.641
0.287
0
0
0
0.091
0
1.019
b
22
0
0
0.999
0
0
0
0
0.999
0
0
0.999
0
0
0
0
0.999Table 12
First-order and total sensitivity indices, ARX to RC transformation.
n
11
n
12
n
22
n
31
n
32
d
1
d
2
Sum
R
c
o
0
0
0.836
0
0
0.061
0.063
0.960
0
0
0.870
0
0
0.082
0.087
1.040
R
c
i
0.716
0
0
0.095
0
0
0
0.811
0.902
0
0
0.286
0
0
0
1.189
R
v
0.646
0
0
0.180
0
0
0
0.826
0.819
0
0
0.355
0
0
0
1.173
R
w
0.715
0
0.002
0.094
0
0
0
0.811
0.901
0
0.002
0.286
0
0
0
1.188
C
w
0.180
0
0
0.577
0
0.014
0.012
0.783
0.380
0
0
0.788
0
0.029
0.034
1.231
C
a
0
0
0
1.003
0
0
0
1.003
0
0
0
1.003
0
0
0
1.003It has been highlighted in Section2.3 that the transformation from ARX to RC is highly sensitive to parameter deviations but, thanks to the factor prioritization, the influence of each parameter on the output variance is known.It has been shown in Section2.3 that the transformation from ARX to SS is very sensitive: good results are obtained only for Rv and Ca because they depend, respectively, on two and three parameters (Table 12). The parameters Rco and Cw present the worst results; by looking at Table 12, they are the only input parameters which are function of the poles d1 and d2. Even if the impact of d1 and d2 is very small, the difference of orders of magnitude between the parameters is unfavourable. The parameters Rci and Rw are function of two and three input parameters, respectively (Table 12); however almost 90% of the output variance is due to the ARX parameter n11. If a priori information is available for these parameters, boundaries on the ARX parameters can be determined such that the physical meaning is not violated. Another example is the parameter n22 in the RC to ARX transformation with an absolute mean error equal to 22% of its value. Both thermal capacities in this transformation, Ca and Cw, account for approximately 30% of the output variance; constraints can be added in the identification of the state-space to limit their variations and therefore make the transformation into ARX model more robust.


## 4. Conclusion

The possibility of switching from a thermal network to state-space representation and to an autoregressive model is appealing since both advantages can be cumulated. The transfer function is practical and not cumbersome whereas the state-space brings physical information which is related to the thermal network. Both forms of the model can be used for control strategy and diagnosis, but which one is more robust to error transmission? To respond to this question, this study is composed of two steps. First, Monte Carlo simulations are used to determine the transformations characteristics and, then, sensitivity analysis is used to rank parameter influences on the output variance. It appears that the transformation from ARX model to RC parameter is less robust to uncertainties. We conclude that identifying the parameters on the state-space representation is a better choice since in this way the results indicate that the transformations are more robust against parameter deviations.



---

*Source: 5636145-2017-09-06.xml*