Here is one approach:
% generate a rectangular boundry
rectX = [0.1 0.1 0.9 0.9];
rectY = [0.1 0.5 0.5 0.1];
rect = polyshape(rectX,rectY);
% initialize a figure for display purposes
figure(1);clf
titles = [...
"All Lines";...
"Partially Within";...
"Exclusively Within"];
for i = 1:3
subplot(1,3,i)
hold on
title(titles(i))
plot(rect)
xlim([0 1])
ylim([0 1])
end
% define the number of line segments to generate
nLines = 25;
% initialize a structure to store the lines (if you need the data later)
lines = struct();
% generate each line, calculate it's intersection conditions, and plot
for i=1:nLines
% generate points for a random line segment
lines(i).points = rand(2);
% section the line into parts at intersection with the rectangle
[in, out] = intersect(rect, lines(i).points);
% plot all the lines that were generated
subplot(1,3,1)
line(lines(i).points(:,1),lines(i).points(:,2))
% plot lines that are partially within the boundary
if isempty(in)==0
subplot(1,3,2)
line(lines(i).points(:,1),lines(i).points(:,2))
end
% plot lines that are exclusively within the boundary
if isempty(out)==1
subplot(1,3,3)
line(lines(i).points(:,1),lines(i).points(:,2))
end
end
Here is the output:
