问ODE15s矩阵是奇异的，接近奇异或严重缩放。结果可能不准确。RCOND = NaN

EN

我试图解一个DAE系统( ODE和代数方程的混合物)。所以，我的系统是这样的

/dt= F(t，y)

矩阵M基本上是对角线上的常数质量矩阵1，所有的零点都在最后几行，最后一列很少。

举个例子，质量矩阵如下所示

M= [1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0]

如果我们以矩阵为例，我的向量Fx是维数(4x1)。

如果我得到了这个错误，我该如何进行？

矩阵是奇异的，接近奇异或严重缩放。结果可能不准确。RCOND = NaN

我的系统解析代码：

Mass = Mass_Matrix( Nc,NT );

opt = odeset('RelTol', 10.0^(-3), 'AbsTol' , 10.0^(-3), 'Mass', Mass ,'MassSingular', 'yes', 'OutputFcn',@odeprint);

[T,Y] = ode15s(@(t,y)Function_Fx( t, y, Resolution, y01),[0 1],y0,opt);

我的向量Fx是由Function_Fx创建的

function Fx = Function_Fx( t, y, Resolution, y01)

分辨率是一个结构变量，y01是存储的值。

我的问题是问题在哪里，是F(t，y)还是M？给定M是常数质量矩阵。

谢谢你的帮助