C++实现欧拉法(Euler's Method)求解微分方程是否正确？

核心逻辑错误修正 #

你代码里的欧拉迭代逻辑存在致命问题，当前代码行：

approx_x[i+1] = approx_x[i+1] + dt*v;

中approx_x[i+1]在赋值前是未初始化的栈垃圾值，完全不符合欧拉法的计算规则。
针对你求解的微分方程 $\frac{dx}{dt}=v$，欧拉法的正确更新公式为：下一时刻位置 = 当前时刻位置 + 步长 * 导数，所以需要把这行改成：

approx_x[i+1] = approx_x[i] + dt * v;

其他优化建议 #

• 输入合法性校验缺失：你定义了数组最大长度max_number_steps = 100000，但没有校验用户输入的步数是否超出这个范围，也没有校验终止时间T是否为正数，一旦输入非法值会直接导致数组越界、计算无意义。建议在用户输入完成后添加校验逻辑：

  if (number_steps <= 0 || number_steps > max_number_steps) {
      cout << "输入步数不合法，有效范围为1~" << max_number_steps << endl;
      return 1;
  }
  if (T <= 0) {
      cout << "终止时间T必须大于0" << endl;
      return 1;
  }
  

• 冗余头文件清理：当前引入的<string>、<stdio.h>、<unistd.h>、<math.h>、<stdlib.h>、<stdarg.h>、<assert.h>、<array>头文件均未在代码中使用，可以直接删除，减少编译冗余。
• 绘图参数修正：你调用gnuplot绘图时的图例名称写的是Velocity，但你绘制的是位置随时间变化的数值解，建议改成Euler Numerical Solution更符合实际含义。

修改后完整可运行代码 #

#include <iostream>
#include "gnuplot.cxx"
using namespace std;

int main()
{
    double v, x0, T, dt;
    int number_steps;
    const int max_number_steps = 100000;
    int i=0;
    
    double value_t [max_number_steps];
    double approx_x [max_number_steps];
 
    cout<<"Enter Initial Condition"<< endl;
    cout<<"Velocity(v) = ";
    cin>> v;
    cout<<"x0 = ";
    cin >> x0;
    cout<<"Final time(T) = ";
    cin >> T;   
    cout << "number steps = ";
    cin >> number_steps;
    
    // 输入合法性校验
    if (number_steps <= 0 || number_steps > max_number_steps) {
        cout << "输入步数不合法，有效范围为1~" << max_number_steps << endl;
        return 1;
    }
    if (T <= 0) {
        cout << "终止时间T必须大于0" << endl;
        return 1;
    }

    dt = T/number_steps;
    value_t[0]= 0.0;
    approx_x[0] = x0;

    for (i= 0; i < number_steps; i++)
    {
        value_t [i+1] = dt*(i+1);
        approx_x[i+1] = approx_x[i] + dt*v;
    }
    
    gnuplot_one_function("Euler Numerical Solution", "linespoints", "Time(t)", "Position(x)", value_t, approx_x, number_steps);

    return 0;
}

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