MATLAB中subs函数对原方程代入正常但对导数代入失效问题求助

问题根因 #

你遇到的不是subs函数失效，也不是替换操作对求导后的表达式无效，本质是两个问题导致的异常表现：

• 方程参数设置完全不符合物理规则，目标方程在气态可行域V > b内没有实根，牛顿迭代直接发散
• 发散后迭代值进入非物理区间，后续计算出现无穷值、非数（NaN），被误判为subs没有返回有效值

原因拆解 #

你当前写的是范德华气体状态方程求体积的逻辑，参数校验如下：
你用到的R=0.4615、a=1.703、b=0.00169是水的范德华常数，而你设置的温度T=50K远低于水的三相点温度273.16K，此时水为固态，根本不存在气态解，迭代第一步就会得到负的体积值，后续计算要么数值溢出，要么出现无穷值，最终触发除以零的报错。

实际测试你的原代码，第一次迭代中subs(F,V,oldguess)和subs(DF,V,oldguess)都能正常返回数值，不存在替换失效的问题。

修复方案 #

1. 修正参数到合理范围，比如将温度调整到水的临界温度以上（如600K），保证方程存在气态实根
2. 初始猜测值改用理想气体状态方程的计算结果，更贴近真实根，降低发散概率
3. 增加迭代次数上限、导数绝对值保护、误差绝对值判断逻辑，避免异常情况卡死或提前终止

修复后可运行代码 #

close all; 
clear all; 
clc;
max_iter = 100; % 新增迭代上限，防止发散死循环
iter_cnt = 0;
root = 0;
T = 600; % 修正温度到合理范围，原T=50K无气态解
R = 0.4615;
b = 0.00169;
a = 1.703;
P = 10000;
syms V
F = (R*T)/(V-b)- a/(V^2)- P;
DF = diff(F,V);
% 初始猜测值改用理想气体状态方程计算，更合理
oldguess = R*T/P; 
newguess = oldguess;
realerror = 100;
acceptable_error = 0.0001;

while realerror > acceptable_error && iter_cnt < max_iter
    iter_cnt = iter_cnt +1;
    f_V = vpa(subs(F,V,oldguess));
    fd_V= vpa(subs(DF,V,oldguess));
    % 导数接近0时提前终止，避免除以零报错
    if abs(fd_V) < 1e-9
        error('迭代点导数接近0，无法继续更新，请调整初始猜测值');
    end
    newguess = oldguess - f_V/fd_V;
    % 误差加绝对值，避免负误差提前终止迭代
    realerror = abs((newguess-oldguess)/newguess)*100;
    oldguess = newguess;
end

if iter_cnt >= max_iter
    disp('迭代达到上限未收敛，请检查参数或初始值');
else
    root = newguess;
    disp(['收敛得到体积根为：',char(vpa(root,6)),' m³/kg']);
end

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