所谓软件修正算法,就是找到从 X到 Y 的映射。由解析几何知识我们知道X到 Y 的映射可以通过一个线性变换C来表示:
这里再多说几句,我个人喜欢使用 gnuplot,gnuplot中有个 fit 命令,也可以完成多元函数拟合。Gnuplot的 fit 命令采用非线性最小二乘 (NLLS) Marquardt-Levenberg算法,拟合能力非常强悍,只要拟合参数的初值给的别太离谱,基本都能够收敛到最优解上。fit命令具体的用法可以参考gnuplot的帮助文件。这里只讲讲如何处理 y是多个值的问题,方法很简单,就是通过多引入一个自变量,将多值函数改造成普通的函数: