使用for循环试错绘制球体重力异常图时绘图异常求助

你的代码核心问题是每次循环只计算了单个x点的重力异常值，而非整条x序列对应的完整曲线——这就导致plot(x, gal1(j))实际在绘制单个点（肉眼几乎看不到），而plot(x, gal1)也只能读取最后一次循环留下的零散值，自然只能出一条线。下面是修正后的代码，以及关键改进点的说明：

修正后的代码 #

x = [-3; -2.5; -2; -1.5; -1; -0.5; 0; 0.5; 1.5]; % 保留你的原始x序列
hold on % 提前开启图窗保留，避免循环内重复调用

% 遍历质量m的所有取值（0.1到0.9，步长0.1）
for m_val = 0.1:0.1:0.9
    % 遍历重心深度h的所有取值（0.1到0.9，步长0.1）
    for h_val = 0.1:0.1:0.9
        % 计算当前(m,h)组合下，整条x对应的重力异常曲线
        gal1 = (6.67 * m_val * h_val) ./ (x.^2 + h_val^2).^(3/2);
        % 绘制曲线并添加图例标记，方便区分不同参数组合
        plot(x, gal1, 'DisplayName', ['m=',num2str(m_val),', h=',num2str(h_val)]);
    end
end

hold off
legend('Location','best') % 自动放置最优位置的图例
xlabel('水平位置x')
ylabel('重力异常值')
title('不同质量m、深度h对应的球体重力异常曲线')

关键改进点说明 #

• 重构循环逻辑：把三重循环简化为m和h的双重循环，对每个参数组合直接计算整条x对应的异常向量，而非单个点的值，确保每次plot都能画出完整曲线。
• 使用点运算保证向量计算：用./、.^这类MATLAB元素级运算，确保对x的每个元素独立计算，得到和x长度一致的异常曲线向量。
• 简化变量存储：直接用m_val和h_val存储当前循环的参数值，不需要额外创建h(i)、m(k)数组，代码更简洁高效。
• 优化图窗控制：提前开启hold on放在循环外，避免重复调用；添加图例和坐标轴标注，让结果更直观可读。

原代码的核心错误解析 #

你原本的gal1(j) = ...仅计算了x(j)这一个点的异常值，每次循环只更新gal1的一个元素——这就导致gal1始终是只有单个有效值的向量，自然无法画出完整曲线。修正后我们直接计算整个x序列的异常值，得到完整的曲线向量，才能实现每次循环绘制一条对应曲线的需求。

内容的提问来源于stack exchange，提问作者Miro Feliciano