成都理工大学地磁学实验报告.doc

地磁学实验报告 指导老师： 姓名： 学号： 实验一 实验目的：测量地磁场，绘制成图，并观测其变化。初步了解质子旋进磁力仪的工作原理和方法；初步掌握仪器的使用；了解仪器的使用时的注意事项； 实验要求：对实验场地进行测量，并记录数据导入电脑，做出成果图； 实验仪器：质子旋进磁力仪 实验场地：砚胡旁马路 实验点数：共68点 实验要求：在做实验测量数据时，实验员身上不能有带磁性或者可被磁化的金属制品。测量时周围尽量不要聚集人，不能有较多的人员走动或车辆的往来。 实验数据: 测线一（宽）： 分析：从图中，可以看出曲线大体上，还是较为平缓的，但是在测点4,和测点34的时候有较大的异常，分析出现该种异常的原因是测线旁有路灯，路灯是金属材料，因埋地时间较久，可能出来漏电现象，产生较大的磁场，从而影响实验数据；实验同学身上的钥匙手机等也会影像数据的准确性。 测线二（长）： 分析：因地下有管道，管道同样也是金属的，影响实验数据，所以变化差异较大 实验二 实验目的：掌握球体的磁场分布特征 实验要求：绘制时对应的磁场等值线图和剖面折线图#includestdio.h #includemath.h void main() { FILE *fp1; fp1=fopen(t10.txt,w); FILE *fp2; fp2=fopen(t20.txt,w); int x,y; double i,a; double t; double M=10; double r=30; // scanf(%f,i); //scanf(%f,a); i=0; a=45; i=double(3.1415*i/180); a=double(3.1415*a/180); for(y=-50;y51;y++) { for(x=-50;x51;x++) { t=M*( (2*r*r-x*x-y*y)*sin(i)*sin(i)+ (2*x*x-y*y-r*r)*cos(i)*cos(i)*cos(a)*cos(a)+ (2*y*y-x*x-r*r)*cos(i)*cos(i)*sin(a)*sin(a)- 3*x*r*sin(2*i)*cos(a)+ 3*x*y*cos(i)*cos(i)*sin(2*a)- 3*y*r*sin(2*i)*sin(a) ) /sqrt(pow(x*x+y*y+r*r,5)); fprintf(fp1,%d\t%d\t%f\n,x,y,t); if(y==0) fprintf(fp2,%f\t%f\n,x,t); } } } 实验结果： 图1 I=0°，a=45° 图2 I=45°；a=45° 图3 I=90°；a=45° 图4 I=120°；a=45° 图5 I=180°；a=45° 实验三 实验目的： 实验内容：要求： #include stdio.h #include math.h #define N 48 //数据个数 #define pi 3.141596 int main() { //从文件中读取数据保存到数组data中 FILE *fj1=fopen(实验数据.dat,r); FILE *fj2=fopen(插值后的数据.dat,w); FILE *fj3=fopen(圆滑后的数据.dat,w); FILE *fj4=fopen(向上延拓后的数据.dat,w); int i,j; double k,t,p,data[N],data1[N];//data[N]是圆滑数据，data1[N]是延拓数据 p=0; if(fj1==NULL) { printf(\n无法打开此文件！确保文件“实验数据.dat”在目录下\n\n); exit(0); } for(i=0;iN;i++) fscanf(fj1,%lf,data[i]); for(i=0;iN;i++) printf(%f\n,data[i]); printf(\n已导入数组!); //线性插值 k=(data[41]-data[37])/3; for(i=38;i=40;i++) data[i]=data[i-1]+k; for(i=0;iN;i++) fprintf(fj2,%f\n,data[i]); printf(\n已插值!\n); //圆滑 五点二次圆滑 for(i=2;iN-2;i++) { data[i]=(1