请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y^2 =e^(-ty) ,y(0)=2 ,y1(0)=y2(0)=0并且绘制出y(t)的曲线.【上式中 ,y(3)中(3)的位置我从书上看处于幂的位置.但是为什么要加一个括号啊?所以我下
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/28 19:29:21
![请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y^2 =e^(-ty) ,y(0)=2 ,y1(0)=y2(0)=0并且绘制出y(t)的曲线.【上式中 ,y(3)中(3)的位置我从书上看处于幂的位置.但是为什么要加一个括号啊?所以我下](/uploads/image/z/5245104-48-4.jpg?t=%E8%AF%B7%E7%BB%99%E5%87%BA%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%9A%84MATLAB%E5%91%BD%E4%BB%A4y%5E%283%29+%2B+t%2Ay1%2Ay2+%2B+t%5E2%2A%28y1%29%2Ay%5E2+%3De%5E%28-ty%29+%2Cy%280%29%3D2+%2Cy1%280%29%3Dy2%280%29%3D0%E5%B9%B6%E4%B8%94%E7%BB%98%E5%88%B6%E5%87%BAy%28t%29%E7%9A%84%E6%9B%B2%E7%BA%BF.%E3%80%90%E4%B8%8A%E5%BC%8F%E4%B8%AD+%2Cy%283%29%E4%B8%AD%283%29%E7%9A%84%E4%BD%8D%E7%BD%AE%E6%88%91%E4%BB%8E%E4%B9%A6%E4%B8%8A%E7%9C%8B%E5%A4%84%E4%BA%8E%E5%B9%82%E7%9A%84%E4%BD%8D%E7%BD%AE.%E4%BD%86%E6%98%AF%E4%B8%BA%E4%BB%80%E4%B9%88%E8%A6%81%E5%8A%A0%E4%B8%80%E4%B8%AA%E6%8B%AC%E5%8F%B7%E5%95%8A%3F%E6%89%80%E4%BB%A5%E6%88%91%E4%B8%8B)
请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y^2 =e^(-ty) ,y(0)=2 ,y1(0)=y2(0)=0并且绘制出y(t)的曲线.【上式中 ,y(3)中(3)的位置我从书上看处于幂的位置.但是为什么要加一个括号啊?所以我下
请给出求解微分方程的MATLAB命令
y^(3) + t*y1*y2 + t^2*(y1)*y^2 =e^(-ty) ,y(0)=2 ,y1(0)=y2(0)=0
并且绘制出y(t)的曲线.
【上式中 ,y(3)中(3)的位置我从书上看处于幂的位置.但是为什么要加一个括号啊?所以我下不了手】
麻烦大家帮我个忙..
谢谢wacs5耐心解释.
y1 ,y2 书上没加括号的.只有y(3)加了括号的.
syms t y; %声明符号变量
dsolve('D3y+t*y1*y2+t^2*y1*y2^2=exp(-t*y)','y(0)=2','Dy(0)=0','D2y(0)=0')
figure;plot3(y(:,1),y(:,2),y(:,3))
由于我是新手,我是模仿书上一步一步来写的.望您指正.
请给出求解微分方程的MATLAB命令y^(3) + t*y1*y2 + t^2*(y1)*y^2 =e^(-ty) ,y(0)=2 ,y1(0)=y2(0)=0并且绘制出y(t)的曲线.【上式中 ,y(3)中(3)的位置我从书上看处于幂的位置.但是为什么要加一个括号啊?所以我下
那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');