求常微分方程组的解。 [图]A、建立函数文件ty.m。 functi...
求常微分方程组的解。
A、建立函数文件ty.m。 function dy=ty(t, y) dy=[ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; 调用函数文件: >> [t, y]=ode45(@ty, [0, 12], [0, 1, 1]); >> plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+')
B、建立函数文件ty.m。 function dy=ty(t, y) dy=[ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; 调用函数文件: >> clear >> h=@ty; >> [t, y]=ode45(h, [0, 12], [0, 1, 1]); >> plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+')
C、ty=@(t, y) [ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; [t, y]=ode45(ty, [0, 12], [0, 1, 1]); plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+')
D、[t, y]=ode45(@(t, y) [ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)], [0, 12], [0, 1, 1]); plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+')