Description Usage Arguments Author(s) See Also Examples
Combines equations of different sources
into a single system. During this combination, the polynomial
maximal degree can be either imposed or optimized to reduce
the model size. All the input have to follow
the convention defined by poLabs
.
1 2 3 4 5 6 7 8 |
inK |
A list of models, each provided as a matrix. A single matrix can also be provided, it will be transformed into a list containing a single matrix. |
inXnote |
A list of vectors with the names of the input variables for each model. If not provided, default notation is used: "X1", "X2", etc. A single matrix can also be provided, it will be transformed into a list containing a single matrix. |
eqNum |
A list of vector, providing each the equations number (relating to the input models) to be kept in the output equation system. If not provided, all the equations are kept. A single matrix can also be provided, it will be transformed into a list containing a single matrix. |
XnoteOut |
A vector with the names of the output variables. If not provided, default notation is used considering that the variables of the input models are all different |
nVarOut |
The dimension of the output equation system (if not provided, this degree is deduced from the input models) |
dMaxOut |
The maximal polynomial degree of the output equation system (if not provided, this degree is deduced from the input models) |
Sylvain Mangiarotti
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 | # Load models
data("allMod_nVar3_dMax2")
# Display equations of system 1
visuEq(nVar = 3, dMax = 2, K = allMod_nVar3_dMax2$NH86, substit = 1)
# Display equations of system 2
visuEq(nVar = 3, dMax = 2, K = allMod_nVar3_dMax2$R76, substit = 1)
# put the two systems in a list
allK <- list()
allK[[1]] <- allMod_nVar3_dMax2$NH86
allK[[2]] <- allMod_nVar3_dMax2$R76
# Example 1: reformulate two autonomous system in a single matrix
visuEq(K = allK[[1]], substit = c('u', 'v', 'w'))
visuEq(K = allK[[2]], substit = c('X', 'Y', 'Z'))
Knew <- combiEq(allK)
visuEq(K = Knew, substit = c('u', 'v', 'w', 'X', 'Y', 'Z'))
# Example 2
inXnote = list()
inXnote[[1]] <- c('u', 'v', 'w')
inXnote[[2]] <- c('X', 'Y', 'Z')
visuEq(K = allK[[1]], substit = inXnote[[1]])
visuEq(K = allK[[2]], substit = inXnote[[2]])
XnoteOut = c('X', 'Y', 'Z', 'u', 'v', 'w')
Knew2 <- combiEq(allK, inXnote = inXnote, XnoteOut = XnoteOut)
visuEq(K = Knew2, substit = XnoteOut)
# Example 3
inXnote = list()
inXnote[[1]] <- c('u', 'v', 'w')
inXnote[[2]] <- c('X', 'Y', 'Z')
visuEq(K = allK[[1]], substit = inXnote[[1]])
visuEq(K = allK[[2]], substit = inXnote[[2]])
XnoteOut = c('u', 'X', 'v', 'Y', 'w', 'Z')
Knew3 <- combiEq(allK, inXnote = inXnote, XnoteOut = XnoteOut, dMaxOut = 3)
visuEq(K = Knew3, substit = XnoteOut)
# Example 4
dim(Knew3)
inXnote = c('x', 'X', 'y', 'Y', 'z', 'Z')
visuEq(K = Knew3, substit = inXnote)
XnoteOut = c('X', 'Y', 'Z')
Knew4 <- combiEq(Knew3, inXnote = inXnote, XnoteOut = XnoteOut)
dim(Knew4)
visuEq(K = Knew4, substit = XnoteOut)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.