没有合适的资源?快使用搜索试试~ 我知道了~
毕业论文-地震数据处理-.doc
0 下载量 76 浏览量
2023-06-30
05:04:46
上传
评论
收藏 481KB DOC 举报
温馨提示
试读
37页
毕业论文-地震数据处理-.doc
资源推荐
资源详情
资源评论
地震资料数据处理方法》课程设计报告
1
《地震资料数据处理》课程设计
总结报告
专业班级:
姓 名:
学 号:
设计时间:
指导老师:
地震资料数据处理方法》课程设计报告
2
目 录
一 、 设计内容………………………………………………………………
(1)褶积滤波………………………………………………
(2)快变滤波………………………………………………
(3)褶积滤波与快变滤波的比较…………………………
(4)设计高通滤波因子……………………………………
(5)频谱分析………………………………………………
(6)分析补零对振幅谱的影响……………………………
(7)线性褶积与循环褶积…………………………………
(8)最小平方反滤波………………………………………
(9)零相位转换……………………………………………
(10)最小相位转换…………………………………………
(11)静校正…………………………………………………
二、附录…………………………………………………………………………
(1)附录 1:相关程序……………………………………
(2)附录 2:相关图件……………………………………
地震资料数据处理方法》课程设计报告
3
【附录 1:有关程序】
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')
地震资料数据处理方法》课程设计报告
4
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=1,I
DO 2 I2=1,J
II=I1+I2-1
2 C(II)=C(II)+A(I1)*B(I2)*0.002
RETURN
END
2. 快变滤波
CCCCCCCCCCCCCCC 频率滤波 CCCCCCCCCCCCCCCCCCCC
PROGRAM MAIN
PARAMETER (PI=3.141592654)
REAL H_LOW(1:200),H_BAND(1:200),Y_LOW(1:200),Y_BAND(1:200)
REAL X(1:200)
INTEGER I,NFFT,K,NP
COMPLEX C(1:200),TEMP(1:200)
REAL DT,DF,F,F1,F2
F=70.0
F1=10.0
F2=80.0
DT=0.002
OPEN(1,FILE='INPUT1.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
READ(1,*)(X(I), I=1,100)
NP=100
K=LOG(100.0)/LOG(2.0)
IF(2**K.LT.100)THEN
K=K+1
NFFT=2**K
END IF
DF=1.0/(NFFT*DT)
DO 10 I=101,128
X(I)=0.0
10 CONTINUE
地震资料数据处理方法》课程设计报告
5
CCCCCCCCCCC 数据变换成复数形式 CCCCCCCCCCC
DO 20 I=1,NFFT
20 C(I)=CMPLX(X(I),0.0)
CCCCCCCCCCC 向频率域转换 CCCCCCCCCCCCCCCCC
CALL FFT(NFFT,C,1.0)
CCCCCCCCCCC 低通滤波因子的设计 CCCCCCCCCCC
DO 30 I=1,NFFT/2
IF(I*DF.LE.F)THEN
H_LOW(I)=1.0
ELSE
H_LOW(I)=0.0
END IF
30 CONTINUE
CCCCCCCCC 使滤波器理想化(对称)CCCCCCCCCC
DO 40 I=NFFT/2+1,NFFT
F=I*DF
H_LOW(I)=H_LOW(NFFT-I+1)
40 CONTINUE
CCCCCCCCCCCCCCC 实现滤波 CCCCCCCCCCCCCCCCC
DO 50 I=1 , NFFT
50 TEMP(I)=C(I)*H_LOW(I)
CCCCCCCCCCCCCCC 向时间域变换 CCCCCCCCCCCCC
CALL FFT(NFFT,TEMP,-1.0)
DO 60 I=1 , NFFT
60 Y_LOW(I)=REAL(TEMP(I))
OPEN(2,FILE='LOWPASS.DAT',FORM='FORMATTED',STATUS='UNKNOWN')
WRITE(2,*)(Y_LOW(I),I=1,NFFT)
CLOSE(2)
CCCCCCCCCCC 带通滤波因子的设计 CCCCCCCCCCC
DO 70 I=1,NFFT/2
IF((I*DF.GE.F1).AND.(I*DF.LE.F2))THEN
H_BAND(I)=1.0
ELSE
H_BAND(I)=0.0
END IF
70 CONTINUE
CCCCCCCCC 使滤波器理想化(对称)CCCCCCCCCC
DO 80 I=NFFT/2+1,NFFT
F=I*DF
H_LOW(I)=H_BAND(NFFT-I+1)
80 CONTINUE
CCCCCCCCCCCCCCC 实现滤波 CCCCCCCCCCCCCCCCC
DO 90 I=1 , NFFT
90 TEMP(I)=C(I)*H_BAND(I)
剩余36页未读,继续阅读
资源评论
Mmnnnbb123
- 粉丝: 686
- 资源: 8万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功