问下当设置tspan为0-1.7的时候，图中显示黑线和蓝线在1.4到1.6秒之间相交，但当设置tspan为1.4-1.6的时候，新的图里两条线又不相交了，这是怎么回事啊
matlab用ode45求解多元微分方程
问下当设置tspan为0-1.7的时候，图中显示黑线和蓝线在1.4到1.6秒之间相交，但当设置tspan为1.4-1.6的时候，新的图里两条线又不相交了，这是怎么回事啊
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
分享- 邀请回答
-
1条回答 默认 最新
joel_1993 2022-04-12 21:58关注只能说你犯了“刻舟求剑”的错误，tspan在
0~1.7时，初值选取是在t=0处的，而tspan在1.4到1.6时，你给的初值时t=1.4处的，而1.4处你根本不知道初值是多少，你只知道0处的啊，所以建议就是tspan还是在0~1.7，但是可以这样写tspan = [0:0.1:1.4, 1.4001:0.0001:1.6];
最后plot结束，加个语句axis([1.4, 1.6, -inf, inf])有帮助望采纳哟谢谢啦
本回答被题主选为最佳回答 , 对您是否有帮助呢?评论 打赏解决 1无用举报 编辑记录 分享
- 2022-04-12 21:14回答 1 已采纳 只能说你犯了“刻舟求剑”的错误，tspan在0~1.7时，初值选取是在t=0处的，而tspan在1.4到1.6时，你给的初值时t=1.4处的，而1.4处你根本不知道初值是多少，你只知道0处的啊，所以建
- 2022-03-08 16:12回答 2 已采纳 你的L是向量，这是不允许的，你可以把L定义成和时间相关的量，或者你的L就是一个数，这样的ODE才是可以求解的，比方说 L = @(t) t; K0 = @(t) pi/2*L(t); K1 = @(t
- 2021-10-13 20:32回答 1 已采纳 你好，两个方程一样求解：代码供参考，该注释的都注释了 % 设置物块质量和弹簧刚度 m1 = 1; m2 = 1; k1 = 1; k2 = 1; k3 = 1; % 设x(1)对应物块1的位移;x(
- 2021-04-24 15:58辟谣的大舌头LONG的博客 简单的微分方程边值问题(BVP)，1stOpt很容易求解：Variable t=[0,2*pi],y=[1,1],y'=[0,0];Plot t[x],y, y'[y2];ODEFunction y'''=y''-y'+y;复制代码目标函数: 1.04603504250813E-13边值估算:y''(t=0): -1....
- 2021-10-07 12:24回答 1 已采纳 你好，既然你知道了y和x那么数值求到了dy/dt和dx/dt，自然就可以通过插值方法得到与x和y等长的dy/dt和dx/dt，然后直接进行拟合就好了啊，lsqcurvefit也好，最小二乘法加上fmi
- 2021-09-11 21:57回答 1 已采纳 同学，你好！我已经找到出问题的原因了，并且帮你更正了。 clear tspan=[0,100]; y0=[9000,100,0,0,10000,100]; beta=0.375; gamma=0.07
- 2021-11-27 21:21回答 1 已采纳 你举个例子啊，不要老提概念，可能你的概念在别人眼里又是一种解释，徒增工作量
- 2022-07-22 19:17studyer_domi的博客 % 调用语句 [T1,Y1] = ode45( @odefun, tspan, y0); % 绘图 figure plot(T1,Y1(:,1)) xlabel 时间 ylabel S title '\Lambda=0.2' figure plot(T1,Y1(:,2)) xlabel 时间 ylabel E title '\Lambda=0.2' figure plot...
- 2021-08-05 14:20回答 1 已采纳 dfun = @(t,y) -y^3 + y + t; [T1,X1] = ode45(dfun,[0 10],1); [T2,X2] = ode23(dfun,[0 10],1); figure p
- 2021-06-18 16:13回答 1 已采纳 你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答 本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务
- 2021-05-19 20:43回答 2 已采纳 odefun = @(t,y) [y(2); -y(1)*cos(t)]; [T,Y] = ode23(odefun,[0 10],[1 0]); figure plot(T,Y)
- 2020-04-30 22:01枇杷鹭的博客 使用 matlab 自带的 ode45 ，方程组用句柄表示。
- 2021-08-16 16:08回答 1 已采纳 同学你好！图片中的第四行代码改成: plot(tSol, XSol(:,1), 'g', tSol, XSol(:,2), 'r'); 即可！
- 2022-06-03 17:57Allen953的博客 ode45就是一种常微分方程求解器，这种求解器采用的是Runge-Kutta解法的 中阶 解法；ode45即Nonstiff（非刚性问题）微分方程式。假设要解下面这个微分防方程：我们把这个Second order differential equation(二阶微分...
- 2021-04-18 17:30小蓝紫蓝莓的博客 matlab关于ode45解二阶微分方程的困惑matlab关于ode45解二阶微分方程的困惑一个二阶微分方程：y''+y'+y=sin(t)初始条件为y(0)=5,y'(0)=6。过程：先降阶为一阶微分方程组y'=zz'=-z-y+sin(t)编制如下函数m文件function...
- 2018-07-08 12:29Junbo_Just_Jia的博客 在工程和计算中常常需要求解微分方程，但微分方程的解析解通常很难得到，所以一种退而求其次的方法就是求其数值解。这里是使用Matlab中的ode45求解常微分方程数值解的简单说明。具体可使用"help ode45"...
- 2022-07-15 16:58在MATLAB中，`ode45`是求解常微分方程（ODE）初值问题的主要函数。这个工具能够高效、稳定地处理各种...这些示例可以帮助初学者和有经验的用户加深对MATLAB数值求解微分方程的理解，从而更好地应用于科研和工程实践中。
- 2021-04-18 03:18weixin_39550172的博客 原问题为利用ode45求解下面微分方程：初值为可以转为二阶微分方程利用Matlab求解，下面为源代码functionMyode45tspan = [2,4]; %求解区间y0 = [1 0]; %y0的初值形式为行向量[t... %使用ode45求解微分方程plot(t,x(:,...
- 2021-04-23 17:43福路娃的博客 (2)).*y(1)+0.0001*(1-y(2)).^2; dy2dx = -1e4*dy1dx + ... 高阶微分方程odefile的编写求解: ......Runge-Kutta法是常微分方程的一种经典解法 ? MATLAB 对应命令:ode45 四阶Runge-Kutta公式 yn?1 ? yn ? h 6 (k1 ? ...
- 2022-03-28 14:12Matlab 求解微分方程（ODE）
- 没有解决我的问题, 去提问
问题事件
系统已结题
4月20日
已采纳回答
4月12日-
创建了问题
4月12日
悬赏问题
- ¥15 Stata链式中介效应代码修改
- ¥15 latex投稿显示click download
- ¥15 请问读取环境变量文件失败是什么原因？
- ¥15 在若依框架下实现人脸识别
- ¥15 网络科学导论，网络控制
- ¥100 安卓tv程序连接SQLSERVER2008问题
- ¥15 利用Sentinel-2和Landsat8做一个水库的长时序NDVI的对比，为什么Snetinel-2计算的结果最小值特别小，而Lansat8就很平均
- ¥15 metadata提取的PDF元数据，如何转换为一个Excel
- ¥15 关于arduino编程toCharArray()函数的使用
- ¥100 vc++混合CEF采用CLR方式编译报错