地震数据处理课程设计.doc

《地震资料数据处理》课程设计 总结报告 专业班级： 姓 名： 学 号： 设计时间： 指导老师： 目 录 一 、 设计内容……………………………………………………………… （1）褶积滤波……………………………………………… （2）快变滤波……………………………………………… （3）褶积滤波与快变滤波的比较………………………… （4）设计高通滤波因子…………………………………… （5）频谱分析……………………………………………… （6）分析补零对振幅谱的影响…………………………… （7）线性褶积与循环褶积………………………………… （8）最小平方反滤波……………………………………… （9）零相位转换…………………………………………… （10）最小相位转换………………………………………… （11）静校正………………………………………………… 二、附录………………………………………………………………………… （1）附录1：相关程序…………………………………… （2）附录2：相关图件…………………………………… 【附录1：有关程序】 褶积滤波 CCCCCCCCCCCCCCCCC 褶积滤波 CCCCCCCCCCCCCCCCC PROGRAM MAIN DIMENSION X(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200) PARAMETER (PI=3.141592654) CCCCCCCC H1是低通滤波因子，H2为带通滤波因子CCCCCC REAL X,H1,H2,Y_LOW,Y_BAND REAL dt,F,F1,F2 INTEGER I dt=0.002 F=70.0 F1=10.0 F2=80.0 OPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=UNKNOWN) READ(1,*)(X(I),I=1,100) CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCC DO 10 I=-50,50 IF ()THEN H1(I)=2*F*PI/PI ELSE H1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt) END IF 10 CONTINUE CCCCCCCCCCCCCCCC输出低通滤波因子CCCCCCCCCCCCCCCC OPEN(2,FILE=H1_LOW.DAT,FORM=FORMATTED,STATUS=UNKNOWN) WRITE(2,*)(H1(I),I=-50,50) CLOSE(2) CALL CON(X,H1,Y_LOW,100,101,200) CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCC OPEN(3,FILE=Y_LOW.DAT,FORM=FORMATTED,STATUS=UNKNOWN) WRITE(3,*)(Y_LOW(I),I=51,150) CLOSE(3) CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCC DO 20 I=-50,50 IF()THEN H2(I)=140 ELSE H2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt) END IF 20 CONTINUE CCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCC OPEN(4,FILE=H2_BAND.DAT,FORM=FORMATTED,STATUS=UNKNOWN) WRITE(4,*)(H2(I),I=-50,50) CLOSE(4) CALL CON(X,H2,Y_BAND,100,101,200) CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCC OPEN(5,FILE=Y_BAND.DAT,FORM=FORMATTED,STATUS=UNKNOWN) WRITE(5,*)(Y_BAND(I), I=51,150) CLOSE(5) END CCCCCCCCCCCCCCCCCCCCC褶积函数CCCCCCCCCCCCCCCCCCCC SUBROUTINE CON(A,B,C,I,J,K) DIMENSION A(I),B(J),C(K) DO 1 K1=1,K 1 C(K1)=0.0 DO 2 I1