gusucode.com > DSISoft是由加拿大地质调查局发布的用于垂直地震剖面(VSP)数据处理的免费软件包 > dsisoftv3/dsisoftv3/dsisoftv3/demos/demoproc.m
%matlab script for processing normetal VSP data dsi_start load nm_bal %spectrally balanced data %take spectral balanced file and use median filter to subtract pwave nm_flat=flat(nm_bal,0.12,15); nm_flat=ener(nm_flat,0.10,0.14); nm_flat=trim(nm_flat,0.10,0.14,0.01); nm_flat=trim(nm_flat,0.10,0.14,0.01); nm_flat=trim(nm_flat,0.10,0.14,0.01); nm_medi=medi_filt(nm_flat,13); nm_sub=subr(nm_medi,nm_flat); nm_unflat=unflat(nm_sub,0.12,15); nm_unflat=mute(nm_unflat,1,15); %mute everything before first breaks load norpoly.mat %polygon for fk filtering nm_fk=nm_unflat; for rec=1:3 nm_fk=fkfilt(nm_fk,norpoly,225,10,0,5,rec); %pass all downgoing end %for nm_final=subr(nm_fk,nm_unflat); %subtract downgoing nm_final=ener(nm_final,0.1,0.5); dsi_end %$Id: demoproc.m,v 3.0 2000/06/13 19:18:09 gilles Exp $ %$Log: demoproc.m,v $ %Revision 3.0 2000/06/13 19:18:09 gilles %Release 3 % %Revision 2.0 1999/05/21 18:41:30 mah %Release 2 % %Revision 1.1 1999/01/06 19:07:57 kay %Initial revision % %