让移动最小二乘法拟合曲线首末端通过定点，写成matlab代码

好的，移动最小二乘法拟合曲线的代码如下：

% 定义首末端点的坐标
x0 = 0; y0 = 1;
xn = 10; yn = 2;

% 定义拟合曲线的阶数
n = 3;

% 定义移动窗口的大小和步长
w = 5;
step = 1;

% 初始化拟合曲线的系数
C = zeros(n+1, 1);

% 逐步移动窗口，并进行拟合
for i = x0+w/2:step:xn-w/2
    % 构造系数矩阵A和常数矩阵B
    A = zeros(n+1, n+1);
    B = zeros(n+1, 1);
    for j = 1:n+1
        for k = 1:n+1
            A(j, k) = sum((i-w/2:i+w/2).^(j+k-2));
        end
        B(j) = sum(((i-w/2:i+w/2)-x0).^(j-1)) - y0*((n==j-1)*1);
    end

    % 求解线性方程组Ax=B
    Ci = A\B;

    % 更新拟合曲线的系数
    C = (i-x0)^2/(i-x0+w/2)^2 * C + (i-x0+w/2)^2/(i-x0)^2 * Ci;

    % 绘制拟合曲线
    x = i-w/2:0.1:i+w/2;
    y = 0;
    for j = 0:n
        y = y + C(j+1)*(x-x0).^j;
    end
    y = y + y0;
    plot(x, y);
    hold on;
end

% 绘制首末端点
plot([x0, xn], [y0, yn], 'o');

其中，首末端点的坐标分别为`(x0, y0) = (0, 1)`和`(xn, yn) = (10, 2)`，拟合曲线的阶数为`n=3`，移动窗口的大小为`w=5`，步长为`step=1`。代码中逐步移动窗口，并在每个窗口中进行最小二乘法拟合，得到拟合曲线的系数C。为了使拟合结果更加平滑，采用了加权移动平均的方法来更新拟合曲线的系数C。最后，根据系数C绘制拟合曲线，并在图中绘制首末端点。