scipy.interpolate.splprep函数系数推导逻辑及结果正确性验证方法

一、系数推导逻辑 #

你得到的系数为2×3结构是参数化B样条的正常输出，推导流程如下：

• splprep 是针对参数化B样条曲线的插值工具，x、y两个坐标分量各自作为参数u的独立B样条函数求解，因此会分别生成两个分量的系数数组，最终拼接为2×3的结构。你参考的资料中提到的3×3矩阵是单个坐标分量求解时用到的基函数矩阵，并非最终系数的维度。
• 具体求解步骤：
  1. 首先完成你已经验证的前置步骤：按弦长参数化得到参数数组u，再生成钳位节点向量，k=2（二次）的钳位样条节点首尾重复k+1=3次，最终节点长度为n+k+1=3+2+1=6，和你得到的[0, 0, 0, 1, 1, 1]完全一致。
  2. 对x、y分量分别求解线性方程组：构造3×3的基函数矩阵N，其中第i行第j列的元素为第j个二次B样条基函数在第i个参数u[i]处的取值，分别求解N * c_x = x原始数组、N * c_y = y原始数组，得到的c_x、c_y就是你输出的系数数组。

二、系数正确性验证方法 #

你已经完成了最核心的验证，除此之外还有多种验证维度：

• 端点特性验证：二次钳位B样条曲线的首尾点严格等于首尾控制点，你得到的c数组第一个元素对应x=150、y=149，最后一个元素对应x=148、y=79，正好和你的原始输入首尾点完全匹配，符合二次钳位B样条的特性。
• 手动基函数计算验证：你可以手动写出当前节点向量下的3个二次B样条基函数表达式，任意取参数u的数值，代入基函数计算得到的x、y值和splev的输出对比，结果完全一致即可验证系数正确。
• 切线方向验证：二次B样条曲线在u=0处的切线方向和第一个控制点到第二个控制点的连线方向完全一致，在u=1处的切线方向和第二个控制点到第三个控制点的连线方向完全一致，你可以计算曲线在首尾点的导数，和控制点连线的方向对比，符合即可验证系数正确。
• 连续性验证：二次B样条曲线在非重复节点处是C1连续的，你可以计算u在节点处的左右导数，结果相等即可验证系数正确。

内容的提问来源于stack exchange，提问作者jani-c