如何在MATLAB中正确测试无穷级数敛散性并求收敛和？

原代码错误原因 #

你之前的代码误用了vpa()转换symsum的结果，导致误判：

• symsum(1/n^0.5, n, 1, inf)本身会返回符号型的inf（表示级数发散），但用vpa()转换后，符号无穷大会被转为一个极大的浮点数（如1.797693134862316e+308），此时isinf()无法识别这个浮点数为无穷大，因此错误输出"convergent"。

通用解决方案 #

直接基于符号计算结果判断敛散性，步骤如下：

1. 定义符号变量与级数通项
   用syms声明求和变量，写出级数的通项表达式。

2. 计算无穷级数的符号和
   使用symsum(term, n, start, inf)计算级数和，其中start是求和起始项（通常为1）。

3. 判断敛散性并输出结果

   • 若结果为inf或-inf，说明级数发散；
   • 若结果为有限的符号表达式或数值，说明级数收敛，可同时输出无穷和。

修正后的示例代码（针对你的原问题）

syms n
term = 1/n^0.5;
y = symsum(term, n, 1, inf);
if isinf(y)
    fprintf('divergent\n');
else
    fprintf('convergent\n');
    fprintf('无穷和为: %s\n', char(y));
end

运行后会正确输出divergent。

更多场景示例 #

收敛的正项级数（p级数，p=2）

syms n
term = 1/n^2;
y = symsum(term, n, 1, inf);
if isinf(y)
    fprintf('divergent\n');
else
    fprintf('convergent\n');
    fprintf('无穷和为: %s\n', char(y)); % 输出 pi^2/6
end

条件收敛的交错级数

syms n
term = (-1)^(n+1)/n;
y = symsum(term, n, 1, inf);
if isinf(y) || isnan(y)
    fprintf('divergent\n');
else
    fprintf('convergent\n');
    fprintf('无穷和为: %s\n', char(y)); % 输出 log(2)
end

无法直接得到解析解的收敛级数

若symsum返回特殊函数形式（如黎曼ζ函数），可通过vpa()获取数值近似：

syms n
term = 1/n^3;
y = symsum(term, n, 1, inf); % 返回 zeta(3)
if isinf(y)
    fprintf('divergent\n');
else
    fprintf('convergent\n');
    fprintf('无穷和解析形式: %s\n', char(y));
    fprintf('无穷和数值近似: %.6f\n', vpa(y));
end

补充：部分和极限验证法 #

对于少数symsum无法直接判断的级数，可计算部分和的极限：

syms n N
term = 1/n^0.5;
partial_sum = symsum(term, n, 1, N);
limit_result = limit(partial_sum, N, inf);
if isinf(limit_result)
    fprintf('divergent\n');
else
    fprintf('convergent\n');
end

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