matlab ode45 求解常微分方程组
dx/dt = c-b*x-w*x*y;
dy/dt = w*x*y-(b+m)*y;
其中c,b,w,m为已知参数,x,y为未知的,分别对t求导,
请问这个常微分方程组要怎么用ode45解呢?
人气:217 ℃ 时间:2020-05-13 16:18:55
解答
这个方程要解出来还要初始条件吧
不妨设(x(0),y(0))=(A,B);
还有用ode45只能得到数值解,不能给出解析式表达
那么先建立一个函数文件
test_fun.m
function dx=tsst_fun(x,t)
dx(1)= c-b*x(1)-w*x(1)*x(2);
dx(2)=w*x(1)*x(2)-(b+m)*x(2);
那么调用ode45
x0=(0,0);
t0=0:0.01:10;
;%根据你所想要的求得值设定t0,间隔是任意的,与求解所用的步长无关,
[x,t]=ode45(@test_fun,t0,x0);
得到了
推荐
- 一个底边长8厘米,高5厘米的三角形,按1:2画在纸上,画在纸上的面积是()
- 变成一般疑问句(没把握的别答) We must finish the work now.____ you_____ the work now?
- 儿童公园里,梨树比桃树少百分之十,桃树比柳树多百分之十,三种树木中谁最多?谁最少?(列出算式)
- 为什么太阳直射180度经线全球分属两天
- 三氧化二铁、氧化铜、二氧化锰怎样用一种试剂区别?
- 1-丁烯,顺2-丁烯,反2-丁烯,异丁烯哪个更稳定
- 在分别写有0、1、2、3、4、5、6、7、8、9的十张卡片中随意抽取一张,则抽到数字小于5的概率是 _.
- 把49分拆成几个自然数的和,这几个自然数的连乘积最大是多少
猜你喜欢