Info

此问题已关闭。 请重新打开它进行编辑或回答。

differential Equations with ode45

1 次查看（过去 30 天）
Alex Pak
Alex Pak 2014-11-2
关闭： MATLAB Answer Bot 2021-8-20
Hello, please help me, I'm trying to solve differential Equations, but I can't understand how I can implement m1 ( m1 is derivative m)
Code
% function OdeCaller2
N=3;
J=12800;
q=800;
S=50;
l=5;
ks=-0.005;
ko=-0.01;
ka=-0.01;
ky=-0.01;
a=3;
kz=-0.08671;
[T,Y]=ode45(@odefun10, [0 4], [0 0 0 0]);
plot(T,Y)
grid on
function dy=odefun10(t,y)
omegaz=y(1);
omegay=y(2);
gamma=y(3);
V=y(4);
m=2*(N-V)+tand(V)*(omegay*cosd(gamma)-omegaz*sind(gamma));
m1=????????????; *How can I implement this?*
dy=zeros(4,1);
delta=(J/(q*S*l*ks))*(3*(m-omegaz)+m1)-(ko+ka*a+kz*omega)/ks;
dy(1)=((q*S*l)/J)*(ko+ko*a+ks*delta+kz*omegaz);
dy(2)=((q*S*l)/J)*(ko+ky*a+ky*delta+ky*omegay);
dy(3)=omegaz*cosd(gamma)+omegay*sind(gamma);
dy(4)=emegay-tand(V)*(omegay*cosd(gamma)-omegaz*sind(gamma));
end
end
  1 个评论
Zoltán Csáti
Zoltán Csáti 2014-11-2
Since the second equation contains delta (and within it, m1) which is in third and fourth equations, the system is implicit. Try to put it into explicit form, or use the ode15i solver.

此问题已关闭。

回答（0 个）

此问题已关闭。

标签

Community Treasure Hunt

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

Start Hunting!

Translated by