matlab求微分解析解function dy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);[T,Y]=ode45('rigid',[0 12],[0 1 1])为何说没有定义y,

来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/27 07:14:51
matlab求微分解析解functiondy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2

matlab求微分解析解function dy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);[T,Y]=ode45('rigid',[0 12],[0 1 1])为何说没有定义y,
matlab求微分解析解
function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
[T,Y]=ode45('rigid',[0 12],[0 1 1])
为何说没有定义y,

matlab求微分解析解function dy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);[T,Y]=ode45('rigid',[0 12],[0 1 1])为何说没有定义y,
这下子就可以了:
function hh
[T,Y]=ode45(@rigid,[0 12],[0 1 1])
plot(T,Y)
function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);