




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、课程设计程序及成果图件报告姓名:学号:专业:地球物理学指导老师:前言:重磁异常正演是在给定地质体的形状大小、空间位置及物性参数的条件下,求它在它外部空间任意点上的物理场值。电子计算机的应用使正演方法进入了一个新的阶段,它可以对任意形状物体进行计算;对物性“非均匀”的地质体,也可以通过将其分解成若干个“均匀”物性的物体的方法和分区块均匀的数值方法来计算,使正演方法可以达到相当高的精度。以下是我对截面为多边形的水平柱体模型的重磁异常的正演程序。C语言代码: #include stdio.h#include math.h#define pi 3.14159 #define n 101 /*测点数*/
2、#define N 360 /*拟合多边形边数*/double fun(double a,double b) /*关于atan(b/a)的讨论*/ doublethet; if(fabs(a)1.0e-15) if(b/a)0.0) if(a0.0&b0.0)thet=atan(b/a);else if(a0.0&b0.0)thet=atan(b/a)-pi; else if(b/a)0.0)thet=atan(b/a)+pi; else if(fabs(a)=1.0e-15) if(b=0.0)thet=pi/2; return(thet); void main() FILE *fp;fp=
3、fopen(输出.txt,w+); double xN=0.0,zN=0.0,Pn=0.0,Qn=0.0; deltagn=0.0,deltaxn=0.0,deltazn=0.0; double M=2000.0,G=6.67e-5,sigma=2.67,r=200.0,x0=1000.0,z0=1000.0; inti,k; doublexk101,x1N,z1N; double gn=0.0,alpha=(2*pi)/N; doubleMx=M*sin(pi/4),Mz=M*cos(pi/4);for(i=0;iN;i+) x1i=x0+r*sin(alpha*i);z1i=z0-r*co
4、s(alpha*i); for(k=0;kn;k+) xkk=20.0*k;for(i=0;iN-1;i+) xi=x1i-xkk;xi+1=x1i+1-xkk;zi=z1i;zi+1=z1i+1;Qk+=(zi+1-zi)*(xi-xi+1)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1)-(1/2)*pow(zi+1-zi),2)*log(pow(xi+1,2)+pow(zi+1,2)/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2); Pk+=(pow(zi+
5、1-zi),2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1)-(1/2)*(zi+1-zi)*(xi-xi+1)*log(pow(xi+1,2)+pow(zi+1,2)/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2);deltaxk=(Mx*Pk+Mz*Qk)/(2*pi);deltazk=(Mx*Qk-Mz*Pk)/(2*pi);gk+=(1/2)*(zi+1-zi)*(xi*zi+1-xi+1*zi)*log(pow(xi+1,2)+pow(zi+1,2)
6、/(pow(xi,2)+pow(zi,2)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)-(xi+1-xi)*(xi*zi+1-xi+1*zi)/(pow(zi+1-zi),2)+pow(xi-xi+1),2)*(fun(xi,zi)-fun(xi+1,zi+1);deltagk=-2*G*sigma*gk; for(i=0;in;i+) printf(%lf %lf %lfn,deltagi,deltaxi,deltazi); fprintf(fp,%lf %lf %lfn,deltagi,deltaxi,deltazi); fclose(fp);数据结果:G X Z0.
7、011686 7.067300 7.0673500.011924 7.062859 7.3543370.012167 7.049233 7.6497860.012415 7.025691 7.9536400.012667 6.991471 8.2658010.012925 6.945778 8.5861210.013186 6.887785 8.9144030.013453 6.816638 9.2503890.013724 6.731454 9.5937600.014000 6.631327 9.9441260.014279 6.515330 10.3010240.014563 6.3825
8、23 10.6639110.014851 6.231952 11.0321550.015142 6.062661 11.4050340.015437 5.873697 11.7817270.015736 5.664118 12.1613120.016037 5.433000 12.5427580.016340 5.179451 12.9249210.016646 4.902621 13.3065440.016954 4.601714 13.6862490.017263 4.276001 14.0625410.017572 3.924832 14.4338060.017882 3.547657
9、14.7983110.018192 3.144038 15.1542070.018501 2.713665 15.4995370.018808 2.256373 15.8322420.019113 1.772162 16.1501690.019416 1.261207 16.4510860.019714 0.723882 16.7326920.020008 0.160769 16.9926400.020297 -0.427326 17.2285520.020579 -1.039362 17.4380450.020855 -1.674053 17.6187530.021122 -2.329861
10、 17.7683580.021381 -3.004993 17.8846130.021629 -3.697400 17.9653790.021867 -4.404783 18.0086530.022094 -5.124601 18.0126010.022308 -5.854081 17.9755920.022508 -6.590239 17.8962250.022694 -7.329901 17.7733630.022865 -8.070195 17.6056940.023020 -8.806829 17.3935070.023159 -9.536516 17.1363240.023280 -
11、10.255658 16.8342670.023384 -10.960629 16.4878160.023469 -11.647843 16.0978720.023536 -12.313736 15.6656380.023584 -12.954885 15.1927260.023612 -13.568017 14.6810740.023621 -14.150072 14.1329370.023611 -14.698194 13.5508990.023581 -15.209814 12.9377990.023532 -15.682675 12.2967040.023464 -16.114835
12、11.6308850.023377 -16.504704 10.9437840.023272 -16.851025 10.2389090.023150 -17.152950 9.5198990.023010 -17.409984 8.7903630.022854 -17.622006 8.0538950.022682 -17.789493 7.3137840.022494 -17.912159 6.5743190.022293 -17.991316 5.8383710.022078 -18.028104 5.1091140.021851 -18.023924 4.3895290.021612
13、-17.980408 3.6823880.021362 -17.899392 2.9902320.021103 -17.782879 2.3153580.020835 -17.633012 1.6598150.020559 -17.452035 1.0253920.020276 -17.242270 0.4136290.019987 -17.006083 -0.1741910.019692 -16.745858 -0.7370270.019393 -16.463975 -1.2740740.019091 -16.162781 -1.7847510.018785 -15.844576 -2.26
14、86850.018478 -15.511596 -2.7257010.018168 -15.165991 -3.1558000.017858 -14.809824 -3.5591480.017548 -14.445052 -3.9360540.017238 -14.073522 -4.2869580.016929 -13.696968 -4.6124110.016622 -13.317006 -4.9130620.016316 -12.935131 -5.1896390.016012 -12.552720 -5.4429400.015711 -12.171032 -5.6738150.015413 -11.791209 -5.8831570.015118 -11.414283 -6.0718890.014826 -11.041176 -6.2409520.014539 -10.672710 -6.3913010.014255 -10.309607 -6.5238920.013975 -9.952497 -6.6396770.013700 -9.601925 -6.7395990.013429 -9.258353 -6.8245820.013162 -8.922171 -6.8955340.012901 -8.593699 -6.9533370.012643 -
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 自助美甲店合作合同范本
- 高空作业安全打协议合同
- 消毒用品捐献协议书模板
- 浴场会所托管合同协议书
- 离婚前三年的财产协议书
- 物业零星工程施工协议书
- 自媒体运营团队合同范本
- 第三方协议护理网签合同
- 续签的合同上没竞业协议
- 糖果批发转让协议书模板
- GA/T 1323-2016基于荧光聚合物传感技术的痕量炸药探测仪通用技术要求
- 2023年苏州国发创业投资控股有限公司招聘笔试题库及答案解析
- 护士注册健康体检表下载【可直接打印版本】
- 高中历史《第一次工业革命》说课课件
- 学生集体外出活动备案表
- SH3904-2022年石油化工建设工程项目竣工验收规定
- 叉车检验检测报告
- DNF装备代码大全
- 基于Qt的俄罗斯方块的设计(共25页)
- 古建筑木构件油漆彩绘地仗施工技术分析
- 食堂投诉处理方案
评论
0/150
提交评论