如何用matlab解二元微分方程组,分别做出x-t,y-t的坐标图? 10
1个回答
展开全部
% 求解微分方程(时间范围0-2)
dx=inline('[(2-3*x(1)+x(2)).*x(1); (4-x(2)+2*x(1)).*x(2)]','t','x');
[t,x]=ode45(dx,[0 2],[0;3])
y=x(:,2);
x=x(:,1);
% 绘制x-t, y-t曲线
% 注意:由于初始条件x(0)=0,计算出导数dx/dt=0,求出的x一直为0
figure(1)
subplot(2,1,1)
plot(t,x)
xlabel('Time (sec)')
ylabel('x')
subplot(2,1,2)
plot(t,y)
xlabel('Time (sec)')
ylabel('y')
% 绘制x-y曲线
% 由于x一直为0,x-y图为一条竖线
figure(2)
plot(x,y)
xlabel('x')
ylabel('y')
dx=inline('[(2-3*x(1)+x(2)).*x(1); (4-x(2)+2*x(1)).*x(2)]','t','x');
[t,x]=ode45(dx,[0 2],[0;3])
y=x(:,2);
x=x(:,1);
% 绘制x-t, y-t曲线
% 注意:由于初始条件x(0)=0,计算出导数dx/dt=0,求出的x一直为0
figure(1)
subplot(2,1,1)
plot(t,x)
xlabel('Time (sec)')
ylabel('x')
subplot(2,1,2)
plot(t,y)
xlabel('Time (sec)')
ylabel('y')
% 绘制x-y曲线
% 由于x一直为0,x-y图为一条竖线
figure(2)
plot(x,y)
xlabel('x')
ylabel('y')
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询