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


رتبه موضوع:
  • 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
پاسخ
 سپاس شده توسط شماره مجازی امارات ، تلگرام ضد فیلتر 2023


پرش به انجمن:


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