作业帮大家帮我看看这段MATLAB解微分方程的程序m-文件rigid.m的第二行“dy=zeros(3,是说“y1(0)=0,y2(0)=1,y3(0)=1”里头的自变量为0么?解微分方程y1'=y2y3y2'=-y1y3y3'=-0.5y1y2y1(0)=0,y2(0)=1,y3(0)=11、建立m-文件rigid.m
来源：学生作业帮助网 编辑：作业帮 时间：2024/11/17 21:21:49
大家帮我看看这段MATLAB解微分方程的程序m-文件rigid.m的第二行“dy=zeros(3,是说“y1(0)=0,y2(0)=1,y3(0)=1”里头的自变量为0么?解微分方程y1'=y2y3y2'=-y1y3y3'=-0.5y1y2y1(0)=0,y2(0)=1,y3(0)=11、建立m-文件rigid.m
大家帮我看看这段MATLAB解微分方程的程序
m-文件rigid.m的第二行“dy=zeros(3,
是说“y1(0)=0,y2(0)=1,y3(0)=1”里头的自变量为0么?
解微分方程
y1'=y2y3
y2'=-y1y3
y3'=-0.5y1y2
y1(0)=0,y2(0)=1,y3(0)=1
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),'+')
m-文件rigid.m的第二行“dy=zeros(3,
是说“y1(0)=0,y2(0)=1,y3(0)=1”里头的自变量为0么?
大家帮我看看这段MATLAB解微分方程的程序m-文件rigid.m的第二行“dy=zeros(3,是说“y1(0)=0,y2(0)=1,y3(0)=1”里头的自变量为0么?解微分方程y1'=y2y3y2'=-y1y3y3'=-0.5y1y2y1(0)=0,y2(0)=1,y3(0)=11、建立m-文件rigid.m
dy=zeros(3,1); 是先要生成方程组的因变量的位置,zeros(m,n)是生成一个m行n列的零矩阵.
y1(0)=0,y2(0)=1,y3(0)=1这是给的初值 ,ode45(‘函数名’,‘取值范围’,‘初值’）为4、5阶的龙格库塔法求微分方程组的解.本题的为[T,Y]=ode45('rigid',[0 12],[0 1 1]);