代码

function GalaxyRotationSimulation()
    % 参数设置
    num_stars = 1000; % 恒星数量
    galaxy_radius = 1; % 银河系半径
    rotation_speed = 0.05; % 旋转速度

    % 生成银河系中的恒星分布
    theta = 2 * pi * rand(num_stars, 1); % 角度
    r = galaxy_radius * sqrt(rand(num_stars, 1)); % 半径
    z = 0.1 * (rand(num_stars, 1) - 0.5); % 添加一些z轴方向的随机性
    x = r .* cos(theta);
    y = r .* sin(theta);

    % 初始化图形
    figure;
    h = scatter3(x, y, z, 1, 'w', 'filled');
    axis equal;
    axis([-1.5 1.5 -1.5 1.5 -0.5 0.5]);
    set(gca, 'Color', 'k');
    xlabel('X');
    ylabel('Y');
    zlabel('Z');
    title('Galaxy Rotation Simulation');

    % 旋转动画
    dt = 0.01; % 时间步长
    for t = 1:1000
        theta = theta + rotation_speed * dt; % 更新角度
        x = r .* cos(theta); % 计算新的x坐标
        y = r .* sin(theta); % 计算新的y坐标

        set(h, 'XData', x, 'YData', y, 'ZData', z); % 更新图形数据
        drawnow; % 刷新图形
        pause(dt); % 暂停以创建动画效果
    end
end

说明

参数设置：
num_stars: 设置银河系中恒星的数量。
galaxy_radius: 设置银河系的半径。
rotation_speed: 设置银河系旋转的速度。

生成银河系中的恒星分布：
使用极坐标生成恒星的随机分布。
将极坐标转换为直角坐标，并添加一些z轴方向的随机性以模拟银河系的厚度。

初始化图形：
使用scatter3函数绘制三维散点图。
设置坐标轴范围和背景颜色。

旋转动画：
定义时间步长dt。
在循环中更新恒星的位置，以模拟银河系的旋转。
使用drawnow函数更新图形

效果