为什么两个方法给不同的结果,哪些是正确的?
1视图(30天)
显示旧的评论
我使用ode15i和数值求解隐式函数,但是两种方法给不同的结果,这是正确的?
清晰的
clc
关闭所有
% % ode15i
f = @ (t, y, yp) yp。^ 2 - 5 / (0.5 * yp - 2. / y);
(t0, y0, yp0 tspan] =交易(0,1.924,2.079,1 [0]);
[y0, yp0] = decic (f, t0, y0, 1 yp0 0);
索尔= ode15i (f, tspan, y0, yp0);
次要情节(1、2、1)
情节(sol.x sol.y)
标题(“ode15i”)
% %数值
tspan = [0 1];%的时间间隔
y0 = 1.924;%初始值
[t、y] =数值(@ (t, y) odefcn (t, y) tspan, y0);%数值
次要情节(1、2、2)
情节(t, y)
标题(“数值”)
%定义函数
函数Dy = odefcn (t, y)
有趣= @ Dy (Dy)。^ 2 - 5 / (0.5 * Dy-2 / y);%隐函数
Dy = fzero(有趣,0);% fzero函数解决对内隐函数
结束
并绘制结果如下:
0评论
接受的答案
大卫。马塞罗说道
2022年4月1日
编辑:大卫。马塞罗说道
2022年4月1日
因为你有一个隐含的颂歌,给出正确的解决方案
ode15i。
的问题
数值
在最初的想用来找到的价值吗
Dy
。我使用下面的东西5而不是0,得到相同的结果
ode15i
。
清晰的
clc
关闭所有
% % ode15i
f = @ (t, y, yp) yp。^ 2 - 5 / (0.5 * yp - 2. / y);
(t0, y0, yp0 tspan] =交易(0,1.924,2.079,1 [0]);
[y0, yp0] = decic (f, t0, y0, 1 yp0 0);
索尔= ode15i (f, tspan, y0, yp0);
次要情节(1、2、1)
情节(sol.x sol.y)
标题(“ode15i”)
% %数值
tspan = [0 1];%的时间间隔
y0 = 1.924;%初始值
[t、y] =数值(@ (t, y) odefcn (t, y) tspan, y0);%数值
次要情节(1、2、2)
情节(t, y)
标题(“数值”)
%定义函数
函数Dy = odefcn (t, y)
有趣= @ Dy (Dy)。^ 2 - 5 / (0.5 * Dy-2 / y);%隐函数
Dy = fzero(有趣,5);% fzero函数解决对内隐函数
结束