York同位素等时线拟合方法的Matlab程序设计

赵玉岩 郝立波 陆继龙 潘贵 侯立宾

等时线拟合方法的选择是提高同位素测年精度的关键因素之一。这里利用Matlab语言编写了基于York双误差回归模型的等时线拟合程序,包括isoacc.m、samw.m、errs.m和resi.m等模块,实现了读入数据,输入衰变常数和斜率迭代初始值,选择权重方法,以及计算截距、年龄、MSWD和绘图的功能。选用Brooks1972的数据,分别作三种权重等时线拟合,结果与Faure1977[9]的计算结果一致。程序简练,易使用,运算速度快,运行稳定。