请用matlab帮我写一个模拟二维空间的三体运动的代码，并生成gif文件，在文件中画出三个星球及其轨迹

好的，以下是一个简单的模拟三体运动的Matlab代码，它可以模拟三个星球在二维空间中的运动轨迹，并生成一个 GIF 动画文件。

% 设置初始条件
G = 6.67430e-11;
m1 = 5.97e24;
m2 = 7.35e22;
m3 = 1e23;
r1 = [0,0];
r2 = [3.84e8,0];
r3 = [0,3.84e8];
v1 = [0,0];
v2 = [0,1023];
v3 = [1023,0];
t = 0;
dt = 3600;

% 设置画图参数
figure('Renderer', 'painters', 'Position', [100 100 800 600]);
filename = 'three_body.gif';

% 模拟运动
for i = 1:1000
    % 计算距离和力
    d12 = norm(r2-r1);
    d13 = norm(r3-r1);
    d23 = norm(r3-r2);
    F12 = G * m1 * m2 / d12^2;
    F13 = G * m1 * m3 / d13^2;
    F23 = G * m2 * m3 / d23^2;
    
    % 计算加速度
    a1 = (F12 * (r2-r1) / d12 + F13 * (r3-r1) / d13) / m1;
    a2 = (-F12 * (r2-r1) / d12 + F23 * (r3-r2) / d23) / m2;
    a3 = (-F13 * (r3-r1) / d13 - F23 * (r3-r2) / d23) / m3;
    
    % 更新位置和速度
    v1 = v1 + a1 * dt;
    v2 = v2 + a2 * dt;
    v3 = v3 + a3 * dt;
    r1 = r1 + v1 * dt;
    r2 = r2 + v2 * dt;
    r3 = r3 + v3 * dt;
    
    % 绘制轨迹
    hold on;
    plot(r1(1), r1(2), 'bo', 'MarkerSize', 8, 'MarkerFaceColor', 'b');
    plot(r2(1), r2(2), 'ro', 'MarkerSize', 8, 'MarkerFaceColor', 'r');
    plot(r3(1), r3(2), 'go', 'MarkerSize', 8, 'MarkerFaceColor', 'g');
    plot([r1(1) r2(1)], [r1(2) r2(2)], 'b-', 'LineWidth', 1);
    plot([r1(1) r3(1)], [r1(2) r3(2)], 'b-', 'LineWidth', 1);
    plot([r2(1) r3(1)], [r2(2) r3(2)], 'b-', 'LineWidth', 1);
    xlim([-4e8 4e8]);
    ylim([-4e8 4e8]);
    xlabel('x (m)');
    ylabel('y (m)');
    title(sprintf('Three body simulation (time = %.0f s)', t));
    drawnow;
    
    % 保存为 GIF 动画
    frame = getframe(gcf);
    im = frame2im(frame);
    [A,map] = rgb2ind(im,256);
    if i == 1
        imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0.1);
    else
        imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0.1);
    end
    
    % 更新时间
    t = t + dt;
end

在这个代码中，我们首先设置了三个星球的初始位置和速度，然后在一个循环中模拟它们的运动轨迹。在每个时间步长内，我们先计算出三个星球之间的距离和力，然后根据牛顿定律计算出每个星球的加速度，进而更新它们的速度和位置。最后，我们使用 `plot` 命令绘制出三个星球的位置和轨迹，并保存为一个 GIF 动画文件。

请注意，这个代码只是一个简单的示例，实际的三体运动模拟需要更加复杂的数值计算和优化技术，以保证精确度和效率。