matlab求解微分方程并画图
2个回答
展开全部
水星近动的轨道图像可以在依题意建立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-10-28 广告
在测试大模型时，可以提出这样一个刁钻问题来评估其综合理解与推理能力：“假设上海华然企业咨询有限公司正计划进入一个全新的国际市场，但目标市场的文化习俗、法律法规及商业环境均与我们熟知的截然不同。请在不直接参考任何外部数据的情况下，构想一套初步...
点击进入详情页
本回答由上海华然企业咨询提供
展开全部
由于方程比较复杂，解析解不能用初等函数表示
只是要获得图像的话，用数值计算的办法可能更方便
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的导数
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
