ybjc.m

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

function [hh]=ybjc(eh);
%对信号预白，加窗
  for i=1:length(eh)
      if i==1
         eh44(i)=eh(i);
      end
      if i==4096
         eh44(i)=eh(i);
      end
      if (i~=4096)&(i~=1)
         eh44(i)=-0.243*eh(i-1)+0.514*eh(i)-0.243*eh(i+1);
      end
      
   hh(i)=eh44(i)*(1-cos(1-2*pi*i/length(eh)));
   end
end