C++循环中变量x无法更新？引力计算程序故障排查

核心问题：x不更新的根源 #

你在初始化initialX时，x的初始值是0.0，后续用户输入x的有效数值后，initialX并没有同步更新，始终保持0。所以循环里的x = x - initialX / 100.0等价于x = x - 0，x自然不会产生任何变化。

解决方法：把initialX = x;的赋值语句移到用户输入x的代码之后，让initialX真正保存用户输入的初始边缘间距。

其他关键错误修正 #

• 引力公式错误：万有引力的正确公式是F = G*m1*m2/(d²)，你写成了除以sqrt(d)，这会导致计算结果完全偏离预期，甚至在d极小时触发数值异常。
• 循环计算顺序错误：当前循环先计算f再更新x和d，第一次循环时d的初始值是0，会触发除以0的未定义行为。应调整顺序：先计算当前的d，再计算引力f，最后更新x。
• 输入验证不严谨：原代码用if判断负数输入，若用户第二次仍输入负数，程序不会再次检查。改用while循环可以强制用户输入合法的正数，同时处理输入字符等非法情况。
• 提示文字错误：输入半径时的提示误写为“Kilograms（千克）”，应改为“Meters（米）”，避免误导用户。

修复后的完整代码 #

#include <iostream>
#include <iomanip>
#include <cmath>
#include <limits>

using namespace std;

const double G = 6.67408e-11;

int main(){
    double m1 = 0.0;
    double m2 = 0.0;
    double r1 = 0.0;
    double r2 = 0.0;
    double x = 0.0;
    double d = 0.0;
    double initialX = 0.0;
    double f = 0.0;
    
    // 输入第一个球体质量，循环验证合法性
    cout << "Enter the Mass of the First Sphere in Kilograms: ";
    while (!(cin >> m1) || m1 < 0) {
        cout << "Error! Enter a Positive Number!" << endl;
        cout << "Enter the Mass of the First Sphere in Kilograms: ";
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(), '\n');
    }
    
    cout << endl << "Enter the Radius of the First Sphere in Meters: ";
    while (!(cin >> r1) || r1 < 0) {
        cout << "Error! Enter a Positive Number!" << endl;
        cout << "Enter the Radius of the First Sphere in Meters: ";
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(), '\n');
    }
    
    cout << endl << "Enter the Mass of the Second Sphere in Kilograms: ";
    while (!(cin >> m2) || m2 < 0) {
        cout << "Error! Enter a Positive Number!" << endl;
        cout << "Enter the Mass of the Second Sphere in Kilograms: ";
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(), '\n');
    }
    
    cout << endl << "Enter the Radius of the Second Sphere in Meters: ";
    while (!(cin >> r2) || r2 < 0) {
        cout << "Error! Enter a Positive Number!" << endl;
        cout << "Enter the Radius of the Second Sphere in Meters: ";
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(), '\n');
    }
    
    cout << endl << "Enter the Distance Between the Edges of the Two Spheres in Meters: ";
    while (!(cin >> x) || x < 0) {
        cout << "Error! Enter a Positive Number!" << endl;
        cout << "Enter the Distance Between the Edges of the Two Spheres in Meters: ";
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(), '\n');
    }
    initialX = x; // 保存用户输入的初始边缘间距
    cout << endl;
    
    cout << setw(15) << "Distance" << setw(15) << "Force" << endl;
    cout << "-------------------------------------" << endl;
    
    for(int i = 0; i <= 10; i++){
        d = x + r1 + r2; // 计算当前球心距
        f = (G * m1 * m2) / (d * d); // 正确的引力计算公式
        cout<< setw(15) << fixed << setprecision(2) << d << setw(15) << scientific << setprecision(5) << f << endl;
        x = x - initialX / 100.0; // 更新边缘间距，可添加判断避免x为负
    }
    return 0;
}

额外提示 #

如果需要避免边缘间距x变为负数（比如两球不能重叠），可以在更新x后添加判断：

x = x - initialX / 100.0;
if(x < 0) x = 0;

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