【matlab】：matlab如何实现计算三个点组成的角的角度值

如题：给出三个点1,2,3，如何计算以点2为角心的角的角度值

方法：
1，利用余弦定理，求出某个角的余弦值
2，利用acos求出该余弦值对应的弧度值
3，利用弧度值与角度值的对应关系求出角度值

图片：

代码如下：
各点的坐标是我随机出来的

x1 = 1;y1 = 1;
x2 = 1;y2 =90;
x3 = 90;y3 = 90;
a2 = (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
b2 = (x3-x2)*(x3-x2)+(y3-y2)*(y3-y2);
c2 = (x1-x3)*(x1-x3)+(y1-y3)*(y1-y3);
a = sqrt(a2);
b = sqrt(b2);
c = sqrt(c2);
pos = (a2+b2-c2)/(2*a*b);    //求出余弦值
angle = acos(pos);         //余弦值装换为弧度值
realangle = angle*180/pi;   //弧度值转换为角度值
disp(realangle);