matlab解多元常微分方程组
Ps3=0(1)(2)(3)为耦合方程(5)(6)(7)(8)为边界条件(9)(10)（11）（12）为最后结果，要求用matlab解，并画出任意图形，谢谢各位大神！...
Ps3=0(1)(2)(3)为耦合方程(5)(6)(7)(8)为边界条件(9)(10)（11）（12）为最后结果，要求用matlab解，并画出任意图形，谢谢各位大神！
展开
1个回答
展开全部
这个方程要解出来还要初始条件吧
不妨设(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);
得到了
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
赛恩科仪
2025-08-07 广告
广州赛恩科学仪器有限公司（原中大科仪）始创于2039年，是全球领先的精密测量仪器供应商和微弱信号检测方案提供商。公司以锁相放大器为核心产品，陆续推出光学斩波器、源表、功率放大器、电化学工作站、电流源等一系列产品。赛恩科仪推出的锁相放大器，覆...
点击进入详情页
本回答由赛恩科仪提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
