需要看一下变量的对应关系
更多代码请参考
clc; clear all; close all;
figure(1)
Rp=100;Cph=85;Cpl=0;Cp=10;Bt=40;Fp=40;Mp=20;Ct=10;Ft=20;Mt=15;Cg=15;Tg=40;
% set(0,'defaultfigurecolor','w');
[t,y]=ode45(@(t,y)xuexi(t,y,Cph,Cpl,Cp,Bt,Rp,Ft,Ct,Cg,Tg,Fp,Mp),[0 50],[0.2 0.2 0.2]);
plot3(y(:,1),y(:,2),y(:,3),'r+','linewidth',1);
hold on;
function dydt=xuexi(t,y,Cph,Cpl,Cp,Bt,Rp,Ft,Ct,Cg,Tg,Fp,Mp)
dydt=zeros(3,1);
dydt(1)=y(1)*(y(1)-1)*(Cph-Cpl-Cp-Bt-y(2)*(Rp-Bt)-y(3)*(Fp+Mp));
end
