Question

用matlab解方程

Answer 1

【1】可以用dsolve（）函数求得，其解析值。然后再用plot（）函数绘制t～x(t)，t～y(t)的图像。实现代码

syms x(t) y(t)

Dx=diff(x,1);Dy=diff(y,1);

[x,y]=dsolve(Dx==2*x+3*y,Dy==2*x+y,x(0)==-2.7,y(0)==2.8);

t=0:0.1:10;x=eval(x);y=eval(y);

plot(t,x,t,y),hold on

xlabel('t'),ylabel('x(t),y(t)')

【2】可以用ode45（）函数求得，其数值值。然后再用plot（）函数绘制t～x(t)，t～y(t)，t～z(t)的图像。实现代码

x0=[0 1 1];

[t,x]=ode45(@func,[0 12],x0); 

disp([t x])

plot(t,x),grid on

xlabel('t'),ylabel('x(t),y(t),z(t)')

运行结果