作业帮matlab求微分解析解
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/09 12:26:25
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);