fftxzhijie.m

基于FFT的大地电磁信号资料处理MATLAB方法实现

function xp=fftxzhijie(duan,cy)
%得到各频点的复信号
%得到EH-4的X-FILE，duan为Y文件中得到的4096*4信号，T为采样时间，XP为得到的X文件。
for i=1:size(duan,2)
    x=duan(:,i);
    [yh(:,i),yw(i,:)]=fftt(x,cy);%得到每一道各频点的HIBERT 谱；
    [h(:,i),w(i,:)]=findxchazhi(yh(:,i),yw(i,:));%得到每一道与S所对应频点的HIBERT 谱；
end
[m,n]=size(duan);

%求得各频率点的自HIBERT 自谱或互谱后，再对M点平均
ZHY=(h(:,1).*conj(h(:,1)))/m;
%ZEX=(h(:,2).*conj(h(:,2)))/m;
ZHX=(h(:,3).*conj(h(:,3)))/m;
%ZEY=(h(:,4).*conj(h(:,4)))/m;

HEXHX=(h(:,2).*conj(h(:,3)))/m;
HEXHY=(h(:,2).*conj(h(:,1)))/m;
HEYHX=(h(:,4).*conj(h(:,3)))/m;
HHYHX=(h(:,1).*conj(h(:,3)))/m;
HHXHY=(h(:,3).*conj(h(:,1)))/m;
HEYHY=(h(:,4).*conj(h(:,1)))/m;

xp=zeros(length(w),9);
xp(:,1)=w(1,:)';
xp(:,2)=HEXHX;
xp(:,3)=HEXHY;
xp(:,4)=HEYHX;
xp(:,5)=HEYHY;
xp(:,6)=ZHX;
xp(:,7)=HHYHX;
xp(:,8)=HHXHY;
xp(:,9)=ZHY;
end