Matlab中多条直线求交点

clear;clc;
x = [-10:1:15];
cx1 = 2;
cy1 = -5;
r1 = 3;

cx2 = 1;
cy2 = 10;
r2 = -5;

cx3 = 3;
cy3 = 14;
r3 = 10;

y1 = (r1 - cx1.*x)/cy1;
y2 = (r2 - cx2.*x)/cy2;
y3 = (r3 - cx3.*x)/cy3;

p1 = [cx1 cy1; cx2 cy2];
b1 = [r1 r2]’;
xy1 = p1\b1;

p2 = [cx1 cy1; cx3 cy3];
b2 = [r1 r3]’;
xy2 = p2\b2;

p3 = [cx2 cy2; cx3 cy3];
b3 = [r2 r3]’;
xy3 = p3\b3;
figure(1);
h=plot(x, y1, ‘r’,x,y2, ‘k’,x, y3, ‘g’,xy1(1),xy1(2), ‘og’,xy2(1),xy2(2), ‘og’,xy3(1),xy3(2), ‘og’);
legend([h(1),h(2),h(3),h(4),h(5),h(6)],‘eq 7’,‘eq 8’, ‘eq 9’,‘solution 1’,‘solution 2’, ‘solution 3’);
t = title(‘Linear Equations and Their Solution’)
ylabel(‘y’)
xlabel(‘x’)
运行结果

重点是求三个交点，分三次计算交点，然后做到一个图中
还有图例显示多个，用句柄方式处理