二元二階常微分方程組求解

發(fā)布時(shí)間：2025-11-03 | 來(lái)源：互聯(lián)網(wǎng)轉(zhuǎn)載和整理

二元二階常微分方程組的求解，可以用龍格-庫(kù)塔法求解其數(shù)值解。

求解方法：

1、自定義二元二階常微分方程組降價(jià)函數(shù)

2、確定初始條件，x1(0)=0，dx1(0)/dt=0，x2(0)=0，dx1(0)/dt=0

3、確定時(shí)間t的范圍，t【0，10】

4、確定時(shí)間t的步長(zhǎng)，h=0.1

5、使用runge_kutta龍格-庫(kù)塔法函數(shù)或ode45函數(shù)，求解其數(shù)值解

6、繪制x1(t)和x2(t)曲線圖

x0=[0;0;0;0];

a=0;b=10;h=0.1;

[t,x]=runge_kutta(@func,x0,h,a,b);

plot(t,x(:,1),t,x(:,3));

gridon

xlabel('t'),ylabel({'x1';'x2'})

h=legend('x1(t)','x2(t)','Location','northwest');

運(yùn)行結(jié)果