Matlab% 太阳系模拟 G = 6.67430e-11; % 万有引力常数 M_sun = 1.989e30; % 太阳质量 M_mercury = 3.3e23; % 水星质量 M_venus = 4.87e24; % 金星质量 M_earth = 5.97e24; % 地球质量 M_mars = 6.39e23; % 火星质量 M_jupiter = 1.898e27; % 木星质量 M_saturn = 5.68e26; % 土星质量 M_uranus = 8.68e25; % 天王星质量 M_neptune = 1.02e26; % 海王星质量 M_pluto = 1.31e22; % 冥王星质量 % 初始位置和速度 P_sun = [0; 0; 0]; P_mercury = [0; 5.7e10; 0]; P_venus = [0; 1.1e11; 0]; P_earth = [0; 1.5e11; 0]; P_mars = [0; 2.2e11; 0]; P_jupiter = [0; 7.8e11; 0]; P_saturn = [0; 1.4e12; 0]; P_uranus = [0; 2.9e12; 0]; P_neptune = [0; 4.5e12; 0]; P_pluto = [0; 5.9e12; 0]; V_sun = [0; 0; 0]; V_mercury = [4.8e4; 0; 0]; V_venus = [3.5e4; 0; 0]; V_earth = [2.98e4; 0; 0]; V_mars = [2.41e4; 0; 0]; V_jupiter = [1.3e4; 0; 0]; V_saturn = [9.7e3; 0; 0]; V_uranus = [6.8e3; 0; 0]; V_neptune = [5.4e3; 0; 0]; V_pluto = [4.7e3; 0; 0]; % 模拟时间和时间步长 t = 0:3600*24*365:3600*24*365*10; dt = 3600*24; % 数值积分 P = [P_sun, P_mercury, P_venus, P_earth, P_mars, P_jupiter, P_saturn, P_uranus, P_neptune, P_pluto]; V = [V_sun, V_mercury, V_venus, V_earth, V_mars, V_jupiter, V_saturn, V_uranus, V_neptune, V_pluto]; M = [M_sun, M_mercury, M_venus, M_earth, M_mars, M_jupiter, M_saturn, M_uranus, M_neptune, M_pluto]; for i = 1:length(t)-1 F = zeros(3, size(P, 2)); for j = 1:size(P, 2) for k = 1:size(P, 2) if j ~= k r = norm(P(:,j)-P(:,k)); F(:,j) = F(:,j) + G*M(j)*M(k)/r^2*(P(:,k)-P(:,j))/r; end end end A = F./M; V = V + A*dt; P = P + V*dt; end % 绘制行星轨道 figure; hold on; plot3(P(1,:), P(2,:), P(3,:), 'k'); grid on; axis equal; view(45, 45)； 代码报错 错误使用：./ 2.矩阵维度必须一致

这个错误通常是由于矩阵的维度不匹配导致的。在这个代码中，错误可能出现在以下几行：

F = zeros(3, size(P, 2));
...
A = F./M;

这里的 `size(P, 2)` 表示矩阵 `P` 的第二个维度的大小，也就是行星数量。因此，`F` 的第二个维度大小应该和行星数量相同。而 `M` 是行星质量的向量，它的大小应该和行星数量相同。如果这两个向量的大小不同，那么在进行除法运算时就会出现维度不匹配的错误。

你可以通过在这些语句之前添加 `disp` 语句来检查这些变量的大小：

disp(size(P, 2));
disp(size(M));

这样就可以找到错误所在，并修复代码。