MATLAB B样条曲线拟合与求导可视化

以下是使用MATLAB进行B样条分别拟合四组数据并求导后画在同一副图中的代码：

% 数据输入
x1 = [0.1, 0.5, 0.9, 1.3, 1.7, 2.1, 2.5, 2.9, 3.3, 3.7, 4.1, 4.5, 4.9, 5.3, 5.7, 6.1, 6.5, 6.9, 7.3, 7.7];
y1 = [0.39, 0.35, 0.4, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42, 0.38, 0.36, 0.36, 0.38, 0.42, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42];
x2 = [0.1, 0.5, 0.9, 1.3, 1.7, 2.1, 2.5, 2.9, 3.3, 3.7, 4.1, 4.5, 4.9, 5.3, 5.7, 6.1, 6.5, 6.9, 7.3, 7.7];
y2 = [0.39, 0.35, 0.4, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42, 0.38, 0.36, 0.36, 0.38, 0.42, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42];
x3 = [0.1, 0.5, 0.9, 1.3, 1.7, 2.1, 2.5, 2.9, 3.3, 3.7, 4.1, 4.5, 4.9, 5.3, 5.7, 6.1, 6.5, 6.9, 7.3, 7.7];
y3 = [0.39, 0.35, 0.4, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42, 0.38, 0.36, 0.36, 0.38, 0.42, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42];
x4 = [0.1, 0.5, 0.9, 1.3, 1.7, 2.1, 2.5, 2.9, 3.3, 3.7, 4.1, 4.5, 4.9, 5.3, 5.7, 6.1, 6.5, 6.9, 7.3, 7.7];
y4 = [0.39, 0.35, 0.4, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42, 0.38, 0.36, 0.36, 0.38, 0.42, 0.48, 0.54, 0.56, 0.54, 0.48, 0.42];

% B样条拟合
t = linspace(0, 1, 200); % 参数t
p1 = spline(x1, y1, t); % 第一组数据的B样条拟合
p2 = spline(x2, y2, t); % 第二组数据的B样条拟合
p3 = spline(x3, y3, t); % 第三组数据的B样条拟合
p4 = spline(x4, y4, t); % 第四组数据的B样条拟合

% 一阶导数
dp1 = gradient(p1) ./ gradient(t); % 第一组数据的一阶导数
dp2 = gradient(p2) ./ gradient(t); % 第二组数据的一阶导数
dp3 = gradient(p3) ./ gradient(t); % 第三组数据的一阶导数
dp4 = gradient(p4) ./ gradient(t); % 第四组数据的一阶导数

% 画图
figure;
plot(x1, y1, 'o', 'MarkerFaceColor', 'b', 'MarkerEdgeColor', 'b'); % 第一组数据的散点图
hold on;
plot(x2, y2, 'o', 'MarkerFaceColor', 'r', 'MarkerEdgeColor', 'r'); % 第二组数据的散点图
plot(x3, y3, 'o', 'MarkerFaceColor', 'g', 'MarkerEdgeColor', 'g'); % 第三组数据的散点图
plot(x4, y4, 'o', 'MarkerFaceColor', 'm', 'MarkerEdgeColor', 'm'); % 第四组数据的散点图

plot(p1, 'b', 'LineWidth', 2); % 第一组数据的B样条曲线
plot(p2, 'r', 'LineWidth', 2); % 第二组数据的B样条曲线
plot(p3, 'g', 'LineWidth', 2); % 第三组数据的B样条曲线
plot(p4, 'm', 'LineWidth', 2); % 第四组数据的B样条曲线

plot(dp1, 'b--', 'LineWidth', 2); % 第一组数据的B样条曲线的一阶导数
plot(dp2, 'r--', 'LineWidth', 2); % 第二组数据的B样条曲线的一阶导数
plot(dp3, 'g--', 'LineWidth', 2); % 第三组数据的B样条曲线的一阶导数
plot(dp4, 'm--', 'LineWidth', 2); % 第四组数据的B样条曲线的一阶导数

legend('Data 1', 'Data 2', 'Data 3', 'Data 4', 'Curve 1', 'Curve 2', 'Curve 3', 'Curve 4', 'Derivative 1', 'Derivative 2', 'Derivative 3', 'Derivative 4');
xlabel('X');
ylabel('Y');
title('B-Spline Fitting and Derivative');

运行以上代码，会得到如下图所示的结果：

本代码演示了如何在MATLAB中使用B样条对四组数据进行拟合，并绘制了拟合后的曲线和一阶导数。代码中使用spline函数进行B样条拟合，gradient函数计算一阶导数，plot函数绘制图形。最后，使用legend函数添加图例，xlabel和ylabel函数添加坐标轴标签，title函数添加标题。

注意:

• 代码中的x1, x2, x3, x4和y1, y2, y3, y4分别代表四组数据的横坐标和纵坐标值。您可以根据自己的需要修改这些数据。
• 代码中的t参数控制了B样条曲线的平滑度。您可以根据自己的需要调整linspace函数的第三个参数。
• 代码中使用了一阶导数，您可以根据需要计算更高阶的导数。

希望本代码和解释对您有所帮助。