资源描述：
《seismiclab地震数据处理包帮助文档归纳翻译》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、SeismicLab地震数据处理包帮助文档归纳翻译一．源代码Codes（一）BP-Filter——Applyaband-passfiltertoagroupoftraces.freq.domainfiltering（带通-频率域滤波）function[o]=bp_filter(d,dt,f1,f2,f3,f4);函数注释：INd:data(columnsaretraces)%数据dt:samplingintervalinsec%采样间隔（秒）f1:freq.inHz%最小频率f2:freq.inHz%
2、f2-f3表示频带宽度f3:freq.inHz%f2-f3表示频带宽度f4:freq.inHz%最大频率
3、___________
4、/Amplitudespectrum振幅谱
5、/%镶边法用于避免吉布斯现象造成的影响
6、/
7、------------------------------->f1f2f3f4OUTo:output(columnsaretraces)%输出数据（地震道）Example:d=linear_events;dout=bp_filter(d,0.004,1,3,30,40);wigb
8、([d,dout]);%wigb画图的函数，[d,dout]表示横向画图在一张figure上。（二）decon——提供的反褶积方法有8种timedomaindeconvolutionandwaveletestimation（时域反卷积与子波评价）【一】kolmog.mKOLMOG:Kolmogoroffspectralfactorization。%柯尔莫哥洛夫（Kolmogoroff）谱分解。Givenawavelet,thisfunctionretrievestheminimumphasewavel
9、etusingKolmogorofffactorization.Iftheinputisatrace,thespectralfactorizationisappliedtotheautocorrelationaftersmoothing.%给定一个小波，该函数检索采用Kolmogoroff分解的最小相位子波。如果输入是一个地震道，在平滑后谱分解再应用于自相关。Function[w_min]=kolmog(w,type_of_input,L);函数注释：INw:awaveletofarbitraryph
10、aseiftype_of_input='w'%任意相位的子波oraseismictraceiftype_of_input='t'%一个地震道L:lenghtofwaveletiftype_of_input='t'%子波长度OUTw_min:aminphasewavelet%输出一个最小相位子波Example:w=[1;2;-1;0;0;0;0]wmin=kolmog(w,'w');subplot(221);stem(w);subplot(222);stem(wmin);【二】ls_inv_filter
11、.mLS_INV_FILTER:Least-squaresinversefilterofawavelet.%一种子波的最小二乘法反滤波。function[f,o]=ls_inv_filter(w,NF,Lag,mu);函数注释：INw:thewavelet%子波NF:lenghtoftheinversefilter%反滤波的长度Lag:thepositionofthespikeinthedesiredoutput%在期望输出中尖峰的位置Lag=1forminimumphasewavelets%最小相位
12、子波mu:prewhiteningaspercentageofthezerolagautocorrelationcoefficientofthewavelet。%依据小波的零延迟自相关系数的百分比白噪化。OUTf:thefilter%滤波器o:theouputorconvolutionofthefilterwiththewavelet%输出或小波与滤波器的卷积。Example:w=[4;-3;1;0;0;0];[f,o]=ls_inv_filter(w,10,1,0.1);subplot(331);s
13、tem(w);title('wavelet');subplot(332);stem(f);title('filter');subplot(333);stem(o);title('convoffilterwithwavelet')%subplot是固定图的位置和大小的函数，stem(data,'fill','b-*');Matlab使用stem和stairs函数绘制离散数据，分别生成火柴棍图像和二维阶梯图像。【三】polar_plot.mPOLAR_PLO
