深圳网站设计 建设科技有限公司,网站建设的常见技术有哪些方面,上饶做网站多少钱,免费下载app应用软件对于你给出的图片#xff0c;我根据之前提到的思路编写了以下程序。感觉还蛮准的。这个程序只能用于两圆相交的情况。当小圆在大圆内部时#xff0c;需要在其中加以判定。仅供参考。clear;clcc imread(1.jpg);cc c;tic[m n] size(c);for i1:mindex find(~c(i,:));if leng…对于你给出的图片我根据之前提到的思路编写了以下程序。感觉还蛮准的。这个程序只能用于两圆相交的情况。当小圆在大圆内部时需要在其中加以判定。仅供参考。clear;clcc imread(1.jpg);cc c;tic[m n] size(c);for i1:mindex find(~c(i,:));if length(index) 2c(i,index(1)1:index(end)-1) true;endend[x y] find(~c);maxdist 0;for i1:length(x)for ji1:length(x)tmpdist (x(i)-x(j))^2(y(i)-y(j))^2;if tmpdist maxdistmaxdist tmpdist;pair [i j];endendendx1 x(pair(1));y1 y(pair(1));x2 x(pair(2));y2 y(pair(2));dist (y2-y1)*(x-x1)-(x2-x1)*(y-y1);index find(dist max(dist));x3 x(index(1));y3y(index(1));index find(dist min(dist));x4 x(index(1));y4y(index(1));dist1 (x3-x1)^2(y3-y1)^2(x4-x1)^2(y4-y1)^2;dist2 (x3-x2)^2(y3-y2)^2(x4-x2)^2(y4-y2)^2;if dist2 dist1tmpx x2; tmpy y2;x2 x1; y2 y1;x1 tmpx; y1 tmpy;endcenterx1 ((y3-y1)*(x4^2y4^2-x1^2-y1^2)-(y4-y1)*(x3^2y3^2-x1^2-y1^2))/((y3-y1)*(x4-x1)-(y4-y1)*(x3-x1))/2;centery1 ((x3-x1)*(x4^2y4^2-x1^2-y1^2)-(x4-x1)*(x3^2y3^2-x1^2-y1^2))/((x3-x1)*(y4-y1)-(x4-x1)*(y3-y1))/2;r1 sqrt((x1-centerx1)^2(y1-centery1)^2);dist (x-centerx1).^2(y-centery1).^2;index find(dist (r13)^2);x(index)[];y(index)[];dist (y2-y1)*(x-x1)-(x2-x1)*(y-y1);index find(dist max(dist));x5 x(index(1));y5y(index(1));index find(dist min(dist));x6 x(index(1));y6y(index(1));centerx2 ((y5-y2)*(x6^2y6^2-x2^2-y2^2)-(y6-y2)*(x5^2y5^2-x2^2-y2^2))/((y5-y2)*(x6-x2)-(y6-y2)*(x5-x2))/2;centery2 ((x5-x2)*(x6^2y6^2-x2^2-y2^2)-(x6-x2)*(x5^2y5^2-x2^2-y2^2))/((x5-x2)*(y6-y2)-(x6-x2)*(y5-y2))/2;r2 sqrt((x2-centerx2)^2(y2-centery2)^2);tocimshow(cc);hold on;plot(centery1,centerx1,r);plot(centery1r1*cos(pi*(0:100)/50),centerx1r1*sin(pi*(0:100)/50),b-)plot(centery2,centerx2,r);plot(centery2r2*cos(pi*(0:100)/50),centerx2r2*sin(pi*(0:100)/50),b-)