求助：(exp(x)-exp(-x))/(2*x)泰勒级数计算函数陷入无限循环

嘿，我来帮你搞定这个无限循环的问题！你的代码卡壳的核心原因是错误地用级数和与真实函数值的差值来判断停止条件——当级数和逼近真实值到浮点数精度极限时，这个差值会在eps附近来回震荡，永远没法严格小于eps，自然就停不下来了。咱们来一步步修正它：

问题拆解 #

先回忆一下目标函数的泰勒展开式：

(exp(x)-exp(-x))/(2x) = Σ(x^(2k))/( (2k+1)! )，其中k从0到∞

你的初始项sum=1是对的（对应k=0时的项），但停止条件选得不对：

1. 当sum和真实函数值f(xo)在浮点数精度下几乎相等时，abs(sum-f(xo))会在eps级别波动，永远满足不了>eps的反向条件，导致循环无限跑。
2. 直接计算xo^n/factorial(n+1)效率低，还容易出现阶乘溢出的问题。

修正后的代码 #

我调整了停止条件，还优化了项的计算方式（用递推代替重复计算阶乘和幂次），同时处理了x=0的特殊情况（避免除以0）：

function k = tayser(xo)
    % 计算(exp(x)-exp(-x))/(2*x)的泰勒级数，直到项的绝对值小于机器epsilon
    if xo == 0
        k = 1; % x=0时函数极限为1，单独处理避免除以0
        return;
    end
    
    sum_val = 1; % 初始项：k=0时，x^0/(1!) = 1
    term = xo^2 / factorial(3); % 第一个新增项：k=1时，x²/3!
    n = 2; % 当前项对应的x的幂次
    
    % 用当前项的绝对值作为停止条件：当项小到机器精度以下，再添加也不会改变sum
    while abs(term) > eps
        sum_val = sum_val + term;
        n = n + 2;
        % 递推计算下一项：利用前一项的结果，避免重复计算阶乘和高次幂
        term = term * xo^2 / (n * (n + 1));
    end
    
    k = sum_val;
end

关键修改点 #

• 停止条件优化：不再用级数和与真实值的差值，而是判断当前要添加的项的绝对值是否小于eps。当项小到这个程度时，它对sum_val的浮点数表示已经没有影响了，循环可以安全停止。
• 递推计算项：每一项都基于前一项计算，比如第k项到k+1项，只需要乘以x²/[(2k+2)*(2k+3)]，既提升了效率，又减少了数值误差。
• 特殊情况处理：x=0时直接返回极限值1，避免原函数的除以0错误。

测试验证 #

比如运行tayser(1)，结果会和(exp(1)-exp(-1))/(2*1)几乎完全一致，而且不会陷入无限循环。运行tayser(0)也能正确返回1。

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