使用origin软件拟合非球面透镜曲线方程实用方1.docVIP

使用origin软件拟合非球面透镜曲线方程的实用方法 赵润2010-10-3 问题来源 我们采用费马原理，很容易求出某个透镜的曲面的各点坐标（见作者较早写的文章），但要把这一系列坐标拟合成非球面透镜的标准方程（下式）似乎没有现成的方法。 y=x2R+R2-(1+K)x2+B4x4+B6x6+B8x8+B10x10+… origin软件中有自定义非线性拟合的功能，但用过的朋友一定知道，如果方程中的参数很多（如上式），直接拟合很多参数，拟合曲线与实际曲线之间的误差会很大，甚至因为计算量太大，最终得不到结果。下面将详细给出手动拟合的具体方法。 二次曲线拟合 首先使用二次曲线拟合，拟合公式为： y=x2R+R2-(1+K)x2 当然也是使用origin软件中的 “自定义非线性拟合”功能，因为参数只有R,K两项，拟合速度和拟合效果都不错。并且我们知道很多时候曲线是接近椭圆的，所以K的取值在区间（-1，0），所以随意给个在此区间的初始值，如-0.5，就可以了。 对二次曲线拟合后的误差进行多项式拟合 我发现直接用下式拟合第一次拟合产生的误差，效果很差。 y=B4x4+B6x6+B8x8+B10x10+… 而用更一般的多项式形式（下式），拟合效果会很好，所以我们要采用下式拟合： y=B2x2+B4x4+B6x6+B8x8+B10x10+… 这样得到的拟合参数为：R，K，B2，B4，B6，B8，B10，... 。与最前面给出的非球面方程的标准形式相比多了一项B2，将B2参数消掉是需要的。 4、消掉参数B2 因为在x值较小的区域有： x2R+R2-(1+K)x2≈x22R 所以我们定义一个R’，求出数值： 12R+B2=12R 或R=(R-1+2B2)-1 其中R和B2为步骤2，3中拟合出的参数，这样我们就得到了R’这个参数的数值。 5、重新计算二次曲线拟合后的误差 用R’的数值代替R，并用步骤2拟合出的参数K，使用公式： y=x2R+R2-(1+K)x2 用“set column value”的方法，在原始数值表中增加两列（第一列为x坐标col(a)，第二列为y坐标col(b)，增加的为第三列col(c)为y’和第四列col(d)为误差δy），设置数值如下： col(c)=col(a)^2/(R’+sqrt(R’^2-(1+K)*col(a)^2)) col(d)=col(b)-col(c) 6、从4次项开始对δy进行多项式拟合 使用公式y=B4x4+B6x6+B8x8+B10x10+… 对δy进行拟合。 小结：经过步骤2-6步，我们拟合出了非球面透镜曲线标准方程的所有参数：R，K，B4，B6，B8，B10...。其中R为第4步计算出的R’，K为第2步拟合出的，而B4，B6，B8，B10...为第6步拟合出的。