matlab矩阵求求解微分方程组
发布网友
发布时间:2023-01-01 13:08
我来回答
共3个回答
热心网友
时间:2023-09-08 21:05
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
热心网友
时间:2023-09-08 21:05
没有解析解
用数值解吧
热心网友
时间:2023-09-08 21:06
你还是去买一个书看吧