问下当设置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(
- 2024-07-24 22:48计算小屋的博客 ode45 是 MATLAB 中用于求解非刚性常微分方程（ODE）的数值方法。它基于 Runge-Kutta 方法，并具有自适应步长调整机制，能够在一定误差控制范围内高效地计算 ODE 的数值解。
- 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 已采纳 你举个例子啊，不要老提概念，可能你的概念在别人眼里又是一种解释，徒增工作量
- 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....
- 2022-08-22 18:08回答 3 已采纳 （1）matlab基于四阶龙格库塔算法求解常微分方程组的函数ode45介绍ode45是matlab中被广泛使用的基于四阶龙格库塔算法求解非线性常微分方程组的函数同时，matlab中ode45的代码是开
- 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-16 12:07回答 2 已采纳 这个错误是使用了未定义的变量.
- 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-11-05 19:44回答 1 已采纳 此函数需要ode23( )调用。例如： [t,x]=ode23(@)odefunction,[0,10], [ 0,1,1 ] )
- 2020-04-30 22:01枇杷鹭的博客 使用 matlab 自带的 ode45 ，方程组用句柄表示。
- 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 ? ...
- 没有解决我的问题, 去提问
问题事件
系统已结题
4月20日
已采纳回答
4月12日-
创建了问题
4月12日
悬赏问题
- ¥15 系统 24h2 专业工作站版，浏览文件夹的图库，视频，图片之类的怎样删除？
- ¥15 怎么把512还原为520格式
- ¥15 MATLAB的动态模态分解出现错误，以CFX非定常模拟结果为快照
- ¥15 求高通平台Softsim调试经验
- ¥15 canal如何实现将mysql多张表(月表)采集入库到目标表中(一张表)？
- ¥15 wpf ScrollViewer实现冻结左侧宽度w范围内的视图
- ¥15 栅极驱动低侧烧毁MOSFET
- ¥30 写segy数据时出错3
- ¥100 linux下qt运行QCefView demo报错
- ¥50 F1C100S下的红外解码IR_RX驱动问题