神马作文网matlab求微分解析解
来源：学生作业帮 编辑：神马作文网作业帮 分类：综合作业 时间：2024/11/11 18:57:49
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 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);
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);