【请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y】
答案:2 悬赏:10 手机版
解决时间 2021-03-01 19:53
- 提问者网友：献世佛
- 2021-02-28 22:20
【请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y】
最佳答案
- 五星知识达人网友：不如潦草
- 2021-02-28 23:39
那y1,y2呢?那个在y(3)是y的三次导的意思.上面所有的y1,y2我认为也是求导的意思（你有一些没有加括号）Dsolve('D3y+t*Dy*D2y+t^2*Dy*y^2=exp(-t*y)','y(0)=2','Dy(0)=0','D2y(0)=0')但是没有精确解.你可以做数值解,下面做数值在Matlab下输入：edit zhidao_rk4_4.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dydt=zhidao_rk4_4(t,y) %y,y的一阶导,二阶导变量分别用y(1),y(2),y(3)表示 dydt=[y(2);y(3);exp(-t*y(1))-t*y(2)*y(3)-t^2*y(2)*y(1)^2]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 在Matlab下面输入：t_end=10; x0=[0;0;0]; [t,y]=ode45('zhidao_rk4_4',[0,t_end],x0); plot(t,y(:,1)); xlabel('t'); ylabel('y');
全部回答
- 1楼网友：笑迎怀羞
- 2021-02-28 23:59
我检查一下我的答案
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
大家都在看
推荐资讯