实验三、化向地磁极

2016年2月2日 305 次阅读 0 条评论 0 人点赞

 实验三、化向地磁极

一、试验要求

根据给定斜交磁化剖面异常数据,计算其化极磁异常

二、实验课时:2课时

三、实验成果提交:化极磁异常图件

四、化极磁异常公式

1234567
a i0.42860.17490.11030.08130.06450.05360.0458

五、化极参数

加权次数

磁倾角(I

磁偏角(A

N=5

45

45

六、化极磁场图

七、源代码

#include
"stdafx.h"

#define
_CRT_SECURE_NO_WARNINGS

#include<stdio.h>

#include<stdlib.h>

#include<math.h>

#define
file2
" file .dat"

int main( )

{

    FILE *fp1,*fp2; float s1[200]; float s2[200];

     int i;

    double Z, I=3.14/4, A=3.14/4, ss, H;

    if ((fopen_s(&fp1,"file1", "rb+")) == NULL) {

        printf("Cannot open file strike any key exit!");

    }//非标准C库函数


    for (i = 0; i<200; i++)

    {

        s1[i] = 0;


    }

    for (i = 0; i<176; i++)

    {

        fscanf(fp1, "%f %f\n", &s2[i], &s1[i + 5]);

        s1[i] = s1[i] + 1;

        printf("%f\t%f\n", s2[i], s1[i + 5]);

    }

    fclose(fp1);


    fopen_s(&fp2, file2, "a+");


    for (i = 5; i <= 180; i++)

    {

        ss = 0.4286*(s1[i + 1] - s1[i - 1]) + 0.1749*(s1[i + 2] - s1[i - 2]) + 0.1103*(s1[i + 3] - s1[i - 3]) + 0.0813*(s1[i + 4] - s1[i - 4]) + 0.0645*(s1[i + 5] - s1[i - 5]);

        Z = 1 / (sin(I)*sin(I) + cos(I)*cos(I)*cos(A)*cos(A))*(s1[i] * (sin(I)*sin(I) - cos(I)*cos(I)*cos(A)*cos(A)) - sin(I)*sin(I)*cos(A) / 3.141592653*ss*3.14);//-sin(I)*sin(I)*cos(A)/3.141592653*ss);

        H = 1 / (sin(I)*sin(I) + cos(I)*cos(I)*cos(A)*cos(A))*(s1[i] * sin(I)*sin(I)*cos(A) + (sin(I)*sin(I) - cos(I)*cos(I)*cos(A)*cos(A)) / 3.141592653*ss*3.14);

        fprintf(fp2, "%f\t%f\t%f\n", s2[i - 5], Z, H);

    }//计算并输出

    fclose(fp2);

    return 0;

}

学生实验 心得

这三次试验分别做了三度体球体和二度体水平圆柱体和斜交磁化厚板的磁异常正演和化向地磁极试验。了解并掌握了简单三度体和二度体简单模型的基本正演,并绘制了平面等值线图和主剖面图等正演图件。进一步提高了c语言编程能力,了解了C语言编程在专业课中的应用,了解了基本的正演方法。了解了简单二度三度体正演磁异常图线的基本形状,加深了对理论知识的理解。

学生(签名):

年 月 日

指导

教师

评语

成绩评定:

指导教师(签名):

年 月 日


菜鸟

文章评论(0)