b=4;c=3;a=1;d=2;
odefun=@(t,x)[x(1)*((1-x(1))*(b-c)+(b+d+a)*x(1)*x(3)-a*x(3));
x(2)*((c-b)*x(1)+(b+d+a)*x(1)*x(3)-a*x(3));
x(1)+x(2)+x(3)-1];%微分方程组
M=[1 0 0;0 1 0;0 0 0];%质量矩阵
options=odeset('mass',M);%对以DAE问题,mass属性必须设置
x0=[0.1;0.4;0.5];%初值
[t,x]=ode15s(odefun,[0 20],x0,options);%这里好像不能使用ode45
figure('numbertitle','off','name','DAE demo—by Matlabsky')
plot(t,x)
legend('x1(t)','x2(t)','x3(t)')
http://scholar.google.com/schhp?hl=zh-CN
重要通知:郭老师组的开题答辩时间安排在5月8日(周三)下午1:30,地点另行通知。请在5月5日之前将PPT发给郭老师。邮箱:guojj@swufe.edu.cn
收到请回复!!!
- 1
- 2
前往页