地震实习报告（共5篇） .docxVIP

地震实习报告（共5篇） 第1篇：地震勘探实_报告目录第一章 绪 论 fclose(fp5); conv(x,m,h1,n,y_band,l); fp3=fopen(\synthesisdata_bandpa.DAT\,\w\); for(i=50;ifpar=fopen(\lowpafilter.par\,\r\); fscanf(fpar,\%f%f\, np=100; k=log(np)/log(2); if(nppow(2,k))k=k+1; nfft=pow(2,k); dt=0.002; df=1.0/(nfft*dt); xr=(double*)calloc(nfft,sizeof(double)); xi=(double*)calloc(nfft,sizeof(double)); H=(float*)calloc(nfft,sizeof(float)); // read x(n) fp1=fopen(\INPUT1.DAT\,\r\); for(i=0;i16fscanf(fp1,\%f\, xr[i]=z; } fclose(fp1); //补零至 128 位 for(i=100;i=fc1)H[i]=1.0; else H[i]=0.0; } //滤波器对称 for(i=nfft/2;i17//获取高通截至频率fpar=fopen(\bandpa.par\,\r\); fscanf(fpar,\%f%f\, fp1=fopen(\INPUT1.DAT\,\r\); for(i=0;i=fc1)H[i]=1.0; else H[i]=0.0; } for(i=nfft/2+1;i18 }3. 褶积滤波与递归滤波 褶积滤波# include \stdio.h\ # include \math.h\ # include \stdlib.h\ # include \conv.c\ # include \fft.c\ # define PI 3.1415926 main() { void conv(); float x[50],h[20],y[69],hreverse[20],hzero[39],yreverse[69]; float dt=0.002; int i,m,n,l,p,q; FILE *fp1,*fp2,*fp3,*fp4,*fp5,*fp6; m=50;n=20; l=m+n-1; //read x(n) fp1=fopen(\INPUT3.DAT\,\r\); for(i=0;iconv(x,m,h,n,y,l);//非零相位褶积滤波 fp3=fopen(\con_filter.dat\,\w\); for(i=0;i19 p=n+n-1; q=m+p-1; //构造零相位滤波因子 conv(h,n,hreverse,n,hzero,p); fp6=fopen(\zerophasefilterfactor.dat\,\w\); for(i=0;ipow(2,k))k=k+1; nfft=pow(2,k); h=(double*)malloc(nfft*sizeof(double)); hi=(double*)malloc(nfft*sizeof(double)); for(i=-50;i22 { if(i==0) h[i+50]=1.0/dt-60; elseh[i+50]=-sin(2*PI*30.0*i*dt)/(PI*i*dt); } for(i=100;ifp1=fopen(\timedomain.dat\,\w\); for(i=0;i# include \stdio.h\ # include \math.h\ # include \dft.c\ # define N 60 # define PI 3.1415926 # define dt 0.004 main() { float x[N],xr[N],xi[N],w[N],wr[N],wi[N],z; int i,k;FILE *fp,*fp1,*fp2,*fp3;23 fp3=fopen(\sin.dat\,\w\); for(i=0;ifp=fopen(\WAVE.DAT\,\r\); for(i=0;i# include \stdio.h\ # include \math.h\ # include \stdlib.h\ # include \fft.c\ # define N 60 # define PI 3.1415926 # define dt 0.004 main() { void fft(); double *x1,*xi1,*x2,*xi2; float z; int i,k,nfft;FI