matlab微分方程1、建立m-文件rigid.m如下： function dy=rigid(t,y) dy=zeros

来源：学生作业帮编辑：灵鹊做题网作业帮分类：综合作业时间：2024/05/02 02:48:39

matlab微分方程
1、建立m-文件rigid.m如下：
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);
2、取t0=0,tf=12,输入命令：
[T,Y]=ode45('rigid',[0 12],[0 1 1]);
plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')
我 1.想问下在句中解微分方程m文件中的dy=zeros(3,1);是什么意思?
2.另外在句子 plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')中 Y（：,1）,Y（：,2）.是为什么这么分?
3.tf是随便找的值么,原题中是没有.

dy=zeros(3,1);是占个空位置.
以下程序和你的程序是等效的,比较一下就可以体会出来.
function hh
[T,Y]=ode45(@rigid,[0 12],[0 1 1]);
plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')
figure
plot(T,Y)
function dy=rigid(t,y)
%dy=zeros(3,1);
dy1=y(2)*y(3);
dy2=-y(1)*y(3);
dy3=-0.51*y(1)*y(2);
dy=[dy1;dy2;dy3]

matlab微分方程1、建立m-文件rigid.m如下： function dy=rigid(t,y) dy=zeros matlab中对于含有导数的函数怎样建立M文件.例如dy=-y+t+1, [T,Y]=ode45(@rigid,[0 1],[rand(15,1)]); 用matlab解M*DY/DT=mg-kv微分方程请帮忙用Matlab求解微分方程dy/dt=[k^(t-1)*y-d]y rigid 微分方程,tanx dy/dx=1+y matlab ode45求解微分方程 D2y-0.01*Dy.^2+2*y=sin(t), 我想用matlab中的ode45 求解常微分方程(Dy)^2-3*Dy+2y=1,y(0)=1,Dy(0)=0. 怎么用MATLAB求解如Dy = y+1/y 的微分方程用MATLAB求解微分方程dy/dx-2y/(x+1)=(x+1)^5/2 微分方程(x+y)(dx-dy)=dx+dy的通解