How to solve the following System of first order differential equations using ode45?please help

1 回表示 (過去 30 日間)
q1=[q11; q21; q13];
q10=[7.1892; 11.4189; 0.5826];
t1 = pi:0.1:2*pi;
dq11dt1= -k1*(12.7734-q11)*cos(q31);
dq21dt1=-k1*(12.7734-q11)*sin(q31);
dq31dt1= -k2*(0.3097+(13.6167-q21)*sin(t1));

回答 (1 件)

Aquatris
Aquatris 2018 年 9 月 7 日
編集済み: Aquatris 2018 年 9 月 7 日
First you create your function that outputs derivative of q when time and q are given to it.
function qd = asd(t,q)
k1 = 1.2;
k2 = 1.3;
k3 = 1.5;
qd = [-k1*(12.7734-q(1))*cos(q(3));
-k1*(12.7734-q(1))*sin(q(3));
-k2*(0.3097+(13.6167-q(2))*sin(t))];
end
Then in the main script, you call
q10=[7.1892; 11.4189; 0.5826];
t1 = pi:0.1:2*pi;
[t,q] = ode45(@asd,t1,q10);
You did not specify k1 k2 k3 so I randomly selected them.
  1 件のコメント
sangita kamat
sangita kamat 2018 年 9 月 8 日
編集済み: sangita kamat 2018 年 9 月 8 日
@Aquatris, thank you very much for your help. Shall do.its perfectly working . thanks once again

サインインしてコメントする。

カテゴリ

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

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by