qq_58208978
2021-06-20 21:18
采纳率: 100%
浏览 11

为什么这个程序是一条直线而不是曲线

function dx=star(t,x) %建立函数G=1;m1=1;m2=1;m3=1

r12=sqrt((x(5)^2-x(1)^2)^2+(x(7)^2-x(3)^2)^2);

r23=sqrt((x(9)^2-x(5)^2)^2+(x(11)^2-x(7)^2)^2);

r13=sqrt((x(9)^2-x(1)^2)^2+(x(11)^2-x(3)^2)^2);

dx=[x(2);(-G*m2*(x(5)-x(1)))/((r12)^3)-(G*m3*(x(9)-x(1)))/(r13)^3; x(4);(-G*m2*(x(7)-x(3)))/((r12)^3)-(G*m3*(x(11)-x(3)))/(r13)^3; x(6);(-G*m1*(x(1)-x(5)))/((r12)^3)-(G*m3*(x(9)-x(5)))/(r23)^3; x(8);(-G*m1*(x(3)-x(7)))/((r12)^3)-(G*m3*(x(11)-x(7)))/(r23)^3; x(10);(-G*m1*(x(1)-x(9)))/((r13)^3)-(G*m2*(x(5)-x(9)))/(r23)^3; x(12);(-G*m1*(x(3)-x(11)))/((r13)^3)-(G*m2*(x(7)-x(11)))/(r23)^3;]end

figure('name','三体运动');

x0=[0.3;-0.466;0.243;-0.433;0.970;-0.466;-0.243;-0.433;0;0.932;0;0.866] %调用函数

[t,x]=ode45(@star,[-1000,1000],x0)

plot(x(:,1),x(:,3),x(:,5),x(:,7),x(:,9),x(:,11))

  • 点赞
  • 写回答
  • 关注问题
  • 收藏
  • 邀请回答

1条回答 默认 最新

  • 有问必答小助手 2021-06-23 16:40
    已采纳

    你好,我是有问必答小助手,非常抱歉,本次您提出的有问必答问题,技术专家团超时未为您做出解答

    本次提问扣除的有问必答次数,已经为您补发到账户,我们后续会持续优化,扩大我们的服务范围,为您带来更好地服务。

     点赞

    点赞 评论

相关推荐 更多相似问题