发布网友 发布时间:2024-10-08 08:53
共1个回答
热心网友 时间:2024-11-20 06:20
给出的微分方程组,可以用matlab的ode45()函数来求解。求解方法:
1、建立自定义微分方程组函数,odefun(t,y),其主要内容
dz(1)=y(1)*y(2)-0.3*y(1); %f(t)=y(1),g(t)=y(2)
dz(2)=-y(1)*y(2);
2、建立求解微分方程组函数命令文件,ode45_main,其主要内容
t0=0;tf=10;
tspan=t0:0.1:tf;
x0=[0.02,0.9]; %初值
options = odeset('RelTol',1e-4,'AbsTol',1e-4); %误差设定1e-4
[t,y]= ode45(@odefun,tspan,x0,options)
figure(1),plot(t,y) %t—f(t),t—g(t)关系图
figure(2),plot(y) %f(t)—g(t)关系图
3、编程后运行,可以得到如下结果
追问可否将程序一次写齐 我总是有点问题 谢谢啊