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.

답변 (2개)

Walter Roberson
Walter Roberson 2021년 6월 9일

1 개 추천

format long g
syms T(t) T8
dTdt = diff(T(t),t)
dTdt = 
eqn = dTdt == -0.0004*(T - 5)
eqn(t) = 
init = T(8) == T8
init = 
sol = simplify(dsolve([eqn, init]))
sol = 
Specific_T8 = 6
Specific_T8 =
6
Tsol = subs(sol, T8, Specific_T8)
Tsol = 
fplot(Tsol, [8 88])
Jan
Jan 2021년 6월 9일
편집: Jan 2021년 6월 9일

0 개 추천

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

카테고리

태그

질문:

2021년 6월 9일

댓글:

2021년 6월 9일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by