matlab阿卡曼公式,阿克曼函数--一个计算方法

在数学上有一个著名的“阿克曼函数”，它是二元函数，其定义式为：

(１)ACK(0, N)＝1＋Ｎ

(２)ACK(M,0)＝ ACK(M-1, 1)(Ｍ>０)

(３)ACK(M, N) ＝ ACK(M-1，ACK(M, N-1)) (Ｍ>0, Ｎ>０)

试用手工求解ACK(3, 7) 的值。

因为这个函数是用递归方式定义的，如果使用递归算法编程求解并不困难。但是，要解这个具体问题，还必须经过将近 70万次 (693964次) 递归调用 !

显然，局限于递归算法，要想用手工求解这个问题是极其困难的，甚至是不可能的！

要想能用手工求解这个困难问题，唯有另想办法，另辟蹊径。

本文给出手工求解这个难题的一种高效算法。

这种高效率算法的基本思路是：先将这个二元函数 ACK(M,N) 的函数值排成一个对应的M行N列的二维数组，并将此数组每一行视为一个数列。只要能求得第M- 1行对应数列的首项值和通项公式：

ACK(M-1，0)和 ACK(M-1，N)= f(N)，

然后，根据式 (2)和式 (3)，即可以得到第M行对应数列的首项值

ACK(M，0)＝ ACK(M-1, 1)

和数列中相邻项的关系公式：

ACK(M，N)＝ ACK(M-1，ACK(M, N-1))= f [ACK(M