基于插值法的煤矿储量估计模型实验 (2).pdf

《基于插值法的煤矿储量估计模型实验 (2).pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西安交通大学数学实验报告实验名称Matlab期末实验报告学院电信学院专业班级计算机26班姓名学号朱增梅2120505129石磊2120505140尹畅畅21205051452013年6月1实验一：一、【实验目的】了解MATLAB软件，熟悉MATLAB的操作与运用，编写程序，解决一些给出的实际问题或数学问题。二、【实验题目】一个对称的地下油库，内部设计如下图所示：横截面为圆，中心位置处的截面半径为3m，上下底处的半径为2m,高为12m，纵截面的两侧是顶点在中心位置的抛物线，试求：（1）油库内油面的深

2、度为h（从底部算起）时，库内的油量容积V（h）（2）设计测量油库油量的标尺，即当油量容积V已知时，算出油的深度h，刻出油量大小，试给出当V=10,20,30……时油的深度。ZyX三、【实验程序】2xaybyc设抛物线方程为：根据题设条件知：y0,x2;y6,x3;y12,x2.2121xyy2解得：363h2y0,12y,ydy2Vhxdy,任取，在上，dVxdy则：0122V12xdy.0程序为：>>symsx>>int('

3、pi*(x^2/(-36)+1/3*x+2)^2',x)得到用求定积分的数值解trapz命令求出油库的容积：>>h0=12;y=0:0.01:h0;>>fy=pi*(-y.^2/36+y/3+2).^2;>>V0=trapz(y,fy);>>fprintf('ThevolumnoftheoilbankisV=%.8f',V0);结果：最大储油量约为271.4335立方米fori=10:10:270c=roots([1/6480,-1/216,0,2/3,4,-i/pi])end得到结果为：10

4、20304050607080900.71171.30921.83722.31822.76563.18783.59083.97904.35591001101201301401501601701804.72405.08575.44315.79786.15156.50606.86287.22387.59091902002102202302402502602707.96638.35278.75349.17259.615810.091410.611611.197111.8880作图显示：312108h642

5、0050100150200250300v四、【实验结果】1.最大储油量约为271.4335立方米。2.1020304050607080900.71171.30921.83722.31822.76563.18783.59083.97904.35591001101201301401501601701804.72405.08575.44315.79786.15156.50606.86287.22387.59091902002102202302402502602707.96638.35278.75349.

6、17259.615810.091410.611611.197111.8880五、【实验总结】通过软件的使用，将作业题通过matlab计算出，使得运算更加迅速和精确，更好的实现一些比较复杂数学问题的解决。数量掌握软件的使用，并适当编程就可以让一些看似运算复杂无法解决的问题交给电脑去处理，能成功实现要求。4实验二：一、【实验目的】本实验的目的是通过对矿区面积的计算,掌握定积分的近似计算方法,对有关数值积分的有关理论和数值计算方法有所了解二、【实验题目】煤矿的储量估计,下表给出了某露天煤矿在平面矩形区域

7、(800m600m)上,在纵横均匀的网格交点处测得的煤层厚度(单位:m)(由于客观原因,有些点无法测量煤层厚度,这里用/标出),其中的每个网格都为(10m8m)的小矩形,试根据这些数据,来估算出该矩形区域煤矿的储藏量(体积)ABCDEFGHIJK1//12.513.517.2/8.814.78.013.0/2///15.618.2136.48.99.211.7/3/1213.513.517.816.913.2////47.512.614.918.717.717.514.713//6.558.9

8、7.812.413.515.717.611.79.69.29.58.66///13.713.616.512.58.79.7//7//8.611.812.511.313.4////三、【实验程序】(1)分行数据拟合:第一行:y=[12.513.517.28.814.78.013.0];x=[160240320480560640720];plot(x,y,'k.','markersize',15)axis([0800020])grid;holdont=0:80:800;p=polyfi