【老生谈算法】Matlab解微分方程组源码.doc
版权申诉
52 浏览量
2023-05-11
21:03:31
上传
评论
收藏 120KB DOC 举报
求解微分方程
2
'
(0) 1
x x
x
ì
= -
í
=
î
可以按如下程序得到数值解 1:
clear;clc;
[t,x]=ode45('xprim1',[0 1],1);
plot(t,x,'o-');
xlabel('time t0=0,tt=1');
ylabel('x values x(0)=1');
function xprim=xprim1(t,x)
xprim=-x.^2;
输出
求解微分方程组
1 1 1 2
2 2 1 2
1
2
' 0.1 0.01
' 0.02 0.04
(0) 30
(0) 20
x x x x t
x x x x t
x
x
= - +
ì
ï
= - + +
ï
í
=
ï
ï
=
î
可以按如下程序得到数值解
clear;clc;
[t,x]=ode45('xprim1',[0 20],[30;20]);
plot(t,x);
xlabel('time t0=0,tt=20');
ylabel('x values x1(0)=30,x2(0)=20');
function xprim=xprim1(t,x)
xprim=[x(1)-0.1*x(1)*x(2)+0.01*t;
-x(2)+0.02*x(1)*x(2)+0.04*t];
资源评论
阿里matlab建模师
- 粉丝: 3374
- 资源: 2786
