当前位置:   article > 正文

matlab ode45_matlab ode45初值

matlab ode45初值

fun.m  文件里面定义微分方程组

function dx = fun(t, x)

dx = zeros(2, 1);

dx(1) = -x(1).^3 - x(2);

dx(2) = x(1) - x(2).^3;

end

 

运行命令:

[t,x] = ode45(@fun, [0,30], [1;0.5]);

[0, 30] 是 t 的区间,[1; 0.5] 是初值。

得到:一列 t, 两列 x。

 

精度参数:

options = odeset('RelTol', 1e-4, 'AbsTol', [1e-4, 1e-5]);

x(1), x(2) 的总的相对误差(relative error)是 1e-4,绝对误差(absolute error)分别是1e-4, 1e-5。

[t,x] = ode45(@fun, [0,30], [1;0.5], options);

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/一键难忘520/article/detail/814150
推荐阅读
相关标签
  

闽ICP备14008679号