# runge kutta for 2 order ODE

3 ビュー (過去 30 日間)
Rio Bratasena 2022 年 11 月 3 日

how can I solve this equation with runge kutta method

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

### 回答 (2 件)

Hiro Yoshino 2022 年 11 月 3 日
Let us assume that and re-write the problem as follows: Let's transform it so it can be used in an ode setting: These are used in the function (vdp1) as follows:
solve this the time interval of [0 10] with initial values of [0.2 0] for x_1, x_2 respectively.
[t,x] = ode45(@vdp1,[0 10],[0.2 0])
t = 89×1
0 0.0001 0.0003 0.0004 0.0005 0.0011 0.0018 0.0024 0.0030 0.0062
x = 89×2
0.2000 0 0.2000 -0.0001 0.2000 -0.0001 0.2000 -0.0002 0.2000 -0.0002 0.2000 -0.0005 0.2000 -0.0007 0.2000 -0.0010 0.2000 -0.0012 0.2000 -0.0024
plot(t,x(:,1),'-o',t,x(:,2),'-o')
title('Solution with ODE45');
xlabel('Time t');
ylabel('Solution x');
legend('x_1','x_2'); function dxdt = vdp1(t,x)
% ODE
dxdt = [x(2); -2*(x(2)+x(1))];
end
##### 0 件のコメント表示非表示 -1 件の古いコメント

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

Sam Chak 2022 年 11 月 3 日   Assuming that the symbol is the Dirac impulse. please check if the following responses are expected when you run it with the Runge-Kutta solver.
A = [0 1;
-2 -2];
B = [0; 0.2];
C = eye(2);
D = [0; 0];
sys = ss(A, B, C, D);
impulse(sys) サインインしてコメントする。

### カテゴリ

Find more on Runge Kutta Methods in Help Center and File Exchange

### Community Treasure Hunt

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

Start Hunting!