solving a second order linear differential equation

7 ビュー (過去 30 日間)
Eliraz Nahum
Eliraz Nahum 2018 年 9 月 23 日
回答済み: Eliraz Nahum 2018 年 9 月 23 日
hello everybody, I was trying to solve a simple pendulum second order linear differential equation of the form y''=-(g/l)*sin(y) while using the ode45 function. since it's a second order equation I understood that I have to manipulate the problem, so it will fit the ode45.
The mathematical manipulation I did is described in the attached picture.
I can't understand why my code fails...please help
g=9.8;
l=0.5;
t_span=[0 30];
teta0=deg2rad(60);
F=@(t,q) [q(2);-(g/l)*sin(q(2))];
[t,q]=ode45(@(t,q) F(t,q),t_span,teta0)

採用された回答

Mischa Kim
Mischa Kim 2018 年 9 月 23 日
編集済み: Mischa Kim 2018 年 9 月 23 日
Almost there:
g = 9.8;
l = 0.5;
t_span = [0 30];
teta0 = deg2rad(60);
tetad0 = 0;
F = @(t,q) [q(2);-(g/l)*sin(q(1))]; % check your derivation
[t,q_t] = ode45(@(t,q) F(t,q),t_span,[teta0; tetad0])
This is a second order DE so you need two initial conditions, one for teta and one for tetad.

その他の回答 (1 件)

Eliraz Nahum
Eliraz Nahum 2018 年 9 月 23 日
thank you very much... sometimes it's not about the code, but the mathematics :-)

カテゴリ

Help Center および File ExchangeOrdinary Differential Equations についてさらに検索

製品


リリース

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by