| 查看: 1482 | 回复: 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 |
回复此楼» 猜你喜欢
求调剂院校信息
已经有4人回复
-
085600材料与化工306
已经有4人回复
-
286求调剂
已经有10人回复
-
328求调剂，英语六级551，有科研经历
已经有9人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有4人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有8人回复
-
资源与环境 调剂申请(333分)
已经有5人回复
-
280求调剂
已经有12人回复
-
269专硕求调剂
已经有5人回复
-
材料学硕301分求调剂
已经有7人回复
2楼2019-05-04 22:49:42
