Can anyone help to see what is wrong with this function to solve a DAE problem
조회 수: 1 (최근 30일)
이전 댓글 표시
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.
댓글 수: 2
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!