搜搜做题作业网常微分方程组matlab程序设计
来源：学生作业帮 编辑：搜搜做题作业网作业帮 分类：综合作业 时间：2024/06/25 03:00:31
常微分方程组matlab程序设计
初始值x=0.7,y1=1.2,y2=0.8,y3=0.7
要求画出x,y1,y2,y3的图像,这个数值可能不对,告诉我程序就行了
第二个方程里等号的右边的d应该是一个系数,可以是0.9
初始值x=0.7,y1=1.2,y2=0.8,y3=0.7要求画出x,y1,y2,y3的图像,这个数值可能不对,告诉我程序就行了
第二个方程里等号的右边的d应该是一个系数,可以是0.9
第二个方程右边的dy1是什么意思?
再问： 右边的d应该是一个系数，可以是0.9
再答： clc
f=@(t,y)([y(1)*(1-y(1))-.5*y(1)*(y(2)+y(3)+y(4))/(1+2*y(1));
2.3*y(1)*(y(2)+1.2*y(3)+.9*y(4))/(1+2*y(1))-.9*y(2)-2.7*y(2)*y(4)/(y(2)+y(3)+y(4));
2.7*y(2)*y(4)/(y(2)+y(3)+y(4))-1.3*y(3);
1.6*y(3)-.7*y(4)]);
[t,y]=ode45(f,[0 100],[0.7 1.2 .8 .7]);
plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4)),xlabel('t')
legend('x','y1','y2','y3')
再问： 右边的d应该是一个系数，可以是0.9
再答： clc
f=@(t,y)([y(1)*(1-y(1))-.5*y(1)*(y(2)+y(3)+y(4))/(1+2*y(1));
2.3*y(1)*(y(2)+1.2*y(3)+.9*y(4))/(1+2*y(1))-.9*y(2)-2.7*y(2)*y(4)/(y(2)+y(3)+y(4));
2.7*y(2)*y(4)/(y(2)+y(3)+y(4))-1.3*y(3);
1.6*y(3)-.7*y(4)]);
[t,y]=ode45(f,[0 100],[0.7 1.2 .8 .7]);
plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4)),xlabel('t')
legend('x','y1','y2','y3')