subs函数返回0而非指定项的MATLAB符号计算问题

问题描述 #

学习MathWorks教程时，执行以下代码发现：x可正常替换为THE_X，但diff(x,t)被替换成0而非预期的THE_XD：

clear
syms m k b F
syms t x(t)
syms THE_X THE_XD THE_XDD
HOLDER_list = [THE_X,THE_XD,THE_XDD];
actual_list = [x, diff(x,t), diff(x,t,2)];

v =diff(x,t)
v_new = subs(v, actual_list, HOLDER_list)

原因分析 #

当x被定义为符号函数x(t)时，diff(x,t)会生成符号函数形式的表达式，直接用数组传入替换对时，MATLAB的符号匹配机制无法正确识别diff(x,t)作为替换目标，最终返回默认未匹配的结果0。

修复方案 #

将替换对象以单元格数组形式传入subs，同时提前把导数赋值给变量明确匹配关系，修改后的代码如下：

clear
syms m k b F
syms t x(t)
syms THE_X THE_XD THE_XDD

% 提前定义导数变量，明确匹配对象
xd = diff(x,t);
xdd = diff(x,t,2);

v = xd;
% 使用单元格数组传入替换对
v_new = subs(v, {x, xd, xdd}, {THE_X, THE_XD, THE_XDD})

补充说明 #

如果替换关系较多，也可以用结构体形式实现更直观的批量替换：

replace_map = struct('x', THE_X, 'diff(x,t)', THE_XD, 'diff(x,t,2)', THE_XDD);
v_new = subs(v, replace_map);

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