陶瓷材料
第26卷 第3期 2006年8月
热 带 地 理
TROPICALGEOGRAPHY
Vol126,No13
Aug.,2006
基于MATLAB的图解粒度参数计算
王 为,吴 正
(华南师范大学地理系,广州510631)
摘 要:粒度特征是沉积物的基本特征之一。,其中图解法必须通过手工作图求累积曲线,是一项相当繁杂的劳动,的函数功能求解百分位数和计算粒度参数。计算结果和验算表明,。该方法原
理简单,,既,有很大的应用价值。关键词:中图分类号:.文献标识码:A 文章编号:1001-5221(2006)03-0239-04 粒度特征受搬运和沉积过程的动力条件控制,
是沉积物的基本特征,与沉积环境关系密切。因此,在沉积分析中往往应用沉积物的粒度参数来判断和恢复沉积环境。计算粒度参数的方法很多,目前应用最广泛的主要是Friedment的矩法和Folk的
[2-3]
图解法两种。有人认为矩法在理论计算和地质
[4-7]
应用效果上都比图解法好,但也有人认为矩法
[8]
不适用于第四纪的沉积物粒度数据分析。仅从粒
[9]
度分析的效果看,很难说那一种方法更好。矩法有一个非常明显的优点,就是它能用于计算机计算,通过简单的编程,就能计算大量样品的粒度参数,省时省工。图解法很复杂,需要手工作图,当要计算大量样品的粒度参数时,要消耗大量的人力和时间。另外,手工操作有一定的不确定性,同样的样品,不同的人作图很难作出完全一样的结果。目前,国外已有用VB语言借助Excel开发的专用粒度分析软
[10]
件GRANDISTAT,可以用于图解法计算,但是,还没有人能利用通用软件进行图解法计算。
MATLAB被称作第四代计算机语言,与BASIC、FORTRAN、C等计算机语言相比较,其编程效率高
[11]
10倍。本文阐述的方法利用MATLAB语言的函数功能,无须作图以图解法计算沉积样品粒度累积曲线的百分位数和粒度参数,其过程并不比矩法复杂,有极大的应用价值。
[1]
1 计算原理
图解法首先要在以纵轴为累积百分比%、横轴为粒径φ值的米格纸上点取样品的各级累积百分含量,然后将各点连成累积曲线,再从累积曲线上点取各个百分位数(图1a),最后由百分位数计算粒度参数。将图解求累积曲线的百分位数的过程转化为数学问题,就是一个不等区间的一维函数内插值问题。如图1,将(a)绕沿夹角为45
度过坐标原点的轴旋
图1 将求累积曲线的百分位数转化为函数插值
Fig.1 Thetranaformofpercentilevaluecalculationinto
functionalinterpolation
转180度,使(a)中粒级累积值的纵轴变为横轴(X轴),粒径φ值的横轴变为纵轴(Y轴)。累积曲线因此可以看成是由已知点构成的函数y=f(x),求百
收稿日期:2005-12-19;修订日期:2006-02-14 基金项目:广东省自然科学基金资助项目(031529)
作者简介:王为(1956-)男,广东人,博士,教授,主要从事海岸地貌和环境研究与教学,(E2mail)wangw@http://doc.xuehai.net。