Question

如何用matlab解微分方程并作图

方程是dy/dt=(qin*Cain-q*y)/V-K*y 然后qin=50；q=50;Cain=1.3;V=1000;K=0.05; 这个要用ODE45怎么作图?每次都有地方不对，画不出图来。 求详细代码的写法，取点是x轴（t值）从0开始 然后是方程组dy1/dt=(qin*Cain-q*y1)/V-K*y1y2 dy2/dt=(qin*Cbin-q*y2)/V-K*y1y2 qin=50；q=50;Cain=1.3;Cbin=3;V=1000;K=0.05 两个东西是连着的，第一个没有头绪，所以第二个也做不出来，在线等跪求大神帮忙

Answer 1

第一个是线性方程，可以直接求解析解
用dsolve
求数值解：
%%%%%%%%%%%%%%%%M-file
以下要单存一个文件在工作路径下%%%%%%%%%%%%%function
dotx=myfun(t,y)
qin=50；q=50;Cain=1.3;V=1000;K=0.05;
dotx=(qin*Cain-q*y)/V-K*y;
%%%%%%%%%%%%%%%%M-file
以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command
窗口
[t,y]=ode45('myfun'，[0,t_final],y_int)%y_int
初值
第二个相同：
求数值解：
%%%%%%%%%%%%%%%%M-file
以下要单存一个文件在工作路径下%%%%%%%%%%%%%function
dotx=myfun2(t,y)
qin=50；q=50;Cain=1.3;Cbin=3;V=1000;K=0.05
dotx=[(qin*Cain-q*y(1))/V-K*y(1)*y(2);(qin*Cbin-q*y(2))/V-K*y(1)*y(2)];
%%%%%%%%%%%%%%%%M-file
以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command
窗口
[t,y]=ode45('myfun2'，[0,t_final],[y1_int,y2_int])%yi_int是
初值