| 查看: 3898 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
文刀星木虫 (小有名气)
|
[求助]
ode45积分出现NaN 已有1人参与
|
|
|
对一个函数积分，给定初始值，积分到某一时刻，出现了NaN的情况，最后一项是一个常数，不用管；然后以出现NaN之前的一个数据重新积分，却可以积分到最后，什么情况？ options=odeset('events',@eventRe); x0=[0.9922,-0.1396,0,1.0157,0,0,0.1667]; [t,x,te,xe,ie]=ode45(@myGAzrfun,[0,2],x0,options) %上面的积分积分到t=0.0324（t（17））的下一个节点0.0357时，出现NaN情况，见下面部分数据； x1=x(17,  ;[t,x,te,xe,ie]=ode45(@myGAzrfun,[t(17),2],x1,options) %然后将t=0.0324的时刻所对应的的初值及时间重新积分，却又能够继续积分下去； 0.9939 -0.1396 0 1.0139 0.0125 0 0.1667 0.9948 -0.1396 0 1.0130 0.0188 0 0.1667 0.9957 -0.1396 0 1.0121 0.0252 0 0.1667 0.9965 -0.1395 0 1.0112 0.0315 0 0.1667 NaN NaN NaN NaN NaN NaN 0.1667 NaN NaN NaN NaN NaN NaN 0.1667 NaN NaN NaN NaN NaN NaN 0.1667 NaN NaN NaN NaN NaN NaN 0.1667 NaN NaN NaN NaN NaN NaN 0.1667 %上面是积分遇到问题时的部分数据； %下面是积分到最后的部分数据，而且也是我想要的数据，说明应当是正确的，可是为啥会出现上面的情况？ 0.4418 -0.1594 0 1.0001 0.1092 0 0.1667 0.4221 -0.1612 0 1.0001 0.1096 0 0.1667 0.4030 -0.1631 0 1.0000 0.1100 0 0.1667 te = 0.1213 xe = 0.4030 -0.1631 0 1.0000 0.1100 0 0.1667 ie = 1  360截图20150516124651022.jpg |
回复此楼» 猜你喜欢
压国家一区线，求导师收留，有恩必谢！
已经有4人回复
-
081200-314
已经有4人回复
-
289求调剂
已经有9人回复
-
085405 考的11408求各位老师带走
已经有3人回复
-
339求调剂
已经有6人回复
-
材料与化工272求调剂
已经有4人回复
-
085701环境工程，267求调剂
已经有15人回复
-
材料求调剂一志愿哈工大324
已经有9人回复
-
292求调剂
已经有15人回复
-
304求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab编程 ODE45
已经有4人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- matlab 计算后 显示NAN
已经有5人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2015-05-16 17:03:21
2楼2015-05-16 14:37:27
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
3楼2015-05-16 14:45:41
5楼2015-05-16 18:56:59
20