基于贝塞尔曲线拟合的GPS路径真实位置预测及不确定性疑问

我有一组代表多次沿固定路径行驶的GPS坐标点（x-y数据），想用统计方法预测路径的真实位置并评估预测精度。我已经对数据排序，并用三次贝塞尔曲线拟合，MATLAB代码如下：

% fit the curve
[px, r_x, J_x, CovB_x, mse_x] = nlinfit(t,xOrdData, @bezier, x0);
[py, r_y, J_y, CovB_y, mse_y] = nlinfit(t,yOrdData, @bezier, y0);

% calculate the predicted confidence interval @ 99% certainty
[xpred,xdelta] = nlpredci(@bezier,t, px, r_x, "Covar", CovB_x, "MSE", mse_x, "Alpha", 1-0.99);
[ypred,ydelta] = nlpredci(@bezier,t, py, r_y, "Covar", CovB_y, "MSE", mse_y, "Alpha", 1-0.99);
x_lower = xpred - xdelta;
x_upper = xpred + xdelta;
y_lower = ypred - ydelta;
y_upper = ypred + ydelta;

% plot the results
figure;
hold on;
plot(xOrdData, yOrdData, '.r');
plot(bezier(px,t), bezier(py,t),'-g');
plot(px,py,"*k");
plot(x_lower,y_lower,'k');
plot(x_upper,y_upper,'k');


function retval = bezier(p, t)
    if size(t,1)>size(t,2)
        t=t';
    end
    TT=[(1-t).^3; 3*(1-t).^2.*t; 3*(1-t).*t.^2; t.^3];
    retval = p*TT;
end

拟合后生成了路径数据图和放大图，但出现疑问：为什么没有显示更大的不确定性区域？看起来95%的数据点都不在当前显示的区间内。是不是因为采用参数曲线x=F(t)、y=G(t)的形式，误用了MATLAB的nlinfit和nlpredci工具？还是有其他问题？

核心问题：参数曲线的独立拟合导致置信区间偏差 #

你分别对x、y分量独立拟合贝塞尔曲线并计算置信区间，这种方式存在三个关键问题：

• 忽略x-y相关性：GPS数据的x、y误差通常是相关的，独立拟合会割裂这种关联，导致计算出的区间无法反映路径的整体不确定性。
• 参数t定义模糊：如果t只是简单的排序索引，而非与路径进度（如行驶时间、累计长度占比）严格绑定，拟合残差会违背nlinfit假设的独立同分布特性，进而扭曲置信区间的计算结果。
• 可视化方式错误：直接将x、y的上下界分别连线，得到的是单分量的区间边界，而非路径点的联合置信区域。实际的不确定性应该是每个t对应的椭圆或不规则区域，而非两条独立曲线。

修正方案 #

1. 改用二维联合拟合模型：将x、y作为整体拟合，同时纳入x-y的协方差信息，示例框架如下：

% 定义二维贝塞尔拟合函数
function [x,y] = bezier2d(p, t)
    if size(t,1)>size(t,2)
        t=t';
    end
    TT = [(1-t).^3; 3*(1-t).^2.*t; 3*(1-t).*t.^2; t.^3];
    x = p(1:4)*TT;
    y = p(5:8)*TT;
end

% 合并响应变量与初始参数
response = [xOrdData; yOrdData];
p0 = [x0; y0];

% 联合拟合
[p, r, J, CovB, mse] = nlinfit(t, response, @(p,t) bezier2d(p,t), p0);

2. 明确t的物理意义：将t定义为每个GPS点对应的累计路径长度占总长度的比例，确保t与路径进度严格对应，让残差符合独立同分布假设。
3. 可视化联合置信区域：对每个t值，根据x、y的协方差计算置信椭圆，通过绘制一系列椭圆展示路径的不确定性（可使用MATLAB File Exchange的ellipse函数或自定义代码实现）。
4. 验证残差假设：检查拟合残差的正态性、独立性，若存在自相关或异方差，需调整模型（如加入时间序列成分）或使用稳健拟合方法。

内容的提问来源于stack exchange，提问作者Eanema