[½»Á÷] ¡¾ÇóÖú¡¿¼±£¡£¡£¡£¡£¡¸ßÈËÇë´Í½Ì

Originally posted by sobigstar11 at 2010-10-13 08:47:15:
ÎÒÏëÇë½ÌÄúһϠÁ½¸öÔ²ÒÑÖªÔ²ÐÄ×ø±êºÍ°ë¾¶Çó½»µã×ø±ê£¬MATLABÔõôʵÏÖ£¬ÎÒÓõÄSOLVEº¯Êý£¬¿ÉÊǽá¹ûµÃµ½µÄÊǸ´ÊýÖµ£¬²»ÖªµÀÔõô»ØÊ£¬ÁíÍ⸴ÊýÄÜÔÚ¶þά×ø±êϵÖÐ×÷ͼÂð

function IntersectPoints( Cx1, Cy1, r1, Cx2, Cy2, r2)
% ÇóÁ½¸öÔ²µÄ½»µã

clc

figure('name','Intersect point(s) of two circle');
theta = 0 : pi / 1000: 2 * pi;
x1 = r1 * cos( theta ) + Cx1;
y1 = r1 * sin( theta ) + Cy1;
plot( x1, y1, 'k');
hold on

x2 = r2 * cos( theta ) + Cx2;
y2 = r2 * sin( theta ) + Cy2;
plot( x2, y2, 'k');
hold off

% ÇóÁ½Ô²µÄÖÐÐľà
a = sqrt( ( Cx2 - Cx1 ) ^ 2 + ( Cy2 - Cy1 ) ^ 2 );

if  a > ( r1 + r2 )
   
    disp( ' There are no intersection points ! ' );
   
elseif a == ( r1 + r2 )
   
    disp( ' There is one intersection point !' );
   
else
   
    A = solve('( x1 - Cx1 ) ^ 2 + (y1 - Cy1 ) ^ 2 - r1 ^ 2','( x1 - Cx2 ) ^ 2 + ( y1 - Cy2 ) ^ 2 - r2 ^ 2');
    Px = subs( A.x1 );
    Py = subs( A.y1);
    disp( ' There are two intersection points ! ' );
    hold on
    plot( Px, Py, 'ro' )
    axis equal
    hold off

end