实验三、化向地磁极

一、试验要求

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

二、实验课时：2课时

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

四、化极磁异常公式

五、化极参数

六、化极磁场图

七、源代码

#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;

}

学生（签名）：

年 月 日

指导

教师

评语

成绩评定：

指导教师（签名）：

年 月 日