Can anyone help to see what is wrong with this function to solve a DAE problem

1 次查看（过去 30 天）
haohaoxuexi1
haohaoxuexi1 2022-1-19
评论： haohaoxuexi1 2022-7-10
dx(1)=x(2);
dx(2)=k1*(x(3)-x(1))-(knl_1*x(1)+knl_3*x(1)^3);
dx(3)=x(4);
dx(4)=(-c*x(4)-k1*(x(3)-x(1))+kAmp*cos(w*t)+theta_p*x(5))/m;
dx(5)=(-theta_p*(ks/(kp+ks))*(x(4)-x(2))*R_s-x(5))/R_s/C_p;
The mass matrix
M=diag([1 0 1 1 1]);
options = odeset('mass', M, 'RelTol',1.e-6);
I am trying to solve the equations below. The y values and V values are unknown and the rest of them are constant.
But, it gave out the error
"This DAE appears to be of index greater than 1."
Can anyone help to solve the problem?
Thank you very much.

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by