主程序

a=input('第一个球的x坐标');

b=input('第一个球的y坐标');

r1=input('第一个球的半径');

c=input('第二个球的x坐标');

d=input('第二个球的y坐标');

r2=input('第二个球的半径');

w=input('矩形宽');

h=input('矩形高');

R=[0 w w 0;

   0 0 h h];

for I=1:4

A(I)=lineangle(a,b,R(1,I),R(2,I));

%A存放第一个球到矩形四个角的角度

end

Z=mirror(c,d,w,h);%z(2,5)存放了射击球的5个射击点

q=[0];

for I=1:5

    [small,big]=fwei(Z(1,I),Z(2,I),a,b,r1+r2);

    for J=1:4

        if (A(J)>=small&&A(J)<=big)||(A(J)-2*pi>=small&&A(J)-2*pi<=big)

            p=answer(A(J),Z(1,I),Z(2,I),c,d,r1,r2);

            q=[q,finalanswer(p,I)];

        end

    end

end

function angle=lineangle(a,b,c,d)
%以(a,b)基点，连接2点直线的角度
if c-a==0
    angle=pi/2;
end
angle=atan((d-b)/(c-a));
if (d-b)/(c-a)<0
    if d-b>0
        angle=angle+pi;
    else
        angle=angle+2*pi;
    end
else
    if d-b<0
        angle=angle+pi;
    end
end
angle=angle+pi;
if(angle>2*pi) angle=angle-2*pi;
end

function [small,big]=fwei(a,b,c,d,r)

%点(a,b)做圆(c,d),半径r两条切线所夹的圆心角范围

angle=lineangle(c,d,a,b);

long=sqrt((a-c)^2+(b-d)^2);

o=acos(r/long);

small=angle-o;

big=angle+o;

function q=finalanswer(p,I)

%在第I种边映射下的p角还原

if I==1

    q=p;

else if (I==2)|(I==4)

        q=-p+2*pi;

    else if (I==3)|(I==5)

            q=pi-p;

        end

    end

end

function p=answer(o,a,b,c,d,r1,r2)

c=c+(r1+r2)*cos(o);

d=d+(r1+r2)*sin(o);

p=lineangle(a,b,c,d);

function zs=mirror(a,b,w,h)

%点(a,b)关于长w,宽h的矩形四条边的映射

zs(:,1)=[a;b];

X=[a -a a 2*w-a];

Y=[-b b 2*h-b b];

for I=2:5

   zs(:,I)=[X(I-1);Y(I-1)];

end