| 查看: 3901 | 回复: 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 |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
材料277求调剂
已经有7人回复
-
数一英一271专硕（085401）求调剂，可跨
已经有5人回复
-
291求调剂
已经有21人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有8人回复
-
药学105500求调剂
已经有3人回复
-
考研调剂
已经有4人回复
-
0703化学
已经有7人回复
-
学硕274求调剂
已经有6人回复
-
321求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab编程 ODE45
已经有4人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- matlab 计算后 显示NAN
已经有5人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
5楼2015-05-16 18:56:59
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
