采矿系统工程上机实验报告.docx
《采矿系统工程上机实验报告.docx》由会员分享，可在线阅读，更多相关《采矿系统工程上机实验报告.docx（13页珍藏版）》请在冰豆网上搜索。
采矿系统工程上机实验报告
中国矿业大学矿业工程学院
采矿系统工程
实验报告
姓名
学号0110
班级采矿工程
指导教师周伟
日期2013年12月
实验报告一
1.实验名称：
用半变异函数解决地质学统计问题
2.题目
设Z（x）为二维区域化变量，满足本征假设，在图所示的正方形网格节点处有一批已知数据，小正方形的边长为a,试计算出图所示α1（东西），α2（南北），α3（北东），α4（北西）四个方向的实验变差函数值，只计算前三个的h值，并做出α1（东），α3（北）两个方向的实验变图案来。
35
35
33
33
34
31
35
37
41
41
35
35
35
35
33
41
37
35
37
35
37
37
39
39
41
37
40
42
34
35
41
34
37
41
33
35
42
33
39
31
30
注：
以上每个数字的位置本应在其所在表格的左上角点处
3.
程序流程图
4.实验思路
1、建立源数据文本文档xtgc.txt，作为输入数据源；
2、在vb窗体中绘制所需图框，其中text1,text2作为接收输入的计算起始和终止步距；4个复选框为选择计算的方向，text3用于输出结果，picture1用于输出绘图内容；
3、编写代码，先编写输入数据的代码，再编写东西、南北、北东、北西四个方向计算方法的子函数，其中子函数的设计是根据方向和输入数据维数的不同运用数组编写，再编写调用子函数的过程，最后根据计算结果编写绘图程序。
4、运行程序，输入步距分别为100,300，再选择计算方向，点击输入数据按钮，再点击计算，文本框中出现显示结果，再点击绘图，在绘图框中出现r（h）的图像。
5.程序代码
OptionBase1
Dima（8,10）AsInteger'数组a用于存放输入数据
Dimb（3）AsSingle
Dimc（3）AsSingle
Dimd（3）AsSingle
Dime（3）AsSingle
PrivateSubCommand1_Click（）
DimiAsInteger,jAsInteger
Open"c:
\ckxtgc1\xtgc.txt"ForInputAs#1
DoWhileNotEOF
（1）
Fori=1To8
Forj=1To10
Input#1,a（i,j）
Nextj
Nexti
Loop
Close#1
EndSub
PrivateSubcommand2_click（）
DimhAsInteger,r1AsSingle,r2AsSingle,r3AsSingle,r4AsSingle,b1AsInteger,b2AsInteger,kAsInteger
b1=Val（h1.Text）
b2=Val（hn.Text）
Fork=1To（b2-b1+100）/100
IfCheck1.Value=1Then
Callewr（k,r1）
b（k）=r1
xtr.Text=xtr.Text&"东西方向r（"&k*100&"）:
"&Str（r1）&vbCrLf
EndIf
IfCheck2.Value=1Then
Callsnr（k,r2）
c（k）=r2
xtr.Text=xtr.Text&"南北方向r（"&k*100&"）:
"&Str（r2）&vbCrLf
EndIf
IfCheck3.Value=1Then
Callner（k,r3）
d（k）=r3
xtr.Text=xtr.Text&"北东方向r（"&k*100&"）:
"&Str（r3）&vbCrLf
EndIf
IfCheck4.Value=1Then
Callnwr（k,r4）
e（k）=r4
xtr.Text=xtr.Text&"北西方向r（"&k*100&"）:
"&Str（r4）&vbCrLf
EndIf
Nextk
EndSub
PrivateSubewr（ByValbAsInteger,rAsSingle）'东西方向变异函数
DimiAsInteger,jAsInteger,a（8,10）AsInteger,jsAsInteger
r=0
Open"c:
\ckxtgc1\xtgc.txt"ForInputAs#1
DoWhileNotEOF
（1）
Fori=1To8
Forj=1To10
Input#1,a（i,j）
Nextj
Nexti
Loop
Close#1
Fori=1To8
Forj=1To10-b
If（a（i,j）<>0Anda（i,j+b）<>0）Then
r=r+（a（i,j）-a（i,j+b））*（a（i,j）-a（i,j+b））
js=js+1
EndIf
Nextj
Nexti
r=r/（2*js）
EndSub
PrivateSubsnr（ByValbAsInteger,rAsSingle）'南北方向半变异函数
DimiAsInteger,jAsInteger,jsAsInteger
r=0
Open"c:
\ckxtgc1\xtgc.txt"ForInputAs#1
DoWhileNotEOF
（1）
Fori=1To8
Forj=1To10
Input#1,a（i,j）
Nextj
Nexti
Loop
Close#1
Fori=1To8-b
Forj=1To10
If（a（i,j）<>0Anda（i+b,j）<>0）Then
r=r+（a（i,j）-a（i+b,j））*（a（i,j）-a（i+b,j））
js=js+1
EndIf
Nextj
Nexti
r=r/（2*js）
EndSub
PrivateSubner（ByValbAsInteger,rAsSingle）'北东方向半变异函数
DimiAsInteger,jAsInteger,jsAsInteger
r=0
Open"c:
\ckxtgc1\xtgc.txt"ForInputAs#1
DoWhileNotEOF
（1）
Fori=1To8
Forj=1To10
Input#1,a（i,j）
Nextj
Nexti
Loop
Close#1
Fori=1+bTo8
Forj=1To10-b
If（a（i,j）<>0Anda（i-b,j+b）<>0）Then
r=r+（a（i,j）-a（i-b,j+b））*（a（i,j）-a（i-b,j+b））
js=js+1
EndIf
Nextj
Nexti
r=r/（2*js）
EndSub
PrivateSubnwr（ByValbAsInteger,rAsSingle）'北西方向半变异函数
DimiAsInteger,jAsInteger,jsAsInteger
r=0
Open"c:
\ckxtgc1\xtgc.txt"ForInputAs#1
DoWhileNotEOF
（1）
Fori=1To8
Forj=1To10
Input#1,a（i,j）
Nextj
Nexti
Loop
Close#1
Fori=1+bTo8
Forj=1+bTo10
If（a（i,j）<>0Anda（i-b,j-b）<>0）Then
r=r+（a（i,j）-a（i-b,j-b））*（a（i,j）-a（i-b,j-b））
js=js+1
EndIf
Nextj
Nexti
r=r/（2*js）
EndSub
PrivateSubcommand3_Click（）'绘图
DimxAsSingle,yAsSingle,iAsInteger,jAsInteger,kAsInteger,hAsInteger
Picture1.Scale（-5,25）-（70,-5）
DrawWidth=1
Picture1.Line（0,0）-（40,0）
Picture1.Print"h"
Picture1.Line（0,0）-（0,20）
Picture1.Print"r"
Picture1.CurrentY=-0.1:
Picture1.CurrentX=-0.4
Picture1.Print"0"
Picture1.CurrentY=0:
Picture1.CurrentX=0
Fori=1To3
Picture1.Line-（10*i,b（i））,RGB（0,0,255）
Nexti
Picture1.CurrentY=0:
Picture1.CurrentX=0
Forj=1To3
Picture1.Line-（10*j,c（j））,RGB（0,255,0）
Nextj
Picture1.CurrentY=0:
Picture1.CurrentX=0
Fork=1To3
Picture1.Line-（10*k,d（k））,RGB（0,255,255）
Nextk
Picture1.CurrentY=0:
Picture1.CurrentX=0
Forh=1To3
Picture1.Line-（10*h,e（h））,RGB（255,0,0）
Nexth
Picture1.Line（45,10）-（50,10）,RGB（0,0,255）
Picture1.Print"东西方向r（h）"
Picture1.Line（45,12）-（50,12）,RGB（0,255,0）
Picture1.Print"南北方向r（h）"
Picture1.Line（45,14）-（50,14）,RGB（0,255,255）
Picture1.Print"北东方向r（h）"
Picture1.Line（45,16）-（50,16）,RGB（255,0,0）
Picture1.Print"北西方向r（h）"
EndSub
6.输入界面及运行结果
图1.输入界面
图2.计算结果及绘图
7.实验小结
实验报告二
1.实验名称
2.题目
某矿原煤生产系统由两部分构成：
破碎筛分系统和储运装车系统。
根据该系统6个月的运行情况，对其进行统计分析并完成如下任务：
（1）分析破碎筛分环节的故障类型，如更换、修补、保养等，给出，各类型的累积故障时间；
（2）分析原煤系统两大环节的各月累积故障时间，绘制柱状图；
（3）假设该原煤系统由破碎筛分系统、储运装车系统和电务班三个环节串联而成，建立串联系统可靠性模型，计算各环节和系统可靠度。
3.
实验步骤
1）各故障累计时间记录
2）
破碎、储运故障累积时间记录
3）建立可靠性模型
假设原煤生产系统年工作时间为330d，每天24小时连续工作。
根据系统累积故障时间和累积故障次数（见下图）
系统累积故障次数统计
根据以上数据，计算破碎和储运环节的可靠度：
（1）破碎环节
R（t）=
==
其中，
为平均故障率，次/小时；m为平均故障间隔时间，h/次，m=
其中，
=
=0.0539次/小时
故，R（t）=
（2）储运环节
=
=0.0125次/小时
R（t）=
4.
实验小结
升级会员 