（a）、在二维坐标系中

已知坐标系中任意四个点，计算任意两点连线的夹角

（x1,y1）（x2,y2）

（x3,y3）（x4,y4）

theta=acosd(dot([x1-x2,y1-y2],[x3-x4,y3-y4])/(norm([x1-x2,y1-y2])*norm([x3-x4,y3-y4])))

（b）、在三维坐标系中

（x1,y1,z1）（x2,y2,z2）

（x3,y3,z3）（x4,y4,z4）

theta=acosd(dot([x1-x2,y1-y2,z2-z1],[x3-x4,y3-y4,z4-z3])/(norm([x1-x2,y1-y2,z2-z1])*norm([x3-x4,y3-y4,z4-z3])))

*解释：

http://s8/mw690/002O3Vnlzy6IiXGlF0bc7&690

程序中，dot为点乘的；norm为计算模值，具体可参考help中的解释