matlab ode45函数设置终点停止条件不起作用的问题?

来源：蜘蛛抓取(WebSpider) 时间：2018-04-08 03:09 标签： matlab ode45
24小时热门版块排行榜&&&&
【悬赏金币】回答本帖问题，作者swuli将赠送您 5 个金币
(小有名气)
在线: 143.3小时
虫号: 960195
注册: 专业: 环境污染化学
matlab的ode45的初值问题
ode45 要求输入初值，但是不知道初值是多少，应该怎样做呢？
& 猜你喜欢
已经有5人回复
已经有16人回复
已经有5人回复
已经有28人回复
已经有63人回复
已经有56人回复
已经有5人回复
已经有17人回复
已经有3人回复
已经有3人回复
& 本主题相关价值贴推荐，对您同样有帮助:
已经有4人回复
已经有8人回复
已经有15人回复
已经有29人回复
已经有5人回复
已经有18人回复
已经有16人回复
已经有11人回复
已经有11人回复
已经有27人回复
已经有8人回复
已经有20人回复
已经有4人回复
已经有11人回复
已经有15人回复
已经有4人回复
已经有3人回复
(著名写手)
散金: 2144沙发: 4
在线: 1236.1小时
虫号: 653607
注册: 性别: GG专业: 机械动力学管辖:
★ 臭水沟(金币+1): 谢谢应助~~
引用回帖:Originally posted by swuli at
ode45 要求输入初值，但是不知道初值是多少，应该怎样做呢？ 什么模型啊？
用不确定的初值，求解出的结果也是没意义的。
明德厚学、求是创新
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请发送邮件至：
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright &
MuChong.com, All Rights Reserved. 小木虫 版权所有matlab中关于ode45的问题
[问题点数：40分，结帖人fengluchaying]
matlab中关于ode45的问题
[问题点数：40分，结帖人fengluchaying]
不显示删除回复
显示所有回复
显示星级回复
显示得分回复
只显示楼主
匿名用户不能发表回复！|查看: 2447|回复: 3|关注: 0
ode45算法中syms 变量与double变量相互转换问题
<h1 style="color:# 麦片财富积分
新手, 积分 16, 距离下一级还需 34 积分
用ode算法解微分方程过程中，若已知微分方程的符号函数形式，如何将其转化为能够参与计算的double类型？具体例子如下：
...(省略)...
[t,y]=ode45('zhong',[0,10],y0);
...(省略)...
_________________________
子m文件1：zhong.m
...此处应该添加什么语句？...
function ydot=[y(1);
& && && && && && && && &V1;
& && && && && && && && &y(3);
& && && && && && && && &V2;]
_________________________
子m文件2：guo.m& && &%（用符号函数推导出V1，V2的表达式）
function [V1,V2]=f(a,b,c,d)
syms a b c d
V1=a+b+c+d;& && &%（随便假设的一种关系式)
V2=a-b-c-d;& && && &%（随便假设的一种关系式）
V1=subs(V1,{'a','b','c','d'},{'y(1)','y(2)','y(3)','y(4)'});
V2=subs(V2,{'a','b','c','d'},{'y(1)','y(2)','y(3)','y(4)'});
小弟在论坛中也看到过这类问题的帖子，说是用eval 函数可以实现此目的。但无奈天资愚钝，对帖子中的eval 函数用法没能看懂。在matlab help中关于该函数的解释也没看太明白。想请教各位高手大侠：
（1）对于小弟这个问题，应该如何解决？
（2）eval 函数的作用
不胜感激！
<h1 style="color:# 麦片财富积分
额，没人看到么？
论坛优秀回答者
<h1 style="color:#99 麦片财富积分
关注者： 260
额，没人看到么？
我就知道eval是将字符串变为函数
你尝试一下
eval('x = 2*3+4');
论坛优秀回答者
关注者： 296
eval的用处是执行表达式，例如
eval('a = 1')
就相当于你写了语句：
你这里用eval不大好（或者说大部分时候用eval都不好），如果你的版本不算太老，可以借助matlabFunction将V1和V2分别转化为关于y的函数，之后在计算ydot时调用他们，不过需要注意的是符号引擎对数组的引用用的是y[1]这种，而非y(1)
另外也可以不用到符号计算推导表达式，而当成隐式微分方程用ode15i求解
站长推荐 /3
车辆动力学工具箱介绍：一款集成3D虚拟环境的汽车仿真产品
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区Matlab中龙格-库塔(Runge-Kutta)方法原理及实现_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
Matlab中龙格-库塔(Runge-Kutta)方法原理及实现
&&卫星定位的算法及matlable应用
阅读已结束，下载本文需要
想免费下载更多文档？
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
还剩2页未读，继续阅读
定制HR最喜欢的简历
你可能喜欢

我要回帖

更多关于 matlab ode45 的文章

 

随机推荐