Solve first ode using ode45
1 回表示 (過去 30 日間)
古いコメントを表示
Hello, my question is i have first ode dT/dt = -0.0004*[T(t) - 5], may i know how to use ode45 to solve t in the range T [8 88]? Can i know how about the script? Thank you for your help.
0 件のコメント
回答 (2 件)
Walter Roberson
2021 年 6 月 9 日
format long g
syms T(t) T8
dTdt = diff(T(t),t)
eqn = dTdt == -0.0004*(T - 5)
init = T(8) == T8
sol = simplify(dsolve([eqn, init]))
Specific_T8 = 6
Tsol = subs(sol, T8, Specific_T8)
fplot(Tsol, [8 88])
Jan
2021 年 6 月 9 日
編集済み: Jan
2021 年 6 月 9 日
You can modify the examples of the documentation easily, see: doc ode45
Temperature0 = 6.2; % Your initial value
[Time, Temperatur] = ode45(@fcn, [8, 88], Temperature0)
plot(Time, Temperature);
function dT = fcn(Time, Temperature)
dT = -0.0004 * (Temperatur - 5);
end
参考
カテゴリ
Help Center および File Exchange で Ordinary Differential Equations についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!