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)
追答
那个速度与时间的曲线等等再弄
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
Sievers分析仪
2025-02-09 广告
是的。传统上，对于符合要求的内毒素检测，最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线；必须有重复的阴性控制；每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪，这些步骤可以通过使用预嵌入的内毒素标准...
点击进入详情页
本回答由Sievers分析仪提供
展开全部
由于方程比较复杂，解析解不能用初等函数表示
只是要获得图像的话，用数值计算的办法可能更方便
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的导数
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
|
广告 您可能关注的内容 |
