Найдите матрицу падения системы уравнений
A = incidenceMatrix(eqs,vars)
Найдите матрицу падения системы пяти уравнений в пяти переменных.
Создайте следующий символьный векторный eqs
, содержащий пять символьных дифференциальных уравнений.
syms y1(t) y2(t) y3(t) y4(t) y5(t) c1 c3 eqs = [diff(y1(t),t) == y2(t),... diff(y2(t),t) == c1*y1(t) + c3*y3(t),... diff(y3(t),t) == y2(t) + y4(t),... diff(y4(t),t) == y3(t) + y5(t),... diff(y5(t),t) == y4(t)];
Создайте вектор переменных. Здесь, c1
и c3
являются символьными параметрами (не переменные) системы.
vars = [y1(t), y2(t), y3(t), y4(t), y5(t)];
Найдите матрицу падения A
для уравнений eqs
и относительно переменных vars
.
A = incidenceMatrix(eqs, vars)
A = 1 1 0 0 0 1 1 1 0 0 0 1 1 1 0 0 0 1 1 1 0 0 0 1 1
daeFunction
| decic
| findDecoupledBlocks
| isLowIndexDAE
| massMatrixForm
| odeFunction
| reduceDAEIndex
| reduceDAEToODE
| reduceDifferentialOrder
| reduceRedundancies
| spy