| 查看: 4158 | 回复: 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 |
回复此楼» 猜你喜欢
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab编程 ODE45
已经有4人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- matlab 计算后 显示NAN
已经有5人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
9楼2016-02-18 13:06:43
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
文刀星
木虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 3651.2
- 散金: 500
- 红花: 5
- 帖子: 60
- 在线: 100.9小时
- 虫号: 3693767
- 注册: 2015-02-22
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2015-05-16 17:03:21
