| 查看: 1489 | 回复: 1 | ||
allengjn新虫 (正式写手)
|
[求助]
为什么我使用ode45计算一个二阶微分方程，计算结果与时间终点会有关系？
|
|
我使用ode45计算一个二阶微分方程 计算程序d1=0;d2=1; [t,y]=ode45('equation',0:0.000001:0.000152,[d1,d2]); 调用函数function dy=equation(t,y) E=2.07e11; v=0.3; r1=0.01;r2=0.0095;m=1;c=0.9;yy=-1; t dy=[y(2);%这里是y的一阶导 -(4*E/(3*2*(1-v^2))*(r2*r1/(r1-r2))^0.5*y(1)^1.5)/m;]; end 计算中ode45的计算时间如果小于0.000151s，计算出的结果都是实数，如果大于0.000151在第128步的时候计算结果就会变成复数，这是为什么？  1.png |
回复此楼» 猜你喜欢
0854求调剂
已经有6人回复
-
求调剂推荐
已经有6人回复
-
材料299专硕求调剂
已经有17人回复
-
一志愿沪9，生物学326求调剂
已经有5人回复
-
294求调剂
已经有4人回复
-
售SCI一区文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
297，工科调剂?河南农业大学本科
已经有11人回复
-
22专硕求调剂
已经有11人回复
-
284求调剂
已经有20人回复
-
291求调剂
已经有7人回复
2楼2019-05-04 22:49:42
