OCC中的参数方程简介

OCC模块概览

参数几何（曲线与曲面类型）

参数

直线

圆

球

Bezier曲线

在Draw Test Harness中绘制Bezier曲线命令：

beziercurve bc 3 0 1 0 4 1 -2 0 5 2 3 0 6
fit

圆的参数方程与隐式方程式的转换

在解析几何中，圆的参数方程和隐式方程之间的转换可以通过代数操作实现。以下我们将详细讲解如何将圆的参数方程转换为隐式方程，以及如何将隐式方程转换回参数方程。

---

1. 圆的参数方程

圆的参数方程通常表示为：

$$\{\begin{array}{l}x=h+r\cos \theta \\ y=k+r\sin \theta \end{array}$$

其中：

• $(h,k)$ 是圆心坐标，
• $r$ 是圆的半径，
• $\theta$ 是参数，取值范围通常是 $0\le \theta <2\pi$。

在原点处的圆，即：
$$\begin{gathered} h=0 \\ k=0 \end{gathered}$$
则：
$$\{\begin{array}{l}x=r\cos \theta \\ y=r\sin \theta \end{array}$$

理解：
参数方程使用一个变量（如 $\theta$）来描述x和y的变化。当 $\theta$ 从0变化到 $2\pi$ 时，点 $(x,y)$ 沿着圆移动一周。

---

2. 将参数方程转换为隐式方程

要将参数方程转换为隐式方程（即标准形式的圆的方程），我们可以利用三角恒等式 ${\cos }^2\theta +{\sin }^2\theta =1$。具体步骤如下：

步骤 1：分离参数

从参数方程中解出 $\cos \theta$ 和 $\sin \theta$：
$$\begin{gathered} \cos \theta =\frac{x-h}{r} \\ \sin \theta =\frac{y-k}{r} \end{gathered}$$

步骤 2：代入三角恒等式

将 $\cos \theta$ 和 $\sin \theta$ 代入 ${\cos }^2\theta +{\sin }^2\theta =1$：
$${\left(\frac{x-h}{r}\right)}^2+{\left(\frac{y-k}{r}\right)}^2=1$$

步骤 3：化简为标准形式

两边同时乘以 $r^2$，得到：
$$(x-h{)}^2+(y-k{)}^2=r^2$$
这就是圆的隐式方程。

---

3. 将隐式方程转换为参数方程

将隐式方程 $(x-h{)}^2+(y-k{)}^2=r^2$ 转换为参数方程，可以通过引入参数 $\theta$ 来表示圆上的点。具体步骤如下：

步骤 1：引入参数

假设 $x$ 和 $y$ 可以用角度 $\theta$ 表示：
$$\{\begin{array}{l}x=h+r\cos \theta \\ y=k+r\sin \theta \end{array}$$
其中，$\theta$ 是从圆心到圆上某点的连线与 x 轴正方向之间的夹角。

步骤 2：验证参数方程满足隐式方程

将 $x=h+r\cos \theta$ 和 $y=k+r\sin \theta$ 代入隐式方程：
$$(x-h{)}^2+(y-k{)}^2=(r\cos \theta {)}^2+(r\sin \theta {)}^2=r^2({\cos }^2\theta +{\sin }^2\theta )=r^2$$
满足隐式方程。

---

4. 总结

• 参数方程：通过角度 $\theta$ 描述圆上点的位置，形式为：
  $$\{\begin{array}{l}x=h+r\cos \theta \\ y=k+r\sin \theta \end{array}$$
• 隐式方程：不涉及参数，直接描述圆上点的坐标关系，形式为：
  $$(x-h{)}^2+(y-k{)}^2=r^2$$

通过上述方法，可以实现圆的参数方程和隐式方程之间的相互转换。

球面的参数方程与隐式方程式的转换

在解析几何中，球面的参数方程和隐式方程之间可以通过坐标变换进行相互转换。以下是一步步的解释：

---

1. 球面的隐式方程

球面在三维空间中的隐式方程通常表示为：
$$(x-h{)}^2+(y-k{)}^2+(z-l{)}^2=r^2$$
其中：

• $(h,k,l)$ 是球心坐标，
• $r$ 是球的半径。

这个方程描述了三维空间中所有到中心点 $(h,k,l)$ 的距离为 $r$ 的点。

---

2. 球面的参数方程

球面可以用两个角度参数 $\theta$ 和 $\varphi$ 来表示，通常称为极角（polar angle）和方位角（azimuthal angle）。参数方程的形式如下：
$$\{\begin{array}{l}x=h+r\sin \theta \cos \varphi \\ y=k+r\sin \theta \sin \varphi \\ z=l+r\cos \theta \end{array}$$
其中：

• $0\le \theta \le \pi$ 是极角，表示点与球心的连线与 $z$-轴之间的夹角，
• $0\le \varphi <2\pi$ 是方位角，表示在水平面（$x$-$y$ 平面）上的旋转角度。

在原点处的球面，即：
$$\begin{gathered} h=0 \\ k=0 \\ l=0 \end{gathered}$$
则：
$$\{\begin{array}{l}x=r\sin \theta \cos \varphi \\ y=r\sin \theta \sin \varphi \\ z=r\cos \theta \end{array}$$

---

3. 隐式方程到参数方程的转换

假设我们有球面的隐式方程：
$$(x-h{)}^2+(y-k{)}^2+(z-l{)}^2=r^2$$
我们可以引入球坐标系来将其转化为参数方程。

步骤 1：定义参数化

在三维空间中，球面对称性非常适合使用球坐标系来表示。我们可以用两个角度 $\theta$ 和 $\varphi$ 来描述球面上的每一点：

• $\theta$ 是从 $z$-轴向下测量的角度，
• $\varphi$ 是绕 $z$-轴旋转的角度。

步骤 2：引入球坐标变换

将隐式方程中的 $x,y,z$ 用球坐标参数化表示：
$$\{\begin{array}{l}x=h+r\sin \theta \cos \varphi \\ y=k+r\sin \theta \sin \varphi \\ z=l+r\cos \theta \end{array}$$
其中 $r$ 是半径，固定不变。

步骤 3：验证参数方程满足隐式方程

将参数方程代入隐式方程：
$$(x-h{)}^2+(y-k{)}^2+(z-l{)}^2=r^2$$
代入 $x,y,z$ 的表达式：
$$\begin{array}{rl}& [r\sin \theta \cos \varphi {]}^2+[r\sin \theta \sin \varphi {]}^2+[r\cos \theta {]}^2\\ & =r^2({\sin }^2\theta {\cos }^2\varphi )+r^2({\sin }^2\theta {\sin }^2\varphi )+r^2{\cos }^2\theta \\ & =r^2{\sin }^2\theta ({\cos }^2\varphi +{\sin }^2\varphi )+r^2{\cos }^2\theta \\ & =r^2{\sin }^2\theta (1)+r^2{\cos }^2\theta \\ & =r^2({\sin }^2\theta +{\cos }^2\theta )\\ & =r^2\end{array}$$
这验证了参数方程确实满足隐式方程。

---

4. 参数方程到隐式方程的转换

假设我们有球面的参数方程：
$$\{\begin{array}{l}x=h+r\sin \theta \cos \varphi \\ y=k+r\sin \theta \sin \varphi \\ z=l+r\cos \theta \end{array}$$
我们需要将其转换为隐式方程。

步骤 1：消去参数 $\theta$ 和 $\varphi$

从参数方程中解出 $\sin \theta \cos \varphi$、$\sin \theta \sin \varphi$ 和 $\cos \theta$：
$$\{\begin{array}{l}x-h=r\sin \theta \cos \varphi \\ y-k=r\sin \theta \sin \varphi \\ z-l=r\cos \theta \end{array}$$

步骤 2：平方并相加

将前两个方程分别平方：
$$(x-h{)}^2=r^2{\sin }^2\theta {\cos }^2\varphi ,(y-k{)}^2=r^2{\sin }^2\theta {\sin }^2\varphi$$
将它们相加：
$$(x-h{)}^2+(y-k{)}^2=r^2{\sin }^2\theta ({\cos }^2\varphi +{\sin }^2\varphi )=r^2{\sin }^2\theta$$

同时，对第三个方程平方：
$$(z-l{)}^2=r^2{\cos }^2\theta$$

步骤 3：相加得到隐式方程

将上述两个结果相加：
$$(x-h{)}^2+(y-k{)}^2+(z-l{)}^2=r^2{\sin }^2\theta +r^2{\cos }^2\theta =r^2({\sin }^2\theta +{\cos }^2\theta )=r^2$$
这就得到了球面的隐式方程：
$$(x-h{)}^2+(y-k{)}^2+(z-l{)}^2=r^2$$

---

总结

1. 球面的隐式方程是关于 $x,y,z$ 的二次方程，描述了所有到球心距离为 $r$ 的点。
2. 球面的参数方程使用两个角度 $\theta$ 和 $\varphi$ 来表示球面上的任意一点。
3. 通过将参数方程代入隐式方程并利用三角恒等式，可以验证两者的一致性。
4. 反过来，通过消去参数 $\theta$ 和 $\varphi$ 并平方相加，可以从参数方程推导出隐式方程。

参数曲线与离散函数表

源码文件：OCCT-7_9_0\src\ElCLib\ElCLib.cxx

曲线类型	离散函数
直线	gp_Pnt2d ElCLib::LineValue(const Standard_Real U, const gp_Ax2d& Pos)
圆	gp_Pnt2d ElCLib::CircleValue(const Standard_Real U, const gp_Ax22d& Pos, const Standard_Real Radius)
椭圆	gp_Pnt2d ElCLib::EllipseValue(const Standard_Real U, const gp_Ax22d& Pos, const Standard_Real MajorRadius, const Standard_Real MinorRadius)
抛物线	gp_Pnt2d ElCLib::ParabolaValue(const Standard_Real U, const gp_Ax22d& Pos, const Standard_Real Focal)
双曲线	gp_Pnt2d ElCLib::HyperbolaValue(const Standard_Real U, const gp_Ax22d& Pos, const Standard_Real MajorRadius, const Standard_Real MinorRadius)
样条曲线	BSplCLib

参数曲面与离散函数表

源码文件：OCCT-7_9_0\src\ElSLib\ElSLib.cxx

曲面类型	离散函数
平面	void ElSLib::PlaneD0(const Standard_Real U, const Standard_Real V, const gp_Ax3& Pos, gp_Pnt& P)
圆柱面	void ElSLib::CylinderD0(const Standard_Real U, const Standard_Real V, const gp_Ax3& Pos, const Standard_Real Radius, gp_Pnt& P)
圆柱面	void ElSLib::ConeD0(const Standard_Real U, const Standard_Real V, const gp_Ax3& Pos, const Standard_Real Radius, const Standard_Real SAngle, gp_Pnt& P)
球面	void ElSLib::SphereD0(const Standard_Real U, const Standard_Real V, const gp_Ax3& Pos, const Standard_Real Radius, gp_Pnt& P)
圆环面	void ElSLib::TorusD0(const Standard_Real U, const Standard_Real V, const gp_Ax3& Pos, const Standard_Real MajorRadius, const Standard_Real MinorRadius, gp_Pnt& P)
样条曲面	BSplSLib

参考文档：

• 刘星《基于Open Cascade的开发培训》(bilibili.com)
• OCC官网（https://dev.opencascade.org/release）
• A Primer on Bézier Curves
• https://zhuanlan.zhihu.com/p/421129261
• https://zhuanlan.zhihu.com/p/55878540