除法器的实现（恢复余数、不恢复余数、级数展开、Newton-Raphson）

参考《基于FPGA的数字信号处理（第二版）》，对其中的除法器算法进行仿真。

1 基于恢复余数（Restoring）算法的除法器

Restoring 算法无法直接用于与有符号数，对于有符号数需要先转换为无符号数。然后根据除数与被除数的符号判断商与余数的符号。

1.1 原理

被除数（Dividend）5621 / 除数（Divisor）3 = 商（Quotient）1873 … 余数（Remainder）2

       1873
     _______
  3 |  5621
    /  3000
       ____
       2621
       2400  
       ____
        221
        210
        ___
         11
          9
          _
          2

$$r_i=r_{i+1}-q_i\cdot D\cdot 10^i$$

$q_i$通过部分余数的正负确定，如果余数为负表示$q_i$较大，需要将其恢复到部分余数为正以获取正确的$q_i$，这便是“恢复”的含义。

用恢复余数算法表示上式计算过程：

得 $q_3{q}_2{q}_1{q}_0$ = 1873，R=2。

以上是十进制算法，对于二进制：

$$r_i=r_{i+1}-q_i\cdot D\cdot 2^i$$

$q_i$只能取0或1，所以每位只需一次运算即可判断是否恢复。

此外还需确定各参数的字长：假定除数D为 nbit，商Q为 nbit，则