faquad.m

最新的HHT时频分析程序By Huang and Wu

% % t = 0:1024;
% % x = exp(-t/256).*cos(pi/64*(t.^2/512+32)+0.3*sin(pi/32*(t.^2/512+32)));
% % xn = cos(pi/64*(t.^2/512+32)+0.3*sin(pi/32*(t.^2/512+32))); 
% % am = exp(-t/256);
% % 
% % %% Instantaneous Frequency
% % trueIF = (pi*t/(64*256)+0.3*pi*t/(32*256).*cos(pi*(t.^2/512+32)/32))/(2*pi);
% % 
% % %%
% % [xnn,an1] = splinenormalize(xn);
% % [xnn,an2] = splinenormalize(xnn);
% % [xnn,an3] = splinenormalize(xnn);
% % [xnn,an4] = splinenormalize(xnn);
% % ann = an1.*an2.*an3.*an4;
% % 
% % r1 = find(abs(xnn) > 1);
% % xnn(xnn<-1) = -1; xnn(xnn>1) = 1;

%% quadrature
function f = faquad(x,dt,nfilter)