检查方程组的差动指数低于2
检查,如果一阶半线性的DAE的系统具有低的差动指数(0
要么1
)。
创建两个差分代数方程组的下列系统。这里,X(t)的
和Y(t)的
是系统的状态变量。指定方程和变量作为两个符号载体:方程式符号方程的矢量和变量为符号函数调用的向量。
SYMS X(t)Y(t)的等式= [DIFF(X(t)中,t)== X(t)+ Y(t)的,X(T)^ 2 + Y(t)的^ 2 == 1];瓦尔= [X(t),Y(t)的];
用isLowIndexDAE
检查系统的差分顺序。这个系统的差分顺序是1
。对于指数系统0
和1
,isLowIndexDAE
回报1
(真正
)。
isLowIndexDAE(均衡器,VARS)
ANS =逻辑1
检查以下DAE系统具有低或高差指数。如果该指数高于1
,然后使用reduceDAEIndex
减少它。
创建两个差分代数方程组的下列系统。这里,X(t)的
,Y(t)的
和Z(t)的
是系统的状态变量。指定方程和变量作为两个符号载体:方程式符号方程的矢量和变量为符号函数调用的向量。
SYMS X(t)Y(t)的Z(t)的F(T)等式= [DIFF(X(t)中,t)== X(t)+ Z(T),...的diff(Y(t)的,T)== F(T),X(T)== Y(吨)];瓦尔= [X(t),Y(t)的,Z(T)];
用isLowIndexDAE
检查系统的差分索引。对于这个系统isLowIndexDAE
回报0
(假
)。这意味着,系统的差分索引是2
或更高。
isLowIndexDAE(均衡器,VARS)
ANS = 0逻辑
用reduceDAEIndex
重写系统,使得所述差分索引是1
。调用带有四个输出参数此功能也显示了原系统的差异指数。新系统有一个额外的状态变量,DYT(t)的
。
[newEqs,newVars,〜,oldIndex] = reduceDAEIndex(方程,乏)
newEqs = DIFF(X(t),T) - Z(吨) - X(t)的DYT(T) - F(T)X(t) - Y(t)的差异(X(t),T) - DYT(吨)newVars = X(t)Y(t)的Z(t)的DYT(t)的oldIndex = 2
检查新系统的微分阶低于2
。
isLowIndexDAE(newEqs,newVars)
ANS =逻辑1
daeFunction
|decic
|findDecoupledBlocks
|incidenceMatrix
|massMatrixForm
|odeFunction
|reduceDAEIndex
|reduceDAEToODE
|reduceDifferentialOrder
|reduceRedundancies