使用Maple识别自然数序列对应有理函数的技术问询

嘿，我明白你的困扰——明明是个很简单的数列，用gfun却没识别出来对吧？其实大概率是你用的命令或者输入方式不对，咱们一步步来搞定它，顺便也给你推荐几个替代工具。

用gfun正确识别通项 #

gfun的guess命令其实完全可以处理这个数列，但需要正确传入数列列表。试试下面的代码：

with(gfun):
# 把你的数列做成列表
my_seq := [4/3, 6/5, 8/7, 10/9];
# 调用guess命令，指定变量n
guess(my_seq, n);

执行后应该会直接返回通项公式(2*n+2)/(2*n+1)。如果之前没成功，可能是你没把数列包装成列表，或者没明确指定变量。

如果guess直接出不来，还可以分两步走：先找递推关系，再解递推得到通项：

# 先猜递推式
rec := guessrec(my_seq, a(n));
# 用rsolve解递推，加上初始条件
rsolve({rec, a(1)=4/3}, a(n));

这一步会先得到递推式a(n+1) = (2*n+4)/(2*n+3)*a(n)，再解出通项，结果和直接猜是一样的。

其他可用的Maple工具 #

如果gfun还是不合你心意，试试这些替代方案：

• ListTools[GuessFormula]：这个工具专门用来猜数列的通项公式，用法也很简单：

  with(ListTools):
  GuessFormula(my_seq, n);
  

• numapprox[Fit]：如果数列是有理函数形式，也可以用拟合的方式，指定有理函数的分子分母次数来求解：

  with(numapprox):
  Fit( (a*n + b)/(c*n + d), [1,2,3,4], my_seq, n );
  

  这里指定分子分母都是一次多项式，拟合后就能得到正确的系数。

总结 #

你的数列完全可以用Maple的工具识别出来，核心是用对命令和输入格式。如果gfun不顺手，上面提到的几个工具也能轻松搞定这个问题。

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