微分方程matlab绘图,用matlab解微分方程组并作图

微分方程matlab绘图,用matlab解微分方程组并作图共回答了19个问题采纳率:89.5%functiondx=appollo(t,x)mu=1/82.45;mustar=1-mu;r1=sqrt((x(1)+mu)^2+x(3)^2);r2=sqrt((x(1)-mustar)^2+x(3)^2);dx=[x(2)2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3x(4)-2*x(2)+…

共回答了19个问题采纳率:89.5%

function dx=appollo(t,x)

mu=1/82.45;

mustar=1-mu;

r1=sqrt((x(1)+mu)^2+x(3)^2);

r2=sqrt((x(1)-mustar)^2+x(3)^2);

dx=[x(2)

2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3

x(4)

-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3];

————————————————————————————–

x0=[1.2;0;0;-1.04935751];%x0(i)对应与xi的初值

options=odeset(‘reltol’,1e-8);%该命令的另一种写法是options=odeset;options.reltol=1e-8;

tic

[t,y]=ode45(@appollo,[0,20],x0,options);%t是时间点,y的第i列对应xi的值,t和y的行数相同

toc

plot(y(:,1),y(:,3))%绘制x1和x3,也就是x和y的图形

title(‘Appollo卫星运动轨迹’)

xlabel(‘X’)

ylabel(‘Y’)

1年前

10

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。
如需转载请保留出处:https://bianchenghao.cn/38099.html

(0)
编程小号编程小号

相关推荐

发表回复

您的电子邮箱地址不会被公开。 必填项已用*标注