地震数据处理课程设计-浙江大学-第一文库 | 海量文档资源下载与分享平台

免费试读已结束，剩余 28 页请下载文档后查看

10 金币

下载此文档

/ 38

下载提示文本预览

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

地震资料数据处理方法》课程设计报告《地震资料数据处理》课程设计总结报告专业班级：姓名：学号：设计时间：指导老师：目录一、设计内容………………………………………………………………（1）褶积滤波………………………………………………（2）快变滤波………………………………………………（3）褶积滤波与快变滤波的比较…………………………（4）设计高通滤波因子……………………………………（5）频谱分析………………………………………………（6）分析补零对振幅谱的影响……………………………（7）线性褶积与循环褶积…………………………………（8）最小平方反滤波………………………………………（9）零相位转换……………………………………………（10）最小相位转换…………………………………………（11）静校正…………………………………………………二、附录…………………………………………………………………………（1）附录1：相关程序……………………………………（2）附录2：相关图件……………………………………【附录1：有关程序】褶积滤波CCCCCCCCCCCCCCCCC褶积滤波CCCCCCCCCCCCCCCCCPROGRAMMAINDIMENSIONX(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200)PARAMETER(PI=3.141592654)CCCCCCCCH1是低通滤波因子，H2为带通滤波因子CCCCCCREALX,H1,H2,Y_LOW,Y_BANDREALdt,F,F1,F2INTEGERIdt=0.002F=70.0F1=10.0F2=80.0OPEN(1,FILE='INPUT1.DAT',FORM='FORMATTED',STATUS='UNKNOWN')READ(1,*)(X(I),I=1,100)CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCCDO10I=-50,50IF(I.EQ.0)THENH1(I)=2*F*PI/PIELSEH1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt)ENDIF10CONTINUECCCCCCCCCCCCCCCC输出低通滤波因子CCCCCCCCCCCCCCCCOPEN(2,FILE='H1_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(2,*)(H1(I),I=-50,50)CLOSE(2)CALLCON(X,H1,Y_LOW,100,101,200)CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCOPEN(3,FILE='Y_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(3,*)(Y_LOW(I),I=51,150)CLOSE(3)CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCCDO20I=-50,50IF(I.EQ.0)THENH2(I)=140ELSEH2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt)ENDIF20CONTINUECCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCCOPEN(4,FILE='H2_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(4,*)(H2(I),I=-50,50)CLOSE(4)CALLCON(X,H2,Y_BAND,100,101,200)CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCCCOPEN(5,FILE='Y_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(5,*)(Y_BAND(I),I=51,150)CLOSE(5)ENDCCCCCCCCCCCCCCCCCCCCC褶积函数CCCCCCCCCCCCCCCCCCCCSUBROUTINECON(A,B,C,I,J,K)DIMENSIONA(I),B(J),C(K)DO1K1=1,K1C(K1)=0.0DO2I1=1,IDO2I2=1,JII=I1+I2-12C(II)=C(I