۱۳۹۶/۱۰/۱۶, ۰۲:۱۸ ب.ظ
سلام . من باید یه دستگاه معادلات دیفرانسیل رو حل کنم (از روش ode) بعد دو تا از معادلات یک ترمی دارند وابسته به زمان است.. ولی ارور میده! متوجه نمیشم باید چیکار کنم!؟
syms t;
a=7600;
mio=398600;
n=sqrt(mio/(a^3));
T=2*pi/n;
J1=165;
J2=150;
J3=55;
k2=(165-55)/150;
teta1 =exp(t*((90801841*2279642322115997427392329^(1/2))));
teta2=(0.000039/(n*n*(165-55)*3))*(1-cos(sqrt(3*n*n*k2)*t));
Mg=-3*n*n*[(J2-J3)*teta1/J1;(J1-J3)*teta2/J2;0];
f0=[0 0 0];
options = odeset('RelTol',1e-8);
[t,f]=ode45(@fn,[0 T],f0,options);
function dydt=fn(~,f,t,Mg)
w1=f(1);
w2=f(2);
w3=f(3);
Dw1=(150-55)*w3*w2/165+Mg(1);
Dw2=(55-165)*w3*w1/150+Mg(2);
Dw3=(165-150)*w1*w2/55;
dydt=[Dw1 Dw2 Dw3]';
end
syms t;
a=7600;
mio=398600;
n=sqrt(mio/(a^3));
T=2*pi/n;
J1=165;
J2=150;
J3=55;
k2=(165-55)/150;
teta1 =exp(t*((90801841*2279642322115997427392329^(1/2))));
teta2=(0.000039/(n*n*(165-55)*3))*(1-cos(sqrt(3*n*n*k2)*t));
Mg=-3*n*n*[(J2-J3)*teta1/J1;(J1-J3)*teta2/J2;0];
f0=[0 0 0];
options = odeset('RelTol',1e-8);
[t,f]=ode45(@fn,[0 T],f0,options);
function dydt=fn(~,f,t,Mg)
w1=f(1);
w2=f(2);
w3=f(3);
Dw1=(150-55)*w3*w2/165+Mg(1);
Dw2=(55-165)*w3*w1/150+Mg(2);
Dw3=(165-150)*w1*w2/55;
dydt=[Dw1 Dw2 Dw3]';
end