单片机C语言求平方根.pdfVIP

. 在单片机 中要开平方 .可以用到下面算法 : 算法 1: 本算法只采用移位、加减法、判断和循环实现，因为它不需要浮点运算，也不需要乘除运算，因此 可以很方便地运用到各种芯片上去。 我们先来看看 10 进制下是如何手工计算开方的。 先看下面两个算式， x = 10*p + q (1) 公式 (1) 左右平方之后得： x^2 = 100*p^2 + 20pq + q^2 (2) 现在假设我们知道 x^2 和 p ，希望求出 q 来，求出了 q 也就求出了 x^2 的开方 x 了。 我们把公式 (2) 改写为如下格式： q = (x^2 - 100*p^2)/(20*p+q) (3) 这个算式左右都有 q ，因此无法直接计算出 q 来，因此手工的开方算法和手工除法算法一样有一步需 要猜值。 我们来一个手工计算的例子：计算 1234567890 的开方 首先我们把这个数两位两位一组分开，计算出最高位为 3 。也就是 (3) 中的 p，最下面一行的 334 为余 数，也就是公式 (3) 中的(x^2 - 100*p^2) 近似值 3 | 12 34 56 78 90 9 | 3 34 下面我们要找到一个 0-9 的数 q 使它最接近满足公式 (3) 。我们先把 p 乘以 20 写在 334 左边： 3 q | 12 34 56 78 90 9 6q| 3 34 我们看到 q 为 5 时(60+q*q) 的值最接近 334 ，而且不超过 334 。于是我们得到： 3 5 | 12 34 56 78 90 9 65| 3 34 | 3 25 9 56 接下来就是重复上面的步骤了，这里就不再啰嗦了。 这个手工算法其实和 10 进制关系不大，因此我们可以很容易的把它改为二进制，改为二进制之后， 公式 (3) 就变成了： q = (x^2 - 4*p^2)/(4*p+q) (4) 我们来看一个例子，计算 100( 二进制 1100100) 的开方： 1 0 1 0 | 1 10 01 00 1 100| 0 10 | 0 00 | 10 011001| 10 01 0 00 这里每一步不再是把 p 乘以 20 了，而是把 p 乘以 4 ，也就是把 p 右移两位，而由于 q 的值只能为 0 ;. . 或者 1，所以我们只需要判断余数 (x^2 - 4*p^2) 和(4*p+1) 的大小关系，如果余数大于等于 (4*p+q) 那么该上一个 1，否则该上一个 0。 下面给出完成的 C 语言程序， 其中 root 表示 p ，rem 表示每步计算之后的余数， divisor 表示 (4*p+1) ， 通过 a30 取 a 的最高 2 位，通过 a=2 将计算后的最高 2 位剔除。其中 root 的两次 1 相当 于 4*p 。程序完全是按照手工