انجمن سایت کلیدستان



حل دستگاه معادلهزمان کنونی: ۱۳۹۷/۰۵/۰۱، ۰۴:۵۰ ق.ظ
کاربران در حال بازدید این موضوع: 2 مهمان
نویسنده: m.mjb
آخرین ارسال: m.mjb
پاسخ 1
بازدید 475

رتبه موضوع:
  • 0 رای - 0 میانگین
  • 1
  • 2
  • 3
  • 4
  • 5
حل دستگاه معادله

حل دستگاه معادله

#1
سلام . من باید یه دستگاه معادلات دیفرانسیل رو حل کنم (از روش 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
پاسخ


پرش به انجمن:


کاربران در حال بازدید این موضوع: 2 مهمان