`massMatrixForm`

Escreve um sistema DAE na forma de matriz de massa M*y' = F(t, y).

Categoria: Mathematics > Equation Solving

Sintaxe

[M, F] = massMatrixForm(eqs, vars)

Retorna a matriz de massa M e o lado direito F.

Descrição

Vários solvers (ode15s, ode23t) aceitam problemas na forma M(t,y) y' = f(t,y). massMatrixForm extrai essa estrutura das equações. Só funciona se o sistema for linear nas derivadas, e o resultado pode alimentar odeFunction e a opção Mass de odeset.

Argumentos de entrada

eqs: Vetor de equações simbólicas.
vars: Vetor das funções simbólicas incógnitas, na ordem das colunas de M.

Argumentos de saída

M: Matriz de massa simbólica. Linhas de zeros indicam equações algébricas.
F: Vetor coluna de expressões em t e nas variáveis.

Exemplo

Uma equação diferencial e uma algébrica

syms x(t) y(t)
eqs = [diff(x,t) == y, x + y == 1];
[M, F] = massMatrixForm(eqs, [x y])

Saída

M = 
[1, 0]
[0, 0]
F = 
           y(t)
1 - x(t) - y(t)

Veja também

daeFunction odeFunction isLowIndexDAE odeset