[已答复] 求助ode45的问题

[复制链接]
糖糖糖1235 发表于 2021-11-23 20:56:10
本帖最后由 糖糖糖1235 于 2021-11-23 20:57 编辑

%该函数解析解应该是2*exp(x)+c,保存为jm.m文件
function dy = odefun(t,y)
dy = zeros(1,1);    % a column vector
dy=2*y(1);
end


tspan = [0 1];
% 初始值;
[T,Y] = ode45(@jm,tspan,2); %该函数解析解应该是2*exp(x)+1
正确答案应该是6.4  为什么我的Y算出来14多,这误差有这么大吗?


1 条回复


TouAkira 发表于 2021-11-24 05:15:47
:L
因为你把解析解算错了
  1. syms y( t ) t
  2. Sol_y = dsolve( [ diff( y, t, 1 ) == 2 * y ], [ y( 0 ) == 2 ] )
  3. vpa( subs( Sol_y, t, 1 ) )
复制代码

遇到这种问题请多检查检查自己算的对不对
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表