关于直线镜像矩阵生成

最近项目需求，在做关于直线镜像变换，网上看了一些，留作备份，以备追溯，

http://wenku.baidu.com/view/9bbcf9252af90242a895e5bf.html

http://blog.youkuaiyun.com/rabbitguiming/article/details/3963213

以上两个是参考的两个帖子，由于我要用到关于任意直线Ax+By+C=0的镜像，采用了第一个说道的方法，并将其代码转换为矩阵表示，更加直观

对于直线ax+by+c=0，有

      M1=[b*b -a*b -a*c;-a*b a*a -b*c;0 0 0]/(a*a+b*b);
      M2=[-1 0 0;0 -1 0; 0 0 1];
      M=M2*(E-M1)+M1;

其中M就是最后所求关于直线对称的镜像矩阵了