matlab练习程序（多线段交点）

很简单的算法，这里是把每对线段都进行比较了。

还有一种似乎先通过x和y排序再进行交点判断的，不过那种方法我还没看太明白。

这里的方法如下：

1.根据线段的端点求两条直线的交点。

2.判断直线的交点是否在两条线段上。

结果如下：

matlab代码如下：

clear all;close all;clc;
 
n=;
p=rand(n,);        %（x1,y1,x2,y2）线段两端点
 
for i=:n
    pbar=p(i,:);
    pbar=reshape(pbar,[,]);
    line(pbar(,:),pbar(,:));
end
 
hold on;
for i=:n-
    p1=p(i,:);
    k1=(p1()-p1())/(p1()-p1());
    b1=p1()-k1*p1();
    for j=i+:n
        p2=p(j,:);
        k2=(p2()-p2())/(p2()-p2());
        b2=p2()-k2*p2();
 
        x=-(b1-b2)/(k1-k2);             %求两直线交点
        y=-(-b2*k1+b1*k2)/(k1-k2);
                                        %判断交点是否在两线段上
        if min(p1(),p1())<=x && x<=max(p1(),p1()) && ...
           min(p1(),p1())<=y && y<=max(p1(),p1()) && ...
           min(p2(),p2())<=x && x<=max(p2(),p2()) && ...
           min(p2(),p2())<=y && y<=max(p2(),p2())
            plot(x,y,'.');
        end
    end
end