- 1、本文档共37页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
地震数据处理
《地震资料数据处理》课程设计
总结报告
专业班级:
姓 名:
学 号:
设计时间:
指导老师:
目 录
一 、 设计内容………………………………………………………………
(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 (I.EQ.0)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(I.EQ.0)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)
文档评论(0)