MSD算扩散系数的几种方法

• 2018-07-19 22:10:28

利用均方位移MSD算扩散系数时, 最重要的是确定线性拟合的区间. 这在以前的博文中有说明. 在这里我总结一下确定拟合区间的几种方法, 并给出matlab代码供参考.

确定拟合区间的方法主要有下面几种:

1. 动态扩散系数(RDC), 也就是查看MSD/6t的变化
2. 稳健拟合
3. 最小协方差行列式方法(MCD)检测异常点, 然后去除异常点后拟合
4. 随机抽样一致性算法(RANSAC)
5. 全局优化

所有以上方法都不是确定性的方法, 所得结果都依赖于一些选择和参数. 其中最后两种方法本质上是随机方法, 无法保证每次所得结果相同.

下图是一个示例, 可以看到各种方法还是有所区别的, 其中MCD方法和RANSAC方法比较接近, 应该是比较好的方法.