13  CPE modeling results

In this section, results to the CPECentroid Position Evolution modeling explained in subsection 10.0.1, shall be presented and assessed . First, a selection of equations, which defines the CPECentroid Position Evolution are given for one model configuration. Next, representative plots of the CPECentroid Position Evolution for different models are analyzed. Finally, the predicted centroid position is compared with the actual clustered centroid position.

Modeling the CPE returns, among other results, analytical equations. These equations describe the behavior of the centroid positions across the range β and can also be used for making predictions for βunseen. The model configuration for which they are be presented is SLS, i.e. Lorenz system 7.1, K=10,βtr=[β0=28;βend=33],nβ,tr=7. The analytical CPECentroid Position Evolution expressions are listed in 13.1 to 13.3 for the centroids [1,2,7], respectively. Recalling that the behavior across the 3 different axes (x, y, z) can vary greatly, each axis has its own regression model (x~,y~,z~). Thus, for each label, 3 different analytical expressions are provided.

(13.1)x~=0.1661cos(3β)y~=0.1375cos(3β)z~=0.8326β

Figure 13.1— SLS, CPE model for centroid: 1

(13.2)x~=0.1543sin(3β)+0.2446cos(3β)y~=0.2638sin(3β)+0.4225cos(3β)z~=0.4877β

Figure 13.2— SLS, CPE model for centroid: 2

(13.3)x~=0.1866β+0.133sin(3β)+0.1411cos(3β)y~=0.3βz~=1.0483+0.6358β

Figure 13.3— SLS, CPE model for centroid: 7

Right to the equations the corresponding plots are depicted in figures 13.1 to 13.3 . Here, the blue and green curves indicate true and modeled CPE, respectively. Each of the figures supports the choice of allowing each axis to be modeled separately. The z-axis appears to undergo less alteration or to be more linear than the x- and y-axis. If a model is supposed to be valid for all 3 axes, a more complex model, i.e., a higher of terms, is required. Although more flexible models fit training data increasingly better, they tend to overfit. In other words, complex models capture the trained data well but could exhibit oscillations for βunseen. The latter is even more severe when the model is employed for extrapolation. The difference between interpolation and extrapolation is that for extrapolation the prediction is made with βunseen which are not in the range of the trained βtr. Therefore, less complexity is preferred.

Next, the performance of predicting the centroid for βunseen is elaborated. For this purpose, figures 13.4 to 13.7 shall be examined. All figures depict the original centroid positions, which are obtained through the clustering step in green and the predicted centroid positions in blue. On closer inspection, orange lines connecting the true and predicted centroid positions can be identified. Note, that they will only be visible if the deviation between the true and predicted state is high enough. Figures 13.4 (a) an 13.4 (b) show the outcome for SLS with βunseen=28.5 and βunseen=32.5, respectively. Visually, both predictions are very close to the true centroid positions. Because of this high performance in showed in figures 13.5 (a) and 13.5 (b) two examples for extrapolation are given for βunseen=26.5 and βunseen=37, respectively. For the first one, the outcome is very applicable. In contrast, βunseen=37 returns some deviations which are notably high.

(a) βunseen=28.5

(b) βunseen=32.5

Figure 13.4— SLS, original vs. modeled centroid position, βunseen=28.5 and βunseen=32.5

(a) βunseen=26.5

(b) βunseen=37

Figure 13.5— SLS, original vs. modeled centroid position, extrapolated βunseen=26.5 and βunseen=37

βunseen=31.75

βunseen=51.75

Figure 13.6— LS20, original vs. modeled centroid position, βunseen=31.75 and βunseen=51.75

βunseen=8.7

βunseen=10.1

Figure 13.7— FW15, original vs. modeled centroid position, βunseen=8.7 and βunseen=10.1

Quantitative measurements are performed by applying the Mean Square Error (MSE) following equation 13.4 . The variables are denoted as the number of samples n, which in this case is equal to the number of centroids n=K, the known f(xk) and the predicted yk centroid position. (13.4)MSE=1ni=1n(f(xk)yk)2

The measured MSE errors for all displayed results are summarized in table 13.1 . The MSE for results of βunseen=28.5 and βunseen=32.5 in figures 13.4 is 0.622 and 0.677, respectively. Consequently, the performance of CNMccontrol-oriented Cluster-based Network Modeling is also confirmed quantitatively. Figures in 13.6 illustrate the outcome for LS20 for βunseen=31.75 and βunseen=51.75. In section 12 it is explained that for LS20 cluster network deformations appear . Nevertheless, the outcome visually and quantitatively endorses the CPE modeling capabilities. Figures in 13.7 depict the outcome for FW15 for βunseen=8.7 and βunseen=10.1. A few orange lines are visible, however overall the outcome is very satisfactory.

Table 13.1— MSE for different Model configurations and βunseen
Figure Model βunseen MSE
13.4 SLS 28.5 0.622
13.4 SLS 32.5 0.677
13.5 SLS 26.5 1.193
13.5 SLS 37 5.452
13.6 LS20 31.75 1.857
13.6 LS20 51.75 2.536
13.7 FW15 8.7 1.617
13.7 FW15 10.1 1.5

It can be concluded that the CPE modeling performance is satisfying. In the case of a few cluster network deformations, CNMccontrol-oriented Cluster-based Network Modeling is capable of providing acceptable results. However, as shown with SLS, if the model’s training range βtr and the number of K was selected appropriately, the MSE can be minimized.