数学地质实验指导书（教材）(5)

三、程序设计框图

建议以实验一的原始数据格式输入，首先要对数据进行预处理，处理完成后，形成实验一的原始数据格式的新文件（详见图1－1）。以下步骤如流程图（图3－1）。

开始打开标准格式的原始数据文件数据预处理和统计模块输出变换后的标准格式数据文件读入n对数据（xij,yi）（i＝1，2，?，n;j＝1，2，?，p）,形成原始矩阵（见式3－3）?1计算X?,(XX?)?1,?XX??X?Y(见式3?7）请参考线性代数矩阵?1根据B?XX??X?Y,求得b0,b1,?,bp运算。?i?b0?b1xi1?b2xi2???bpxip求y（i?1,2,?,n）n?S?yi2?ny2?总?i?1n??i2?ny2?S回?y计算?i?1?S?S?S剩总回?S回p?F???S剩n?p?1???计算r?S回S总输入F?,r?（或者建立数据库，输 入?，系统自动查询F?,r? ）F>Fα且r>rαY输出F,r,S，F?,r???b0?b1x1?b2x2???bpxp读入待预测样品x1,x2,?,xp值，求出y计算剩余标准差???S剩/(n?p?1)N??2??L21＝y??2??L22＝y??3??L31＝y??3??L32＝y输出95％的置信区间（L22,L21)和99.7%的置信区间（L32,L31)结束

图3-1 多元线性回归分析模块流程图

13

四、算例

（一）、某矿区从18个矿样中测得Cu,Pd,Pt的含量如下所示，试求Pt对Cu和Pd的回归方程，并检验其有无实用价值。 编号 1 2 3 4 5 6 7 8 9 x1 Cu(%) 0.36 0.05 0.24 0.05 0.22 0.13 0.10 0.07 0.16 x2 Pd(g/t) 0.169 0.035 0.248 0.015 0.216 0.086 1.300 0.009 0.160 x3 Pt(g/t) 0.176 0.020 0.262 0.035 0.234 0.092 1.669 0.006 0.150 编号 10 11 12 13 14 15 16 17 18 x1 Cu(%) 0.06 0.53 0.06 0.21 0.03 0.33 0.01 0.21 0.08 x2 Pd(g/t) 0.023 0.329 0.065 0.283 0.030 0.21 0.030 0.252 0.053 x3 Pt(g/t) 0.045 0.292 0.080 0.235 0.030 0.237 0.050 0.24 0.103 （二）、由于碳、氢、氧是煤燃烧过程中产生热量的主要元素，对某矿（褐煤、常烟煤、肥煤、

r焦煤和无烟煤）中取12块煤样，经分析化验后，其发热量QDT（焦耳/克）与Cr(%)、Hr（%）、Or（%）

元素的数据如下。问：今后能否不再对该矿煤进行发热量测定，而由元素分析结果对其进行预测？

编号 1 2 3 4 5 6 x1 62 70 75 75 78 80 x2 5.0 6.0 6.5 5.0 5.5 5.2 x3 15 20 25 10 15 4.0 y rDT rrC(%) H（%） Or（%） Q（J/g） 6.5 7.0 7.2 7.5 7.7 8.0 编号 7 8 9 10 11 12 rrr C(%) H（%） Or（%） QDT（J/g）x1 x2 x3 y 85 88 90 90 92 95 6.0 3.0 5.0 2.5 3.0 3.5 6.0 3.0 5.0 2.5 3.0 3.5 8.4 8.5 8.8 8.0 8.5 8.7

14

试验四 趋势面分析

一、目的：

通过趋势面分析程序设计及完成算例，掌握趋势面分析方法原理及工作步骤。

二、方法概要

1、建立原始数据矩阵，并进行必要的预处理

（1）选取应变量

选取什么样的变量进行趋势分析，取决于研究对象和研究目的。 （2）确定控制点

在确定控制点时，一般要考虑以下几个方面：①控制点的代表性、真实性和可靠性。②控制点要具有面性、均匀、随机分布的特点。

（3）整理原始数据统计表

收集因变量观测值及相应坐标，填入趋势分析原始数据统计表，其格式可参考表4.1 。

表4.1 趋势分析原始数据统计表格式

编 号

顺序号

原编号

控制点坐标

横坐标(u)