MATLAB实现Euler法求解微分方程报错：变量x未定义

你遇到的这个错误，本质是在迭代计算导数时没有正确传递当前的x值，导致MATLAB找不到变量x的定义。我帮你修正代码，同时拆解每一步的逻辑，确保你能理解：

首先明确我们要解的微分方程：

y'(x) = -8y(x) + 0.5x + 1/16，初始条件y(0)=2

修正后的完整代码 #

function [xar, yar] = Euler(a, b, ybouco, N)
    % Euler法求解一阶常微分方程
    % 输入参数:
    %   a: x的下限（初始点）
    %   b: x的上限（终点）
    %   ybouco: 初始条件y(a)
    %   N: 迭代步数
    % 输出参数:
    %   xar: 存储所有x值的数组
    %   yar: 存储对应y值的数组
    
    % 计算迭代步长
    h = (b - a) / N;
    
    % 初始化x数组，从a到b按步长h生成
    xar = a:h:b;
    % 初始化y数组，长度与xar一致，首元素为初始值
    yar = zeros(size(xar));
    yar(1) = ybouco;
    
    % 执行Euler迭代
    for n = 1:N
        % 取出当前迭代步的x和y值
        x_current = xar(n);
        y_current = yar(n);
        % 计算微分方程右端的导数值
        dy_dx = -8*y_current + 0.5*x_current + 1/16;
        % 用Euler公式更新下一个y值
        yar(n+1) = y_current + h*dy_dx;
    end
end

错误原因分析 #

你之前的代码应该是直接用了未定义的全局变量x来计算导数，但Euler法需要用到当前迭代步的具体x值（也就是xar(n)），必须在循环里明确取出并传递给导数计算式，不能直接用一个未赋值的x。

测试验证 #

你可以用这段代码测试函数是否正常工作：

% 设定参数：x从0到1，初始y(0)=2，迭代100步
[xar, yar] = Euler(0, 1, 2, 100);
% 输出前5组x和y值查看结果
disp([xar(1:5)' yar(1:5)'])

运行后就能得到正确的xar和yar数组，不会再出现变量未定义的错误。

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