空间解析几何之直线与平面:推导直线与直线、直线与平面交点

2024-03-15 21:36

本文主要是介绍空间解析几何之直线与平面:推导直线与直线、直线与平面交点,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

空间解析几何——直线与平面

三维空间中的直线和平面与二维空间中的性质有一定的类似之处,但是其相交关系的求解方式有所差异。本文回顾了三维空间中直线和平面的解析表达,然后推导线-线线-面交点。

平面

空间平面的表达式为:
A x + B y + C z + D = 0 (1) Ax+By+Cz+D=0\tag{1} Ax+By+Cz+D=0(1)
包含了4个参数 A , B , C , D A,B,C,D A,B,C,D ( A , B , C ) (A,B,C) (A,B,C)是平面的一个法向量。但是它们并非独立的,即法向量的长度可以是任意的。
若限定 A , B , C A,B,C A,B,C三个参数满足:
A 2 + B 2 + C 2 = 1 (2) A^2+B^2+C^2=1\tag{2} A2+B2+C2=1(2)
此时 ( A , B , C ) (A,B,C) (A,B,C)是平面的单位法向量, D D D表示坐标原点到该平面的距离。方便起见,下面的讨论默认平面的A,B,C参数满足(2)式。

本质上,空间平面仅需要三个参数确定: ( n x , n y , 1 ) (n_x,n_y,1) (nx,ny,1)描述其法向量(仅需2参数), d d d描述其到原点的距离,表达式为 n x x + n y y + z + d = 0 n_x x+n_y y+z+d=0 nxx+nyy+z+d=0,与(1)是等价的。

直线

空间直线的一种表达是:两个平面的交。用 π 1 , π 2 \pi_1,\pi_2 π1,π2表示两个不平行、不共面的平面,则一条空间直线可以表达为
l = π 1 ∩ π 2 → l = { A 1 x + B 1 y + C 1 z + D 1 = 0 A 2 x + B 2 y + C 2 z + D 2 = 0 l=\pi_1\cap\pi_2\rightarrow l=\begin{cases}A_1 x+B_1 y+C_1 z+D_1=0\\A_2 x+B_2 y+C_2 z+D_2=0\end{cases} l=π1π2l={A1x+B1y+C1z+D1=0A2x+B2y+C2z+D2=0

显然,这种表达方式不是唯一的。一条空间直线可能是无数对平面的交线。
空间直线更常用的表达式为:
x − x 0 a = y − y 0 b = z − z 0 c \dfrac{x-x_0}{a}=\dfrac{y-y_0}{b}=\dfrac{z-z_0}{c} axx0=byy0=czz0

包含6个参数 x 0 , y 0 , z 0 , a , b , c x_0,y_0,z_0, a,b,c x0,y0,z0,a,b,c ( a , b , c ) (a,b,c) (a,b,c)是直线的方向向量,其长度任意,因而相当于只有两个独立参数; ( x 0 , y 0 , z 0 ) (x_0,y_0,z_0) (x0,y0,z0)是直线上一点。因此,一条空间直线只需要5个独立参数即可描述。

点到平面的距离

( x 1 , y 1 , z 1 ) (x_1,y_1,z_1) (x1,y1,z1)到平面 A x + B y + C z + D = 0 Ax+By+Cz+D=0 Ax+By+Cz+D=0的距离为:

d = ∣ A x 1 + B y 1 + C z 1 + D ∣ A 2 + B 2 + C 2 d=\dfrac{|Ax_1+By_1+Cz_1+D|}{\sqrt{A^2+B^2+C^2}} d=A2+B2+C2 Ax1+By1+Cz1+D

点到直线的距离

( x 1 , y 1 , z 1 ) (x_1,y_1,z_1) (x1,y1,z1)到直线 x − x 0 a = y − y 0 b = z − z 0 c \dfrac{x-x_0}{a}=\dfrac{y-y_0}{b}=\dfrac{z-z_0}{c} axx0=byy0=czz0的距离为:

d = ∣ ( x 1 − x 0 , y 1 − y 0 , z 1 − z 0 ) × ( a , b , c ) ∣ a 2 + b 2 + c 2 d=\dfrac{|(x_1-x_0,y_1-y_0,z_1-z_0)\bold{\times} (a,b,c)|}{\sqrt{a^2+b^2+c^2}} d=a2+b2+c2 (x1x0,y1y0,z1z0)×(a,b,c)

表达式中为何出现叉积?这是因为, ∣ a × b ∣ = ∣ a ∣ ∣ b ∣ ∣ sin ⁡ θ ∣ |a\times b |=|a||b||\sin\theta| a×b=a∣∣b∣∣sinθ,而 d d d恰恰就等于 ∣ ( x 1 − x 0 , y 1 − y 0 , z 1 − z 0 ) ∣ ∗ ∣ sin ⁡ θ ∣ |(x_1-x_0,y_1-y_0,z_1-z_0)|*|\sin\theta| (x1x0,y1y0,z1z0)sinθ

在这里插入图片描述

直线之间的距离

两条直线 x − x 0 a = y − y 0 b = z − z 0 c \dfrac{x-x_0}{a}=\dfrac{y-y_0}{b}=\dfrac{z-z_0}{c} axx0=byy0=czz0 x − x 1 a ′ = y − y 1 b ′ = z − z 1 c ′ \dfrac{x-x_1}{a'}=\dfrac{y-y_1}{b'}=\dfrac{z-z_1}{c'} axx1=byy1=czz1,分别通过点P1 ( x 0 , y 0 , z 0 ) (x_0,y_0,z_0) (x0,y0,z0)和P2 ( x 1 , y 1 , z 1 ) (x_1,y_1,z_1) (x1,y1,z1)
两条直线方向向量的叉积即为待求距离所在的方向,我们记作 N = ( a , b , c ) × ( a ′ , b ′ , c ′ ) N=(a,b,c)\times(a',b',c') N=(a,b,c)×(a,b,c)
则它们之间的距离为:

d = ∣ ( x 1 − x 0 , y 1 − y 0 , z 1 − z 0 ) ⋅ N ∣ ∣ ∣ N ∣ ∣ d = \dfrac{|(x_1-x_0,y_1-y_0,z_1-z_0)\cdot N|}{||N||} d=∣∣N∣∣(x1x0,y1y0,z1z0)N
其中, ∣ ∣ ∗ ∣ ∣ ||*|| ∣∣∣∣表示取模。
在这里插入图片描述

直线与直线相交

两条直线 x − x 0 a = y − y 0 b = z − z 0 c \dfrac{x-x_0}{a}=\dfrac{y-y_0}{b}=\dfrac{z-z_0}{c} axx0=byy0=czz0 x − x 1 a ′ = y − y 1 b ′ = z − z 1 c ′ \dfrac{x-x_1}{a'}=\dfrac{y-y_1}{b'}=\dfrac{z-z_1}{c'} axx1=byy1=czz1,分别通过点P1 ( x 0 , y 0 , z 0 ) (x_0,y_0,z_0) (x0,y0,z0)和P2 ( x 1 , y 1 , z 1 ) (x_1,y_1,z_1) (x1,y1,z1)
若两条空间直线相交,则它们必共面。与求解直线的距离类似,为了求两条空间直线的交点,首先计算两条直线方向向量的叉积 N N N

解法一

由于两条直线的方向向量已知,所以求交点 C C C的问题,可以转化为求解CP1长度的问题。将CP1长度记为L,如图所示,显然有:

L = d / ∣ cos ⁡ θ ∣ L=d/|\cos\theta| L=d/∣cosθ

d d d即P1点到另一直线的距离。即有:
d = ∣ ( x 1 − x 0 , y 1 − y 0 , z 1 − z 0 ) × ( a ′ , b ′ , c ′ ) ∣ a ′ 2 + b ′ 2 + c ′ 2 d=\dfrac{|(x_1-x_0,y_1-y_0,z_1-z_0)\bold{\times} (a',b',c')|}{\sqrt{a'^2+b'^2+c'^2}} d=a′2+b′2+c′2 (x1x0,y1y0,z1z0)×(a,b,c)

两直线夹角的余弦为:
∣ cos ⁡ θ ∣ = a a ′ + b b ′ + c c ′ a 2 + b 2 + c 2 a ′ 2 + b ′ 2 + c ′ 2 |\cos\theta|=\dfrac{aa'+bb'+cc'}{\sqrt{a^2+b^2+c^2}\sqrt{a'^2+b'^2+c'^2}} cosθ=a2+b2+c2 a′2+b′2+c′2 aa+bb+cc
在这里插入图片描述
那么,由于距离d没有符号,C点坐标为:
C = P 1 + L ⋅ ( a , b , c ) 或 C = P 1 − L ⋅ ( a , b , c ) C=P1+L\cdot (a,b,c)\text{或}C=P1-L\cdot (a,b,c) C=P1+L(a,b,c)C=P1L(a,b,c)
为确定C坐标为上面两个结果中的哪一个,将它们代入直线 l 2 l_2 l2的表达式,误差较小的那个就是真正的交点。

解法二

设未知数t,由于C在直线 l 1 l_1 l1上,其坐标可以表达为 ( x 0 + t a , y 0 + t b , z 0 + t c ) (x_0+ta,y_0+tb,z_0+tc) (x0+ta,y0+tb,z0+tc);代入直线 l 2 l_2 l2的表达式,有:

{ x 0 + t a − x 1 a ′ = y 0 + t b − y 1 b ′ x 0 + t a − x 1 a ′ = z 0 + t c − z 1 c ′ \begin{cases}\dfrac{x_0+ta-x_1}{a'}=\dfrac{y_0+tb-y_1}{b'}\\\dfrac{x_0+ta-x_1}{a'}=\dfrac{z_0+tc-z_1}{c'}\\\end{cases} ax0+tax1=by0+tby1ax0+tax1=cz0+tcz1

由上面两式,可分别求出两个t值:

t 1 = ( x 1 − x 0 ) b ′ − ( y 1 − y 0 ) a ′ a b ′ − a ′ b t 2 = ( x 1 − x 0 ) c ′ − ( z 1 − z 0 ) a ′ a c ′ − a ′ c t_1=\dfrac{(x_1-x_0)b'-(y_1-y_0)a'}{ab'-a'b}\\ t_2=\dfrac{(x_1-x_0)c'-(z_1-z_0)a'}{ac'-a'c}\\ t1=abab(x1x0)b(y1y0)at2=acac(x1x0)c(z1z0)a
如果 t 1 , t 2 t_1,t_2 t1,t2非常接近,则说明两直线确实交于一点;如果它们差异较大,说明两直线可能不相交(即不共面或平行)。

多条直线相交时,应按照最小二乘法求交点。

直线与平面相交

已知直线 x − x 0 a = y − y 0 b = z − z 0 c \dfrac{x-x_0}{a}=\dfrac{y-y_0}{b}=\dfrac{z-z_0}{c} axx0=byy0=czz0,直线上一点P ( x 0 , y 0 , z 0 ) (x_0,y_0,z_0) (x0,y0,z0),平面方程为: A x + B y + C z + D = 0 Ax+By+Cz+D=0 Ax+By+Cz+D=0,如图所示求交点坐标:

在这里插入图片描述
设未知数t,由于C在直线上,其坐标可以表达为 ( x 0 + t a , y 0 + t b , z 0 + t c ) (x_0+ta,y_0+tb,z_0+tc) (x0+ta,y0+tb,z0+tc)

因为C在平面上,所以满足
A ( x 0 + t a ) + B ( y 0 + t b ) + C ( z 0 + t c ) + D = 0 A(x_0+ta)+B(y_0+tb)+C(z_0+tc)+D=0 A(x0+ta)+B(y0+tb)+C(z0+tc)+D=0

整理可得
t = − A x 0 + B y 0 + C z 0 + D a A + b B + c C t=-\dfrac{Ax_0+By_0+Cz_0+D}{aA+bB+cC} t=aA+bB+cCAx0+By0+Cz0+D

那么,直线与平面的交点坐标为:
C = P + t ∗ ( a , b , c ) C=P+t*(a,b,c) C=P+t(a,b,c)

这篇关于空间解析几何之直线与平面:推导直线与直线、直线与平面交点的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/813334

相关文章

MySQL8.0临时表空间的使用及解读

《MySQL8.0临时表空间的使用及解读》MySQL8.0+引入会话级(temp_N.ibt)和全局(ibtmp1)InnoDB临时表空间,用于存储临时数据及事务日志,自动创建与回收,重启释放,管理高... 目录一、核心概念:为什么需要“临时表空间”?二、InnoDB 临时表空间的两种类型1. 会话级临时表

Three.js构建一个 3D 商品展示空间完整实战项目

《Three.js构建一个3D商品展示空间完整实战项目》Three.js是一个强大的JavaScript库,专用于在Web浏览器中创建3D图形,:本文主要介绍Three.js构建一个3D商品展... 目录引言项目核心技术1. 项目架构与资源组织2. 多模型切换、交互热点绑定3. 移动端适配与帧率优化4. 可

Python进阶之列表推导式的10个核心技巧

《Python进阶之列表推导式的10个核心技巧》在Python编程中,列表推导式(ListComprehension)是提升代码效率的瑞士军刀,本文将通过真实场景案例,揭示列表推导式的进阶用法,希望对... 目录一、基础语法重构:理解推导式的底层逻辑二、嵌套循环:破解多维数据处理难题三、条件表达式:实现分支

Ubuntu如何分配​​未使用的空间

《Ubuntu如何分配​​未使用的空间》Ubuntu磁盘空间不足,实际未分配空间8.2G因LVM卷组名称格式差异(双破折号误写)导致无法扩展,确认正确卷组名后,使用lvextend和resize2fs... 目录1:原因2:操作3:报错5:解决问题:确认卷组名称​6:再次操作7:验证扩展是否成功8:问题已解

MySQL之InnoDB存储页的独立表空间解读

《MySQL之InnoDB存储页的独立表空间解读》:本文主要介绍MySQL之InnoDB存储页的独立表空间,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、背景2、独立表空间【1】表空间大小【2】区【3】组【4】段【5】区的类型【6】XDES Entry区结构【

CSS实现元素撑满剩余空间的五种方法

《CSS实现元素撑满剩余空间的五种方法》在日常开发中,我们经常需要让某个元素占据容器的剩余空间,本文将介绍5种不同的方法来实现这个需求,并分析各种方法的优缺点,感兴趣的朋友一起看看吧... css实现元素撑满剩余空间的5种方法 在日常开发中,我们经常需要让某个元素占据容器的剩余空间。这是一个常见的布局需求

MySQL启动报错:InnoDB表空间丢失问题及解决方法

《MySQL启动报错:InnoDB表空间丢失问题及解决方法》在启动MySQL时,遇到了InnoDB:Tablespace5975wasnotfound,该错误表明MySQL在启动过程中无法找到指定的s... 目录mysql 启动报错:InnoDB 表空间丢失问题及解决方法错误分析解决方案1. 启用 inno

在Java中基于Geotools对PostGIS数据库的空间查询实践教程

《在Java中基于Geotools对PostGIS数据库的空间查询实践教程》本文将深入探讨这一实践,从连接配置到复杂空间查询操作,包括点查询、区域范围查询以及空间关系判断等,全方位展示如何在Java环... 目录前言一、相关技术背景介绍1、评价对象AOI2、数据处理流程二、对AOI空间范围查询实践1、空间查

MySQL表空间结构详解表空间到段页操作

《MySQL表空间结构详解表空间到段页操作》在MySQL架构和存储引擎专题中介绍了使用不同存储引擎创建表时生成的表空间数据文件,在本章节主要介绍使用InnoDB存储引擎创建表时生成的表空间数据文件,对... 目录️‍一、什么是表空间结构1.1 表空间与表空间文件的关系是什么?️‍二、用户数据在表空间中是怎么

CentOS7增加Swap空间的两种方法

《CentOS7增加Swap空间的两种方法》当服务器物理内存不足时,增加Swap空间可以作为虚拟内存使用,帮助系统处理内存压力,本文给大家介绍了CentOS7增加Swap空间的两种方法:创建新的Swa... 目录在Centos 7上增加Swap空间的方法方法一:创建新的Swap文件(推荐)方法二:调整Sww