Octave ODE程序出现未明错误：lsode调用时't'未定义

Got it，我帮你分析下这个Octave 4.2.x版本里的ODE求解报错问题～

这个error: 't' undefined near line 2 column 14错误本质是变量作用域+旧版本Octave的参数传递行为导致的：

• 你在调用lsode时，把t的定义嵌在了参数里：(t=linspace(0,5,50)')，在Octave 4.2.1/4.2.2这类旧版本中，这种赋值表达式的作用域仅限于参数传递的上下文，并没有把t真正定义到当前工作区中。
• 当lsode调用你的f函数时，f需要访问t变量，但此时t不在工作区，也没有被正确传递进来，所以就触发了未定义错误。

有两种简单可靠的修复方式，任选其一即可：

方案1：先单独定义t变量，再调用lsode #

把t的赋值和lsode调用分开，确保t存在于当前工作区：

# 先定义时间向量，确保它在工作区中可见
t = linspace(0,5,50)';
# 再调用lsode求解ODE
x = lsode("f", 2, t);
# 绘图
plot(t,x)

方案2：改用函数句柄（更推荐，兼容性更好） #

在旧版本Octave中，使用函数句柄@f替代字符串"f"来指定ODE函数，参数传递更稳定，也更符合现代语法：

# 先定义时间向量
t = linspace(0,5,50)';
# 使用函数句柄@f调用lsode
x = lsode(@f, 2, t);
plot(t,x)

如果你不想单独维护f.m文件，还可以直接用匿名函数：

# 定义匿名形式的ODE函数
f = @(x,t) -exp(-t)*x^2;
t = linspace(0,5,50)';
x = lsode(f, 2, t);
plot(t,x)

这两种方案在Octave 4.2.1/4.2.2版本中都能正常运行，不会再出现t未定义的错误，并且能正确求解并绘制出ODE的数值解曲线。

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