余弦相似度在很多地方都有介绍，这里只是用IML简单的验证。
数据引自：http://www.ruanyifeng.com/blog/2013/03/cosine_similarity.html
关键步骤：

对应代码：
note：这样只是为了和公式更好对应。事实上可以用a×b(转置)的方法写，更加简洁。

proc iml;

a={1 2 2 1 1 1 0};

b={1 2 2 1 1 2 1};

ab=a#b;

cos=ab[+]/(sqrt(a*a`)*sqrt(b*b`));

aa=a*a`;

bb=b*b`;

print cos aa bb;

quit;