anglec.m

在matlab环境下

function y = angle(a,b)
        aa = sqrt(a'*a);
        bb = sqrt(b'*b);
        ab = a'*b;
        cab = ab/(aa*bb);
        y1 = acos(cab)*180/pi;
%-------------------------------------------------------------        
        yrad = y1*pi/180;
         m = [cos(yrad) -sin(yrad); sin(yrad) cos(yrad)];
         nb = m*b;
         nnb = sqrt(nb'*nb);
         nip = a'*nb;
         ncos = nip/(aa*nnb);
         y2 = acos(ncos)*180/pi;
         if y2 == 0 & y1 ~= 180
             y = -y1+360;
         else    
             y = y1;
         end;