如何让MATLAB的symsum识别发散级数并给出不收敛提示？

MATLAB的符号数学工具箱symsum函数可以计算收敛级数的解析结果，比如执行symsum(2^k/factorial(k),k,1,Inf)会返回exp(2)-1；但对于像log(k+1)/k这类发散级数，symsum只会返回未求值的原表达式，不会像Mathematica的Sum函数那样明确提示级数发散。要实现类似的提示功能，需要自定义封装函数，结合收敛性判别逻辑来实现。

实现思路 #

核心是在调用symsum后，判断结果是否为未求值表达式，再通过级数收敛性判别方法（通项极限判别、积分判别、比较判别等）判定级数是否发散，最后给出对应提示：

• 先调用symsum计算级数，若返回未求值的symsum表达式，进入收敛性判断流程；
• 先检查通项的极限：若当k→∞时通项极限不为0，直接判定级数发散；
• 对正项级数，使用积分判别法：若通项的无穷积分发散，则级数发散；
• 针对特定通项（如log(k+1)/k），可以补充比较判别法，和已知的发散级数（如1/k）对比判定；
• 若无法判定收敛性，返回原表达式并给出提示。

自定义函数示例 #

function result = symsum_with_check(a_k, k, lower, upper)
    syms k;
    result = symsum(a_k, k, lower, upper);
    
    % 判断是否返回未求值的symsum表达式
    if isa(result, 'sym') && strcmp(char(result), ...
            ['symsum(' char(a_k) ', ' char(k) ', ' char(lower) ', ' char(upper) ')'])
        % 第一步：通项极限判别
        term_limit = limit(a_k, k, Inf);
        if term_limit ~= 0
            warning('级数发散：通项当k→∞时极限不为0');
            return;
        end
        
        % 第二步：正项级数积分判别法
        if isAlways(a_k > 0, [k >= lower])
            integral_result = int(a_k, k, lower, Inf);
            if isinf(integral_result)
                warning('级数发散：积分判别法判定无穷积分发散');
                return;
            end
        end
        
        % 第三步：针对log(k+1)/k的比较判别示例
        if strcmp(char(a_k), 'log(k+1)/k')
            comp_term = 1/k;
            comp_ratio_limit = limit(a_k / comp_term, k, Inf);
            % 1/k是调和级数，已知发散，若比值极限为正有限值则原级数发散
            if comp_ratio_limit > 0 && isinf(symsum(comp_term, k, lower, upper))
                warning('级数发散：与调和级数1/k比较判别，比值极限为正，原级数发散');
                return;
            end
        end
        
        % 若以上判别都无法确定，给出提示并返回原结果
        warning('无法判定级数收敛性，返回未求值表达式');
    end
end

使用方法 #

在MATLAB中定义上述函数后，调用方式和symsum一致：

syms k
symsum_with_check(log(k+1)/k, k, 1, Inf)

执行后会弹出对应的发散警告。

需要注意的是，这个自定义函数的判别逻辑可以根据需求扩展，比如添加更多级数类型的判别规则，覆盖更多发散或收敛的场景。

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