地震數(shù)據(jù)處理課程設(shè)計(jì)(報(bào)告)

《地震數(shù)據(jù)處理課程設(shè)計(jì)(報(bào)告)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、中國石油大學(xué)（華東）《地震資料數(shù)據(jù)處理方法》課程設(shè)計(jì)報(bào)告《地震資料數(shù)據(jù)處理》課程設(shè)計(jì)總結(jié)報(bào)告專業(yè)班級(jí)：姓名：學(xué)號(hào)：設(shè)計(jì)時(shí)間：指導(dǎo)老師：2011年5月30日37中國石油大學(xué)（華東）《地震資料數(shù)據(jù)處理方法》課程設(shè)計(jì)報(bào)告目錄一、設(shè)計(jì)內(nèi)容………………………………………………………………（1）褶積濾波………………………………………………（2）快變?yōu)V波………………………………………………（3）褶積濾波與快變?yōu)V波的比較…………………………（4）設(shè)計(jì)高通濾波因子……………………………………（5）頻譜分析……………………

2、…………………………（6）分析補(bǔ)零對(duì)振幅譜的影響……………………………（7）線性褶積與循環(huán)褶積…………………………………（8）最小平方反濾波………………………………………（9）零相位轉(zhuǎn)換……………………………………………（10）最小相位轉(zhuǎn)換…………………………………………（11）靜校正…………………………………………………二、附錄…………………………………………………………………………（1）附錄1：相關(guān)程序……………………………………（2）附錄2：相關(guān)圖件……………………………………37中國石油大學(xué)（華

3、東）《地震資料數(shù)據(jù)處理方法》課程設(shè)計(jì)報(bào)告【附錄1：有關(guān)程序】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.

4、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輸

5、出低通濾波因子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輸出濾波后的數(shù)據(jù)CCCCCCCCCCCCCCCCOPEN(3,FILE='Y_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(3,*

6、)(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',FO

7、RM='FORMATTED',STATUS='UNKNOWN')37中國石油大學(xué)（華東）《地震資料數(shù)據(jù)處理方法》課程設(shè)計(jì)報(bào)告WRITE(4,*)(H2(I),I=-50,50)CLOSE(4)CALLCON(X,H2,Y_BAND,100,101,200)CCCCCCCCCCCCCCCC輸出濾波后的數(shù)據(jù)CCCCCCCCCCCCCCCCCOPEN(5,FILE='Y_BAND.DAT',FORM='FORMATTED',STATUS='UNKNOWN')WRITE(5,*)(Y_BAND(I),I=51,

8、150)CLOSE(5)ENDCCCCCCCCCCCCCCCCCCCCC褶積函數(shù)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(II)+A(I1)*B(I2)*0.002RETURNEND1.快變?yōu)V波CCCCCCCCCCCCCCC頻率濾波CC