Maple中伊藤微积分的微分乘积替换方法及赋值报错疑问

为什么直接赋值会报错？ #

Maple的赋值运算符:=要求左侧必须是合法的单个变量名（比如X__t、B__t这类标识符），而dB__t*dt、dB__t^2是包含乘法、幂次运算符的表达式，不属于可赋值的变量名称范畴，因此直接执行这类赋值会触发Error, illegal use of an object as a name错误。这些伊藤微积分的规则是表达式替换规则，而非变量赋值操作。

如何替换伊藤微积分的规则？ #

你可以使用Maple的subs函数对展开后的dY__t进行批量替换，直接传入所有需要替换的规则即可：

# 先确保dY__t已完成展开化简
dY__t := simplify(expand(dY__t));
# 执行伊藤规则替换
dY__t := subs(
    dB__t*dt = 0,
    dt*dB__t = 0,
    dt^2 = 0,
    dB__t^2 = dt,
    dY__t
);

如果需要更精准的模式匹配（比如避免意外替换复杂乘积项），可以使用subsindets结合模式匹配函数，针对特定表达式结构进行替换：

# 替换dB__t的平方项为dt
dY__t := subsindets(dY__t, power(dB__t, 2), () -> dt);
# 替换dt的平方项为0
dY__t := subsindets(dY__t, power(dt, 2), () -> 0);
# 替换dB__t与dt的乘积项（两种顺序）为0
dY__t := subsindets(dY__t, product(identical(dB__t), identical(dt)), () -> 0);
dY__t := subsindets(dY__t, product(identical(dt), identical(dB__t)), () -> 0);

执行上述操作后，dY__t就会按照伊藤微积分的规则完成化简。

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