简单版matlab实现两个速度不同的球体的转动动画

以下是一个简单版的 Matlab 代码实现两个速度不同的球体的转动动画：

% 设定球体半径和旋转速度
r1 = 1; % 球体1半径
r2 = 0.5; % 球体2半径
w1 = 0.1; % 球体1旋转速度
w2 = 0.2; % 球体2旋转速度

% 初始化画布
figure;
axis equal;
axis([-1.5 1.5 -1.5 1.5 -1.5 1.5]);
xlabel('x');
ylabel('y');
zlabel('z');
hold on;

% 生成球体表面坐标
[x1, y1, z1] = sphere(10);
[x2, y2, z2] = sphere(10);
x1 = x1 * r1;
y1 = y1 * r1;
z1 = z1 * r1;
x2 = x2 * r2;
y2 = y2 * r2;
z2 = z2 * r2;

% 绘制球体
h1 = surf(x1, y1, z1, 'FaceColor', 'b', 'EdgeColor', 'none');
h2 = surf(x2, y2, z2, 'FaceColor', 'r', 'EdgeColor', 'none');

% 设置视角
view(45, 30);

% 动画循环
t = 0;
while true
    % 计算球体旋转角度
    theta1 = w1 * t;
    theta2 = w2 * t;
    
    % 计算球体旋转矩阵
    R1 = [cos(theta1), -sin(theta1), 0; sin(theta1), cos(theta1), 0; 0, 0, 1];
    R2 = [cos(theta2), -sin(theta2), 0; sin(theta2), cos(theta2), 0; 0, 0, 1];
    
    % 旋转球体并更新绘制
    set(h1, 'XData', x1 * R1(1,1) + y1 * R1(1,2) + z1 * R1(1,3));
    set(h1, 'YData', x1 * R1(2,1) + y1 * R1(2,2) + z1 * R1(2,3));
    set(h1, 'ZData', x1 * R1(3,1) + y1 * R1(3,2) + z1 * R1(3,3));
    set(h2, 'XData', x2 * R2(1,1) + y2 * R2(1,2) + z2 * R2(1,3));
    set(h2, 'YData', x2 * R2(2,1) + y2 * R2(2,2) + z2 * R2(2,3));
    set(h2, 'ZData', x2 * R2(3,1) + y2 * R2(3,2) + z2 * R2(3,3));
    drawnow();
    
    % 更新时间
    t = t + 0.1;
end

这个代码同样会生成一个包含两个不同大小和旋转速度的球体的动画，但是没有使用 `now()` 函数获取当前时间，而是使用固定时间步长 `0.1` 来更新时间。可以使用 `Ctrl+C` 停止动画循环。