MATLAB优化算法分享 http://blog.sciencenet.cn/u/Jasion Optimization Problem

博文

matlab ode 常微分方程组

已有 1759 次阅读 2018-1-23 21:47 |个人分类:所见所闻|系统分类:科研笔记

http://cn.mathworks.com/help/matlab/ref/ode45.html


%main.m

clc
clear
a=0.034;
b=0.005;
r=0.274;
x0=[0.2 0.1 0.01 0.092];
tspan=[0 80.1];
[t,x] = ode45('fun',tspan,x0);
plot(t,x(:,1),'r-o',t,x(:,2),'g-*',t,x(:,3),'b-*',t,x(:,4),'k-*')



%fun.m

function dx=fun(t,x)
a=0.034;
b=0.005;
r=0.274;
dx=zeros(4,1);
dx(1) = -a*x(1)*x(2);
dx(2) = a*x(1)*x(2)-b*x(2)-r*x(2);
dx(3) = b*x(2);
dx(4) = r*x(2);

end



http://blog.sciencenet.cn/blog-734913-1096442.html

上一篇:hisat2
下一篇:笔记

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2020-1-20 02:12

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部