MATLAB B样条曲线拟合散点图 - 代码示例

使用MATLAB B样条曲线拟合散点图 - 代码示例

本文将使用MATLAB代码演示如何使用B样条曲线拟合一组散点图数据。假设您有一组横纵坐标分别为'x'和'y'的散点数据，例如：

x = [0.05, 0.103448, 0.151515, 0.2, 0.25];
y = [-0.447321, -0.277143, -0.101953, 0.009087, 0.10446];

以下代码展示了如何利用B样条曲线拟合这些数据并绘制拟合曲线：

% 构造节点矢量
t = [zeros(1,3), linspace(0,1,length(x)-2), ones(1,3)];

% 计算B样条基函数
B = spcol(t,3,length(x),linspace(0,1,100));

% 构造系数矩阵
A = zeros(length(x),length(x));
for i = 1:length(x)
    A(i,:) = B(i,:);
end

% 求解系数向量
c = A \ y';

% 计算拟合曲线
y_fit = B * c;

% 绘制散点图和拟合曲线
plot(x,y,'o',linspace(min(x),max(x),100),y_fit);

代码说明：

1. 构造节点矢量: 节点矢量t定义了B样条曲线的参数域。本例中，我们使用了一个长度为length(x)+4的节点矢量，其中前三个和后三个节点为0和1，中间部分均匀分布。
2. 计算B样条基函数: 使用spcol函数计算B样条基函数。spcol(t,3,length(x),linspace(0,1,100))表示计算三次B样条基函数，其中t为节点矢量，3为样条曲线阶数，length(x)为控制点个数，linspace(0,1,100)为参数域中的100个等间距点。
3. 构造系数矩阵: 系数矩阵A用于求解B样条曲线的系数向量c。矩阵的每一行对应于一个控制点，每一列对应于一个B样条基函数。
4. 求解系数向量: 使用A \ y'求解系数向量c，其中y'为纵坐标数据矩阵的转置。
5. 计算拟合曲线: 利用系数向量c和B样条基函数B计算拟合曲线y_fit。
6. 绘制散点图和拟合曲线: 使用plot函数绘制原始散点图和拟合曲线。

通过以上步骤，即可使用B样条曲线拟合给定散点数据，并绘制出相应的拟合曲线。