Re: [問題] 非線性ODE求解已回收
clear
clc
Time = [0 500];
a = 0.4;
b = 15055*10^(-6);
c = 523;
m = 120;
e = 0.8;
x_ini = [0; 298];
[t x] = ode15s(@CaseII_ode,Time,x_ini,[],a,b,c,e,m);
plot(t,x)
function dx = CaseII_ode(t,x,a,b,c,e,m)
y = 4422 - (80-10*cos(t)-(70^2-(10*sin(t))^2)^0.5)*201;
z = -2010*sin(t) - 20100/(4900-100*sin(t)^2)^0.5*sin(t)*cos(t);
dx = [x(2); ...
-a*x(1)*z/y];
不曉得是不是Y兄要的結果
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 123.193.30.39
推
03/11 10:33, , 1F
03/11 10:33, 1F
→
03/11 17:40, , 2F
03/11 17:40, 2F
→
03/11 17:47, , 3F
03/11 17:47, 3F
討論串 (同標題文章)
完整討論串 (本文為第 4 之 4 篇):
問題
1
3