MATLAB数值求解微分 f'=fg-0.3f,g'=-fg,f(0)=0.02,g(0)=0.9_百度...
发布网友
发布时间：2024-10-19 14:47
共1个回答
热心网友
时间：2024-11-29 20:33
给出的微分方程组，可以用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、编程后运行，可以得到如下结果
热心网友
时间：2024-11-29 20:33
给出的微分方程组，可以用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、编程后运行，可以得到如下结果
