【建模教程】-Dimine矿体线导入及建模

----宋停云与您分享--------宋停云与您分享----矿体线导入及建模注：以-275米分段为例。1、先在CAD中，把各个分段中的图纸坐标都转换移动到原坐标中。2、在CAD中把矿体线用多段线连接好。（以-275米分段为例）----宋停云与您分享--------宋停云与您分享----3（CAD右键无法截图。4、在迪迈开始下的导入中设置好高程（以-275米分段为例）注：高程设置为-274米，上样标高为本分段上面一米。----宋停云与您分享--------宋停云与您分享----5CAD数据粘贴转换成多段线Continuous线型。----宋停云与您分享--------宋停云与您分享----6、在线编辑中用双线法的相邻连接把矿体线全部连接起来。7、选择矿体线右键直接闭合成一个封闭的闭合多段线。----宋停云与您分享--------宋停云与您分享----8、对矿体界线进行线清理（冗余点检测，检测出来有问题则右键确定清除异常点、线。----宋停云与您分享--------宋停云与您分享----9、对矿体线的多段线需进行加密处理。加密前----宋停云与您分享--------宋停云与您分享----加密后10、用上面相同的步骤方法把其他分段的矿体界线粘贴到迪迈中，并对矿体线进行处理好。注：建模以-235米中段跟-243米分段建模为例。-235米中段跟-243-243-253米分段视图（-253-26、-263与-27-275与-28。----宋停云与您分享--------宋停云与您分享----12进行建模。用同样的操作完成这块矿体的建模9中）----宋停云与您分享--------宋停云与您分享----13、对于上个分段为两个矿体线，而下分段为一个矿体线的，可以在----宋停云与您分享--------宋停云与您分享----12双线连接成一条闭合多段线----宋停云与您分享--------宋停云与您分享----14、用普通连三角网连接法进行建模，再用同样的方法（步骤13、14）对另外的矿体线进行建模。----宋停云与您分享--------宋停云与您分享----15、封闭为面：选择矿体模型最上和最下的矿体线右键封闭为面。16、合并：把矿体线全部隐藏（关闭显示面，然后选中所有面片右键合并（打开显示面。----宋停云与您分享--------宋停云与您分享--------宋停云与您分享--------宋停云与您分享----17、实体有效性检测：选择该实体右键实体有效性检测。注：如果检测有无效边或重合面的话，则找到无效边的位置，把跟其相关的面片删除，右键提取开口线，选择开口线封闭为面，再全部选中面片合并，继续检测修改直至通过有效性检测为止。----宋停云与您分享--------宋停云与您分享----统计建模与r软件第八章答案r软件-4~9章习题】x-c(0.1,0.2,0.9,0.8,0.7,0.7)n-length(x)a1-mean(x);m2-(n-1)/n*var(x)a1-1/(1-a1)-2;a1[1]0.3076923极大似然估计f-function(a2){sum(log(x))+n/(1+a2)}out-uniroot(f,c(0,1));a2-out$root;a2[1]0.2111824.2x-c(rep(5,365),rep(15,245),rep(25,150),rep(35,100),rep(45,70),rep(55,45),rep(65,25)) lamda-length(x)/sum(x);lamda[1]0.054.3取均值即可。x-c(rep(0,17),rep(1,20),rep(2,10),rep(3,2),rep(4,1))mean(x)[1]11个。4.4obj-function(x){f-c(-13+x[1]+((5-x[2])*x[2]-2)*x[2],-29+x[1]+((x[2]+1)*x[2]-14)*x[2]);sum(f^2)}x0-c(0.5,-2)nlm(obj,x0)$minimum[1]48.98425$estimate[1]11.4127791-0.8968052$gradient[1] 1.411401e-08-1.493206e-07$code[1]1$iterations----宋停云与您分享--------宋停云与您分享----[1]164.5x-c(54,67,68,78,70,66,67,70,65,69)t.test(x) #t.test()做单样本正态分布区间估计onesamplet-testdata: xt=35.947,df=9,p-value=4.938e-11alternativehypothesis:truemeanisnotequalto095percentconfidenceinterval:63.158571.6415sampleestimates:meanofx67.4平均脉搏点估计为67.4，95%区间估计为63.158571.6415。t.test(x,alternative=less,mu=72) #t.test()做单样本正态分布单侧区间估计onesamplet-testdata: xt=-2.4534,df=9,p-value=0.01828alternativehypothesis:truemeanislessthan7295percentconfidenceinterval:-inf70.83705sampleestimates:meanofx67.4p0.05，拒绝原假设，平均脉搏低于常人。4.6x-c(140,137,136,140,145,148,140,135,144,141);x[1]140137136140145148140135144141y-c(135,118,115,140,128,131,130,115,131,125);y[1]135118115140128131130115131125t.test(x,y,var.equal=true)twosamplet-testdata: xandyt=4.6287,df=18,p-value=0.0002087alternativehypothesis:truedifferenceinmeansisnotequalto095percentconfidenceinterval:7.5362620.06374----宋停云与您分享--------宋停云与您分享----sampleestimates:meanofxmeanofy140.6 126.8期望差的95%置信区间为7.5362620.06374。4.7x-c(0.143,0.142,0.143,0.137)y-c(0.140,0.142,0.136,0.138,0.140)t.test(x,y,var.equal=true)twosamplet-testdata: xandyt=1.198,df=7,p-value=0.2699alternativehypothesis:truedifferenceinmeansisnotequalto095percentconfidenceinterval:-0.001996351 0.006096351sampleestimates:meanofxmeanofy0.141250.13920期望差的的区间估计为-0.001996351 0.0060963514.84.6var.test(x,y)ftesttocomparetwovariancesdata: xandyf=0.2353,numdf=9,denomdf=9,p-value=0.04229alternativehypothesis:trueratioofvariancesisnotequalto195percentconfidenceinterval:0.058452760.94743902sampleestimates:ratioofvariances0.2353305var.test 可做两样本方差比的估计。此结果可认为方差不等，因此，4.6中，计算期望差时应该采取方差不等的参数。t.test(x,y)welchtwosamplet-testdata: xandyt=4.6287,df=13.014,p-value=0.0004712alternativehypothesis:truedifferenceinmeansisnotequalto0----宋停云与您分享--------宋停云与您分享----95percentconfidenceinterval:7.35971320.240287sampleestimates:meanofxmeanofy140.6 126.8期望差的95%置信区间为7.35971320.240287。4.9x-c(rep(0,7),rep(1,10),rep(2,12),rep(3,8),rep(4,3),rep(5,2))n-length(x)tmp-sd(x)/sqrt(n)*qnorm(1-0.05/2)mean(x)[1]1.904762mean(x)-tmp;mean(x)+tmp[1]1.494041[1]2.315483平均呼唤次数为1.90.95的置信区间为1.49,2,324.10x-c(1067,919,1196,785,1126,936,918,1156,920,948)t.test(x,alternative=greater)onesamplet-testdata: xt=23.9693,df=9,p-value=9.148e-10alternativehypothesis:truemeanisgreaterthan095percentconfidenceinterval:920.8443infsampleestimates:meanofx997.1灯泡平均寿命置信度95%的单侧置信下限为920.84435.1x-c(220,188,162,230,145,160,238,188,247,113,126,245,164,231,256,183,190,158,224,175)t.test(x,mu=225)###双边检验onesamplet-testdata: xt=-3.4783,df=19,p-value=0.002516alternativehypothesis:truemeanisnotequalto22595percentconfidenceinterval:----宋停云与您分享--------宋停云与您分享----172.3827211.9173sampleestimates:meanofx192.15p0.05，拒绝原假设，认为油漆工人的血小板计数与正常成年男子有差异。（原假设：油漆工人的血小板计数与正常成年男子无差异；备择假设：油漆工人的血小板计数与正常成年男子有差异。）t.test(x,mu=225,alternative=less)##单边检验（备择假设）onesamplet-testdata: xt=-3.4783,df=19,p-value=0.001258alternativehypothesis:truemeanislessthan22595percentconfidenceinterval:-inf208.4806sampleestimates:meanofx192.15结论：油漆工人的血小板计数小于正常成年男子。5.2pnorm(1000,mean(x),sd(x))[1]0.5087941x[1]1067 9191196 7851126 936 9181156 920 948pnorm(1000,mean(x),sd(x))[1]0.5087941结论：x=1000的概率为0.509,x大于1000的概率为0.491.5.3a-c(113,120,138,120,100,118,138,123)b-c(138,116,125,136,110,132,130,110)t.test(a,b,paired=true)pairedt-testdata: aandbt=-0.6513,df=7,p-value=0.5357alternativehypothesis:truedifferenceinmeansisnotequalto095percentconfidenceinterval:-15.628898.87889sampleestimates:----宋停云与您分享--------宋停云与您分享----meanofthedifferences-3.375p0.05，接受原假设，即两种治疗方法无差异。5.4（1）正态性w检验：x-c(-0.7,-5.6,2,2.8,0.7,3.5,4,5.8,7.1,-0.5,2.5,-1.6,1.7,3,0.4,4.5,4.6,2.5,6,-1.4)y-c(3.7,6.5,5,5.2,0.8,0.2,0.6,3.4,6.6,-1.1,6,3.8,2,1.6,2,2.2,1.2,3.1,1.7,-2) shapiro.test(x)shapiro-wilknormalitytestdata: xw=0.9699,p-value=0.7527shapiro.test(y)shapiro-wilknormalitytestdata: y【篇二：统计建模与r软件课后习题答案2-5章】用薛毅编的《统计建模与r软件》吧，找不出更好的了……工作环境仍是linux。第二章答案：ex2.1x-c(1,2,3)y-c(4,5,6)e-c(1,1,1)z=2*x+y+ez1=crossprod(x,y)#z1x1x2x%*%yz2=tcrossprod(x,y)#z1x1x2x%o%yz;z1;z2要点：基本的列表赋值方法，内积和外积概念。内积为标量，外积为矩阵。ex2.2a-matrix(1:20,c(4,5));ab-matrix(1:20,ow=4,byrow=true);bc=a+b;c#ab这种写法e=a*b;ef-a[1:3,1:3];fh-matrix(c(1,2,4,5),ow=1);h#h起过渡作用，不规则的数组下标g-b[,h];g要点：矩阵赋值方法。默认是byrow=false,数据按列放置。----宋停云与您分享--------宋停云与您分享----取出部分数据的方法。可以用数组作为数组的下标取出数组元素。ex2.3x-c(rep(1,times=5),rep(2,times=3),rep(3,times=4),rep(4,times=2));x#或者省略times=，如下面的形式x-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));x要点：rep（）的使用方法。rep（a,b）ab次ex2.4n-5;h-array(0,dim=c(n,n))for(iin1:n){for(jin1:n){h[i,j]-1/(i+j-1)}};hgsolve(h);gh的逆矩阵eveigen(h);evh的特征值和特征向量要点：数组初始化；for循环的使用待解决：如何将很长的命令（for循环）用几行打出来再执行？每次想换行的时候一按回车就执行了还没打完的命令...ex2.5studentdata-data.frame(name=c(zhangsan,lisi,wangwu,zhaoliu,dingyi),sex=c(f,m,f,m,f),age=c(14,15,16,14,15),height=c(156,165,157,162,159),weight=c(42,49,41.5,52,45.5));studentdata要点：数据框的使用待解决：ssh登陆linux服务器中文显示乱码。此处用英文代替。ex2.6write.table(studentdata,file=studentdata.txt)#studentdata在工作目录里输出，输出的文件名为studentdata.txt.studentdata_a-read.table(studentdata.txt);studentdata_a#以数据框的形式读取文档studentdata.txt，存入数据框studentdata_a中。write.csv(studentdata_a,studentdata.csv)#把数据框studentdata_a在工作目录里输出，输出的文件名为studentdata.csv,可用excel打开.要点：读写文件。read.table(file)write.table(rdata,file)read.csv(file) write.csv(rdata,file)外部文件，不论是待读入或是要写出的，命令中都得加双引号。ex2.7fun-function(n){----宋停云与您分享--------宋停云与您分享----if(n=0)list(fail=pleaseinputaintegerabove0!)else{repeat{if(n==1)breakelseif(n%%2==0){n-n/2}elsen-3*n+1}list(sucess!)}linuxr2.7.rrsource(2.7.r)，即打开了这个程序脚本。然后就可以执行函数了。输入fun(67)，显示输入fun(-1)，显示$failpleaseinputaintegerabove0!待解决：source(*.r)rr环r文件呢？ok,自己写的第一个r程序~~第二章答案：ex2.1x-c(1,2,3)y-c(4,5,6)e-c(1,1,1)z=2*x+y+ez1=crossprod(x,y)#z1x1x2x%*%yz2=tcrossprod(x,y)#z1x1x2x%o%yz;z1;z2要点：基本的列表赋值方法，内积和外积概念。内积为标量，外积ex2.2a-matrix(1:20,c(4,5));ab-matrix(1:20,ow=4,byrow=true);bc=a+b;c#ab这种写法e=a*b;ef-a[1:3,1:3];fh-matrix(c(1,2,4,5),ow=1);h----宋停云与您分享--------宋停云与您分享----#h起过渡作用，不规则的数组下标g-b[,h];gbyrow=false,ex2.3x-c(rep(1,times=5),rep(2,times=3),rep(3,times=4),rep(4,times=2));xtimes=，如下面的形式x-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));x要点：rep（）的使用方法。rep（a,b）ab次ex2.4n-5;h-array(0,dim=c(n,n))for(iin1:n){for(jin1:n){h[i,j]-1/(i+j-1)}};hgsolve(h);gh的逆矩阵eveigen(h);evh的特征值和特征向量要点：数组初始化；for循环的使用待解决：如何将很长的命令（for循环）用几行打出来再执行？每次想换行的时候一按回车就执行了还没打完的命令...ex2.5studentdata-data.frame(name=c(zhangsan,lisi,wangwu,zhaoliu,dingyi),sex=c(f,m,f,m,f),age=c(14,15,16,14,15),height=c(156,165,157,162,159),weight=c(42,49,41.5,52,45.5));studentdata要点：数据框的使用待解决：ssh登陆linux服务器中文显示乱码。此处用英文代替。ex2.6write.table(studentdata,file=studentdata.txt)#studentdata在工作目录里输出，输出的文件名为studentdata.txt.studentdata_a-read.table(studentdata.txt);studentdata_a#以数据框的形式读取文档studentdata.txt，存入数据框studentdata_a中。write.csv(studentdata_a,studentdata.csv)#把数据框studentdata_a在工作目录里输出，输出的文件名为studentdata.csv,可用excel打开.要点：读写文件。read.table(file)write.table(rdata,file) read.csv(file)write.csv(rdata,file)----宋停云与您分享--------宋停云与您分享----外部文件，不论是待读入或是要写出的，命令中都得加双引号。ex2.7fun-function(n){if(n=0)list(fail=pleaseinputaintegerabove0!)else{rep