平均品位计算方法.docVIP

平均品位计算方法

平均品位计算新方法

张亮

（江西贵溪银矿，江西贵溪流335413）

摘要：本文采用三次多项式曲线拟合方法，对单位工程品位指标值多边形图进行分段拟合处理，使得多边形图成为一条连续光滑的复合曲线，并对拟合的复合曲线进行积分，求出品位平均值，并给出计算编程计算程序及使用说明。

关键词：单工程，三次多项式，拟合，平均品位计算，新方法。

引言

如图1所示，沿某一直线方向品位的直方图和多边形图，设图1中的直方图间距分别为L1、L2、…Ln、且L=L1+L2+…+Ln（L取样品跨度），与其相对应的品位指标什值分别为C1、C2、…Cn。现按常规的加权平均值法计算品位平均值，有

又设Si=Li×Ci，则上式变为

……(1)

Y

P2

Pi

P1

P0

P-1

SHAPE

上述p-1、p0、pn-1、pn+2各点坐标仅供曲线拟合计算之用，不参与平均品位的计算。

3、拟合曲线方程系数求解

如图2，以曲线边pi～pi+1为例，由于三次多项式曲线通过pi点和pi+1点，将pi点、pi+1点的坐标值代入（2）式有

yi=a0+a1xi+a2x2i+a3x3i……(4)

yi+1=a0+a1xi+1+a2x2i+1+a3x3i+1……(5)

又由于曲线在pi点、pi+1处有一阶导数ti和ti+1,则有

ti=al+2a2xi+3a3x2i……(6)

ti+1=a1+2a2xi+1+3a3x2i+1……(7)

由(4)、(5)、(6)、(7)解方程组可得

……(8)

平均值的计算

品位指标值平均也即为图2中复印合曲线函数的平均值，而复

y

合曲线函数的平均平均等于函数曲线与坐标轴所围成的面积除以函数曲线跨度。

如图2所示，复合函数曲线跨度（即取样总长）OG=L，P1，P2，．．．，Pn为边形图转折点；H点为P0~P1段曲线与Y轴的交点，其坐标为H（0，yh）；K点为Pn~Pn+1段曲线与取样边界线（边界品位线）GK的交点，交点坐标为K（L，yk）。

设HOEP1所围成的面积为Sa，EP1PiPnF所围成的面积为Sb，FPnKG所围成的面积为Sc，则复合曲线所围成的总面积为：

S=Sa+Sb+Sc．．．．．．（9）

由区边Pｉ~Pｉ+1所确定的面积为：

dX

=．．．．．．（10）

利用（10）式可以计算每分段的面积，且（10）式中各分段曲线方程系数am（m=0，1，2，3）由（8）式所得。由P1~Pn所围成的总面积Sb为：

Sb=．．．．．．（11）

同理可得：

s=

=．．．．．．（12）

=．．．．．．（13）

由（11）式，（12）式，（13）式可得总面积S，故平均品位C

为：C=S/L．．．．．．（14）

上述解算虽然简单，但计算点较我时，其计算工作量较大，如果采用计算机编程计算，那么上述计算十分方便、快捷。

5、计算程序设计

5．1计算程序

由于BASIC语言是最基本且通用性较强的语言，故此，本文采用BASIC语言编制计算程序，程序如下：

PRINT“取样个数，取样长度分别为：”：

20INPUTN，L

30DIMX（N+2），Y（N+2），T（N+1），A（3）

40FORI=0T0N+2

50READX（I），Y（1）

60NEXTI

65FORI=ITON+1

70T（I）=1/2*（（Y（I）-Y（I-1））/（X（I）-X（I-1））+（Y（I+1）

-Y（1））

X（N+2），Y（N+2）存多边形图中P0~Pn+2各点的坐标值

T（N+1）存多边形图中P0~Pn+1各点的导数值

A（3）存每间分段三次多项式曲线系数

X1，Y1为P-1点的坐标值

（2）运行过程

将多边形图中各点（P0~Pn+2）横纵坐标依次置入DATA语句中，运行该程序，输入N，L入P-1点的坐标X1，Y！的值后，转入子程序500语句计算P0~P1分段三次多项式系数，接着，转入子程序600语句计算H-P1分段的面积Sa。自220语入~380语句，计算P1~K段的面积（Sb+Sc），最后计算并打印出