1、实验变异函数的计算机求解

1、 釆矿系统工程实验报告一、实验名称实验变异函数的计算机求解二、程序流程图三、实验思路1、明确题目中要求实现的功能，即实现数据文件的读取，计算各个方向上不同步距间的半变异函数值并将结束显示并输出到文件里；2、绘制程序流程图，根据程序流程图进行程序的设计；3、建立数据文件data.txt；4、编写读取文件的代码，实现数据文件的读取功能；5、通过用户输入要计算的内容，即要计算哪些方向和步距；编写每个方向上的函数：定义东西方向上半变异函数rew,南北方向上半变异函数isn,北东方向上半变异函数ine和北西方向上半变异函数niw；6、显示并将结果输出到result.txt文件中。四、源代码存includ

2、e存includeusingnamespacestd;#defiiiem8definen10intamn;intmam()ffloatrew(iiit);floatisn(int);floatrne(int);floatniw(iiit);ifstreaminfile(Hdata.txt,)jif(!infile)定义行数定义列数定义东西方向上半变异函数rew定义南北方向上半变异函数rsil定义北东方向上半变异函数me定义北西方向上半变异函数rawcerropenerror?endl;exit(l);for(inti=0;im;i+)for(intj=0;jn;j+)ofstreamoutfi

3、le(Mresult.txfios:out);if(!outfile)cerrHopenerror!Mendl;exit(l);cout请选择要计算的数据：”endl;coutMi.东西方向HH”2南北方向yv”nn3.北东方向yy”4北西方向nnHH5.所有方向1endl;intdir;方向序号inthl,bl,112.b2;步距while(l)coutH请输入序号:,endl;cindir;if(dir5)coutH输入错误!,fendl;contmue;cout请输入起始步距(米)：endl;cinhl;cout请输入终止步距(米)：endl;cinh2;bl=hl/100.b2=h2/

4、100;switch(dir)case1:for(i=bl;i=b2;i+)intj=i*100;coutH东西方向米：,frew(i)endl;outfileM东西方向”vvjvv”米：Hrew(i)endl;break;case2:for(i=bl;i=b2;i+)intj=i*100;coutH南北方向,VjH,rrsn(i)endl;outfileH南北方向,jHt|C:：,rrsn(i)endl;break;case3:for(i=bl;i=b2;i-H-)intj=i*100;coutn北东方向,jH米:,rrne(i)endl;outfileH北东方向njH米:,rrne(i)e

5、ndl;break;case4:for(i=bl;i=b2;i-H-)intj=i*100;coutn北西方向jH米:,rrnw(i)endl;outfileH北西方向,jH米:,rrnw(i)endl;break;case5:for(i=bl;i=b2;i+)intj=i*100;coutn东西方向jH米：nrew(i)endl;coutH南北方向njH,rrsn(i)endl;coutn北东方向jH米:,rrne(i)endl;coutH北西方向HjH米：Hrnwendl;,rrew(i)endl;,rrsn(i)endl;Hrne(i)endl;Hrnw(i)endl;outfileH东

6、西方向HjH米outfile,r南北方向HjH米outfileH北东方向米outfileH北西方向HjH米break;coutM是否继续(y/n)?Mendl;if(getchar()=y)continue;elseif(getchai()=,n,)break;return0;floatrew(iiitk)东西方向上半变异函数lewfloatb=O.t=O;for(inti=0;im;i+)for(intj=0;jn-k;j+)if(aij!=0&ai|j+k!=0)b+;t+=(aij-ai|j+k)*何i|j-ai|j+k);return(t/2/b);floatrsn(intk)南北方向

7、上半变异函数rsnfloatb=O.t=O;Mintj=0Jnj+)for(inti=0;im-k;i+)if(aij!=0&ai+k|J!=0)b+;t+=(aij-ai+k|j)*(ai|j-ai+kIj);return(t/2/b);floatine(intk)北东方向上半变异函数mefloatb=0.t=0;for(inti=0;im-k;i+)for(intj=0;jn-k;j+)if(a1j!=0&ai+kj+k!=0)b+;t+=(aijai+kj+k)qaij卜ai+kj+k);return(t/2/b);floatrnw(mtk)北西方向上半变异函数inwfloatb=O.t=O;for(inti=0;im-k;i+)for(intj=O;ju+k;j+)if(a1j!=0&ai+kj-k!=0)b+;t+=(aij-ai+kj-k)*(aij-ai+kj-k);return(t/2/b);3北东方向4北西方向5所有方向CD|回五、运行结果请逵择票计算的数蜃“亠亠券酉方向2南北方向请输入序号：请输入起始步距（米）：100请输入终止步距（米）：300南北方向佃0米：4.25南:匕方向200米:8.22222商北方向300米：10-