数值微分法（DDA）详解

数值微分法（$DDA$）

一、原理

1. 假定直线的起点、终点分别为：$(x_0,y_0)$，$(x_1,y_1)$，且都为整数。
   
   则过端点$P_0(x_0,y_0)$，$P_1(x_1,y_1)$的直线段$L:y=kx+b$，直线斜率为：
   $$k=\frac{y_1-y_0}{x_1-x_0}(x_1\ne x_0)$$

2. 假设$x$已知，即从$x$的起$x_0$开始，沿$x$方向前进一个像素（步长= 1)，可以计算出相应的y值。因为像素的坐标是整数，所以y值还要进行取整处理：
   $$y=round(y)$$

3. $DDA$算法就是一个增量算法，即在一个迭代算法中，如果每一步的$x$、$y$值是用前一步的值加上一个增量来获得，则称为增量算法。

4. 这种方法直观，但效率太低，因为每一步需要一次浮点乘法和一次舍入运算。

二、计算

$$\begin{array}{rl}{y}_{i+1}& =k{x}_{i+1}+b\\ & =k(x_i+\Delta x)+b\\ & =k{x}_i+b+\Delta x\\ & =y_i+k\Delta x\end{array}$$

• 当$x$每递增$1$，$y$递增$k$(即直线斜率)，即当$x_{i+1}=x_i+1$时，$y_{i+1}=y_i+k$；
• 取整：$y_{i+1}=round(y_i+k+0.5)$；
• 注意上述分析的算法仅适用于$|k|\le 1$的情形。在这种情况下，$x$每增加$1$，$y$最多增加$1$。
• 当$|k|\ge 1$时，$y$每增加$1$，$x$增加$\frac{1}{k}$，之后同理。

三、举例

写出下列直线段$P_0(0,0)\to P_1(5,2)$所经过实际坐标。

1. 判断$k$值，决定使用谁为增量
   $k=\frac{2-0}{5-0}=0.4<1$，使用$x$为增量，步长为$1$。
2. 计算
   $$\begin{array}{rl}{x}_{i+1}& =x_i+1\\ y_{i+1}& =y_i+0.4\end{array}$$
   即：当$x_{i+1}=x_i+1$时，$y_{i+1}=y_i+0.4$

x	y+0.5	int(y+0.5)
0	0+0.5	0
1	0.4+0.5	0
2	0.8+0.5	1
3	1.2+0.5	1
4	1.6+0.5	2
5	2.0+0.5	2

3. 结果
   经过如下5个坐标：
   $P_0(0,0)、P_1(1,0)、P_2(2,1)、P_3(3,1)、P_4(4,2)、P_5(5,2)$