绘制给出的微分方程组的图像解,可以这样来实现。
syms x(t) y(t)
Dx=diff(x,1);Dy=diff(y,1);
[x,y]=dsolve(Dy==x,Dx==y,x(0)==1,y(0)==1);
t=0:0.1:10;
x=eval(x);
y=eval(y);
figure(1)
plot(x,y,'LineWidth',1.5),grid on
xlabel('x'),ylabel('y')
figure(2)
plot(t,x,'r*',t,y,'k-','LineWidth',1.2),grid on
xlabel('t'),ylabel('x(t),y(t)')
legend('x(t)','y(t)',2)
运行上述代码后,得到的图像