关于直线镜像矩阵生成

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

• 关于任意直线的对称变换
• 矩阵运算-矩阵投影，镜像，切变

以上两个是参考的两个帖子，由于我要用到关于任意直线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就是最后所求关于直线对称的镜像矩阵了