xsmooth.m

这是一个关于hht变换很有用的工具箱

function h=xsmooth(nt,k)
 
% The function XSMOOTH returns data that is horizontally smoothed using MATLAB Library FILTFILT function.
% 
% Data x that requires smoothing should be in column form.
%
% Calling sequence-
% h=xsmooth(nt,k)
%
% Input-
%	nt	- 2-D input data nt(in,jn)
%	k	- number of smoothing points
% Output-
%	h	- 2-D smoothed output data h(in,jn)
 
% Z.Shen (JHU)			1997 Initial

%----- Get dimensions
[in,jn]=size(nt);

%----- Initialize
h=zeros(in,jn);

%----- Apply the filter to smooth data
q1=ones(1,k)/k;
for i=1:in
   h(i,:)=filtfilt(q1,1,nt(i,:));
end