Matlab中使用RANSAC拟合直线报错求助：fitFun类型不符问题

问题背景 #

我用RANSAC算法对30×2的double类型数据进行直线拟合，参考Matlab官方示例编写代码后运行报错，无法定位问题。

原始代码 #

load linedata
data = [xm,ym];
N = length(xm); % 数据点数量
sampleSize = 2; % 每次迭代采样的点数
maxDistance = 2; % 内点允许的最大距离
fitLineFcn = polyfit(xm,ym,1); % 用polyfit定义拟合函数
evalLineFcn =@(model) sum(ym - polyval(fitLineFcn, xm).^2,2); % 距离评估函数
[modelRANSAC, inlierIdx] = ransac(data,fitLineFcn,evalLineFcn,sampleSize,maxDistance);

错误提示 #

使用ransac时出错
预期fitFun为以下类型之一：
function_handle
但其实际类型为double。
出错位置：ransac>parseInputs (第202行)
validateattributes(fitFun, {'function_handle'}, {'scalar'}, mfilename, 'fitFun');
出错位置：ransac (第148行)
[params, funcs] = parseInputs(data, fitFun, distFun, sampleSize, ...

问题分析与修复 #

• 错误根源：fitLineFcn的赋值不符合RANSAC要求。直接调用polyfit(xm,ym,1)会得到一个存储直线系数的double数组，但RANSAC要求fitFun必须是函数句柄——这个函数需要接收随机采样的子集数据，返回对应模型的参数。
• 修复步骤：
  1. 重构fitLineFcn为函数句柄，使其能处理采样的子集数据：

     fitLineFcn = @(sampleData) polyfit(sampleData(:,1), sampleData(:,2), 1);
     

  2. 修正evalLineFcn逻辑：原代码错误地使用了全局拟合的参数，正确的做法是接收RANSAC传入的model参数，计算所有点到该模型的距离：

     evalLineFcn = @(model) sum((ym - polyval(model, xm)).^2, 2);
     

• 修复后完整代码：

  load linedata
  data = [xm,ym];
  N = length(xm); % 数据点数量
  sampleSize = 2; % 每次迭代采样的点数
  maxDistance = 2; % 内点允许的最大距离
  fitLineFcn = @(sampleData) polyfit(sampleData(:,1), sampleData(:,2), 1); % 拟合函数句柄
  evalLineFcn = @(model) sum((ym - polyval(model, xm)).^2, 2); % 距离评估函数
  [modelRANSAC, inlierIdx] = ransac(data,fitLineFcn,evalLineFcn,sampleSize,maxDistance);
  

额外说明 #

RANSAC的核心逻辑是：每次随机采样数据子集 → 调用fitFun拟合模型 → 用distFun筛选内点。因此fitFun必须是能动态处理不同子集的函数句柄，而不是提前拟合好的固定参数数组。

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