基于MATLAB的B样条计算与可视化

B样条基函数

令 U = { u 0 , u 1 , ⋯   , u m } \mathbf{U}=\{u_0,u_1,\cdots,u_m\} U={u0​,u1​,⋯,um​}为一个单调增的实数序列，即$u_i\le u_{i+1},i=0,\cdots ,m-1$。那么对于一个$p$次B样条，其第$i$个B样条基函数 $N_{i,p}(u)$可以由p−1阶的两个B样条基函数线性表出：

$$N_{i,p}(u)=\frac{u-u_i}{u_{i+p}-u_i}{N}_{i,p-1}(u)+\frac{u_{i+p+1}-u}{u_{i+p+1}-u_{i+1}}{N}_{i,p}(u)$$

其中$u_i$称为节点(knots)，$\mathbf{U}$称为节点向量。特别的，0阶的B样条基函数定义为:

$$N_{i,0}(u)=\{\begin{array}{ll}1& \text{if }{u}_i\le u<u_{i+1}\\ 0& \text{otherwise}\end{array}$$

B样条曲面

B样条曲面是由多条B样条曲线在u,v两个方向上由$(m+1)\times (n+1)$个控制点构成一张网格：

$$S(u,v)={\Sigma }_{i=0}^n{\Sigma }_{j=0}^m{N}_{i,p}(u)N_{j,p}(v)P_{i,j}$$

其中 $P_{ij}$为曲面的控制点，下图展示了一个基于双3次样条的贝壳曲面

基于MATLAB语言的实现

基于MATLAB平台开发了样条B-曲面的计算和可视化功能，通过3个实例演示了相关开发功能。

主要包括matlab软件图标、虎头、贝壳3个类型的曲面，函数自动读取控制点、节点向量，并构造和可视化样条，适合样条初学者学习

Demo演示

基于matlab的B样条曲面开发与可视化