求直线与平面的交点

求直线与平面的交点今天遇到一个小问题，如何求线与平面的交点，以及如何将其高效实现，下面就从数学的原理以及编程实现来探究一下这个小问题。目录文章目录问题重述与几何模型已知P1(x1,y1,z1),P2(x2,y2,z2),平面α:ax+by+cz+d=0,求直线P1P2与平面α的交点P。数学分析存在性：直线与平面的交点可能有零个，一个，或无数个。可行性：已知直线上不重合两点，可以确定一条直线，已知直线与平面，则一定可以得到两者之间的关系。###常规解法（联立方程）当P1P2不与任一坐标平面平行时：$\dfrac{x-x1}{x1-x2}=\dfrac{y-y1}{y1-y2}=\dfrac{z-z1}{z1-z2}$用xxx表示yyy和zzz:y=y1+(x1−x2)(y1−y2)x−x1y=y1+\dfrac{(x1-x2)(y1-y2)}{x-x1}y=y1+x−x1(x1−x2)(y1−y2)​z=z1+(x1−x2)(z−z1)x−x1z=z1+\dfrac{(x1-x2)(z-z1)}{x-x1}z=z1+x−x1(x1−x2)(z−z1)​与直线方程ax+by+cz+d=0ax+by+cz+d=0ax+by+cz+d=0联立可得：ax(x1−x2)+b((y1−y2)(x−x1)+y1(x1−x2))+c((z1−z2)(x−x1)+z1(x1−x2))+d=0ax(x1-x2)+b((y1-y2)(x-x1)+y1(x1-x2))+c((z1-z2)(x-x1)+z1(x1-x2))+d=0ax(x1−x2)+b((y1−y2)(x−x1)+y1(x1−x2))+c((z1−z2)(x−x1)+z1(x1−x2))+d=0解该方程即可得：x=(−d+bx1(y1−y2)−by1(x1−x2)+cx1(z1−z2)−cz1(x1−x2))a(x1−x2)+b(y1−y2)+c(z1−z2)x=\dfrac{(-d+bx1(y1-y2)-by1(x1-x2)+cx1(z1-z2)-cz1(x1-x2))}{a(x1-x2)+b(y1-y2)+c(z1-z2)}x=a(x1−x2)+b(y1−y2)+c(z1−z2)(−d+bx1(y1−y2)−by1(x1−x2)+cx1(z1−z2)−cz1(x1−x2))​然后即可得到y,zy,zy,z。：该方法原理简单，容易理解，在一般考试解题中使用较多，但是不容易步骤化，编程实现需要很多判断，代码比较繁琐。###向量法当已知平面的一般式方程时(ax+by+cz+d=0ax+by+cz+d=0ax+by+cz+d=0)，n⃗=(a,b,c)′\vecn=(a,b,c)'n=(a,b,c)′就是平面的法矢量，也就能够很容易求出点到平面的距离和一个向量到法矢量的投影。由前面的几何模型，我们容易得出：P1P2→P1P→=P1D2→P1D→=∣P1D2→∣∣P1D→∣=n\dfrac{\overrightarrow{P1P2}}{\overrightarrow{P1P}}=\dfrac{\overrightarrow{P1D2}}{\overrightarrow{P1D}}=\dfrac{\lvert\overrightarrow{P1D2}\rvert}{\lvert\overrightarrow{P1D}\rvert}=nP1PP1P2​=P1DP1D2​=∣P1D∣∣P1D2∣​=n其中：∣P1D→∣=ax1+by1+cz1+da2+b2+c2\lvert\overrightarrow{P1D}\rvert=\dfrac{ax1+by1+cz1+d}{\sqrt{a^{2}+b^{2}+c^{2}}}∣P1D∣=a2+b2+c2​ax1+by1+cz1+d​∣P1D2→∣=a(x2−x1)+b(y2−y1)+c(z2−z1)a2+b2+c2\lvert\overrightarrow{P1D2}\rvert=\dfrac{a(x2-x1)+b(y2-y1)+c(z2-z1)}{\sqrt{a^{2}+b^{2}+c^{2}}}∣P1D2∣=a2+b2+c2​a(x2−x1)+b(y2−y1)+c(z2−z1)​可得：n=∣P1D2→∣∣P1D→∣=a(x2−x1)+b(y2−y1)+c(z2−z1)ax1+by1+cz1+dn=\dfrac{\lvert\overrightarrow{P1D2}\rvert}{\lvert\overrightarrow{P1D}\rvert}=\dfrac{a(x2-x1)+b(y2-y1)+c(z2-z1)}{ax1+by1+cz1+d}n=∣P1D∣∣P1D2∣​=ax1+by1+cz1+da(x2−x1)+b(y2−y1)+c(z2−z1)​易得：OP→=OP1→+P1P→=OP1→+nP1P2→\overrightarrow{OP}=\overrightarrow{OP1}+\overrightarrow{P1P}=\overrightarrow{OP1}+n\overrightarrow{P1P2}OP=OP1+P1P=OP1+nP1P2上面过程主要使用向量运算，巧妙避开了求模运算中的开方运算，减少了很多的运算量，而且，容易步骤化，易于编程实现。##编程实现###常规解法。//c//p[]:XYZforpointP(intersection)//P1[]:XYZforP1//P2[]:XYZforP2//coffi[]:a,b,c,dfortheplaneboolintersectionLinePlane(floatP[],floatP1[],floatP2[],floatcoeffi[]){floatP1P2[3]={0};floattemp_x=P1[0]-P2[0],temp_y=P1[1]-P2[1],temp_z=P1[2]-P2[2];inti=0;for(i=0;i<3;i++){P1P2[i]=P2[i]-P1[i];}//parallel?if(fabs(P1P2[0]*coeffi[0]+P1P2[1]*coeffi[1]+P1P2[2]*coeffi[2])<1e-5){//nointersectionorthelineisontheplanereturnfalse;}//ifitcan'tbesolvebythismethod?if(fabs(temp_x)<1e-5||fabs(temp_y)<1e-5||fabs(temp_z)<1e-5){returnfalse;}*des_x=(-coeffi[3]*temp_x+b*P1[0]*temp_y-coeffi[1]*P1[1]*temp_x+coeffi[2]*p1[0]*temp_z-coeffi[2]*P1[2]*temp_x)/(coeffi[0]*temp_x+coeffi[1]*temp_y+coeffi[2]*temp_z);*des_y=(*des_x-P1[0])*temp_y/temp_x+P1[1];*des_z=(*des_x-P1[0])*temp_z/temp_x+P1[2];returntrue;}。###向量解法//C//p[]:XYZforpointP(intersection)//P1[]:XYZforP1//P2[]:XYZforP2//coffi[]:a,b,c,dfortheplaneboolintersectionLinePlane(floatP[],floatP1[],floatP2[],floatcoeffi[]){floatP1P2[3]={P2[0]-P1[0],P2[1]-P1[1],P2[2]-P1[2]};inti=0;floatnum,den,n;num=coeffi[0]*P1[0]+coeffi[1]*P1[1]+coeffi[2]*P1[2]+coeffi[3];den=coeffi[0]*P1P2[0]+coeffi[1]*P1P2[1]+coeffi[2]*P1P2[2];if(fabs(d