| 查看: 3836 | 回复: 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 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程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
- 专业: 导航、制导与传感技术
3楼2015-05-16 14:45:41
2楼2015-05-16 14:37:27
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2015-05-16 17:03:21
5楼2015-05-16 18:56:59