




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 今生今世的证据教案17 教案教学设计
- 探索推进速度效应以增强坚硬顶板大采高工作面煤壁稳定性
- 脊髓骶段良性肿瘤护理
- B族链球菌的健康宣教
- 胃神经官能症个案护理
- 低空经济发展趋势分析
- 条纹状角膜炎护理课件
- 蛛网膜粘连的护理查房
- 2025届曲靖第一中学物理高二下期末质量跟踪监视模拟试题含解析
- 过敏性皮疹的健康教育
- GB/T 8312-2002茶咖啡碱测定
- 2023年苏州国发创业投资控股有限公司招聘笔试题库及答案解析
- 通信线路工程施工组织设计方案【实用文档】doc
- 护士注册健康体检表下载【可直接打印版本】
- 预计财务报表编制及分析课件
- 学生集体外出活动备案表
- Q∕SY 1347-2010 石油化工蒸汽透平式压缩机组节能监测方法
- 西门子顺序功能图语言S7-Graph的应用
- 中医治疗室工作制度管理办法
- 提花装造工艺技术培训课程
- 直播传媒公司简介PPT课件(参考)
评论
0/150
提交评论