|
2018-07-19 22:10:28
利用均方位移MSD算扩散系数时, 最重要的是确定线性拟合的区间. 这在以前的博文中有说明. 在这里我总结一下确定拟合区间的几种方法, 并给出matlab代码供参考.
确定拟合区间的方法主要有下面几种:
动态扩散系数(RDC), 也就是查看MSD/6t的变化
最小协方差行列式方法(MCD)检测异常点, 然后去除异常点后拟合
全局优化
所有以上方法都不是确定性的方法, 所得结果都依赖于一些选择和参数. 其中最后两种方法本质上是随机方法, 无法保证每次所得结果相同.
下图是一个示例, 可以看到各种方法还是有所区别的, 其中MCD方法和RANSAC方法比较接近, 应该是比较好的方法.