大家好,又见面了,我是你们的朋友全栈君。
上篇博客介绍了Matlab求解常微分方程组解析解的方法:博客地址
微分方程组复杂时,无法求出解析解时,就需要求其数值解,这里来介绍。...以下内容按照Matlab官方文档提供的方程来展开(提议多看官方文档)
介绍一下核心函数ode45()
一般形式:[t,y] = ode45(odefun,tspan,y0) 其中 tspan = [t0...t,y]=ode45(@(t,y) 2*t,tspan,y0); %定义函数y'=2*t,使用ode45求解
plot(t,y,'-o'); %绘制求得的数值曲线
说明:简单的odefun参数就是这个形式...0-20
y0=[2;0];%定义因变量的初值,当x=0时,y1=2,y2=y1'=0;
[t,y]=ode45(@odefun,tspan,y0); %使用ode45求解
%%下面为作图过程,不解释...更多形式
讲到这里,大部分我们用到的微分方程形式都可以求解了,Matlab还支持带有时变项和额外参数的微分方程求解,这里不再赘述,大家可以自行参阅官方文档。