fx_eigen.m

著名的seismiclab的代码 是地震学研究人员必备的工具

function [m] = fx_eigen(d,p,flow,fhigh,dt,mu);
%
% Eigenfiltering 
%

[nt,nh]=size(d);

nf = 2^nextpow2(nt);

D = fft(d,nf,1);
M = zeros(nf,nh);

i = sqrt(-1);

ilow  = floor(flow*dt*nf)+1; if ilow<1; ilow=1;end;
ihigh = floor(fhigh*dt*nf)+1;if ihigh>floor(nf/2)+1; ihigh=floor(nf/2)+1;end

for k=ilow:ihigh

y = D(k,:)';
[s,w,g] = eigen_filtering(y,p,mu);

M(k,:) = s';
end;

for k=nf/2+2:nf;
M(k,:) = conj(M(nf-k+2,:));
end

m = real(ifft(M,[],1));
m = m(1:nt,:);
 
return