特殊微分方程的数值解——微分代数方程

所谓的微分代数方程（DAE），是指在微分方程中，某些变量间满足某些代数方程的约束。

微分代数方程的一般形式为：

其中，M（t，x）矩阵的为奇异矩阵。在求解微分方程时，两边不能同时乘以M矩阵的逆矩阵M‘，所以这类方程并不能变成普通的微分方程取求解。

例题1：

假设我们这里有三个微分方程，第一个和第二个是常规形式的微分方程，第三个方程退化成了一个代数方程，所以这类方程又称为微分代数方程。

给定初值条件，x1(0)=0.8, x2(0)=x3(0)=0.1

用矩阵形式表示该微分代数方程：

f=@(t,x)[-0.2*x(1)+x(2)*x(3)+0.3*x(1)*x(2);
    2*x(1)*x(2)-5*x(2)*x(3)-2*x(2)*x(2);
    x(1)+x(2)+x(3)-1];
M=[1,0,0;0,1,0;0,0,0];
options=odeset;options.Mas