Пояснить реализацию алгоритма отсечения

n = 100;
Xmin = -4; Xmax = 4;
Ymin = -3; Ymax = 3;
 
figure; hold on;
 
P1x = rand(1,n)*20 - 10;  
P2x = rand(1,n)*20 - 10;
P1y = rand(1,n)*20 - 10;  
P2y = rand(1,n)*20 - 10;
P1code = zeros(n,4);
P2code = zeros(n,4);
 
for i = 1:n    
 
    if P1x(i) < Xmin
        P1code(i,1) = 1;
    end
    if P1x(i) > Xmax
        P1code(i,2) = 1;
    end
    if P1y(i) < Ymin
        P1code(i,3) = 1;
    end
    if P1y(i) > Ymax
        P1code(i,4) = 1;
    end    
 
    if P2x(i) < Xmin
        P2code(i,1) = 1;
    end
    if P2x(i) > Xmax
        P2code(i,2) = 1;
    end
    if P2y(i) < Ymin
        P2code(i,3) = 1;
    end
    if P2y(i) > Ymax
        P2code(i,4) = 1;
    end
 
    plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'b-');
end
hold off;
 
P_label = zeros(1,n);
figure; hold on;
for i = 1:n
    P_or = P1code(i,1:4) | P2code(i,1:4);
    if sum(P_or) == 0
        P_label(i) = 1;
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'r-');
    end
    P_and = P1code(i,1:4) & P2code(i,1:4);
    if sum(P_and) > 0
        P_label(i) = 2;
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'g-');
    end
end
hold off;
 
figure; hold on;
for i = 1:n
    if P_label(i) == 0
 
        P_or = P1code(i,1:4) | P2code(i,1:4);
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'g-');
 
        if P_or(1) == 1
            Py = P1y(i) + (Xmin-P1x(i))*(P2y(i)-P1y(i))/(P2x(i)-P1x(i));
            if P1x(i) < Xmin
                P1x(i) = Xmin; P1y(i) = Py;
            elseif P2x(i) < Xmin
                P2x(i) = Xmin; P2y(i) = Py;
            end            
        end
 
        if P_or(2) == 1
            Py = P1y(i) + (Xmax-P1x(i))*(P2y(i)-P1y(i))/(P2x(i)-P1x(i));
            if P1x(i) > Xmax
                P1x(i) = Xmax; P1y(i) = Py;
            elseif P2x(i) > Xmax
                P2x(i) = Xmax; P2y(i) = Py;
            end            
        end
 
        if P_or(3) == 1
            Px = P1x(i) + (Ymin-P1y(i))/(P2y(i)-P1y(i))*(P2x(i)-P1x(i));
            if P1y(i) < Ymin
                P1x(i) = Px; P1y(i) = Ymin;
            elseif P2y(i) < Ymin
                P2x(i) = Px; P2y(i) = Ymin;
            end            
        end
 
        if P_or(4) == 1
            Px = P1x(i) + (Ymax-P1y(i))/(P2y(i)-P1y(i))*(P2x(i)-P1x(i));
            if P1y(i) > Ymax
                P1x(i) = Px; P1y(i) = Ymax;
            elseif P2y(i) > Ymax
                P2x(i) = Px; P2y(i) = Ymax;
            end            
        end
        if P1x(i) >= Xmin & P1x(i) <= Xmax & ...
                P2x(i) >= Xmin & P2x(i) <= Xmax & ...
                P1y(i) >= Ymin & P1y(i) <= Ymax & ...
                P2y(i) >= Ymin & P2y(i) <= Ymax
            plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'m-','LineWidth',2);
        end
        plot([-4 4 4 -4 -4],[-3 -3 3 3 -3],'b-','LineWidth',2);
    end
end
hold off;

clc,clear
 
n = 20;                    % количество отрезков
    % границы области
% t = linspace(0,2*pi,5)+pi/4;
% xv = 8*cos(t);
% yv = 6*sin(t);
xv = [-4 -4 4 4 -4];
yv = [-3 3 3 -3 -3];
 
    % генерируем точки отрезков в границах -10 : 10
a = 20;
x1 = rand(1,n)*a - a/2;  
x2 = rand(1,n)*a - a/2;
y1 = rand(1,n)*a - a/2;  
y2 = rand(1,n)*a - a/2;
 
    % точки пересечения отрезков и границ области
CR = false(1,n);
[xc1,xc2,yc1,yc2] = deal( zeros(1,n) );
for i = 1:n
        % отрезок
    x = [x1(i) x2(i)];
    y = [y1(i) y2(i)];
    [xc,yc] = polyxpoly(x,y,xv,yv);
    if ~isempty(xc)            % если точки пересечения есть
        CR(i) = true;
        xc1(i) = xc(1);        % первая точка пересечения
        yc1(i) = yc(1);
        if length(xc) == 2     % если точек - две
            xc2(i) = xc(2);
            yc2(i) = yc(2);
        else                   % если точка одна, определяем 
                               % какая из точек отрезка 
                               % находится внутри области
            if inpolygon(x1(i),y1(i),xv,yv)
                xc2(i) = x1(i);
                yc2(i) = y1(i);
            else
                xc2(i) = x2(i);
                yc2(i) = y2(i);
            end
        end
    end
end
 
in1 = inpolygon(x1,y1,xv,yv);
in2 = inpolygon(x2,y2,xv,yv);
IN  = in1 & in2;                % отрезки находяющиемя внутри
OUT = ~(IN | CR);            % отрезки находяющиемя снаружи
 
plot(xv, yv, 'b-','LineWidth',2);                   % границы области
hold on
plot([x1; x2], [y1; y2],'g')                        % все отрезки
plot([xc1(CR); xc2(CR)], [yc1(CR); yc2(CR)],'r')    % пересекающие
plot([x1(IN);   x2(IN)], [y1(IN);   y2(IN)],'k')    % внутри
plot([x1(OUT); x2(OUT)], [y1(OUT); y2(OUT)],'b')    % снаружи
hold off

clc,clear
 
n = 50;                    % количество отрезков
    % границы области
% t = linspace(0,2*pi,5)+pi/4;
% xv = 8*cos(t);
% yv = 6*sin(t);
xv = [-4 -4 4 4 -4]';
yv = [-3 3 3 -3 -3]';
 
    % генерируем точки отрезков в границах -10 : 10
a = 20;
x1 = rand(n,1)*a - a/2;  
x2 = rand(n,1)*a - a/2;
y1 = rand(n,1)*a - a/2;  
y2 = rand(n,1)*a - a/2;
 
    % точки пересечения отрезков и границ области
CR = false(n,1);
[xc1,xc2,yc1,yc2] = deal( zeros(n,1) );
for i = 1:n
        % отрезок
    x = [x1(i) x2(i)];
    y = [y1(i) y2(i)];
    [xc,yc] = intersections(x,y,xv,yv);
    if ~isempty(xc)            % если точки пересечения есть
        CR(i) = true;
        xc1(i) = xc(1);        % первая точка пересечения
        yc1(i) = yc(1);
        if length(xc) == 2     % если точек - две
            xc2(i) = xc(2);
            yc2(i) = yc(2);
        else                   % если точка одна, определяем 
                               % какая из точек отрезка 
                               % находится внутри области
            if inpoly2([x1(i) y1(i)],[xv yv])
                xc2(i) = x1(i);
                yc2(i) = y1(i);
            else
                xc2(i) = x2(i);
                yc2(i) = y2(i);
            end
        end
    end
end
 
in1 = inpoly2([x1 y1],[xv yv]);
in2 = inpoly2([x2 y2],[xv yv]);
IN  = in1 & in2;                % отрезки находяющиемя внутри
OUT = ~(IN | CR);            % отрезки находяющиемя снаружи
 
plot(xv, yv, 'b-','LineWidth',2);                   % границы области
hold on
plot([x1 x2]', [y1 y2]','g')                        % все отрезки
plot([xc1(CR) xc2(CR)]', [yc1(CR) yc2(CR)]','r')    % пересекающие
plot([x1(IN)   x2(IN)]', [y1(IN)   y2(IN)]','k')    % внутри
plot([x1(OUT) x2(OUT)]', [y1(OUT) y2(OUT)]','b')    % снаружи
hold off