Question

用matlab解齐次线性微分方程组

du/dt=v
dv/dt=-u
满足t=0时，（u，v）=（1,0）
用ode23解 并能画出这三个图的matlab程序 希望能详细写出来

Answer 1

用不着ode23

首先dsolve即可求出
u=cost，v=-sint

然后画图即可

clc
s=dsolve('Du=v','Dv=-u','u(0)=1','v(0)=0');
figure(1);
ezplot3('t',s.u,s.v,[0 20*pi])
grid off
figure(2);
ezplot(s.u,s.v,[0 2*pi])

追问

可是老师要求用ode解 想让我们学会这种方法 不在于方程组是什么

追答

clc
dy=@(t,y)([y(2);-y(1)]);
[t,y]=ode23(dy,[0 20*pi],[1 0]);
figure(1);plot3(t,y(:,1),y(:,2)),xlabel('t'),ylabel('u'),zlabel('v');
grid on
figure(2);plot(y(:,1),y(:,2)),xlabel('u'),ylabel('v'),grid on

Answer 2

太难了哑哑哑哑哎呀眼