求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图 就剩10分了、、、、、
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/15 20:32:42
求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图 就剩10分了、、、、、
求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图
求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图 就剩10分了、、、、、
求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图求用matlab解微分方程组 [dx1;dx2]=[0.06 -0.0196;40 -9.9]*[x1-900;x2-10800000] 得X跟T的关系图 就剩10分了、、、、、
function dx = myeq(t,x)
dx=[0.06 -0.0196;40 -9.9]*[x(1)-900;x(2)-10800000]
end
主程序:
x0=[0;0];
tspan = [0,500];
[tt,yy]=ode45(@myeq,tspan,x0)
plot(yy(:,1),yy(:,2)) ;
figure
plot(tt,yy(:,1)) ;
figure
plot(tt,yy(:,1)) ;
初值取得都是0
[x1,x2]=dsolve('Dx1=0.06*(x1-900)-0.0196*(x2-10800000)','Dx2=40*(x1-900)-9.9*(x2-10800000)')
求解后可以直接用plot绘图,因为你没有给初值,有无穷多解,所以没有办法直接绘制。
绘图非常简单,参考 help plot
也可用数值解绘图,但一样需要定解条件,参考 help ode45...
全部展开
[x1,x2]=dsolve('Dx1=0.06*(x1-900)-0.0196*(x2-10800000)','Dx2=40*(x1-900)-9.9*(x2-10800000)')
求解后可以直接用plot绘图,因为你没有给初值,有无穷多解,所以没有办法直接绘制。
绘图非常简单,参考 help plot
也可用数值解绘图,但一样需要定解条件,参考 help ode45
收起
你没给初始条件,我用的0,0,但不太理想,我把函数写出你自己带吧。
建立M文件
function dydt=myfun(t,y)
dydt=[0.06 -0.0196;40 -9.9]*[y(1)-900;y(2)-10800000];
在工作空间输入[t y]=ode45(@myfun,[0 5],[10 10 ])
plot(t,y(1),t,y(2),'r')