大家帮我看看这段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/06/23 23:54:02
![大家帮我看看这段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](/uploads/image/z/2105622-54-2.jpg?t=%E5%A4%A7%E5%AE%B6%E5%B8%AE%E6%88%91%E7%9C%8B%E7%9C%8B%E8%BF%99%E6%AE%B5MATLAB%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%9A%84%E7%A8%8B%E5%BA%8Fm-%E6%96%87%E4%BB%B6rigid.m%E7%9A%84%E7%AC%AC%E4%BA%8C%E8%A1%8C%E2%80%9Cdy%3Dzeros%283%2C%E6%98%AF%E8%AF%B4%E2%80%9Cy1%280%29%3D0%2Cy2%280%29%3D1%2Cy3%280%29%3D1%E2%80%9D%E9%87%8C%E5%A4%B4%E7%9A%84%E8%87%AA%E5%8F%98%E9%87%8F%E4%B8%BA0%E4%B9%88%3F%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8By1%27%3Dy2y3y2%27%3D-y1y3y3%27%3D-0.5y1y2y1%280%29%3D0%2Cy2%280%29%3D1%2Cy3%280%29%3D11%E3%80%81%E5%BB%BA%E7%AB%8Bm-%E6%96%87%E4%BB%B6rigid.m)
大家帮我看看这段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]);