使用matlab编写代码要求将形如x1 x2 x3 x4 y1 y2 y3 y4四组xy散点图分别进行B样条拟合并画出 并将拟合后的曲线求导后画出

以下是一个示例代码，其中对于每组散点图，我们使用B样条插值函数fit函数进行拟合，并使用ppval函数计算导数。注意，由于B样条插值是局部插值，因此在插值时需要为每个组合适当选择插值节点数。

% 测试数据
x1 = [0 1 2 3];
y1 = [0 1 4 9];
x2 = [0 1 2 3];
y2 = [0 2 4 6];
x3 = [0 1 2 3];
y3 = [0 -1 2 7];
x4 = [0 1 2 3];
y4 = [0 3 8 15];

% 插值节点数
k = 3;

% B样条拟合
pp1 = fit(x1.', y1.', 'cubicinterp').pp;
pp2 = fit(x2.', y2.', 'cubicinterp').pp;
pp3 = fit(x3.', y3.', 'cubicinterp').pp;
pp4 = fit(x4.', y4.', 'cubicinterp').pp;

% 画出拟合曲线
figure
subplot(2,2,1)
xx = linspace(x1(1), x1(end), 1000);
yy = ppval(pp1, xx);
plot(x1, y1, 'o', xx, yy)
title('Group 1')

subplot(2,2,2)
xx = linspace(x2(1), x2(end), 1000);
yy = ppval(pp2, xx);
plot(x2, y2, 'o', xx, yy)
title('Group 2')

subplot(2,2,3)
xx = linspace(x3(1), x3(end), 1000);
yy = ppval(pp3, xx);
plot(x3, y3, 'o', xx, yy)
title('Group 3')

subplot(2,2,4)
xx = linspace(x4(1), x4(end), 1000);
yy = ppval(pp4, xx);
plot(x4, y4, 'o', xx, yy)
title('Group 4')

% 画出导数曲线
figure
subplot(2,2,1)
dy = ppval(fnder(pp1), xx);
plot(xx, dy)
title('Group 1')

subplot(2,2,2)
dy = ppval(fnder(pp2), xx);
plot(xx, dy)
title('Group 2')

subplot(2,2,3)
dy = ppval(fnder(pp3), xx);
plot(xx, dy)
title('Group 3')

subplot(2,2,4)
dy = ppval(fnder(pp4), xx);
plot(xx, dy)
title('Group 4')