减少高次微分方程系统到一阶微分方程的等效系统
降低含有更高阶的DAE仅含有一阶DAE的一个系统的系统。
创建微分方程的系统,它包括一个二阶表达式。这里,X(t)的
和Y(t)的
是系统的状态变量,C1
和C2
为参数。指定方程和变量作为两个符号载体:方程式符号方程的矢量和变量为符号函数调用的向量。
SYMS X(t)Y(t)的C1 C2等式= [DIFF(X(t)的,T,T)+的sin(x(t))的+ Y(t)的== C1 * COS(T),...的diff(Y(T),T)== C2 * X(t)];瓦尔= [X(t),Y(t)的];
重写这个系统,使所有方程成为一阶微分方程。该reduceDifferentialOrder
功能通过一阶表达式通过引入新的变量取代高阶DAEDXT(t)的
。这也代表了所有的方程式符号表达式。
[newEqs,newVars] = reduceDifferentialOrder(方程,乏)
newEqs = DIFF(DXT(T),T)+的sin(x(t))的+ Y(T) - C1 * cos(t)的差异(Y(T),T) - C2 * X(t)的DXT(吨) - DIFF(X(t),T)newVars = X(t)Y(t)的DXT(t)的
减少含有第二和第三阶表达式仅含有一阶DAE的一个系统的系统。此外,返回表示所产生的变量的矩阵来reduceDifferentialOrder
通过该系统的原始变量。
创建差分方程的系统,其包括第二和第三阶表达式。这里,X(t)的
和Y(t)的
是系统的状态变量。指定方程和变量作为两个符号载体:方程式符号方程的矢量和变量为符号函数调用的向量。
SYMS X(t)Y(t)的F(T)等式= [DIFF(X(t)的,T,T)==差异(F(T),T,T,T),DIFF(Y(t)的吨,T,T)==差异(F(T),T,T)];瓦尔= [X(t),Y(t)的];
呼叫reduceDifferentialOrder
有三个输出参数。这句法收益矩阵[R
有两列:第一列包含新的变量,和第二列表示新的变量作为原始变量的衍生物,X(t)的
和Y(t)的
。
[newEqs,newVars,R] = reduceDifferentialOrder(方程,乏)
newEqs = DIFF(DXT(T),T) - DIFF(F(T),T,T,t)的差异(Dytt(T),T) - DIFF(F(T),T,T)DXT(t)的- DIFF(X(t),T)DYT(T) - DIFF(Y(T),T)Dytt(T) - DIFF(DYT(T),T)newVars = X(t)Y(t)的DXT(吨)DYT(t)的Dytt(t)的R = [DXT(t)的差异(X(t)中,t)] [DYT(t)的差异(Y(T),T)] [Dytt(t)的的diff(Y(t)的,T,T)]
daeFunction
|decic
|findDecoupledBlocks
|incidenceMatrix
|isLowIndexDAE
|massMatrixForm
|odeFunction
|reduceDAEIndex
|reduceDAEToODE
|reduceRedundancies