题目
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解呢?
dx/dt = c-b*x-w*x*y;
dy/dt = w*x*y-(b+m)*y;
其中c,b,w,m为已知参数,x,y为未知的,分别对t求导,
请问这个常微分方程组要怎么用ode45解呢?
提问时间:2021-03-25
答案
这个方程要解出来还要初始条件吧
不妨设(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);
得到了
不妨设(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);
得到了
举一反三
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
奥巴马演讲不用看稿子.为什么中国领导演讲要看?
想找英语初三上学期的首字母填空练习……
英语翻译
1,人们染上烟瘾,最终因吸烟使自己丧命.
最新试题
- 1陶渊明的简介
- 2“你没有告诉我.我知道了.”加上关联词,把两句话合成一句话.
- 3《贝多芬传》
- 4请问y=-x∧2-3的图像顶点?
- 5求苏教版五年级下册口算、混合运算、解方程、和应用题多一点!
- 6已知椭圆的焦点是F1(-1,0),F2(1,0),P为椭圆上一点,且|F1F2|是|PF1|和|PF2|的等差中项,若点P在第三象限,且∠PF1F2=120°,求tan∠F1PF2的值.
- 715.5化成分数是多少
- 8在三角形中 有一角为30 夹这一角的两边比为根号三比二这个三角形是直角三角形么
- 9两个质数的和是15,则积是()
- 10There was asily wolf in the old house怎么改成复数句
热门考点
- 1恬淡如水 什么意思
- 2初三2道一元2次方程的解答题
- 3英语翻译
- 4用180厘米长的铁丝做一个长方体框架,长·宽·高的比是4:3:2.这个长方体框架的体积是多少?
- 5I wish I would do sth.与I wish I +过去式区别是什么?我搞糊涂了
- 6太阳与地球相距离1.5乘以10的八次方,则太阳光传到地球需要多少时间?
- 7带近反义词的成语
- 81.FeCl2溶液呈浅绿色,其中存在着下列平衡:Fe2+ +2H2O===Fe(OH)2+2H+往该溶液滴加稀盐酸,发生的变化是(溶液由浅绿色变为深绿色)
- 9Has somebody been drinking my fat可以修改成Has somebody drunk my 意思一样吗
- 10《项脊轩志》最后一段在表达情感的方式上有什么特点?