| 查看: 4200 | 回复: 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 |
回复此楼» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程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
