matlab实现——DDA画线和Bresenham画线

DDA算法：

步骤：

1、求斜率K

2、x坐标增加

3、y坐标增加K，并进行＋0.5取整

4、重复步骤2、3直到绘制完

例如：

Bresenham算法：

步骤：

1、计算Δx   Δy 

2、e=-Δx  x0=x0   y=y0

3、绘制点（x,y）

4、更新e=e+2*Δy,且

1、e>0 时：
    (x,y)变为（x+1，y+1）
    e变为e-2*Δx
2、e<=0时：
    （x,y）变为（x+1,y）

5、重复步骤3、4直到绘制完

MATLAB代码：

DDA算法：

function DDA(x,y)
    L=max(abs(x(1)-y(1)),abs(x(2)-y(2)));
    k=(x(2)-y(2))/(x(1)-y(1));
    X=[];Y=[];
    for i=1:L+1  
        X(i)=x(1);
        Y(i)=x(2);
        x(1)=x(1)+1;
        x(2)=x(2)+k;     
    end
    Y=fix(Y+0.5);
    figure,
    plot(X,Y);
end

 Bresenham算法：

function Bresenham(x,y)
    detax=abs(x(1)-y(1));
    detay=abs(x(2)-y(2));
    e=-detax;
    X=[];Y=[];
    L=max(abs(x(1)-y(1)),abs(x(2)-y(2)));
    for i=1:L+1  
        X(i)=x(1);
        Y(i)=x(2);
        e=e+2*detay;
        if(e>0)
            x(1)=x(1)+1;
            x(2)=x(2)+1;
            e=e-2*detax;
        else
            x(1)=x(1)+1;
            x(2)=x(2); 
        end
    end
    figure,
    plot(X,Y);
end

演示：

clc;
clear;
close all;
x=zeros(1,2);
string=['Enter starting point Coordinates:'];
x=input(string);
y=zeros(1,2);
string=['Enter the coordinates of the end point:'];
y=input(string);
Bresenham(x,y);
DDA(x,y);
Midpoint(x,y);

---

（0,0）——>(5,3)

感觉这什么东西，哪里是直线

---

（0,0）——>(100,40)

有点直线那意思了

 

---

（0,0）——>(1000,1000)

已经看不出来是弯的

如果有错误欢迎指正。