频率波数域滤波

% CREWES elastic finite-difference seismic modelling % The top boundary is a free surface suitable for propagating surface % (Rayleigh) waves. %% CREWES elastic finite-difference seismic modelling % The top boundary is a free surface suitable for propagating surface % (Rayleigh) waves. %

clear;
clc;

fm=30;
dt=0.001;
t=dt:dt:1;
dist=1:48;
t0=0.05;
tran=48;
data=zeros（100048）;
v（1）=100;v（2）=300;v（3）=500;


s=（1-2*（pi.*fm.*（t-t0））.^2）.*exp（-（pi.*fm.*（t-t0））.^2）;
dt1=1/v（1）;dt2=1/v（2）;dt3=1/v（3）;
for j=1:3
dtt=dist/v（j）;
for luoi=1:tran
ds=round（dtt（luoi）/dt）;
tem=zeros（10001）;
tem（ds:end）=s（1:length（data（ds:endluoi）））;
data（:luoi）=data（:luoi）+tem;
end
end
 seiplot（data48100010.0010.5）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         437  2013-08-23 19:54  F-K\dizhenjilu.m
     文件        4239  2013-08-23 19:54  F-K\fftrl.m
     文件        1523  2013-08-23 19:54  F-K\fk_separate.m
     文件        4145  2013-08-23 19:54  F-K\fktran.m
     文件        4032  2013-08-23 19:54  F-K\ifftrl.m
     文件        2265  2013-08-23 19:54  F-K\zhenyuanmoni.m
     目录           0  2015-09-17 13:01  F-K\