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/05/12 15:16:58
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,

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);