matlab求解微分方程并画图
展开全部
水星近动的轨道图像可以在依题意建立ODE方程组后，使用matlab常微分方程组求解器ode45得出：
代码的图和运行效果如图
clc;clear all;
a=1;b=0.06;
mer_orbit=@(theta,r)[r(2);
b*r(1)^2-r(1)+a];
figure(1);clf;
[theta,r]=ode45(mer_orbit, [0:.01*pi:30*pi],[.1,0]);
axis equal;axis off;hold on;
plot(0,0,'*r');
%%太阳位置
[x, y]=pol2cart(theta,1./r(:,1));
%极坐标转换为直角坐标
plot(x,y)
追答
那个速度与时间的曲线等等再弄
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
东莞大凡
2024-11-14 广告
标定板认准大凡光学科技,专业生产研发厂家，专业从事光学影像测量仪,光学投影测量仪.光学三维测量仪,光学二维测量仪,光学二维测量仪,光学三维测量仪,光学二维测量仪.的研发生产销售。东莞市大凡光学科技有限公司创立于 2018 年，公司总部坐落于...
点击进入详情页
本回答由东莞大凡提供
展开全部
由于方程比较复杂，解析解不能用初等函数表示
只是要获得图像的话，用数值计算的办法可能更方便
fun=@(t,y) 1.44*(10^9)*(1-y).*exp(-109170./(8.314*t));
[T,Y] = ode23t(fun,500:600,0);
DY=fun(T,Y);
AX=plotyy(T,Y,T,DY);
set(get(AX(1),'Xlabel'),'String','T');
set(get(AX(1),'Ylabel'),'String','\alpha');
set(get(AX(2),'Ylabel'),'String','d\alpha/dT');
T是自变量，Y是变量也就是alpha
DY是，Y对T的导数
只是要获得图像的话，用数值计算的办法可能更方便
fun=@(t,y) 1.44*(10^9)*(1-y).*exp(-109170./(8.314*t));
[T,Y] = ode23t(fun,500:600,0);
DY=fun(T,Y);
AX=plotyy(T,Y,T,DY);
set(get(AX(1),'Xlabel'),'String','T');
set(get(AX(1),'Ylabel'),'String','\alpha');
set(get(AX(2),'Ylabel'),'String','d\alpha/dT');
T是自变量，Y是变量也就是alpha
DY是，Y对T的导数
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
