matlab矩阵求求解微分方程组
来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:综合作业 时间:2024/07/02 19:16:42
matlab矩阵求求解微分方程组
{用matlab编程}
数值求解微分方程组![](http://img.wesiedu.com/upload/5/3c/53cc748fb7d3fe53ca4a930aee53cdb5.jpg)
求程序
{用matlab编程}
数值求解微分方程组
![](http://img.wesiedu.com/upload/5/3c/53cc748fb7d3fe53ca4a930aee53cdb5.jpg)
求程序
![matlab矩阵求求解微分方程组](/uploads/image/z/410965-61-5.jpg?t=matlab%E7%9F%A9%E9%98%B5%E6%B1%82%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84)
function myCalc
clear;
clc;
close all;
opt=odeset('RelTol',1E-5,'AbsTol',[1E-6 1E-6]);
[T Y]=ode45(@getdY,[0 100],[0.02 0.98],opt);
figure;
hold on;
grid on;
box on;
plot(T,dY(:,1),'r')
plot(T,dY(:,2),'b')
return
function dY=getdY(t,Y)
dY=zeros(2,1);
dY(1)=Y(1)*Y(2)-0.3*Y(1);
dY(2)=-Y(1)*Y(2);
return
clear;
clc;
close all;
opt=odeset('RelTol',1E-5,'AbsTol',[1E-6 1E-6]);
[T Y]=ode45(@getdY,[0 100],[0.02 0.98],opt);
figure;
hold on;
grid on;
box on;
plot(T,dY(:,1),'r')
plot(T,dY(:,2),'b')
return
function dY=getdY(t,Y)
dY=zeros(2,1);
dY(1)=Y(1)*Y(2)-0.3*Y(1);
dY(2)=-Y(1)*Y(2);
return