matlab练习程序（射线法判断点与多边形关系）

clear all;close all;clc;

polyn=20;
poly=rand(polyn,2);
poly=createSimplyPoly(poly);  %创建简单多边形

ployn=polyn+1;                %连线方便，把第一个点添加到最后一个点后面
poly(polyn,:)=poly(1,:);

pn=500;
p=rand(pn,2);

hold on;
for i=1:pn
   
    flag=0;
    for j=2:polyn
        x1=poly(j-1,1);         %多边形前后两个点
        y1=poly(j-1,2);
        x2=poly(j,1);
        y2=poly(j,2);
       
        k=(y1-y2)/(x1-x2);      %多边形一条边直线
        b=y1-k*x1;
        x=p(i,1);               %过当前点直线和多边形交点
        y=k*x+b;          
               
        if min([x1 x2])<=x && x<=max([x1 x2]) && ...        %点同时在射线和多边形边上
           min([y1 y2])<=y && y<=max([y1 y2]) &&  y>=p(i,2)
               flag=flag+1;
        end
    end
   
    if mod(flag,2)==0               %偶数则在外部
        plot(p(i,1),p(i,2),‘r.‘);
    else                            %奇数则在内部
        plot(p(i,1),p(i,2),‘g.‘);        
    end
end

plot(poly(:,1),poly(:,2));