用matlab6.5模拟方程组

来源:百度知道 编辑:UC知道 时间:2024/06/01 19:20:59
用matlab6.5模拟方程组 dx/dt= α(y-x)- γyn/(θ+yn) ;dy/dt=c-β(y-x)-δyn/(θ+yn) ymin=0.4, ymax=0.8,α=0.3,β=2,c=0.25,γ=2,δ=30,θ=2,n=3。麻烦详细点,新手上路,谢谢。

%函数文件,保存为fun1.m文件
%令y(1)=y
%y(2)=x
function dy=fun1(t,y)

alpha=0.3;beta=2;c=0.25;garma=2;
delta=30;theta=2;n=3;

dy=[alpha*(y(1)-y(2))-garma*y(1)*n/(theta+y(1)*n);c-beta*(y(1)-y(2))-delta*y(1)*n/(theta+y(1)*n)];
%-------------------------------------------------------------------------------------------------------------------
%
%以下是脚本程序:(t取值为[0,1])
clear
t=0:0.1:1;
y0=[0.4,0]
[t,ya]=ode45('fun1',t,y0)